About the Authors:
J. Alberto Gallardo-Cruz
* E-mail: [email protected]
Affiliation: Departamento de Ecología y Recursos Naturales, Facultad de Ciencias, Universidad Nacional Autónoma de México, Mexico City, Mexico
Jorge A. Meave
Affiliation: Departamento de Ecología y Recursos Naturales, Facultad de Ciencias, Universidad Nacional Autónoma de México, Mexico City, Mexico
Edgar J. González
Affiliation: Departamento de Ecología y Recursos Naturales, Facultad de Ciencias, Universidad Nacional Autónoma de México, Mexico City, Mexico
Edwin E. Lebrija-Trejos
Affiliations Forest Ecology and Management Group, Wageningen University and Research Centre, Wageningen, The Netherlands, Department of Plant Sciences, Faculty of Life Sciences, Tel Aviv University, Tel Aviv, Israel
Marco A. Romero-Romero
Affiliation: Departamento de Ecología y Recursos Naturales, Facultad de Ciencias, Universidad Nacional Autónoma de México, Mexico City, Mexico
Eduardo A. Pérez-García
Affiliation: Departamento de Ecología y Recursos Naturales, Facultad de Ciencias, Universidad Nacional Autónoma de México, Mexico City, Mexico
Rodrigo Gallardo-Cruz
Affiliation: Departamento de Ecología y Recursos Naturales, Facultad de Ciencias, Universidad Nacional Autónoma de México, Mexico City, Mexico
José Luis Hernández-Stefanoni
Affiliation: Centro de Investigación Científica de Yucatán A.C., Mérida, Yucatán, Mexico
Carlos Martorell
Affiliation: Departamento de Ecología y Recursos Naturales, Facultad de Ciencias, Universidad Nacional Autónoma de México, Mexico City, Mexico
Introduction
In the dawn of the 21st century the magnitude of the human footprint on the planet's ecological systems has become undeniable [1]–[4]. Although much emphasis has been placed on the effects of industrial activities and their potential contribution to global change through greenhouse gas emissions [5]–[8], the chronic effects of land clearance for the purpose of food production on the Earth's natural vegetation are likely to be among the most long-lasting human legacies [9]–[11]. Ecologists now acknowledge that the majority of the planet's vegetation during the present century will consist of secondary or successional communities: from now on we will co-exist with secondary forests, use them, and entirely depend on them [12], [13]. The maintenance of terrestrial biodiversity will be possible as long as we are capable of keeping expanses of secondary forests [14], [15], and the regulation of the world's ecosystems will be closely linked to their existence [13], [16]. Secondary forests have also been identified as important carbon reservoirs and may play a crucial role in mitigating future global warming [17]–[27].
Vegetation ecologists currently struggle in their attempts to distinguish secondary forests from primary vegetation through remote sensing [28]–[32]. More critical, however, is the difficulty in differentiating the various successional stages that secondary forests normally comprise and measure their extent [20], [33]–[37]. As their structure and functions depend on their succesional status, there is a strong need to efficiently evaluate the extent and complexity of secondary vegetation existing in any region and to discern its attributes.
Our success in achieving these goals will depend largely on our ability to estimate accurately the structure and diversity of secondary communities at broad geographic scales. Efforts to assess the extent of secondary vegetation and to distinguish its successional variants through remote sensing have followed several routes. Most studies estimating forest structure and diversity through satellite imagery have exclusively used image spectral features and their derived vegetation indices [20], [30], [34], [38]–[45]. There are, however, several problems related to this approach; for example, some remotely-sensed vegetation indices face the problem of saturation, i.e. they are unable to discriminate different plant communities beyond a certain biomass or canopy development threshold [22], [30], [39], [46], [47]. Moreover, although some studies have succeeded in discriminating forest successional stages accurately, they have been limited to the recognition of few broad stages, which do not reflect the continuous nature of the successional process [28], [34], [48]–[52].
Recent theoretical developments in Landscape Ecology have established the link between the structural and compositional complexity of vegetation and the spatial variability of its remotely-sensed signal [53]–[55]. This spatial variability is directly related to the heterogeneity of the plant community and can be assessed by analyzing the texture of a remotely-sensed image [33], [56]–[58]. Texture refers to the spatial variation of the elements of which any image is composed [59]. Although measures of texture have been commonly used as image descriptors in remote sensing analyses [46], [60], [61], the resolution of most sensors currently employed for this purpose has prevented the examination of the internal heterogeneity of plant communities, as the size of commonly-used pixels is too large to detect such small-scale variation [37], [62]. This drawback may be overcome by using very-high resolution imagery (VHR; pixels <10 m), currently available for most of the Earth's surface, as it provides a better match between pixel size and the internal variation of vegetation [33], [56]–[58], [63]–[66]. Therefore, we hypothesized that the textural information contained in VHR images has the potential to reflect the variability of secondary vegetation, allowing us to model the successional process.
The goal of this study was to examine the potential of textural properties of a VHR Quickbird image to model secondary vegetation attributes measured in the field, in a seasonally dry tropical region. We wanted to test the power of the texture of remote images to describe and predict vegetation attributes, while identifying those texture attributes with the highest predictive potential. In modeling the relationship between textural and vegetation attributes, we succeeded in producing simple models that can be easily obtained for many regions, and that have a straightforward biological interpretation.
Materials and Methods
Study area
The study was conducted in the dry tropical region of Nizanda, Oaxaca State, Mexico (16°39.49′N, 95°0.66′W; Fig. 1). Mean annual temperature is 26°C and the total average annual rainfall is 900 mm, largely concentrated between June and October. The prevailing vegetation matrix is a low-stature (7–10 m) seasonally dry tropical forest [67], [68]. Traditional slash-and-burn agriculture is practiced in the area. Fields are typically cropped for one or two years before being abandoned [69], which results in a mosaic of differently-aged fallows spread across the area.
[Figure omitted. See PDF.]
Figure 1. Study area (UTM zone 15n) and location of the secondary plots (▪) used for modeling their attributes from the texture derived from a Quickbird satellite image.
https://doi.org/10.1371/journal.pone.0030506.g001
Field data
Fourteen 30×30 m fallows with time since abandonment (age) ranging from 2 to ca. 60 years, and one mature forest site, were selected from field surveys conducted in 2005. Site age was obtained from interviews with landowners and verified through dendrochronology [70]. This set of fallows included a large range of environmental and vegetational heterogeneity [71], [72], from sites with a little dense canopy dominated by shrubs of open foliage, to sites with a dense plant cover and a low percentage of bare soil.
The sites were marked and designated as permanent sites in 2003 [69], [71], and detailed structural and floristic information was collected annually (Table 1 lists the variables used for vegetation description and their abbreviations). In each fallow, four 20×5 m transects (400 m2, subdivided in four 5×5 m quadrats) were established for the sampling of woody and succulent plants. In each transect, all individuals ≥5 cm DBH were sampled in the four quadrats; individuals with DBH≥2.5 cm but <5 cm were sampled in two quadrats, and individuals with ≥1 but <2.5 cm DBH were sampled in one quadrat only. For each individual, DBH and two orthogonal crown diameters (used to calculate crown areas) were measured. Structural variables were obtained by scaling the data to 1 ha. Based on this information, for each plot the T and U sets (i.e. Total and Upper, respectively) of structural and diversity attributes were prepared. The T set included all sampled plants in the plot, whereas the U set included only those plants that are more likely to be remotely sensed [73]. The U set comprised those plants above the median in the frequency distribution of canopy cover; this subset represented between 50% and 75% of the basal area in a sites. For these two sets we calculated Dn (individuals in the sampled area), CC (the sum of the individual crown areas), BA (the sum of individual basal areas), S (number of species), and Simpson's D' and Shannon's H' diversity indices [74]. In addition, Hgt was calculated as the average of the heights of eight trees, each the tallest tree in the zones formed by two adjacent quadrats. Structural and diversity data for the study sites are shown in Table 2; please refer to Table S1 for information on within-site variability for those variables for which the calculation of such variation is feasible and sensible.
[Figure omitted. See PDF.]
Table 1. Vegetational attributes used in the analysis and their abbreviations.
https://doi.org/10.1371/journal.pone.0030506.t001
[Figure omitted. See PDF.]
Table 2. Structural and diversity attribute values for 15 plots.
https://doi.org/10.1371/journal.pone.0030506.t002
Image processing
We used a high-resolution Quickbird satellite image (pixel size = 2.6 m) acquired in early December 2005. This date, which corresponds to the beginning of the dry season, was chosen to minimize cloud cover while ensuring the presence of foliage in the plants. The image was geometrically and atmospherically corrected to surface reflectance following Krause [75].
From the four available bands in this image, we selected the red (RED; 0.63–0.69 µm) and the near infrared (IR; 0.76–0.90 µm) [41], both of which are known to reflect the condition of vegetation functioning and overall condition. We also estimated two commonly used vegetation indices, namely the normalized difference vegetation index (NDVI) and the enhanced vegetation index (EVI) [47], both derived from the combination of the two bands, as follows:EVI incorporates empirical parameters (C1 = 6, C2 = 7.5) and the blue band (BLUE, 0.45–0.52 µm) for atmospheric correction, and sensitivity minimization of soil background reflectance variation (L = 1). The EVI does not saturate under dense canopy conditions as the NDVI does, and it also appears to be more sensitive to canopy structural characteristics [76].
Image texture analysis
Image texture refers to the spatial variation and arrangement of the pixels of which any image is composed [59], [77]. Although this property can be extracted through a wide array of methods [77], [78], we chose to follow a statistical approach for this study. We calculated texture variables known as first-order and second-order measurements. First-order texture measures are statistical properties that do not consider pixel neighbor relationships and are derived from the original image values within a certain window (group of pixels); for this group of texture measures, the spectral variability within the window was assessed by calculating the range and the skewness of the values.
Unlike the textural variables from the first group, second-order measurements consider the spatial relations between groups of two neighboring pixels within the window [59], therefore these measurements were also selected because of their greater potential to reflect the heterogeneity in successional vegetation stands. The calculation of second-order variables involves the construction of Gray-Level Co-occurrence Matrices (GLCMs), which are matrices containing the probabilities of co-occurrence of pixel values for pairs of pixels in a given direction and distance. To construct such matrices we used a spatial distance of one pixel, four directions (0°, 45°, 90°, 135°), and 64 levels of gray. A co-occurrence matrix was constructed for each direction, and from each co-occurrence matrix a specific texture measurement was calculated for the window. The texture measurements of each direction were then averaged to obtain a single spatially-invariant texture value. This procedure was applied for the variables described in Table 3, which correspond to three groups of variables describing the degree of contrast between pixels (homogeneity, contrast, dissimilarity), the regularity in the pixels within a window (entropy, angular second moment), and the statistics derived from the GLCM (mean, variance, correlation).
[Figure omitted. See PDF.]
Table 3. Texture variables derived from the grey-level co-occurrence matrix (GLCM).
https://doi.org/10.1371/journal.pone.0030506.t003
The ten textural variables (two first-order and eight second-order variables) were calculated for the RED and the IR bands, as well as for the two vegetation indices (NDVI and EVI). We used a moving-window approach with a window size of 15 pixels to match the size of the sampling plots in the field; the central pixel value of this window was extracted from each of the 40 texture layers (four layers and ten variables). The entire procedure was programmed in the ENVI+IDL environment [79].
Statistical analysis
We assessed the potential of the 40 texture variables to describe the observed changes in each of the 14 vegetation variables by means of linear models. In all cases the response variables were log-transformed before model fitting. This procedure guaranteed that the estimated values, when transformed to their original scale, would be positive; additionally, such transformation tends to homogenize the residuals, which are usually proportional to the mean in most probability density functions restricted to the positive numbers domain.
Assuming that different texture variables provide supplementary information about the remotely-sensed vegetation, we fitted three types of models, depending on whether they included one (560 models), two (10,920) or three (138,320) texture variables; no interaction was examined due to the limited degrees of freedom. For each of the three model types we selected the one having the largest coefficient of determination (R2). For our analysis it was crucial to be able to compare the potential of textural variables to model different vegetational variables; for this purpose, R2 was appropriate owing to its fixed range (from 0 to 1), unlike the commonly used Akaike's Information Criterion (AIC), whose range depends, among other things, on the sum of squares of the error of the dependent variable [80], making it impossible to compare models produced for different vegetational attributes.
Nevertheless, R2 does not provide a good measure to determine whether, for any given vegetational variable, increasing the number of variables incorporated into the model improves its performance. Therefore, we compared the best-fit models of each type (i.e. with one, two, or three textural variables) through the small-sample-size bias-correction version of AIC (AICc); two models were considered to be equally good when the difference in AICc between them was <2 [80].
Due to the large number of models that were fitted and the small data set, it was expected that large R2 values would be obtained by chance. To minimize this possibility, we produced null models by randomly sorting the texture- and vegetation-attribute data, fitting the same linear models as above, and selecting those with the largest R2-values. This procedure was repeated 1,000 times and an empirical distribution of the largest expected R2 was obtained. We then estimated the P-value associated with each model as the fraction of the simulated R2 values that were greater than the observed ones. The median of the empirical distribution was used as a measure of the expected magnitude of R2 under a completely random scenario.
To assess the predictive power of the models, we used leave-two-out cross-validation, i.e. we used a linear model fitted to the data from 13 plots (the calibration data) to predict the vegetation attribute of the remaining two plots (the validation data). We fitted a model for each possible split of the data set into calibration and validation subsets, and calculated the sum of squares between the estimated and observed values of the vegetation variable for the two validation plots. The model with the highest predictive power would be the one with the smallest average sum of squares () averaged over all the possible splits of the data set. We used a leave-two-out procedure to avoid the problems arising from the more popular leave-one-out approach [81]–[83], while keeping a sample size as large as possible for the estimation of the four parameters of the most complex models (those based on three texture variables).
The actual values of cannot be compared between vegetation attributes as they depend on the scale of measurement. A more informative statistic is given by the proportion of the variation of the vegetation attribute that can be predicted by the model. This statistic, expressed as a leave-d-out cross-validation R2CV, is:where yi is the vegetation attribute for the ith plot, and n is the number of plots. R2CV ranges from −∞ to 1, where a value of 1 means that a model predicts the validation data perfectly, while a negative value means that a model is over-fitted, because it would make worse predictions than those made by a null model. This statistic is dimensionless, so it can be used for comparing vegetation variables.
All statistical analyses were performed using R [84].
Results
Structural and diversity attributes of secondary vegetation
Two sets of structural and diversity attributes (T and U) were analyzed for each sampled plot (see Table 2). For the T set, which included all sampled plants in the plot, most structural variables exhibited clear increasing trends with successional development. However, this was not the case for DnT, which did not show a clear successional pattern, and for DT (dominance), which tended to decrease from young to old fallows. Some variables clearly showed stabilizing trends, particularly Hgt and CCT, or non-monotonic responses, as was the case of ST.
The U set was arbitrarily defined as those trees that were above the median of the canopy cover cumulative distribution of each fallow. The proportion of BAU with respect to BAT was generally high, mostly above 60%, and often higher in mid-age fallows, despite the very low number of individuals included in this community's subset (in all cases less than 2% of stem density of the entire community).
Descriptive models of successional vegetation attributes
We constructed three types of linear models to describe the relationship between the 14 vegetation attributes and the 40 image texture variables. These types corresponded to the number of texture variables used in constructing the model: one, two or three texture variables.
Most of the best fit models had relatively high significant R2 values, and these increased as more variables were included. Among one-variable models, five out of 14 models had R2>0.80, and this number increased to 10 and 12 for two- and three-variable models, respectively (Table 4). According to AICc, the latter were always the best-fit models (Table 4, TV entries in bold typeface). However, it must be noted that for half of the 14 vegetational variables the best two-variable models were equally good.
[Figure omitted. See PDF.]
Table 4. Best descriptive linear models for 14 vegetation attributes (VA) as a function of 1, 2 and 3 textural variables (TV) with corresponding R2, P and R2CV values.
https://doi.org/10.1371/journal.pone.0030506.t004
We measured the expected magnitude of R2 under a completely random scenario through the median of its empirical distribution. The difference between this median and the observed one decreased as more variables were included, from 0.37 for models with one variable, to 0.24 and 0.11 in two- and three-variable models, respectively (Fig. 2). Also, the proportion of non-significant best-fit models increased from 0.14, through 0.21 to 0.43 as more variables were included. There were small, non-significant differences (paired t = 0.190, P = 0.851) between the R2 values calculated using the two data sets (T and U).
[Figure omitted. See PDF.]
Figure 2. Fraction of the variation in vegetation attributes (median and range) explained by the descriptive (—♦—), predictive (- -○- -) and null (—▴—) models using a different number of textural attributes as explanatory variables.
For the descriptive and null models conventional R2 values are reported, while for the predictive models R2CV is reported, so the values are not strictly comparable (see Methods for explanation).
https://doi.org/10.1371/journal.pone.0030506.g002
BAT and BAU were the response variables for which the best descriptive models were obtained (Table 4). For the best two-texture variable models, R2 values were 0.93 and 0.92 for BAT and BAU, respectively. These vegetation attributes had an R2 value of 0.96 for the best three-texture variable models. Conversely, R2 for the best descriptive model with one textural variable was 0.82. Unlike most vegetation variables, the best-fit models for Dn had much lower R2 values, and these did not differ significantly from the values derived from the null model. The same occurred for the best-fit models for H' and D indices, both for the entire community and for the upper canopy, in the case of the three texture-variable models.
Comparison between the three types of models revealed that the most complex models did not necessarily incorporate the same variables as simpler models (Table 4). Most descriptive one-variable models included REDVAR or REDMEAN as the best explanatory variables of the behavior of vegetation variables (See Methods and Table 3 for a full description of texture variables and their abbreviations); NDVICORR and NDVISM were important only for DnU and DnT. For two-variable models, either REDVAR or REDMEAN were retained in six models only, whereas variables incorporating textural information derived from vegetation indices (NDVI and EVI) became prominent. When moving to three-texture variable models, IRMEAN emerged in eight models as capable of making a significant contribution to the descriptive power of the models. Conversely, in this set REDMEAN and other RED-related textural variables became much less important, which indicates their limited descriptive ability in the presence of other textural variables.
The signs of the coefficients of the textural variables in the models changed according to the way in which they relate to the different response variables. Within the group of one-variable models, DT and DU were the only vegetation variables whose models had positive coefficients associated to the textural variables. This relation is less obvious for two- and three-texture variable models, yet the coefficients of the textural variables still have different sign when involved in the descriptive modeling of dominance as opposed to other vegetation variables.
Predictive models of successional vegetation attributes
Most (27 out of 28 models) one- or two-texture variable best-fit models predicting vegetation response variables were identical to the respective descriptive models regarding the identity of the explanatory variables (see Table S2). Moreover, Spearman correlations between all possible R2 and cross-validation R2 (R2CV) pairs of values for each vegetation variable and group of texture variables were very high. These correlations, which were higher for the one-texture variable models, were inversely related to the number of texture variables involved (Table 5). In general, R2CV values were lower than R2 values in descriptive models, as low as 25%, but more often around 10% lower (Fig. 2, Table 4). Despite such reduction, BAU and BAT had R2CV>0.90 in three-texture variable models. As was the case with descriptive models, there were no significant differences between the R2 values of predictive models developed from the upper canopy and total sets (paired t = 0.634, P = 0.534).
[Figure omitted. See PDF.]
Table 5. Spearman's ρ between descriptive R2 and predictive R2CV values calculated for all linear models resulting from modeling each vegetation attribute as a function of one, two and three textural variables (TV).
https://doi.org/10.1371/journal.pone.0030506.t005
Departure of predictive models from descriptive ones occurred mostly within the set of three-texture variable models (Fig. 3, Table S2). For the new set of predictive models, textural variables derived from the RED band (i.e. REDVAR or REDMEAN) became prominent again among models with high R2CV.
[Figure omitted. See PDF.]
Figure 3. Observed (x-axes) vs. estimated (y-axes) values for the best descriptive (▴) and predictive (red +) linear models for vegetation attributes.
See Table 1 for vegetational attributes abbreviations. Digits 1, 2, and 3 refer to the number of textural variables included in the model as explanatory variables.
https://doi.org/10.1371/journal.pone.0030506.g003
Discussion
Predictive potential of satellite image texture
In this study we demonstrate the large potential of image texture for predicting vegetation attributes during tropical forest succession. Texture is an emergent property of satellite images that is related to the neighborhood relationships among pixels [59], and thus it is capable of reflecting the internal organization (i.e. heterogeneity, directionality, entropy) of a region of interest, rather than on its mean properties. This seems to be the reason why the performance of texture-based analyses tends to exceed those based on spectral information in discriminating different successional stages [33], [46], [56]–[58], [62], [85]–[87]. High R2 values comparable to those obtained by us have been reported by some studies, but only after complex image processing and modeling protocols based on mean canopy reflectance [30], [45], [73], [88]–[90]. Apparently this complexity has limited the broad application of these procedures, hence motivating the ongoing search for simpler solutions that are useful in a variety of circumstances. The method proposed in this study contrasts by its simplicity: the analysis was performed with a single image, and the models were linear and included few variables. Moreover, textural information can presently be extracted with ease. This simplicity, which becomes an asset in studying secondary vegetation and its attributes, depends on the basic principle that image texture actually reflects the internal heterogeneity of successional vegetation at the proper scale [62].
The ability to predict characteristics of secondary vegetation accurately depends on a combination of three relevant methodological aspects, all of which synergically contribute to the high predictive value of the models.
The first aspect is high image resolution. Typically, scholars interested in predicting vegetation attributes from space have used images with pixel size ≥30 m [28], [73], [88], [89], [91]. Agricultural fields derived from non-mechanized practices in tropical dry regions often have relatively small sizes (in our study area, the sizes of most successional stands range from 900 to 2,500 m2); thus a single such large pixel covers just one secondary vegetation stand. Therefore, high spatial resolution is required to detect and analyze the internal spatial variation typical of each secondary stand. Proisy et al. [33] came to a similar conclusion while mapping biomass in successional mangrove communities.
The second aspect was the inclusion of stand-level heterogeneity, an essential feature of successional vegetation. This inclusion was achieved by using a range of image textural attributes, some of which proved to have a very high predictive potential, even though our study also shows that many textural attributes do not have such potential (Table 4 shows nine variables that were included in only one model, in addition to 17 variables that were not included in any of them), in agreement with other studies [60], [62].
The third aspect was the decision to assess and contrast two large sets of models derived from alternative modeling procedures: one set included a limited number of descriptive models that included all sampling sites, whereas the other consisted of numerous predictive models constructed through leave-two-out cross-validation. The high degree of consistency between the models selected from either procedure confers increased reliability to the results, and implies that constructing descriptive models may suffice for assessing the secondary vegetation in a region. This is a valuable result as the construction of predictive models may require a large computational capacity as well as ample programming and statistical skills.
It is not uncommon for this kind of studies to face a limitation derived from the high cost of obtaining field information for every vegetation stand; thus having a large sample size, which would increase the accuracy in the predictions of the models, may not be feasible. This limitation was a strong motivation for this investigation. One would expect the prediction of vegetation attributes for new plots using our models to be flawed in two cases, neither of which occurred in our study. The first case would be if the models were used to estimate the attributes of plots with ages beyond those used in model fitting (extrapolation). We did not need to extrapolate because our plots represented the broadest possible successional gradient. The second case would be if the model's estimated coefficients were inaccurate due to a small sample size. We avoided this problem by using the leave-two-out cross-validation procedure [82]. Our high predictive R2CV values confirm that even a model based on a rarified sample was capable of providing reliable estimates for the vegetation attributes of new sites and warrants that our conclusions are not the artificial result of a small sample size.
An unanticipated conclusion from our study is that two-texture variable models should be preferred over three variable ones for describing and predicting vegetational attributes from image texture. This conclusion derives from two different results. On the one hand, AICc indicated that three-texture variable models were better than those with two variables only for half of the vegetational attributes, whereas for the other half both types of models were equally good. On the other, as the number of variables included in the models increased, the departure between observed and null R2 distributions decreased, rendering three-variable models less reliable than two-variable ones. This conclusion is also relevant from a practical perspective, as the construction and validation of three-texture variable models requires much larger computing time and costs.
Despite previous suggestions that forest structure and diversity characteristics are preferably predicted from canopy-reflectance information [20], [73], in our case restricting the analysis to the upper canopy did not necessarily result in a better predictive capacity. In fact, our models predicting BAT had higher R2 values than BAU. In the case of an analysis based on texture of VHR imagery, predicting total community or upper canopy attributes can be done with comparable accuracy.
Even though we were able to demonstrate a high potential of GLCM textural indices to predict successional vegetation attributes, some caution must be exerted in using them. Like other indices, GLCM face potential important limitations that must be acknowledged. A particularly worrisome one is the fact that texture may be sensitive to image sun-view acquisition conditions [56]. Recently, Barbier et al. [92] proposed a mitigation method for FOTO (Fourier Transform Textural Ordination) indices that seems promising, albeit expensive and not totally straightforward. Further research is required aimed to develop a similar procedure to GLCM indices.
The significance of image textural information
Understanding why some textural attributes are more useful than others in predicting vegetation properties, and therefore why they were repeatedly incorporated into the models, is important for a number of reasons. From a practical perspective this knowledge will orient future efforts to assess secondary vegetation by guiding researchers as to which variables they should focus on. Also, this information will provide a firmer ground for theoretical inquiry, as it represents an efficient way to identify relevant biological properties of the vegetation system and its spatially explicit spectral expression.
In this study it became clear that the three textural attributes that excelled in their predictive capabilities were IRMEAN, REDVAR and REDMEAN (Table 4), which indicates that in the context of texture, the predictive potential of the raw information contained in these bands exceeds that of NDVI and EVI, both based on RED and IR [47]. This implies that in calculating these indices the relevant spatial information that reflects the internal heterogeneity is lost. Under an approach centered on the examination of the internal heterogeneity of successional stands such loss of information is crucial at this scale of analysis, as objects that are not well differentiated spectrally can be finely discerned [20].
One finding that deserves particular attention is the inverse relationship between satellite-sensed heterogeneity, in particular the mean and the variance of textural variables, and ground-level vegetation development (i.e. stand age and other vegetation attributes). Large mean values obtained from a GLCM denote high levels of between-pixel spectral heterogeneity. Likewise, large GLCM variances indicate that such changes are highly variable regardless of the mean change. Therefore, although the inverse texture/vegetation relationship may seem counterintuitive and even contradictory to some recent findings reported in the literature [93]–[96], the explanation might lie in comprehending what the satellite actually perceives. One conceivable explanation is that the pixels corresponding to an early successional stand not only contain the reflectance properties of the plants, but also the spectral properties of the substrate on which they grow. If this interpretation is correct, it follows that textural attributes should show decreasing trends as vegetation structure becomes more complex and covers the soil. Thus the internal heterogeneity of a mature successional stand would be mostly related to the differences between less contrasting reflectance properties of the plants.
Our research is in agreement with other studies that have shown texture of satellite imagery to be closely related to the heterogeneity of the vegetation stand [57], [58]. For example, Frazer et al. [97] reported that LiDAR-derived indices such as lacunarity (the degree to which an object departs from a geometric pattern) are sensitive to canopy structure attributes. Despite the obvious ability of GLCM indices to reflect community-level attributes (e.g. total basal area, stand age), it is not yet clear how they relate to individual-level or other finer-scale traits (e.g. crown size). Thus, a promising line of future research will consist in finding out how these textural metrics relate to fine-scale vegetation and overall stand properties.
Potential applications of image-texture-based modeling
The application of the method described here may produce important information related to two of the most relevant threats to biosphere integrity: climate change and biodiversity loss [98]–[101]. Basal area, the vegetation variable that was best predicted from image textural attributes, is strongly correlated with the standing biomass of a forest community [102], and thus to carbon storage [103]. Carbon sequestration rates may also be obtained by considering stand age [36], another variable accurately predicted from image texture. Therefore, by applying this procedure, it should be possible to assess and map with high confidence the spatial distribution of the potential carbon storage and sequestration in regions dominated by secondary vegetation in different stages of development.
Biodiversity conservation is one of the major goals of tropical ecologists nowadays [104]–[106]. Several efforts have been made recently to assess the possibility that local floras and faunas may persist in regions where native vegetation has undergone major transformations [14], [107]–[109]. Therefore, the possibility to predict species richness is of utmost importance. Our results show that species richness can be predicted with a precision close to 80%. This figure implies the existence of a relationship between the occurrence of different species in the terrain and the information sensed by a satellite. At present this topic is receiving much attention from researchers [55], [65], [110], [111], and our analysis opens new avenues to pursue it.
Canopy cover and vegetation height were also well predicted by our models. Again, there are several potential applications of this result. For example, information on canopy cover in a region dominated by secondary vegetation may help in assessing the potential soil erosion due to the kinetic energy of rainfall [112]–[114]. Similarly, it will provide information that can be used to assess habitat quality for a regional fauna, particularly for those animals whose survival depends on a closed canopy [115]–[117].
Concluding remarks
Analyzing the extent and complexity of secondary vegetation by recognizing the spatial variation of its spectral information opens new and attractive research avenues; these differ substantially from previous efforts to study secondary vegetation that have been primarily based on the examination of spectral reflectance properties. Overall, the procedure is potentially usable in any successional plant community whose development involves large changes in heterogeneity through time. The current availability of VHR imagery, together with increasing computing capabilities, make it possible to develop faster and more efficient ways to assess the amount and condition of secondary vegetation in increasingly human-impacted regions world-wide based on this approach.
Supporting Information
[Figure omitted. See PDF.]
Table S1.
Within-site variability for vegetational variables. See Table 1 for vegetational attributes abbreviations and units of measurement. SE: Standard error. M: mature forest.
https://doi.org/10.1371/journal.pone.0030506.s001
(XLS)
Table S2.
Best descriptive and predictive linear models for vegetation attributes as a function of one, two and three textural variables (TV). See Table 1 for vegetational attributes abbreviations. IR: near infra-red band, RED: red band, NDVI: Normalized Difference Vegetation Index, EVI: Enhanced Vegetation Index. See Table 3 for the description of textural variables denoted by subindexed terms MEAN, VAR, ASM, CORR, DR, SKEW, DISS, CONT, ENT, and HOM. Only those best predictive models that differed from the best descriptive ones are presented.
https://doi.org/10.1371/journal.pone.0030506.s002
(XLS)
Acknowledgments
Armando Peralta provided support in image processing. JAGC and EJG thank the Graduate Program in Biological Sciences of the Universidad Nacional Autónoma de México for providing continuous support during their PhD studies. Thanks to Trudy Kavanagh for revision and edition of the manuscript. We warmly thank the people of Nizanda for their hospitality, permission to work on their land and assistance during fieldwork. We thank two anonymous reviewers for commenting on a previous draft of this work.
Author Contributions
Conceived and designed the experiments: JAGC JAM EJG. Performed the experiments: JAGC JAM EJG EEL MAR EAP RGC JLH CM. Analyzed the data: JAGC JAM EJG EEL MAR EAP RGC CM. Contributed reagents/materials/analysis tools: JAGC JAM EJG EEL MAR EAP RGC JLH CM. Wrote the paper: JAGC JAM EJG EEL MAR EAP RGC JLH CM.
Citation: Gallardo-Cruz JA, Meave JA, González EJ, Lebrija-Trejos EE, Romero-Romero MA, Pérez-García EA, et al. (2012) Predicting Tropical Dry Forest Successional Attributes from Space: Is the Key Hidden in Image Texture? PLoS ONE 7(2): e30506. https://doi.org/10.1371/journal.pone.0030506
1. Hassan R, Scholes R, Ash N (2005) Ecosystems and Human Well-being: current state and trends. Washington, D.C.: Island Press. 922 p.R. HassanR. ScholesN. Ash2005Ecosystems and Human Well-being: current state and trendsWashington, D.C.Island Press922
2. Janzen D (1998) Gardenification of wildland nature and the human footprint. Science 279: 1312–1313.D. Janzen1998Gardenification of wildland nature and the human footprint.Science27913121313
3. Lepers E, Lambin E, Janetos A, DeFries R, Achard F, et al. (2005) A synthesis of information on rapid land cover change for the period 1981–2000. Bioscience 55: 115–124.E. LepersE. LambinA. JanetosR. DeFriesF. Achard2005A synthesis of information on rapid land cover change for the period 1981–2000.Bioscience55115124
4. Sanderson EW, Jaiteh M, Levy MA, Redford KH, Wannebo AV, et al. (2002) The human footprint and the last of the wild. Bioscience 52: 891–904.EW SandersonM. JaitehMA LevyKH RedfordAV Wannebo2002The human footprint and the last of the wild.Bioscience52891904
5. Fearnside P, Laurance W (2004) Tropical deforestation and greenhouse-gas emissions. Ecol Appl 14: 982–986.P. FearnsideW. Laurance2004Tropical deforestation and greenhouse-gas emissions.Ecol Appl14982986
6. Bonan G (2008) Forests and climate change: forcings, feedbacks, and the climate benefits of forests. Science 320: 1444–1449.G. Bonan2008Forests and climate change: forcings, feedbacks, and the climate benefits of forests.Science32014441449
7. Satterthwaite D (2008) Cities' contribution to global warming: notes on the allocation of greenhouse gas emissions. Environ Urban 20: 539–549.D. Satterthwaite2008Cities' contribution to global warming: notes on the allocation of greenhouse gas emissions.Environ Urban20539549
8. Thomson A, Calvin K, Chini L, Hurtt G, Edmonds J, et al. (2010) Climate mitigation and the future of tropical landscapes. P Natl Acad Sci USA 107: 19633–19638.A. ThomsonK. CalvinL. ChiniG. HurttJ. Edmonds2010Climate mitigation and the future of tropical landscapes.P Natl Acad Sci USA1071963319638
9. Foley J, DeFries R, Asner G, Barford C, Bonan G, et al. (2005) Global consequences of land use. Science 309: 570–574.J. FoleyR. DeFriesG. AsnerC. BarfordG. Bonan2005Global consequences of land use.Science309570574
10. Foley J, Monfreda C, Ramankutty N, Zaks D (2007) Our share of the planetary pie. P Natl Acad Sci USA 104: 12585–12586.J. FoleyC. MonfredaN. RamankuttyD. Zaks2007Our share of the planetary pie.P Natl Acad Sci USA1041258512586
11. Nelson G (2005) Drivers of ecosystem change: summary chapter. In: Hassan R, Scholes R, Ash N, editors. Ecosystems and human well-being: current state and trends. Washington, D.C.: Island Press. pp. 73–76.G. Nelson2005Drivers of ecosystem change: summary chapter.R. HassanR. ScholesN. AshEcosystems and human well-being: current state and trendsWashington, D.C.Island Press7376
12. Chazdon R, Harvey C, Komar O, Griffith D, Ferguson B, et al. (2009) Beyond reserves: a research agenda for conserving biodiversity in human-modified tropical landscapes. Biotropica 41: 142–153.R. ChazdonC. HarveyO. KomarD. GriffithB. Ferguson2009Beyond reserves: a research agenda for conserving biodiversity in human-modified tropical landscapes.Biotropica41142153
13. Guariguata M, Ostertag R (2001) Neotropical secondary forest succession: changes in structural and functional characteristics. For Ecol Manage 148: 185–206.M. GuariguataR. Ostertag2001Neotropical secondary forest succession: changes in structural and functional characteristics.For Ecol Manage148185206
14. Chazdon RL, Peres CA, Dent D, Sheil D, Lugo AE, et al. (2009) The potential for species conservation in tropical secondary forests. Conserv Biol 23: 1406–1417.RL ChazdonCA PeresD. DentD. SheilAE Lugo2009The potential for species conservation in tropical secondary forests.Conserv Biol2314061417
15. Dent DH, Wright SJ (2009) The future of tropical species in secondary forests: a quantitative review. Biol Conserv 142: 2833–2843.DH DentSJ Wright2009The future of tropical species in secondary forests: a quantitative review.Biol Conserv14228332843
16. Chinea J, Helmer E (2003) Diversity and composition of tropical secondary forests recovering from large-scale clearing: results from the 1990 inventory in Puerto Rico. For Ecol Manage 180: 227–240.J. ChineaE. Helmer2003Diversity and composition of tropical secondary forests recovering from large-scale clearing: results from the 1990 inventory in Puerto Rico.For Ecol Manage180227240
17. Aide T, Zimmerman J, Pascarella J, Rivera L, Marcano H (2000) Forest regeneration in a chronosequence of tropical abandoned pastures: implications for restoration ecology. Restor Ecol 4: 328–338.T. AideJ. ZimmermanJ. PascarellaL. RiveraH. Marcano2000Forest regeneration in a chronosequence of tropical abandoned pastures: implications for restoration ecology.Restor Ecol4328338
18. Brearley F, Prajadinata S, Kidd P, Proctor J, Suriantata (2004) Structure and floristics of an old secondary rain forest in Central Kalimantan, Indonesia, and a comparison with adjacent primary forest. For Ecol Manage 195: 385–397.F. BrearleyS. PrajadinataP. KiddJ. ProctorSuriantata2004Structure and floristics of an old secondary rain forest in Central Kalimantan, Indonesia, and a comparison with adjacent primary forest.For Ecol Manage195385397
19. Canadell J, Pataki D, Gifford R, Houghton R, Luo Y, et al. (2007) Saturation of the terrestrial carbon sink. In: Canadell J, Pataki D, Pitelka L, editors. Terrestrial ecosystems in a changing world. Berlin: Springer-Verlag. pp. 59–78.J. CanadellD. PatakiR. GiffordR. HoughtonY. Luo2007Saturation of the terrestrial carbon sink.J. CanadellD. PatakiL. PitelkaTerrestrial ecosystems in a changing worldBerlinSpringer-Verlag5978
20. Castro K, Sanchez-Azofeifa G, Rivard B (2003) Monitoring secondary tropical forests using space-borne data: implications for Central America. Int J Remote Sens 24: 1853–1894.K. CastroG. Sanchez-AzofeifaB. Rivard2003Monitoring secondary tropical forests using space-borne data: implications for Central America.Int J Remote Sens2418531894
21. Ewel J (1977) Differences between wet and dry successional tropical ecosystems. Geo-Eco-Trop 1: 103–117.J. Ewel1977Differences between wet and dry successional tropical ecosystems.Geo-Eco-Trop1103117
22. Helmer E, Brandeis T, Lugo A, Kennaway T (2008) Factors influencing spatial pattern in tropical forest clearance and stand age: implications for carbon storage and species diversity. J Geophys Res 113: G02S04.E. HelmerT. BrandeisA. LugoT. Kennaway2008Factors influencing spatial pattern in tropical forest clearance and stand age: implications for carbon storage and species diversity.J Geophys Res113G02S04
23. Jepsen M (2006) Above-ground carbon stocks in tropical fallows, Sarawak, Malaysia. For Ecol Manage 225: 287–295.M. Jepsen2006Above-ground carbon stocks in tropical fallows, Sarawak, Malaysia.For Ecol Manage225287295
24. Kenzo T, Ichie T, Hattori D, Kendawang J, Sakurai K, et al. (2010) Changes in above- and belowground biomass in early successional tropical secondary forests after shifting cultivation in Sarawak, Malaysia. For Ecol Manage 260: 875–882.T. KenzoT. IchieD. HattoriJ. KendawangK. Sakurai2010Changes in above- and belowground biomass in early successional tropical secondary forests after shifting cultivation in Sarawak, Malaysia.For Ecol Manage260875882
25. Rhoades CC, Eckert GE, Coleman DC (2000) Soil carbon differences among forest, agriculture, and secondary vegetation in lower montane Ecuador. Ecol Appl 10: 497–505.CC RhoadesGE EckertDC Coleman2000Soil carbon differences among forest, agriculture, and secondary vegetation in lower montane Ecuador.Ecol Appl10497505
26. Shevliakova E, Pacala S, Malyshev S, Hurtt G, Milly P, et al. (2009) Carbon cycling under 300 years of land use change: Importance of the secondary vegetation sink. Global Biogeochem Cy 23: GB2022.E. ShevliakovaS. PacalaS. MalyshevG. HurttP. Milly2009Carbon cycling under 300 years of land use change: Importance of the secondary vegetation sink.Global Biogeochem Cy23GB2022
27. Yang X, Richardson T, Jain A (2010) Contributions of secondary forest and nitrogen dynamics to terrestrial carbon uptake. Biogeosciences 7: 3041–3050.X. YangT. RichardsonA. Jain2010Contributions of secondary forest and nitrogen dynamics to terrestrial carbon uptake.Biogeosciences730413050
28. Arroyo-Mora J, Sanchez-Azofeifa G, Kalacska M, Rivard B, Calvo-Alvarado J, et al. (2005) Secondary forest detection in a neotropical dry forest landscape using Landsat 7 ETM+ and IKONOS imagery. Biotropica 4: 497–507.J. Arroyo-MoraG. Sanchez-AzofeifaM. KalacskaB. RivardJ. Calvo-Alvarado2005Secondary forest detection in a neotropical dry forest landscape using Landsat 7 ETM+ and IKONOS imagery.Biotropica4497507
29. Carreiras J, Pereira J, Campagnolo M, Shimabukuro Y (2006) Assessing the extent of agriculture/pasture and secondary succession forest in the Brazilian Legal Amazon using SPOT VEGETATION data. Remote Sens Environ 101: 283–298.J. CarreirasJ. PereiraM. CampagnoloY. Shimabukuro2006Assessing the extent of agriculture/pasture and secondary succession forest in the Brazilian Legal Amazon using SPOT VEGETATION data.Remote Sens Environ101283298
30. Nelson R, Kimes D, Salas W, Routhier M (2000) Secondary forest age and tropical forest biomass estimation using Thematic Mapper imagery. Bioscience 50: 419–431.R. NelsonD. KimesW. SalasM. Routhier2000Secondary forest age and tropical forest biomass estimation using Thematic Mapper imagery.Bioscience50419431
31. Tottrup C, Rasmussen M, Samek J, Skole D (2007) Towards a generic approach for characterizing and mapping tropical secondary forests in the highlands of mainland Southeast Asia. Int J Remote Sens 28: 1263–1284.C. TottrupM. RasmussenJ. SamekD. Skole2007Towards a generic approach for characterizing and mapping tropical secondary forests in the highlands of mainland Southeast Asia.Int J Remote Sens2812631284
32. Wunderle A, Franklin S, Guo X (2007) Regenerating boreal forest structure estimation using SPOT-5 pan-sharpened imagery. Int J Remote Sens 28: 4351–4364.A. WunderleS. FranklinX. Guo2007Regenerating boreal forest structure estimation using SPOT-5 pan-sharpened imagery.Int J Remote Sens2843514364
33. Proisy C, Couteron P, Fromard F (2007) Predicting and mapping mangrove biomass from canopy grain analysis using Fourier-based textural ordination of IKONOS images. Remote Sens Environ 109: 379–392.C. ProisyP. CouteronF. Fromard2007Predicting and mapping mangrove biomass from canopy grain analysis using Fourier-based textural ordination of IKONOS images.Remote Sens Environ109379392
34. Helmer E, Cohen W, Brown S (2000) Mapping montane tropical forest successional stage and land use with multi-date Landsat imagery. Int J Remote Sens 21: 2163–2183.E. HelmerW. CohenS. Brown2000Mapping montane tropical forest successional stage and land use with multi-date Landsat imagery.Int J Remote Sens2121632183
35. Helmer E, Ramos O, López T, Quiñones M, Diaz W (2002) Mapping forest type and land cover of Puerto Rico, a component of the Caribbean biodiversity hotspot. Caribb J Sci 38: 165–183.E. HelmerO. RamosT. LópezM. QuiñonesW. Diaz2002Mapping forest type and land cover of Puerto Rico, a component of the Caribbean biodiversity hotspot.Caribb J Sci38165183
36. Sánchez-Azofeifa G, Castro-Esau K, Joyce A (2009) Monitoring carbon stocks in the tropics and the remote sensing operational limitations: from local to regional projects. Ecol Appl 19: 480–494.G. Sánchez-AzofeifaK. Castro-EsauA. Joyce2009Monitoring carbon stocks in the tropics and the remote sensing operational limitations: from local to regional projects.Ecol Appl19480494
37. Vieira ICG, de Almeida AS, Davidson EA, Stone TA, Reis de Carvalho CJ, et al. (2003) Classifying successional forest using Landsat spectral properties and ecological characteristics in eastern Amazônia. Remote Sens Environ 87: 470–481.ICG VieiraAS de AlmeidaEA DavidsonTA StoneCJ Reis de Carvalho2003Classifying successional forest using Landsat spectral properties and ecological characteristics in eastern Amazônia.Remote Sens Environ87470481
38. Boyd DS, Foody GM, Curran PJ, Lucas RM, Honzak M (1996) An assessment of radiance in Landsat TM middle and thermal infrared wavebands for the detection of tropical forest regeneration. Int J Remote Sens 17: 249–261.DS BoydGM FoodyPJ CurranRM LucasM. Honzak1996An assessment of radiance in Landsat TM middle and thermal infrared wavebands for the detection of tropical forest regeneration.Int J Remote Sens17249261
39. Foody GM, Cutler ME, McMorow J, Dieter P, Tangki H, et al. (2001) Mapping the biomass of Bornean tropical rain forest from remotely sensed data. Global Ecol Biogeogr 10: 379–387.GM FoodyME CutlerJ. McMorowP. DieterH. Tangki2001Mapping the biomass of Bornean tropical rain forest from remotely sensed data.Global Ecol Biogeogr10379387
40. Lucas RM, Honzak M, Foody GM, Curran PJ, Corves C (1993) Characterizing tropical secondary forests using multi-temporal Landsat sensor imagery. Int J Remote Sens 14: 3061–3067.RM LucasM. HonzakGM FoodyPJ CurranC. Corves1993Characterizing tropical secondary forests using multi-temporal Landsat sensor imagery.Int J Remote Sens1430613067
41. Nagendra H (2001) Using remote sensing to assess biodiversity. Int J Remote Sens 22: 2377–2400.H. Nagendra2001Using remote sensing to assess biodiversity.Int J Remote Sens2223772400
42. Sader SA, Waide RB, Laurence WT, Joyce AT (1989) Tropical forest biomass and successional age class relationships to a vegetation index derived from Landsat TM data. Remote Sens Environ 28: 143–156.SA SaderRB WaideWT LaurenceAT Joyce1989Tropical forest biomass and successional age class relationships to a vegetation index derived from Landsat TM data.Remote Sens Environ28143156
43. Sohn Y, Moran E, Gurri F (1999) Deforestation in North-Central Yucatan (1985–1995): mapping secondary succession of forest and agricultural land use in Sotuta using the cosine of the angle concept. Photogramm Eng Rem S 65: 947–958.Y. SohnE. MoranF. Gurri1999Deforestation in North-Central Yucatan (1985–1995): mapping secondary succession of forest and agricultural land use in Sotuta using the cosine of the angle concept.Photogramm Eng Rem S65947958
44. Steininger MK (1996) Tropical secondary forest-regrowth in the Amazon: age, area and change estimation with Thematic Mapper data. Int J Remote Sens 17: 9–27.MK Steininger1996Tropical secondary forest-regrowth in the Amazon: age, area and change estimation with Thematic Mapper data.Int J Remote Sens17927
45. Steininger MK (2000) Satellite estimation of tropical secondary forest above-ground biomass: data from Brazil and Bolivia. Int J Remote Sens 21: 1139–1157.MK Steininger2000Satellite estimation of tropical secondary forest above-ground biomass: data from Brazil and Bolivia.Int J Remote Sens2111391157
46. Lu D (2005) Aboveground biomass estimation using Landsat TM data in the Brazilian Amazon. Int J Remote Sens 26: 2509–2525.D. Lu2005Aboveground biomass estimation using Landsat TM data in the Brazilian Amazon.Int J Remote Sens2625092525
47. Schowengerdt R (2007) Remote sensing: Models and methods for image processing. Oxford: Elsevier. 515 p.R. Schowengerdt2007Remote sensing: Models and methods for image processingOxfordElsevier515
48. Brondizio E, Moran E, Mausel P, Wu Y (1996) Land cover in the Amazon estuary: linking of the Thematic Mapper with botanical and historical data. Photogramm Eng Rem S 62: 921–929.E. BrondizioE. MoranP. MauselY. Wu1996Land cover in the Amazon estuary: linking of the Thematic Mapper with botanical and historical data.Photogramm Eng Rem S62921929
49. Kalacska M, Sanchez-Azoifefa G, Calvo-Alvarado J, Quesada B, Janzen D (2004) Species composition, similarity and diversity in three successional stages of a seasonally dry tropical forest. For Ecol Manage 200: 227–247.M. KalacskaG. Sanchez-AzoifefaJ. Calvo-AlvaradoB. QuesadaD. Janzen2004Species composition, similarity and diversity in three successional stages of a seasonally dry tropical forest.For Ecol Manage200227247
50. Kimes D, Nelson R, Salas W, Skole D (1999) Mapping secondary tropical forest and forest age from SPOT HRV data. Int J Remote Sens 20: 3625–3640.D. KimesR. NelsonW. SalasD. Skole1999Mapping secondary tropical forest and forest age from SPOT HRV data.Int J Remote Sens2036253640
51. Mausel P, Wu Y, Li Y, Moran EF, Brondizio ES (1993) Spectral identification of successional stages following deforestation in the Amazon. Geocarto Int 4: 61–71.P. MauselY. WuY. LiEF MoranES Brondizio1993Spectral identification of successional stages following deforestation in the Amazon.Geocarto Int46171
52. Rignot E, Salas WA, Skole DL (1997) Mapping deforestation and secondary growth in Rondonia, Brazil using imaging radar and Thematic Mapper data. Remote Sens Environ 59: 167–179.E. RignotWA SalasDL Skole1997Mapping deforestation and secondary growth in Rondonia, Brazil using imaging radar and Thematic Mapper data.Remote Sens Environ59167179
53. Rocchini D, Balkenhol N, Carter G, Foody G, Gillespie T, et al. (2010) Remotely sensed spectral heterogeneity as a proxy of species diversity: recent advances and open challenges. Ecol Inform 5: 318–329.D. RocchiniN. BalkenholG. CarterG. FoodyT. Gillespie2010Remotely sensed spectral heterogeneity as a proxy of species diversity: recent advances and open challenges.Ecol Inform5318329
54. Rocchini D, Chiarucci A, Loiselle SA (2004) Testing the spectral variation hypothesis by using satellite multispectral images. Acta Oecol 26: 117–120.D. RocchiniA. ChiarucciSA Loiselle2004Testing the spectral variation hypothesis by using satellite multispectral images.Acta Oecol26117120
55. Rocchini D, Vannini A (2010) What is up? Testing spectral heterogeneity vs. NDVI relationship by quantile regression. Int J Remote Sens 31: 2745–2756.D. RocchiniA. Vannini2010What is up? Testing spectral heterogeneity vs. NDVI relationship by quantile regression.Int J Remote Sens3127452756
56. Barbier N, Couteron P, Proisy C, Malhi Y, Gastellu-Etchegorry J-P (2010) The variation of apparent crown size and canopy heterogeneity across lowland Amazonian forests. Global Ecol Biogeogr 19: 72–84.N. BarbierP. CouteronC. ProisyY. MalhiJ-P Gastellu-Etchegorry2010The variation of apparent crown size and canopy heterogeneity across lowland Amazonian forests.Global Ecol Biogeogr197284
57. Couteron P, Pelissier R, Nicolini EA, Paget D (2005) Predicting tropical forest stand structure parameters from Fourier transform of very high-resolution remotely sensed canopy images. J Appl Ecol 42: 1121–1128.P. CouteronR. PelissierEA NicoliniD. Paget2005Predicting tropical forest stand structure parameters from Fourier transform of very high-resolution remotely sensed canopy images.J Appl Ecol4211211128
58. Malhi Y, Román-Cuesta RM (2008) Analysis of lacunarity and scales of spatial homogeneity in IKONOS images of Amazonian tropical forest canopies. Remote Sens Environ 112: 2074–2087.Y. MalhiRM Román-Cuesta2008Analysis of lacunarity and scales of spatial homogeneity in IKONOS images of Amazonian tropical forest canopies.Remote Sens Environ11220742087
59. Haralick RM (1979) Statistical and structural approaches to texture. P IEEE 67: 786–804.RM Haralick1979Statistical and structural approaches to texture.P IEEE67786804
60. Fuchs H, Magdon P, Kleinn C, Flessa H (2009) Estimating aboveground carbon in a catchment of the Siberian forest tundra: combining satellite imagery and field inventory. Remote Sens Environ 113: 518–531.H. FuchsP. MagdonC. KleinnH. Flessa2009Estimating aboveground carbon in a catchment of the Siberian forest tundra: combining satellite imagery and field inventory.Remote Sens Environ113518531
61. Lu D (2006) The potential and challenge of remote sensing-based biomass estimation. Int J Remote Sens 27: 1297–1328.D. Lu2006The potential and challenge of remote sensing-based biomass estimation.Int J Remote Sens2712971328
62. Sarker LR, Nichol JE (2011) Improved forest biomass estimates using ALOS AVNIR-2 texture indices. Remote Sens Environ 115: 968–977.LR SarkerJE Nichol2011Improved forest biomass estimates using ALOS AVNIR-2 texture indices.Remote Sens Environ115968977
63. Boyd DS, Danson FM (2005) Satellite remote sensing of forest resources: three decades of research development. Prog Phys Geog 29: 1–26.DS BoydFM Danson2005Satellite remote sensing of forest resources: three decades of research development.Prog Phys Geog29126
64. Hay GJ, Niemann KO, McLeab GF (1996) An object-specific image texture analysis of H-resolution forest imagery. Remote Sens Environ 55: 108–122.GJ HayKO NiemannGF McLeab1996An object-specific image texture analysis of H-resolution forest imagery.Remote Sens Environ55108122
65. Kerr JT, Ostrovsky M (2003) From space to species: ecological applications for remote sensing. Trends Ecol Evol 18: 299–305.JT KerrM. Ostrovsky2003From space to species: ecological applications for remote sensing.Trends Ecol Evol18299305
66. Laliberte AS, Rango A, Havstad KM, Paris JF, Beck RF, et al. (2004) Object-oriented image analysis for mapping shrub encroachment from 1937 to 2003 in southern New Mexico. Remote Sens Environ 93: 198–210.AS LaliberteA. RangoKM HavstadJF ParisRF Beck2004Object-oriented image analysis for mapping shrub encroachment from 1937 to 2003 in southern New Mexico.Remote Sens Environ93198210
67. Gallardo-Cruz J, Pérez-García E, Meave J (2009) β-diversity and vegetation structure as influenced by slope aspect and altitude in a seasonally dry tropical landscape. Landscape Ecol 24: 473–482.J. Gallardo-CruzE. Pérez-GarcíaJ. Meave2009β-diversity and vegetation structure as influenced by slope aspect and altitude in a seasonally dry tropical landscape.Landscape Ecol24473482
68. Pérez-García E, Meave J, Villaseñor J, Gallardo-Cruz J, Lebrija-Trejos E (2010) Vegetation heterogeneity and life-strategy diversity in the flora of the heterogeneous landscape of Nizanda, Oaxaca, Mexico. Folia Geobot 45: 143–161.E. Pérez-GarcíaJ. MeaveJ. VillaseñorJ. Gallardo-CruzE. Lebrija-Trejos2010Vegetation heterogeneity and life-strategy diversity in the flora of the heterogeneous landscape of Nizanda, Oaxaca, Mexico.Folia Geobot45143161
69. Lebrija-Trejos E, Bongers F, Pérez-García E, Meave J (2008) Successional change and resilience of a very dry tropical deciduous forest following shifting agriculture. Biotropica 40: 422–431.E. Lebrija-TrejosF. BongersE. Pérez-GarcíaJ. Meave2008Successional change and resilience of a very dry tropical deciduous forest following shifting agriculture.Biotropica40422431
70. Brienen R, Lebrija-Trejos E, van Breugel M, Bongers F, Meave J, et al. (2009) The potential of tree rings for the study of forest succession in southern Mexico. Biotropica 4: 186–195.R. BrienenE. Lebrija-TrejosM. van BreugelF. BongersJ. Meave2009The potential of tree rings for the study of forest succession in southern Mexico.Biotropica4186195
71. Lebrija-Trejos E, Pérez-García E, Meave J, Bongers F, Poorter L (2010) Functional traits and environmental filtering drive community assembly in a species-rich tropical system. Ecology 91: 386–389.E. Lebrija-TrejosE. Pérez-GarcíaJ. MeaveF. BongersL. Poorter2010Functional traits and environmental filtering drive community assembly in a species-rich tropical system.Ecology91386389
72. Lebrija-Trejos E, Pérez-García EA, Meave JA, Poorter L, Bongers F (In Press) Environmental changes during secondary succession in a tropical dry forest in southern Mexico. J Trop Ecol. E. Lebrija-TrejosEA Pérez-GarcíaJA MeaveL. PoorterF. BongersIn PressEnvironmental changes during secondary succession in a tropical dry forest in southern Mexico.J Trop Ecol
73. Kalacska M, Sanchez-Azoifefa G, Rivard B, Caelli T, White H, et al. (2007) Ecological fingerprinting of ecosystem succession: estimating secondary tropical dry forest structure and diversity using imaging spectroscopy 108: 82–96. Remote Sens Environ 108: 82–96.M. KalacskaG. Sanchez-AzoifefaB. RivardT. CaelliH. White2007Ecological fingerprinting of ecosystem succession: estimating secondary tropical dry forest structure and diversity using imaging spectroscopy 108: 82–96.Remote Sens Environ1088296
74. Magurran A (2004) Measuring biological diversity. London: Blackwell-Science. 256 p.A. Magurran2004Measuring biological diversityLondonBlackwell-Science256
75. Krause K (2005) Radiometric use of Quickbird imagery. K. Krause2005Radiometric use of Quickbird imagery.Technical Note (DigitalGlobe Inc). Technical Note (DigitalGlobe Inc).
76. Gao X, Huete AR, Ni WG, Miura T (2000) Optical–biophysical relationships of vegetation spectra without background contamination. Remote Sens Environ 74: 609–620.X. GaoAR HueteWG NiT. Miura2000Optical–biophysical relationships of vegetation spectra without background contamination.Remote Sens Environ74609620
77. Petrou M, García-Sevilla P (2006) Image processing: dealing with texture. Sussex: John Wiley & Sons. 618 p.M. PetrouP. García-Sevilla2006Image processing: dealing with textureSussexJohn Wiley & Sons618
78. Marquez J (2008) Texture Characterization and Analysis Tutorial I - version 3.0. Mexico City: CCADET-UNAM. 90 p.J. Marquez2008Texture Characterization and Analysis Tutorial I - version 3.0Mexico CityCCADET-UNAM90
79. ITT (2008) ENVI 4.5/IDL Version 7.0 User's Guide. Boulder: ITT Visual Information Solutions. ITT2008ENVI 4.5/IDL Version 7.0 User's GuideBoulderITT Visual Information Solutions
80. Burnham KP, Anderson DR (2002) Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer. 515 p.KP BurnhamDR Anderson2002Model selection and multimodel inference: a practical information-theoretic approachNew YorkSpringer515
81. Baumann K (2003) Cross-validation as the objective function for variable-selection techniques. Trac-Trend Anal Chem 22: 395–406.K. Baumann2003Cross-validation as the objective function for variable-selection techniques.Trac-Trend Anal Chem22395406
82. Browne M (2000) Cross-validation methods. J Math Psychol 44: 108–132.M. Browne2000Cross-validation methods.J Math Psychol44108132
83. Shao J (1993) Linear model selection by cross-validation. J Am Stat Assoc 88: 486–494.J. Shao1993Linear model selection by cross-validation.J Am Stat Assoc88486494
84. R Development Core Team (2010) R: A language and environment for statistical computing (v.2.11.1). Vienna: R Foundation for Statistical Computing. R Development Core Team2010R: A language and environment for statistical computing (v.2.11.1)ViennaR Foundation for Statistical Computing
85. Kuplich TM, Curran PJ, Atkinson PM (2005) Relating SAR image texture to the biomass of regenerating tropical forests. Int J Remote Sens 26: 4829–4854.TM KuplichPJ CurranPM Atkinson2005Relating SAR image texture to the biomass of regenerating tropical forests.Int J Remote Sens2648294854
86. Murray H, Lucieer A, Williams R (2010) Texture-based classification of sub-Antarctic vegetation communities on Heard Island. Int J Appl Earth Obs 12: 138–149.H. MurrayA. LucieerR. Williams2010Texture-based classification of sub-Antarctic vegetation communities on Heard Island.Int J Appl Earth Obs12138149
87. Wijaya A, Liesenberg V, Gloaguen R (2010) Retrieval of forest attributes in complex successional forests of Central Indonesia: Modeling and estimation of bitemporal data. For Ecol Manage 259: 2315–2326.A. WijayaV. LiesenbergR. Gloaguen2010Retrieval of forest attributes in complex successional forests of Central Indonesia: Modeling and estimation of bitemporal data.For Ecol Manage25923152326
88. Helmer E, Lefsky M, Roberts D (2009) Biomass accumulation rates of Amazonian secondary forest and biomass of old-growth forests from Landsat time series and the Geoscience Laser Altimeter System. J Appl Remote Sens 3: 033505.E. HelmerM. LefskyD. Roberts2009Biomass accumulation rates of Amazonian secondary forest and biomass of old-growth forests from Landsat time series and the Geoscience Laser Altimeter System.J Appl Remote Sens3033505
89. Liu W, Song C, Schroeder T, Cohen W (2008) Predicting forest successional stages using multitemporal Landsat imagery with forest inventory and analysis data. Int J Remote Sens 29: 3855–3872.W. LiuC. SongT. SchroederW. Cohen2008Predicting forest successional stages using multitemporal Landsat imagery with forest inventory and analysis data.Int J Remote Sens2938553872
90. Song C, Schroeder TA, Cohen WB (2007) Predicting temperate conifer forest successional stage distributions with multitemporal Landsat Thematic Mapper imagery. Remote Sens Environ 106: 228–237.C. SongTA SchroederWB Cohen2007Predicting temperate conifer forest successional stage distributions with multitemporal Landsat Thematic Mapper imagery.Remote Sens Environ106228237
91. Gillespie T, Zutta B, Early M, Saatchi S (2005) Predicting and quantifying the structure of tropical dry forests in South Florida and the Neotropics using spaceborne imagery. Global Ecol Biogeogr 15: 225–236.T. GillespieB. ZuttaM. EarlyS. Saatchi2005Predicting and quantifying the structure of tropical dry forests in South Florida and the Neotropics using spaceborne imagery.Global Ecol Biogeogr15225236
92. Barbier N, Proisy C, Vega C, Sabatier D, Couteron P (2011) Bidirectional texture function of high resolution optical images of tropical forest: an approach using LiDAR hillshade simulations. Remote Sens Environ 115: 167–179.N. BarbierC. ProisyC. VegaD. SabatierP. Couteron2011Bidirectional texture function of high resolution optical images of tropical forest: an approach using LiDAR hillshade simulations.Remote Sens Environ115167179
93. Gould W (2006) Remote sensing of vegetation, plant species richness, and regional biodiversity hot spots. Ecol Appl 10: 1861–1870.W. Gould2006Remote sensing of vegetation, plant species richness, and regional biodiversity hot spots.Ecol Appl1018611870
94. Foody GM, Cutler MEJ (2006) Mapping the species richness and composition of tropical forests from remotely sensed data with neural networks. Ecol Model 195: 37–42.GM FoodyMEJ Cutler2006Mapping the species richness and composition of tropical forests from remotely sensed data with neural networks.Ecol Model1953742
95. Levin N, Shmida A, Levanoni O, Tamari H, Kark S (2007) Predicting mountain plant richness and rarity from space using satellite-derived vegetation indices. Divers Distrib 13: 692–703.N. LevinA. ShmidaO. LevanoniH. TamariS. Kark2007Predicting mountain plant richness and rarity from space using satellite-derived vegetation indices.Divers Distrib13692703
96. Palmer MW, Earls P, Hoagland BW, White PS, Wohlgemuth T (2002) Quantitative tools for perfecting species lists. Environmetrics 13: 121–137.MW PalmerP. EarlsBW HoaglandPS WhiteT. Wohlgemuth2002Quantitative tools for perfecting species lists.Environmetrics13121137
97. Frazer GW, Wulder MA, Niemann KO (2005) Simulation and quantification of the fine-scale spatial pattern and heterogeneity of forest canopy structure: A lacunarity-based method designed for analysis of continuous canopy heights. For Ecol Manage 214: 65–90.GW FrazerMA WulderKO Niemann2005Simulation and quantification of the fine-scale spatial pattern and heterogeneity of forest canopy structure: A lacunarity-based method designed for analysis of continuous canopy heights.For Ecol Manage2146590
98. Díaz S, Fargione J, Chapin FS III, Tilman D (2006) Biodiversity loss threatens human well-being. PLoS Biol 4: E277.S. DíazJ. FargioneFS Chapin IIID. Tilman2006Biodiversity loss threatens human well-being.PLoS Biol4E277
99. Dirzo R, Raven PH (2003) Global state of biodiversity and loss. Annu Rev Env Resour 28: 137–167.R. DirzoPH Raven2003Global state of biodiversity and loss.Annu Rev Env Resour28137167
100. Thomas CD, Cameron A, Green RE, Bakkenes M, Beaumont LJ, et al. (2004) Extinction risk from climate change. Nature 427: 145–148.CD ThomasA. CameronRE GreenM. BakkenesLJ Beaumont2004Extinction risk from climate change.Nature427145148
101. Walther G-R, Post E, Convey P, Menzel A, Parmesan C, et al. (2002) Ecological responses to recent climate change. Nature 416: 389–395.G-R WaltherE. PostP. ConveyA. MenzelC. Parmesan2002Ecological responses to recent climate change.Nature416389395
102. Slik JWF, Aiba SI, Brearley FQ, Cannon CH, Forshed O, et al. (2010) Environmental correlates of tree biomass, basal area, wood specific gravity and stem density gradients in Borneo's tropical forests. Global Ecol Biogeogr 19: 50–60.JWF SlikSI AibaFQ BrearleyCH CannonO. Forshed2010Environmental correlates of tree biomass, basal area, wood specific gravity and stem density gradients in Borneo's tropical forests.Global Ecol Biogeogr195060
103. Brown S (1997) Estimating biomass and biomass change of tropical forests: a primer. Rome: Food and Agriculture Organization of the United Nations. 57 p.S. Brown1997Estimating biomass and biomass change of tropical forests: a primerRomeFood and Agriculture Organization of the United Nations57
104. Barlow J, Gardner TA, Araujo IS, Ávila-Pirez TC, Bonaldo AB, et al. (2007) Quantifying the biodiversity value of tropical primary, secondary, and plantation forests. P Natl Acad Sci USA 104: 18555–18560.J. BarlowTA GardnerIS AraujoTC Ávila-PirezAB Bonaldo2007Quantifying the biodiversity value of tropical primary, secondary, and plantation forests.P Natl Acad Sci USA1041855518560
105. Brodie J, Post E, Laurance W (2010) How to conserve the tropics as they warm. Nature 468: 634.J. BrodieE. PostW. Laurance2010How to conserve the tropics as they warm.Nature468634
106. Putz FE, Blate GM, Redford KH, Fimber R, Robinson J (2001) Tropical forest management and conservation of biodiversity: an overview. Conserv Biol 15: 7–20.FE PutzGM BlateKH RedfordR. FimberJ. Robinson2001Tropical forest management and conservation of biodiversity: an overview.Conserv Biol15720
107. Daily GC, Ehrlich PR, Sanchez-Azofeifa GA (2001) Countryside biogeography: use of human-dominated habitats by the avifauna of southern Costa Rica. Ecol Appl 11: 1–13.GC DailyPR EhrlichGA Sanchez-Azofeifa2001Countryside biogeography: use of human-dominated habitats by the avifauna of southern Costa Rica.Ecol Appl11113
108. Mayfield MM, Daily GC (2005) Countryside biogeography of neotropical herbaceous and shrubby plants. Ecol Appl 15: 423–439.MM MayfieldGC Daily2005Countryside biogeography of neotropical herbaceous and shrubby plants.Ecol Appl15423439
109. Perfecto I, Vandermeer J (2008) Biodiversity conservation in tropical agroecosystems: a new conservation paradigm. Ann N Y Acad Sci 1134: 173–200.I. PerfectoJ. Vandermeer2008Biodiversity conservation in tropical agroecosystems: a new conservation paradigm.Ann N Y Acad Sci1134173200
110. Hernández-Stefanoni JL, Gallardo-Cruz JA, Meave JA, Dupuy JM (2011) Combining geostatistical models and remotely sensed data to improve tropical tree richness mapping. Ecol Indicators. JL Hernández-StefanoniJA Gallardo-CruzJA MeaveJM Dupuy2011Combining geostatistical models and remotely sensed data to improve tropical tree richness mapping.Ecol Indicators
111. Turner W, Spector S, Gardiner N, Fladeland M, Sterling E, et al. (2003) Remote sensing for biodiversity science and conservation. Trends Ecol Evol 18: 306–314.W. TurnerS. SpectorN. GardinerM. FladelandE. Sterling2003Remote sensing for biodiversity science and conservation.Trends Ecol Evol18306314
112. Brandt J (1988) The transformation of rainfall energy by a tropical rain forest canopy in relation to soil erosion. J Biogeogr 15: 41–48.J. Brandt1988The transformation of rainfall energy by a tropical rain forest canopy in relation to soil erosion.J Biogeogr154148
113. Marques MJ, Bienes R, Jiménez L, Pérez-Rodríguez R (2007) Effect of vegetal cover on runoff and soil erosion under light intensity events. Rainfall simulation over USLE plot. Sci Total Environ 378: 161–165.MJ MarquesR. BienesL. JiménezR. Pérez-Rodríguez2007Effect of vegetal cover on runoff and soil erosion under light intensity events. Rainfall simulation over USLE plot.Sci Total Environ378161165
114. van Dijk AIJM, Bruijnzeel LA, Rosewell CJ (2002) Rainfall intensity-kinetic energy relationships: a critical literatura appraisal. J Hydrol 261: 1–23.AIJM van DijkLA BruijnzeelCJ Rosewell2002Rainfall intensity-kinetic energy relationships: a critical literatura appraisal.J Hydrol261123
115. Bélisle M, Desrochers A, Fortin M-J (2001) Influence of forest cover on the movements of forest birds: a homing experiment. Ecology 82: 1893–1904.M. BélisleA. DesrochersM-J Fortin2001Influence of forest cover on the movements of forest birds: a homing experiment.Ecology8218931904
116. Harvey CA, Medina A, Sánchez DM, Vílchez S, Hernández B, et al. (2006) Patterns of animal diversity in different forms of tree cover in agricultural landscapes. Ecol Appl 16: 1986–1999.CA HarveyA. MedinaDM SánchezS. VílchezB. Hernández2006Patterns of animal diversity in different forms of tree cover in agricultural landscapes.Ecol Appl1619861999
117. Trzcinski MK, Fahrig L, Merriam G (1999) Independent effects of forest cover and fragmentation on the distribution of forest breeding birds. Ecol Appl 9: 586–593.MK TrzcinskiL. FahrigG. Merriam1999Independent effects of forest cover and fragmentation on the distribution of forest breeding birds.Ecol Appl9586593
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2012 Gallardo et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License: https://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Biodiversity conservation and ecosystem-service provision will increasingly depend on the existence of secondary vegetation. Our success in achieving these goals will be determined by our ability to accurately estimate the structure and diversity of such communities at broad geographic scales. We examined whether the texture (the spatial variation of the image elements) of very high-resolution satellite imagery can be used for this purpose. In 14 fallows of different ages and one mature forest stand in a seasonally dry tropical forest landscape, we estimated basal area, canopy cover, stem density, species richness, Shannon index, Simpson index, and canopy height. The first six attributes were also estimated for a subset comprising the tallest plants. We calculated 40 texture variables based on the red and the near infrared bands, and EVI and NDVI, and selected the best-fit linear models describing each vegetation attribute based on them. Basal area (R2 = 0.93), vegetation height and cover (0.89), species richness (0.87), and stand age (0.85) were the best-described attributes by two-variable models. Cross validation showed that these models had a high predictive power, and most estimated vegetation attributes were highly accurate. The success of this simple method (a single image was used and the models were linear and included very few variables) rests on the principle that image texture reflects the internal heterogeneity of successional vegetation at the proper scale. The vegetation attributes best predicted by texture are relevant in the face of two of the gravest threats to biosphere integrity: climate change and biodiversity loss. By providing reliable basal area and fallow-age estimates, image-texture analysis allows for the assessment of carbon sequestration and diversity loss rates. New and exciting research avenues open by simplifying the analysis of the extent and complexity of successional vegetation through the spatial variation of its spectral information.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer