Mapping Land Cover and Estimating the Grassland Structure in a Priority Area of the Chihuahuan Desert

A field characterization of the grassland vegetation structure, represented by the coverage of grass canopy (CGC) and the grass height, was carried out during three years (2009–2011) in a priority area for the conservation of grasslands of North America. Landsat Thematic Mapper (TM5) images were selected and the information of reflectance was obtained based on the geographical location of each field-sampling site. Linear models, constructed with field and satellite data, with high coefficients of determination for CGC (R2 = 0.81, R2 = 0.81 and R2 = 0.72) and grass height (R2 = 0.82, R2 = 0.79 and R2 = 0.73) were obtained. The maps showed a good level of CGC (>25%) and grass height (>25 cm), except for the year 2009, which presented the lowest values of grass height in the area. According to the Kappa Index, a moderate concordance among the three CGC maps was presented (0.49–0.59). Conversely, weak and moderate concordances were found among the grass height maps (0.36–0.59). It was observed that areas with a high CGC do not necessarily correspond to areas with greater grass height values. Based on the data analyzed in this study, the grassland areas are highly dynamic, structurally heterogeneous and the spatial distribution of the variables does not show a definite pattern. From the information generated, it is possible to determine those areas that are the most important for monitoring to then establish effective strategies for the conservation of these grasslands and the protection of threatened migratory bird species.


Introduction
Given their great biodiversity and the environmental services they provide; grasslands are ecosystems of great importance.According to Adams et al. [1] grasslands account for 40.5% of the Earth's surface and therefore play an important role in the global carbon cycle [2].Currently, these ecosystems are undergoing significant deterioration due to overgrazing, urbanization, land use change, the presence of invasive species and habitat fragmentation [3][4][5].Indeed, Hoesktra et al. [6] considered that grasslands are the most threatened ecosystems on Earth.Consequently, many species of migratory birds that use grasslands as habitat in the North American deserts have declined their populations continuously and steadily [7,8].Previous studies on grasslands have focused on herbaceous productivity trends and on the factors that may affect the herbaceous coverage [9,10].For instance, from studies on grasslands of the Chihuahuan Desert, which is distributed in northern Mexico and the southern United States of America, some studies have focused on listing the plant species, on studies of its fauna and on its ecosystem fragmentation [11,12].However, more research is needed on the spatio-temporal variation of the coverage of grass canopy (CGC) and grass height.This topic is particularly relevant because these variables are strongly related to the habitat quality for migratory birds and because the grasslands immersed in the Chihuahuan Desert do not escape from the issues and pressures that characterize this ecosystem [13,14].
One of the tools widely employed for monitoring habitat quality is remote sensing.This tool offers the possibility of studying extensive areas with a wide temporal margin at a low cost [15,16].Data from the thematic mapper (TM) sensor (Landsat 5 satellite) provides the capability to perform synoptic monitoring.The Landsat archive provides detailed and consistent data about the change dynamics experienced in the terrestrial ecosystems [17][18][19][20].The spectral data provided by the sensor have allowed the determination of the structural characteristics and composition of plants in various ecosystems [21][22][23].Some of the parameters that have been correlated to habitat quality are CGC and grass height, parameters that have been employed to estimate the net annual biomass productivity of grasslands in the field.However, these kinds of methodologies require a significant investment of time and economic resources [24].Thus, the use of remote sensing tools represents an alternative for the determination of grassland variables in a more economical and faster manner [25].
Several studies have shown the condition of winter habitat as a determining factor in the decline of some grassland avian species.In fact, the abundance and distribution of the native species of the Chihuahuan Desert, such as Ammodramus bairdii, A. savannarum, Anthus spragueii, Calcarius ornatus and Pooecetes gramineus, are strongly related to the structure of the grasslands, as well as to variables such as CGC and grass height [14,26].Therefore, it is very important to obtain information about the dynamics of the grassland where the structural variables are highly susceptible for drastic changes in relatively short periods [27].The objective of this study was to determine the spatial and temporal variations of CGC and grass height during a period of three years in an area under conservation called Cuchillas de la Zarca, located within the Chihuahuan Desert in Mexico.

