Maintaining a land base that supports safe and realistic training operations is a significant challenge for military land managers which can be informed by frequent monitoring of land condition in relation to management practices. This study explores the relationship between fire and trends in tallgrass prairie vegetation at military and non-military sites in the Kansas Flint Hills. The response variable was the long-term linear trend (2001-2010) of surface greenness measured by MODIS NDVI using BFAST time series trend analysis. Explanatory variables included fire regime (frequency and seasonality) and spatial strata based on existing management unit boundaries. Several non-spatial generalized linear models (GLM) were computed to explain trends by fire regime and/or stratification. Spatialized versions of the GLMs were also constructed. For non-spatial models at the military site, fire regime explained little (4%) of the observed surface greenness trend compared to strata alone (7% - 26%). The non-spatial and spatial models for the non-military site performed better for each explanatory variable and combination tested with fire regime. Existing stratifications contained much of the spatial structure in model residuals. Fire had only a marginal effect on surface greenness trends at the military site despite the use of burning as a grassland management tool. Interestingly, fire explained more of the trend at the non-military site and models including strata improved explanatory power. Analysis of spatial model predictors based on management unit stratification suggested ways to reduce the number of strata while achieving similar performance and may benefit managers of other public areas lacking sound data regarding land usage.
The US Department of Defense is one of the largest land stewards in the nation [
Maintaining a land base that supports safe and realistic training operations is a significant challenge impacting the mission readiness of US military forces. A number of federal laws (e.g., Sikes Act of 1960, National Environmental Policy Act of 1969) and military regulations outline specific requirements for natural resource conservation and management plans and activities, including that sustainable use should be a programmatic goal within the US Army [
As in temperate grasslands and savannas across the world, fire frequency and timing can have significant impacts on the spatial pattern of plant productivity, vegetation structure, and nutrient cycling in the Flint Hills ecoregion of North America [
In this study, promoting sustainable grasslands is considered the primary objective with sustainability being defined using results from a long-term analysis of vegetation trends derived from remotely-sensed estimates of normalized difference vegetation index (NDVI) data. The influence of fire and land use practices on grassland vegetation trends is assessed to provide land managers with site-specific feedback about progress towards their training land sustainability goals.
Spatial statistical tools are used to analyze spatiotemporal trends in vegetation and fire regime using a time series of multisource remote sensing images over the period 2001-2010. This approach has previously been used by [
This research was conducted at Fort Riley, a US Army installation located in northeastern Kansas, within the framework of the Army’s Integrated Training Area Management (ITAM) program. In conjunction with installation ITAM staff, we seek to demonstrate methods for continuously monitoring, assessing and identifying trends in key indicators of military training land sustainability over appropriate timescales, rapidly analyze data, and ensure information is available when needed by military land managers. Though military lands are the focus here, these objectives are equally applicable for effective management of other grasslands and savannas in the western US, South America, Europe, Africa, and Asia.
Fort Riley is a 41,128 ha U.S. Army installation in northeast Kansas (39˚11'N, 96˚48'W), on the Kansas River, between Junction City and Manhattan within Geary, Riley and Clay counties. The installation is located within both the Flint Hill ecoregion [
The installation is subdivided into 103 administrative units within which military training occurs. Disturbances from military training on Fort Riley include a variety of on- and off-road field maneuvers (including tracked and wheeled combat vehicle operations), mortar and artillery fire, small arms fire, and aircraft flights. For the past four decades, the majority of combat vehicle maneuvers take place in the western part of the installation.
The Konza Prairie Biological Station (KPBS) was included in this analysis to contrast vegetation change measured at Fort Riley with that at a relatively natural grassland site. KPBS (3487 ha) is owned by the Nature Conservancy and operated by the Division of Biology at Kansas State University (http://kpbs.konza.ksu.edu). One of the U.S. National Science Foundation’s Long Term Ecological Research Sites, KPBS has similar vegetation, soils and climate due to its close proximity (less than 10 km) to Fort Riley.
It is assumed that any differences observed between the two sites are caused by varying anthropogenic activities and land management practices. Though not impacted by military training, vegetation at KPBS is influenced by other human actions including prescribed burning. The KPBS site is subdivided into 48 sub-watershed areas which serve as experimental units where specific combinations of fire (by both season and frequency) and grazing (including both domestic cattle and native bison) are applied.
Data used in this analysis consists of vegetation activity change and fire regime data. All procedures used to process data for Fort Riley were also applied for KPBS. The study period is limited to 2001-2010 to match the study duration of the prerequisite vegetation trend analysis work previously completed for Fort Riley and KPBS and the lack of data needed to characterize the fire regime (high spatial resolution images) after 2010.
A vegetation change indicator was developed to characterize long-term trends in activity over the study period (
A fire regime indicator was also developed to characterize burning activities according to season and frequency by analyzing a time series of multisource remote sensing data for the period 2001-2010. The primary remote sensing data were high spatial resolution (HSR) images most appropriate for delineating burn scars at local and regional scales [
Due to persistent cloud contamination during the study period, it was impossible to acquire and classify HSR images at a frequency sufficient to adequately identify burn scars. This was especially a problem in 2004-2005,
2007-2008, and 2009-2010 when the November through March period was not well covered by Landsat 5. To fill these gaps, 135 additional MODIS Burned Area Product (MCD45) images were used to complement the HSR time series. Produced by the MODIS Fire Team at the University of Maryland (http://modis-fire.umd.edu/pages/BurnedArea.php), the Burned Area Product provided monthly images of burned pixels at a 500 m spatial resolution. Though the spatial resolution is coarse, preventing detection of small and fragmented burns, the MCD45 images were a practical and inexpensive source of data to ensure a continuous and reliable source of fire information within the time series.
Next, monthly burned area maps were created for each year and used to calculate annual burned area maps. To ensure compatibility with the vegetation change indicator, monthly burned area maps were overlaid with a grid of the same spatial resolution. Only grid cells with a minimum burned area of 80% were retained for further analysis. Summing the annual maps over the study period yielded fire frequency (
defined: Fall (September-October-November), winter (December-January-February), and spring (March-April- May). In determining both fire frequency and seasonality (i.e., fire regime), it was assumed that an individual pixel could burn only once per year with the timing of that burn defining the season. With information on both fire frequency and seasonality, a fire regime indicator consisting of 13 classes was defined.
Modeling the relationship between explanatory and response variables is a fundamental activity encountered in statistics. Multiple linear regression is used to investigate the relationship between a single response variable and several explanatory variables. Often, the response is simply a designation of one or two possible outcomes (i.e.,
a binary response). In this study, the possible outcomes are “degraded” or “non-degraded” vegetation activity and this binary response is examined in relation to fire regime and/or other anthropogenic activities represented by stratification using non-spatial and spatial generalized linear models (GLM).
The grassland vegetation activity change map was first overlaid with the fire regime data. Next, GLM models were computed on a per-pixel basis to study how well fire regime and/or stratification explained degraded vegetation activity. In the first non-spatial model, where independence between observations was assumed, a classical estimation method was used. Due to the nature of the response variable, the model was a binomial GLM with a probability of success depending on the characteristics of each observation site. The link between this probability and the linear predictor based on individual characteristics is the logit and the estimation was done by maximum likelihood [
Because correlation between observations was suspected, a second model accounting for spatial correlation was applied. In this model, an existing spatially random field was assumed which takes into account what is not explained in the independent model. Given this unobserved random field, the model adds to the linear predictor the value of the random field at the observation site. The maximum likelihood method, however, is now problematic because the likelihood cannot be derived analytically and must be approximated. Instead, the Bayesian Markov Chain Monte Carlo method based on Metropolis dynamics was used to estimate model parameters. The estimates are then taken as the posterior mean [
Both types of models (non-spatial and spatial) were evaluated with two metrics. The first was the Akaike Information Criterion (AIC) which summarizes the tradeoff between model accuracy and complexity [
Five models using different explanatory variables were tested and their ability to explain degraded vegetation activity assessed. The initial assumption was that vegetation degradation would be explained adequately by only the fire regime variable (model 1), or only stratification (model 2), or a combination of fire regime and stratification (model 3). However, the original stratification, especially that at Fort Riley, was not specifically designed for use in monitoring grassland vegetation dynamics. To address this, a simplified stratification more adapted to the response variable was also tested. This simplified stratification was developed by analyzing the distribution of the model predictor based only on the original stratification (model 2) as a function of the stratification units. Two additional models were then developed. One only with simplified stratification (model 4) and one with fire regime and simplified stratification (model 5).
Comparing model 1 vs. model 2 and model 2 vs. model 3 provided results sufficient to assess the effect of fire regime on grassland vegetation degradation. The comparison of model 2 to model 4 or model 5 allowed evaluation of the simplified stratification and its ability to better explain observed vegetation degradation as well as its potential for serving as a spatial framework for monitoring grassland condition in the future.
Each GLM model generates a probability map for vegetation degradation. Model performance is then validated using two different reference datasets. The first is the vegetation degradation class resulting from the previously mentioned MODIS NDVI trend analysis. The second is a new HSR vegetation degradation class from a change analysis performed using two Landsat 5 TM images acquired near the beginning (08/05/2001) and end (08/30/ 2010) of the study period. This method was recommended by [
The HSR dataset consists of a raster with the same spatial resolution as the vegetation activity change map (250 m) in which each grid cell corresponds to one of two classes (degraded or non-degraded) based on the interpretation of NDVI change derived from the Landsat 5 TM images. Since the GLM model estimates only a probability of vegetation degradation, only the NDVI change class for degraded vegetation is retained during this validation step.
A contingency matrix was constructed with rows representing the reference data (Landsat 5 NDVI change class for degraded vegetation or MODIS NDVI trend class for degraded vegetation) and columns with two classes of vegetation degradation probabilities (>0.50 and ≤0.50).
The ratio (1 − residual deviance/null deviance) * 100 was used to determine how well explanatory variables in the five GLM models explained variability in the vegetation activity change indicator. Associated AIC values were also calculated to characterize the relative explanatory power of each model (
To assess the effect of fire regime on vegetation activity change, results from GLM models with only the fire regime variable (model 1) were compared with those of GLM models using only stratification (model 2). Fire alone explains very little at Fort Riley (4% - 11%) versus to stratification alone (26% - 31%). The marginal effect of fire on vegetation dynamics is confirmed when also considering the results of the GLM model with fire and stratification (model 3). The combination of fire regime and stratification does not improve the explanatory power of the model (27% - 31%). This same applies to KPBS, except that the effect of fire on vegetation activity change is higher than Fort Riley (14% - 29%).
At Fort Riley, regardless of the GLM model, using a spatial approach in the statistical analysis is needed. Spatial versions of the GLM models provide a 4% - 7% improvement in performance compared to non-spatial models with the added benefit of systematically lower AIC values. The situation at KPBS is slightly different. Model spatialization is not needed if the analysis is done with the original stratification, as the percentage of vegetation activity change explained and AIC values are nearly the same with both non-spatial and spatial models. At KPBS, the original stratification already accounts for the type and distribution of explanatory variables in the models. However, when modifying the original stratification (model 4), spatialization provides a 9% improvement in model performance with a decrease in AIC values from 577 to 511.
By analyzing the range of predictor values for GLM models with only stratification (model 2) for each of the original management units (103 training areas at Fort Riley and 48 sub-watersheds at KPBS), a simplified stratification was proposed. For both study areas, this simplification consists only of five strata obtained after reclassification of training areas or sub-watersheds.
Results from the GLM model with only original stratification (model 2) were compared to those of the GLM model with only simplified stratification (model 4). At Fort Riley and KPBS, simplification explains nearly the same amount of vegetation activity change but with only five strata. As already seen with the original stratification, adding the fire regime variable to simplified stratification in a GLM model (model 5) does not improve the ability to explain vegetation activity change. Therefore, the GLM model with only simplified stratification was considered most appropriate.
GLM Models | % of Vegetation Degradation Explained | |||||
---|---|---|---|---|---|---|
Fort Riley Models | KPBS Models | |||||
n | Non-Spatial (AIC) | Spatial (AIC) | n | Non-Spatial (AIC) | Spatial (AIC) | |
Model with only fire (model 1) | 13 | 4% (8244) | 11% (7614) | 13 | 14% (806) | 29% (673) |
Model with only original stratification (model 2) | 103 | 26% (6537) | 31% (6118) | 48 | 39% (650) | 40% (653) |
Model with fire and original stratification (model 3) | 116 | 27% (6507) | 31% (6092) | 61 | 40% (661) | 43% (646) |
Model with only simplified stratification (model 4) | 5 | 23% (6575) | 27% (6234) | 5 | 37% (577) | 46% (511) |
Model with fire and simplified stratification (model 5) | 18 | 23% (6579) | 27% (6296) | 18 | 37% (587) | 40% (580) |
The vegetation degradation probability maps for Fort Riley and KPBS generated by the GLM model with only simplified stratification are shown in
Evaluation of the contingency matrix helps validate the vegetation degradation probability maps from the selected GLM model (
Results from this analysis provide two major advances. First, from a thematic point of view and contrary to what
Reference Data for Vegetation Degradation Class | Vegetation Degradation Probability Classes from GLM Model 4 | ||
---|---|---|---|
>0.50 | ≤0.50 | ||
Fort Riley | Landsat NDVI change class | 59 | 41 |
MODIS NDVI change class | 76 | 24 | |
KPBS | Landsat NDVI change class | 40 | 60 |
MODIS NDVI change class | 62 | 38 |
was expected, fire regime did not adequately explain patterns of vegetation activity change. At both sites, based on comparison of models with only original stratification (model 2) with original stratification and fire (model 3), vegetation degradation is better explained by the original stratification. This is especially true at Fort Riley and points to disturbances other than fire as the driver of local grassland trends.
Second, from a methodological point of view, the soundest approach for identifying a model that best explained vegetation activity change involved first using a model with original stratification to form a new stratification more adapted to the response variable using a spatial approach. At Fort Riley, where simplified stratification is very informative (almost 3 times more than fire regime only), the five new strata were used to characterize the probability of vegetation degradation provided by the GLM model (
The role of military training activities was also investigated. High intensity training associated with mechanized military maneuvers has been cited as the cause of increased bare soil, reduced plant cover, soil compaction, and compositional shifts in plant communities [
In strata 4 and 5, the percentage area classified as having high military training intensity is almost double that of strata 1, 2, and 3 (
Though fire appears to have only a minimal impact on vegetation dynamics based on GLM model results, fire regime was also analyzed within the new simplified strata. All strata were almost entirely burned (87%) during the 2001-2010 study period at predominantly low or very low frequency (i.e., no more than one fire every three years). In strata 5, the zone with the highest level of vegetation degradation (85%) and high intensity training (43%), fire typically occurred in spring (70%). In strata 1, where vegetation degradation and training intensity were both low (0.4% and 16%, respectively), fire occurred almost equally across the spring, fall, and winter seasons (
Spring burning in the Flint Hills ecoregion removes senescent vegetation from the previous growing season to
promote growth of new grasses. Annual spring fires favor C4 vegetation over C3 species, but also cause long- term increases in bare ground [
Considering KPBS, strata 5 had an average vegetation degradation probability of 0.69% and 61% of the area is classified as degraded based on MODIS NDVI trend analysis. Strata 1 has a very low probability of degradation (0.03) and 99% of the area is classified as non-degraded. Both strata 5 (75%) and 1 (90%) were almost entirely burned during the study period with spring the most common season of burning (96% for strata 5 and 100% for strata 1). The main difference between the strata is fire frequency and other land uses. In strata 1, 56% of the area is affected by annual burns without any grazing pressure and 36% by one fire every four years with bison grazing conducted annually between May and October. In strata 5, fire frequency is variable but well balanced (54% low and 46% high frequency), but 40% of the area is under intensive grazing by cattle (northeast) and 37% is used for trail-based recreation (north central). Spring fires at KPBS seem to be a stabilizing factor for the grassland ecosystem with vegetation degradation due mainly to grazing activities combined with fire frequency and human disturbances.
There appears to be contrasting approaches to fire management at Fort Riley. In areas where training intensities are high, spring fires are perhaps used to promote grassland health where military training disturbances are significant and frequent. In other training areas, no fire regime dominates. While the literature and results here from KPBS support spring burning as a good practice from a grassland resource perspective, when combined with frequent and intense training during the wettest time of the year, it may actually be counterproductive and amplify vegetation degradation. At minimum, the influence of military training on vegetation dynamics is important and serves to minimize any potential benefit realized from a typical Flint Hills fire regime.
To confirm such conclusions, more detailed and temporally consistent training data are required. While these data do exist, they are embedded within complex US Army database management systems and are difficult to extract for use in natural resource studies. Until such a time when training data are routinely available for land managers, the new simplified strata presented here may be used as a spatial guide to modify training schedules and prevent further declines in vegetation condition where degradation has been predicted and measured. Land managers suffering from similar data deficiencies in land usage, such as off-road vehicle use in public recreation areas, may also find utility in defining custom strata for their sites.
This work was supported by the French Agence Nationale de la Recherche through the Moduland project (ANR-11-BSH1-005). Funding to support international research collaboration was provided by Ecole d’Ingénieurs de Purpan and the Geographic Information Systems Spatial Analysis Laboratory (GISSAL) at Kansas State University. Additional travel funding was provided by the National Polytechnic Institute (INP) of Toulouse to support researcher exchange.
Anne Jacquin,Michel Goulard,J. M. Shawn Hutchinson,Thomas Devienne,Stacy L. Hutchinson, (2016) A Statistical Approach for Predicting Grassland Degradation in Disturbance-Driven Landscapes. Journal of Environmental Protection,07,912-925. doi: 10.4236/jep.2016.76081