Materials and Methods
The area of Cuchillas de la Zarca (CUZA) is a priority region for the conservation of grasslands and is located within the Chihuahuan Desert.This desert grassland is the most extensive in North America and is recognized for having high species richness despite being a dry area.In CUZA we can find both transition areas dominated by grasslands and extensive plains covered by them [28].The area of CUZA is located between the coordinates 105 • 04 30.86 W and 24 • 41 41.28 N, as well as 104 • 57 07.48 W and 25 • 27 04.24N. It comprises the northern part of the state of Durango and the southern part of the state of Chihuahua, Mexico.It has an area of 11,600 km 2 (Figure 1).The vegetation consists mainly of induced pasture, natural grassland, microphyll desert scrub, chaparral areas and pine-oak forest.In each of these vegetation types, grass species of the genera Bouteloa, Aristida, Buchloe, Andropogon, Melinis, Muhlenbergia, Sporobolus, Heteropogon and Pleuraphis may be present [29,30].

Samples and Processing Date
The grassland structure, represented by CGC and grass height, was characterized during January and February by 96 transects in 2009, as well as 102 transects in 2010 and 2011.The transects of one kilometer of length were located according to a grid of squares with side lengths of 18 km covering the study area (Figure 2).The grid design was based on information from vector files of land use and roads network [31].From all the 17 quadrants, only those with roads of acceptable condition were selected for easy access to the sampling points (red marks, Figure 1), which were located at every 500 m in each quadrant, as was established in a previous study [14].For vegetation sampling, only the first three points of each quadrant were selected.In them, two sampling transects, perpendicular to the path, were established.Along the transects, circles of 5 m of radius were located each 100 m.In these circles, the area covered by grass canopy was visually assessed and the heights of three plants were measured at the ground level.For the height measurements, a small, a medium and a tall plant were chosen.From the three values, an average was estimated.All the measurements were made by only one person.For a summary of field data, the average CGC and grass height were obtained.

Samples and Processing Date
The grassland structure, represented by CGC and grass height, was characterized during January and February by 96 transects in 2009, as well as 102 transects in 2010 and 2011.The transects of one kilometer of length were located according to a grid of squares with side lengths of 18 km covering the study area (Figure 2).The grid design was based on information from vector files of land use and roads network [31].From all the 17 quadrants, only those with roads of acceptable condition were selected for easy access to the sampling points (red marks, Figure 1), which were located at every 500 m in each quadrant, as was established in a previous study [14].For vegetation sampling, only the first three points of each quadrant were selected.In them, two sampling transects, perpendicular to the path, were established.Along the transects, circles of 5 m of radius were located each 100 m.In these circles, the area covered by grass canopy was visually assessed and the heights of three plants were measured at the ground level.For the height measurements, a small, a medium and a tall plant were chosen.From the three values, an average was estimated.All the measurements were made by only one person.For a summary of field data, the average CGC and grass height were obtained.

Samples and Processing Date
The grassland structure, represented by CGC and grass height, was characterized during January and February by 96 transects in 2009, as well as 102 transects in 2010 and 2011.The transects of one kilometer of length were located according to a grid of squares with side lengths of 18 km covering the study area (Figure 2).The grid design was based on information from vector files of land use and roads network [31].From all the 17 quadrants, only those with roads of acceptable condition were selected for easy access to the sampling points (red marks, Figure 1), which were located at every 500 m in each quadrant, as was established in a previous study [14].For vegetation sampling, only the first three points of each quadrant were selected.In them, two sampling transects, perpendicular to the path, were established.Along the transects, circles of 5 m of radius were located each 100 m.In these circles, the area covered by grass canopy was visually assessed and the heights of three plants were measured at the ground level.For the height measurements, a small, a medium and a tall plant were chosen.From the three values, an average was estimated.All the measurements were made by only one person.For a summary of field data, the average CGC and grass height were obtained.Two scenes (path/row 31/41 and 31/42) of Landsat TM5 were acquired for each of the studied years from the United States Geological Survey (USGS) Global Visualization Viewer (GloVis, https: //www.usgs.gov).The Landsat images correspond to the date of ground measurements, 5th February in 2009, 8th February in 2010, and 11th February 2011.For each year, the scenes were merged into a single mosaic.
To establish the relationship between the variables of vegetation structure and spectral values, the spectral information of the digital number (DN), extracted from the bands 1-5 and 7 of the sensor, the radiometrically-corrected values of reflectance and three synthetic bands, generated for the tasseled cap (TC), were correlated with CGC and grass height.The obtained bands of TC corresponded to the components of brightness (B), greenness (G) and wetness (W) [32].The radiometric calibration was made by converting the gross value of each pixel (DN) into values of absolute spectral radiance.To reduce the variability among scenes, the spectral radiance was converted into exoatmospheric top of the atmosphere (TOA) reflectance.Formulae and coefficients provided by Chander et al. [33], represented by Equations ( 1) and ( 2), were used to perform the radiometric corrections applied to the TM sensor.
where QCAL is DN; QCALmin and QCALmax are the minimum and maximum quantized calibrated pixel values, respectively; Lmin λ is the spectral radiance scales to QCALmin; Lmax λ is the spectral radiance scales to QCALmax; d is the distance from the earth to the sun; ESUN λ is the mean solar exoatmospheric irradiance; and θ s is the solar zenith angle.The atmospheric correction seeks to reduce or eliminate the image distortions resulting from the interaction of the atmosphere components with the sensor [34].Such correction allows the standardization of image data and thus permits a comparison of images from different dates.From each of the circles of vegetation sampling located in the transects, the value of reflectance and TC components for each band were obtained from the corresponding pixel.Although some bushes and forests exist in the area, when sampling of vegetation in areas dominated by pasture, the sensor data were mainly associated to the spectral signatures of pastures.The bands that best explained the variables of CGC and grass height were selected to feed the models.

Statistical Analysis
The normality of the variables was proven with the Shapiro-Wilk test.After checking it, an analysis of variance was applied to determine if the differences among the quadrants were significant.Subsequently, multiple regression models [35] were run to explore the relationships of CGC and grass height with the spectral data.The CGC and grass height were established as response variables while the spectral data were used as independent variables.The spectral data included the DN, reflectance values and the values from the tasseled cap.The models were produced separately by dependent variable and by year.Once the models were obtained, the Gauss-Markov assumptions, which include normality, independence and homogeneity of variance, were verified on the data.To verify normality, the test was applied to the residuals of each model; to check independence, it was determined whether the correlation between the residual and predicted values were zero; and, to verify the homogeneity of variance, it was verified whether the correlation between the absolute value of the residuals and predicted values were close to zero (p > 0.05).

Maps of Cover of Grass Canopy and Grass Height
The estimators of the selected full models were used to construct maps of the variables CGC and grass height.For the case of CGC (%), their continuous values were reclassified in the following five classes: water bodies, 0-25 (low), 25-50 (medium), 50-75 (high) and 75-100 (very high).For the case of grass height (cm), the classes were water bodies, 0-25 (low), 25-50 (medium), 50-75 (high) and >75 (very high).To make a comparison between the resulting maps, a cross-tabulation of images was made.Such a cross-tabulation determines the frequency of pixels that stayed constant during the evaluation period and serves to make inferences about the dynamics of CGC and grass height through time in the study area.Through this operation, the categories of a raster image were compared with those of a second image; both images had to have the same number of categories and the same data type.The operation casts the Kappa Index (KI) as a measure of agreement.The KI ranges from zero, indicating that there is no agreement, and up to one, indicating perfect agreement [36,37].

Meteorological Data
Given that information about precipitation in the study area was scarce, vector files from the National Meteorological Service, providing information on drought events occurring in the country, were used.These data are based on obtaining and interpreting indicators such as the standardized precipitation index (SPI), precipitation percent anomaly, satellite vegetation health index (VHI), NOAA/CPC leaky bucket soil moisture model, the normalized difference vegetation index (NDVI), the average temperature anomaly, the water availability percentage in the country's dams and the contribution of local experts.These indicators are displayed in the form of layers in a geographic information system and, by consensus, the drought-affected regions are determined according to a scale of intensities ranging from abnormally dry, moderate drought, severe drought, drought extreme and exceptional drought.As a result of the consensus, polygons for the intensity of drought were delineated, generating shape files.When the polygons correspond to the mid-month analysis (issued on the 15th day of each month) they are used to quantify the drought on the domestic territory.Likewise, when the polygons correspond to the final evaluation of the month, they complement the regional or continental map of the North American Drought Monitor [38].Weather information linked to drought events was used to explain the different scenarios shown by the maps, considering that these events are related to precipitation and, in turn, rains occurring in a given year can largely explain the grassland conditions [39].

Field Data
The average CGC ranged between 50%-75% and there were no differences between the three years (p > 0.05).Conversely, significant differences in grass height (cm) were found among the years studied (p < 0.05).The lowest values were recorded in 2009 (Table 1).

Spectral Data
The digital numbers and reflectance values followed the same pattern over the three years evaluated; bands 5, 7 (middle-infrared) and 4 (near-infrared) registered the highest values.A decline of reflectance was observed for the quadrants located in the central and southeast parts of the study area, particularly from the bands 1-3 (blue, green and red, respectively).In 2011, the values from the bands 2-4 and 7 showed a slight difference compared to the other two years (Figure 3).

Spectral Data
The digital numbers and reflectance values followed the same pattern over the three years evaluated; bands 5, 7 (middle-infrared) and 4 (near-infrared) registered the highest values.A decline of reflectance was observed for the quadrants located in the central and southeast parts of the study area, particularly from the bands 1-3 (blue, green and red, respectively).In 2011, the values from the bands 2-4 and 7 showed a slight difference compared to the other two years (Figure 3).Forty-two models for CGC and grass height were obtained.They corresponded to three years, two dependent variables, three elements of TC with the original bands, radiometrically-corrected bands and three transformations of dependent with independent variables.Six of the models were selected, one for each variable per year.The selected models met the Gauss-Markov assumptions, and had the lowest statistical deviance information criterion, Akaike information criterion, Bayesian information criterion, (DIC, AIC, BIC) and the Root Mean Square Error (RMSE).In all the cases, the selected models only used images from bands into original digital values or in combination with an Forty-two models for CGC and grass height were obtained.They corresponded to three years, two dependent variables, three elements of TC with the original bands, radiometrically-corrected bands and three transformations of dependent with independent variables.Six of the models were selected, one for each variable per year.The selected models met the Gauss-Markov assumptions, and had the lowest statistical deviance information criterion, Akaike information criterion, Bayesian information criterion, (DIC, AIC, BIC) and the Root Mean Square Error (RMSE).In all the cases, the selected models only used images from bands into original digital values or in combination with an element of the TC.The coefficients of determination for the six models were high and with a significant probability (p < 0.01).In addition, the validation tests allowed us to verify the usefulness of each model, since the values of the mean square error of validation were slightly lower than those of calibration (Table 2).For CGC, the best selected models for 2010 and 2011 used the bands in combination with the B element of the TC.For 2009, the model only used bands (Table 3).Similarly, the grass height models of 2009 and 2011 used only bands.Conversely, the models of 2010 used both the bands and the W element of TC (Table 4).

Distribution Maps of CGC, Grass Height and Kappa Index
The selected models allowed the obtainment of CGC maps for three years (Figures 4a, 5a and 6a).According to the Kappa Index, among all the CGC maps a moderate concordance was presented, which was slightly higher between 2009 and 2010 and lower between 2010 and 2011 (Table 5).This was perhaps due to the fact that in 2009 and 2010 much of the surface of the study area presented a higher CGC (50-75%), while in 2011 the predominant interval was mainly located in low or intermediate levels (25-50%), as can be seen in Figure 4a.In 2009, the areas with higher CGC were in the zones of grassland and in the desert scrub.For 2010, the highest CGC was found in the northwest, a zone in which grassland areas and oak-pine forest are predominant.A similar pattern of distribution of the classes was seen in 2010 and 2011, although in a lower and more dispersed way for the latter.
As in the case of CGC, three maps of grass height were obtained from the selected models (Figures 4b, 5b and 6b).Based on the Kappa index only, a moderate concordance was presented between the maps of 2010 and 2011, while 2009 had a weak concordance for the two years (Table 5).According to the graph of the surface of the grass height ranges (Figure 7), the basic difference among the years with moderate concordance lay mainly in the fact that during 2009 the interval with the highest frequency was 0-25 cm, and for 2010 and 2011 it was 50-75 and 25-50 cm, respectively.Regarding 2009, it was characterized by grasses of low grass height along almost all of CUZA, only 8% of the area had a height of grass in the range of 25-50 cm, located primarily in the oak-pine areas.In 2010, the taller-height range (50-75 cm) was distributed generally almost along the entire area, except for the southeastern region, where grasslands with a lower grass height prevailed.Grasslands of greater grass height (>75) were conserved in both forest-oak and grassland areas.In 2011, the height ranges had a behavior similar to 2010, except that in the former, the range with the greatest presence was 25-50 cm in 57.2% of the area.In terms of distribution, the lowest grasses were in the oak-pine forests, whereas the middle and taller heights were distributed more or less uniformly through the area.In 2010, the taller-height range (50-75 cm) was distributed generally almost along the entire area, except for the southeastern region, where grasslands with a lower grass height prevailed.Grasslands of greater grass height (>75) were conserved in both forest-oak and grassland areas.In 2011, the height ranges had a behavior similar to 2010, except that in the former, the range with the greatest presence was 25-50 cm in 57.2% of the area.In terms of distribution, the lowest grasses were in the oak-pine forests, whereas the middle and taller heights were distributed more or less uniformly through the area.In 2010, the taller-height range (50-75 cm) was distributed generally almost along the entire area, except for the southeastern region, where grasslands with a lower grass height prevailed.Grasslands of greater grass height (>75) were conserved in both forest-oak and grassland areas.In 2011, the height ranges had a behavior similar to 2010, except that in the former, the range with the greatest presence was 25-50 cm in 57.2% of the area.In terms of distribution, the lowest grasses were in the oak-pine forests, whereas the middle and taller heights were distributed more or less uniformly through the area.According to the Kappa index, the CGC images exhibit a moderate agreement among them.In the case of the variable grass height, the highest agreement occurred between 2010 and 2011 (Table 5).In the case of the variable CGC, the greatest agreement was observed the years 2009 and 2010.The lowest agreement between the images was obtained for the variable grass height when comparing 2009 and 2010, as well as 2009 and 2011, which can be explained by the large surface area with a low grass height showed in the map of 2009.The field data showed lower mean and statistical differences for this year, as mentioned above.However, the best agreement between maps of 2010 and 2011 does not necessarily mean a spatial correspondence between the variables.Thus, an area with high CGC does not necessarily correspond to an area with a high grass height of the grassland.Moreover, although the descriptions of the variable distribution were generally completed, various degrees of structural heterogeneity can be appreciated with a greater detail in the maps.According to the Kappa index, the CGC images exhibit a moderate agreement among them.In the case of the variable grass height, the highest agreement occurred between 2010 and 2011 (Table 5).In the case of the variable CGC, the greatest agreement was observed the years 2009 and 2010.The lowest agreement between the images was obtained for the variable grass height when comparing 2009 and 2010, as well as 2009 and 2011, which can be explained by the large surface area with a low grass height showed in the map of 2009.The field data showed lower mean and statistical differences for this year, as mentioned above.However, the best agreement between maps of 2010 and 2011 does not necessarily mean a spatial correspondence between the variables.Thus, an area with high CGC does not necessarily correspond to an area with a high grass height of the grassland.Moreover, although the descriptions of the variable distribution were generally completed, various degrees of structural heterogeneity can be appreciated with a greater detail in the maps.

Discussion
In the literature, there is a diversity of papers similar to this study in which the coefficients of determination are lower or similar, especially in cases where the same sensor is used [23,25].The precision of the results may be partially due to the method of obtaining information in the field, the characteristics of the sensor (i.e., temporal and spatial resolution) and its age, the type of ecosystem, as well as the statistical analyses carried out.The latter appears to substantially affect the fit of the models between the data from the field and the sensor data [40][41][42][43].Although some studies have obtained a good fit with the use of a limited number of bands [25,42], others studies [23] have required the interaction of a greater number of bands, indices or synthetic bands to achieve a better fit, as it was the case in this study.It is worth mentioning that the models selected for one year may not be appropriate to be used for a different year, because the images from different years have different digital values [44].Possibly, a good fit on the models was not achieved with the data from the images due to an unsatisfactory correction of the data, which did not allow the use of a single model for each variable for the three years [45].
To explain the slight quantitative differences between the CGC maps, the study of Khumalo and Holechek [46] was considered.They found that the precipitation of the previous year is one of the major factors that are associated with biomass production in grasslands of the Chihuahuan Desert in southern New Mexico.Their results showed that individually, precipitation in August is the variable most highly correlated to the grassland condition.When analyzing the drought data for the three years, the presence of this phenomenon was not observed in the area during this month [38].This may explain why there were basically no differences in the Kappa index of CGC for the three years.Therefore, for these variations among the years, other aspects such as spatial variation in grazing, forest fires and the phenology of the grass species should be considered.The CGC levels recorded during the three winter seasons studied could have enabled different species of wildlife, such as migratory grassland birds, to find this important habitat requirement [47].
The high reflectance values for the middle-infrared and near-infrared bands correspond to the general pattern of spectral behavior in vegetated surfaces [42].The highest values for bands 2-4, 7 in 2011 may indicate high moisture contents in the soil as well as vegetation vigor.The latter factor may be attributed to the photosynthetic activity of perennial grasses, which grow and develop during the rains of February.However, according to the comparison of the maps, this difference does not seem to have greater repercussions with respect to the marked differences between the levels of CGC and grass height between 2010 and 2011.The decline in reflectance values of the bands 1-3 in some of the quadrants in the central and the southeastern zones, seems to be more associated with the particular characteristics of the composition and structure of the vegetation in some transects, since the resulting maps did not reflect the same pattern.
The CGC and grass height showed quantitative inter-annual fluctuations and spatio-temporal variations without an apparent specific pattern.Jin et al. [48] found similar results for biomass, which is directly related to CGC and grass height.The greatest quantitative variation was registered for the variable height, which was possibly due to the fact that this variable is more susceptible to changes

Discussion
In the literature, there is a diversity of papers similar to this study in which the coefficients of determination are lower or similar, especially in cases where the same sensor is used [23,25].The precision of the results may be partially due to the method of obtaining information in the field, the characteristics of the sensor (i.e., temporal and spatial resolution) and its age, the type of ecosystem, as well as the statistical analyses carried out.The latter appears to substantially affect the fit of the models between the data from the field and the sensor data [40][41][42][43].Although some studies have obtained a good fit with the use of a limited number of bands [25,42], others studies [23] have required the interaction of a greater number of bands, indices or synthetic bands to achieve a better fit, as it was the case in this study.It is worth mentioning that the models selected for one year may not be appropriate to be used for a different year, because the images from different years have different digital values [44].Possibly, a good fit on the models was not achieved with the data from the images due to an unsatisfactory correction of the data, which did not allow the use of a single model for each variable for the three years [45].
To explain the slight quantitative differences between the CGC maps, the study of Khumalo and Holechek [46] was considered.They found that the precipitation of the previous year is one of the major factors that are associated with biomass production in grasslands of the Chihuahuan Desert in southern New Mexico.Their results showed that individually, precipitation in August is the variable most highly correlated to the grassland condition.When analyzing the drought data for the three years, the presence of this phenomenon was not observed in the area during this month [38].This may explain why there were basically no differences in the Kappa index of CGC for the three years.Therefore, for these variations among the years, other aspects such as spatial variation in grazing, forest fires and the phenology of the grass species should be considered.The CGC levels recorded during the three winter seasons studied could have enabled different species of wildlife, such as migratory grassland birds, to find this important habitat requirement [47].
The high reflectance values for the middle-infrared and near-infrared bands correspond to the general pattern of spectral behavior in vegetated surfaces [42].The highest values for bands 2-4, 7 in 2011 may indicate high moisture contents in the soil as well as vegetation vigor.The latter factor may be attributed to the photosynthetic activity of perennial grasses, which grow and develop during the rains of February.However, according to the comparison of the maps, this difference does not seem to have greater repercussions with respect to the marked differences between the levels of CGC and grass height between 2010 and 2011.The decline in reflectance values of the bands 1-3 in some of the quadrants in the central and the southeastern zones, seems to be more associated with the particular characteristics of the composition and structure of the vegetation in some transects, since the resulting maps did not reflect the same pattern.
The CGC and grass height showed quantitative inter-annual fluctuations and spatio-temporal variations without an apparent specific pattern.Jin et al. [48] found similar results for biomass, which is directly related to CGC and grass height.The greatest quantitative variation was registered for the variable height, which was possibly due to the fact that this variable is more susceptible to changes than CGC.This may be related to the grasses' characteristics, such as senescence [49].In addition, when grazing is practiced in an area, the height of the grassland is the first variable affected.However, and similar to the variable CGC, other abiotic factors should be discarded, such as climate variables, which may affect the characteristics of the grassland [50,51].
The spatial heterogeneity observed in all the maps on a larger scale, agreed with the results found in other studies [48,52].Even though in our study we found evidence of spatial gradients, these were not as marked as the ones found at the larger scale.Among other things, this element allows the presence of various species of grassland birds that have different CGC and grass height requirements during the winter season [14].Hence, it has been found that the CUZA region is an area of high diversity of not only bird species, but also of some others, such as grasses and forbs [53,54].According to McGranahan et al. [55], additional factors that could influence this heterogeneity are fires and grazing.
The maps generated for CGC and grass height showed moderately high accuracies that can be representative and used as a reference to help assess current and future land use development, grasslands capacities, animal production potentials and the status of grassland avian habitats [41,[56][57][58].

Conclusions
Through the use of images of the Landsat TM5 sensor, it was possible to predict the spatial distribution of CGC and grass height of a grassland, using linear models in combination with field data.Although a field validation was not performed, the maps contributed to the understanding of the characteristics and dynamics of the Chihuahuan Desert grasslands.The information generated may serve to establish better strategies for monitoring grasslands and determine those areas that are most relevant for conservation.In addition, the methodology implemented in this study can be applied to perform a program for monitoring the grassland condition.
A benefit of the methodology applied in this study is the ability to identify and focus fragmented and degraded areas.Other benefits include the capacity to extrapolate from limited ground monitoring locations with increased confidence to the possibility of monitoring inaccessible areas, as well as to analyze changes across the landscape over time.
The maps obtained with this approach can be used as layers to model the potential distribution of migratory grassland birds given that the structural characteristics of vegetation are among the main factors that determine their distribution in the area.

Figure 1 .
Figure 1.Study area.Sampling sites are marked in red.

Figure 2 .
Figure 2. Spatial location of the samples within the study area.

Figure 1 .
Figure 1.Study area.Sampling sites are marked in red.

Figure 1 .
Figure 1.Study area.Sampling sites are marked in red.

Figure 2 .
Figure 2. Spatial location of the samples within the study area.

Figure 2 .
Figure 2. Spatial location of the samples within the study area.
moderate concordance lay mainly in the fact that during 2009 the interval with the highest frequency was 0-25 cm, and for 2010 and 2011 it was 50-75 and 25-50 cm, respectively.Regarding 2009, it was characterized by grasses of low grass height along almost all of CUZA, only 8% of the area had a height of grass in the range of 25-50 cm, located primarily in the oak-pine areas.

Figure 7 .
Figure 7. Area in percentages occupied by the ranges of (a) Coverage of Grass Canopy (CGC) and (b) grass height in the region Cuchillas de la Zarca.

Figure 7 .
Figure 7. Area in percentages occupied by the ranges of (a) Coverage of Grass Canopy (CGC) and (b) grass height in the region Cuchillas de la Zarca.

Table 1 .
Average (or mean) of field data for Coverage of Grass Canopy (CGC) and grass height in the quadrants measured in Cuchillas de la Zarca.

Table 2 .
Coefficients of determination of linear regression models (p < 0.01).

Table 3 .
Estimators of the regression models for coverage grass canopy (CGC).

Table 4 .
Estimators of the regression models for grass height.

Table 5 .
Kappa index between the maps for Coverage of Grass Canopy (CGC) and grass height during the three years studied.

Table 5 .
Kappa index between the maps for Coverage of Grass Canopy (CGC) and grass height during the three years studied.