Received 21 June 2016; accepted 18 July 2016; published 21 July 2016
Natural selection pressures act as filters for establishment and succession of vegetative species. Via dispersal, environmental, and physiological filters, community composition develops to best befit the condition of any given area   . Presence of anthropogenic contaminants adds an additional filter to the system, reducing fecundity of sensitive species, and increasing that of tolerant ones  . Loss or suppression of dominant/key species can have significant effects on ecosystem function  , but even changes in non-dominant species density may have significant community level impacts  . Community shifts have been previously investigated and documented for a range of anthropogenic contaminants   from mining operations  -  to radiation exposure  . Richness and diversity of plant species and functional diversity may be vital to ecosystem stability and function  -  .
Long term impacts of contaminants on vegetation are studied less than shorter exposures  . This poses a problem as many anthropogenic contaminants are recalcitrant (i.e. heavy metals, explosives, agrochemicals, etc.)  -  when released into the environment and elicit chronic stress for many years in surrounding flora and fauna  . Even short term exposure to anthropogenic contaminants may have long-lasting implications  . The persistent nature of these compounds adds to this complexity of community level responses. Chronic exposure over extended periods of time has significant impacts on established communities    and the seedbank of contaminated areas   .
Once released into the environment anthropogenic contaminants exhibit a range of behaviors   . Some rapidly bind to soil particles, while others are not constrained to a static location  . Compound mobility and structure play large roles in contaminant behavior in the environment  . For example, heavy metals behave differently from pharmaceuticals, agrochemicals, and explosives   . Nanoparticles of a particular element behave very differently than larger particles of the same element  . Just as there are numerous behaviors of contaminants, responses of vegetation to these compounds are greatly varied.
Impacts on vegetation can range from morphological damage to physiological interruptions and exact response depends on the contaminant and plant species in the target area. Plant functional traits (e.g. life history, leaf morphology, root structure, photosynthetic pathway, etc…) can also influence response to stress  -  and ability to persist in contaminated conditions. As plants age, biochemical processes change  -  , which may cause variations in response at different life stages. Mature plants have more biomass relative to earlier life stage which may allow the plant to withstand stress associated with the explosive compounds and then recover  . There are dualities present with tolerance to contaminants; plants capable of surviving may be harmed more as exposure time increases or may pass the contaminant through the system  . High variability in impacts on vegetation based on species, functional group, and even life stage adds a selection pressure to con- taminated areas.
Our study objective was to investigate long term impacts of explosives presence on the current vegetative community of a minefield contaminated 17 years prior with three different explosive compounds/mixtures. We hypothesized that 1) community richness and diversity in contaminated areas would not exhibit change from reference plots as species fill in vacant niches, 2) species composition would differ from reference plots as presence of contaminants filter out intolerant species, and 3) changes in species functional traits would differ in contaminated treatment plots from reference.
2. Materials and Methods
The two most commonly used explosive compounds are hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) and 2-methyl-1,3,5-trinitrobenzene (TNT). These compounds are the most widely distributed organic explosive contaminates found in the environment    and both can have significant impacts on vegetation health. Munitions contain either a single compound (TNT most commonly) or more commonly, a mixture of compounds (TNT and RDX). Composition B (Comp B) a mixture of RDX (60%) and TNT (40%) is a historically common munitions mixture. It was used extensively throughout World War II up through the 1950s, and is still used today  .
2.2. Study Site
Data were collected from a privately owned experimental minefield in South Carolina, USA (33˚54'5.39"N, 82˚2'1.79"W)   (Figure 1). Geologically the site is located on the Carolina terrane (ct) composed of clastic rocks in the upper region and intermediate to felsic pyroclastic rocks (>3 km) in the lower. The soil at the plots was Herndon very fine sandy loam and Kirksey silt loam, both with 2% - 6% slopes  .
The property was owned and managed by Force Protection Inc. at the time of data acquisition. Initially used for a microbial study in 1999 this site was razed to the ground and sectioned off into 5 × 5 m plots. Within a
Figure 1. Location of field sites in Edgefield County, South Carolina. Colored circles show field plots where data was collected. Red marks reference plots, yellow marks RDX plots, blue TNT, and green Composition B.
number of those plots explosive compounds were buried in “dummy landmines” to mimic leaching of explosives from unexploded ordnance (UXOs). These dummy mines were boxes filled with TNT (2-methyl-1,3,5- trinitrobenzene) and RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine), or a mixture of the two compounds, Composition B (Comp B). Comp B is comprised of 60% RDX, 39% TNT, and 1% wax binder. Dummy mines were designed to be similar to antipersonnel landmines and contained between 75 - 125 g of the designated explosive compound   . Reference areas were established upslope of the contaminated plots in the razed area. Once the dummy mines had been planted the area was left to naturally revegetate.
In each of the treatment and reference areas 1x1m plots were delineated (Figure 1) and contaminated plots were checked using site maps from the property owner prior to data collection to ensure presence of a dummy mine. Herbaceous species were identified and percent cover established in each plot. Woody species abundance was also collected after the species had been identified. Minimally destructive metrics were used due to site protocol.
2.3. Species Diversity
Density and frequency were calculated for all species. Species richness (S; Equation (1)), Shannon-Weiner (H′; Equation (2)), and Simpson (D1; Equation (3)) metric values were calculated for total species as well as herbaceous and woody species separately. Both H′ and D1 were used to take into consideration both rare, and common species respectively   . Plants not readily identifiable in the field were transported to the lab and identified using Radford et al. (2010)  . All species metrics were compared across treatment groups via analysis of variance (ANOVA) using a Tukey’s post hoc test.
Equation (1). Equation for calculating species richness (S) where “s” is the number of species present and “N” represents the number of individuals present.
Equation (2). Equation for calculating Shannon-Weiner (H′) values where “pi” is the proportional abundance of species “i” in the plots.
Equation (3). Equation for calculating Simpson (D1) values where “pi” is the proportional abundance of species “i” in the plots.
2.4. Species Composition
Jaccard’s index of similarity was calculated for total, herbaceous, and woody species composition to investigate whether similarities in species composition were present across groups. Cluster analysis (with Sorenson distance) was used to determine how the plots relate to one another based on species composition. A dendrogram was constructed using PC-Ord (MjM Software Design) to visually represent the cluster analysis results. To determine if the resulting groups were significantly different, a multi-response permutation procedure (MRPP) was used (with Sorenson distance). MRPP was followed by pairwise comparison to compare species composition among groups.
2.5. Functional Trait Composition
Due to constraints on site access, functional trait data for all species recorded in the field were aggregated from literature sources  -  . Functional traits included: life history, native status, growth form, habit, maximum height, root structure, leaf morphology, seed size, dispersal mechanism, and photosynthetic pathway. To standardize trait values, community weighted mean (CWM) for each plot were calculated using the FD package in R   . CWM allows categorical and continuous functional traits to be weighted by abundance of representative species, limiting any bias which may come from rare species or traits. This allows for the identification of dominant traits and provides insight into ecosystem processes  . CWM values for each functional trait were compared across treatment groups via ANOVA and Tukey’s post-hoc test. Values were also used in a cluster analysis and MRPP (both with Sorenson distance) to investigate plot grouping. MRPP was followed by a pairwise comparison to compare species composition among groups.
3.1. Species Diversity
Woody species density was variable across treatment groups for most species with Composition B plots distinguishing themselves from both reference and other treatment plots (Figure 2). RDX and TNT showed no significant change in woody species density compared to reference plots. Comp B plots exhibited significantly higher density for Ulmus alata Michx. (F = 6.40, p < 0.01) and Quercus virginiana Mill. (F = 3.92, p = 0.02) relative to reference values (p < 0.01 and p < 0.01, respectively). Ligustrum sinensce Lour, an invasive species, was common across all sites (F = 2.15, p = 0.11) but was significantly more dense in Comp B treated areas when compared to densities in TNT plots (p = 0.02). Diospyros virginiana L. was the only woody species which exhibited a significantly reduced presence in all three contaminated areas relative to reference (Figure 1; F = 6.11, p < 0.01). Density of the other three woody species was similar across treatments; however, only one of those, Rhus copallina L., was present in abundance (Figure 2).
No differences were observed in richness, diversity, or Simpson index for herbaceous and woody species combined (p < 0.05; Table 1). Separating out herbaceous and woody species did provide some differentiation of treatment groups. Species richness remained unchanged for herbaceous species yet woody species showed significant variation (F = 3.37, p = 0.03), between Composition B and both RDX and TNT. Shannon-Weiner diversity, which is sensitive to rare species, was significantly different for herbaceous (F = 4.41, p = 0.01) and woody (F = 3.79, p = 0.02) plant communities (Table 1). Only woody species showed a significant difference between the reference and treatment groups. Simpson’s index, which is sensitive to abundant species, also showed significant differences between treatment groups in herbaceous and woody species (Table 1). This metric was the
Figure 2. Woody species density across treatment types for reference and explosives contaminated plots. Letters denote statistical grouping. The lack of letter codes signifies that no difference was present between treatments.
Table 1. Mean and standard error of diversity metrics (richness [S], Shannon-Weiner [H′], and Simpson index [D1] values) for reference and explosives contaminated plots.
only one to show a significant difference between the reference and a treatment group for both vegetation types. Overall it appeared that abundant species, rather than rare ones, were better able to differentiate reference and treatment communities suggesting shifts in species dominants with presence of explosives contamination.
3.2. Species Composition
Jaccard’s index of similarity revealed that community composition across treatments was between 40% and 60% similar to one another (Table 2). Herbaceous community had similar values. The woody assemblage exhibited the least amount of similarity among treatments groups ranging from 14% to 43%. The reference in particular saw a large reduction for woody species exhibiting only 14% - 28% similarity.
Cluster analysis of species composition across treatments delineated three very distinct groups: 1) reference plots, 2) single contaminant plots (RDX and TNT), and 3) mixed contaminant plots (Composition B; Figure 3). Dendrogram grouping was supported by multi-response permutation procedure (MRPP) of species composition. MRPP showed that communities present in the plots were able to define the preexisting grouping scheme (A = 0.11, p < 0.01). Pairwise comparisons showed that reference communities were significantly different from those in contaminated areas as well as separate single compound contaminant plots (RDX and TNT) from mixed
Table 2. Jaccard’s index of similarity values for reference (Ref) and explosives contaminated plots.
Figure 3. Dendrogram of plots based on species composition for reference and explosives contaminated plots. Axis along the top represent information remaining and the red dashed line represents the cutoff point used here.
compound treatment plots (Comp B). RDX and TNT were not different from one another based on species composition (Table 3).
3.3. Functional Composition
The number of species representing animal dispersed (F = 3.15, p = 0.04), annual (F = 7.91, p<0.01), graminoid (F = 6.19, p < 0.01), invasive (F = 3.14, p = 0.03), monocot (F = 5.45, p < 0.01), perennial (F = 3.41, p = 0.03), trailing (F = 3.88, p = 0.02), tree (F = 4.730, p < 0.01), and vine (F = 5.26, p < 0.01) functional groups were significantly different between groups (Table 4). RDX plots were most often responsible for the detected differences, exhibiting far more annual (p < 0.01), monocot (p = 0.02), graminoid (p < 0.01), trailing (p = 0.02), and vine (p < 0.01) species, as well as fewer animal dispersed species (p = 0.04) relative to reference. RDX also distinguished itself from TNT and Comp B plots with significantly higher counts of vine (Comp B p = 0.03), grass (TNT p < 0.01, Comp B p = 0.01), annual (Comp B p = 0.02), and monocot (TNT p < 0.01, Comp B p = 0.02) species. Additionally Comp B plots exhibited elevated numbers of tree and invasive species but was only significantly different from TNT plots (p < 0.01 and p = 0.02 respectively).
Table 3. P-values from pairwise comparisons of species and functional composition of treatment groups according to multi response permutation procedure (MRPP) analysis for reference and explosives contaminated plots.
Bold text indicates the presence of a significant difference (p < 0.05).
Table 4. Mean and standard error of species count in each treatment for various functional types for reference and explosives contaminated plots.
Letter codes designate statistical grouping. Rows lacking letter codes have no significant differences present. Bold text signifies a significant difference from reference values.
Community weighted means (CWM) revealed significant differences across treatments for several functional traits. CWM for species native status exhibited significant reductions (F = 3.59, p = 0.02) for TNT relative to RDX (p = 0.03) and reference (p = 0.04) plots. A similar trend was observed for life history (F = 3.78, p = 0.02) as values for TNT plots were significantly lower than those of RDX (p = 0.04) and reference (p = 0.03) ones. TNT plots also exhibited significantly higher CWM values for photosynthetic mechanism (F = 4.76, p = 0.01) compared to RDX (p = 0.01) and reference (p = 0.03). Despite these differences cluster analysis was unable to produce any distinctive pattern and the MRPP was unable to define the preexisting grouping scheme (A = −0.004, p = 0.48).
Due to the inherent risk of injury at many locations where unexploded ordnance (UXOs) are present, field monitoring of vegetative responses is rare; legacy community level information is even more so. As such, the objective of our study was to investigate and quantify long term impacts of explosives soil contamination on a vegetative community. Our initial hypotheses were that 1) community richness and diversity in contaminated areas would not exhibit change from reference plots as species would fill in vacant niches, 2) species composition would differ from reference plots as the presence of contaminants filtered out intolerant species, and 3) species functional traits would differ in contaminated plots relative to those of reference. Based on the results presented here all three hypotheses were supported.
Presence or absence of explosives compounds was unable to differentiate treatment groups from one another using standard species diversity metrics (S, H′, D1). Woch et al. (2016)  found a similar response along a heavy metal gradient and suggested that similarity was due to less tolerant species being removed and more tolerant species coming in and establishing in the newly formed gaps and this trend had also been observed repeatedly in systems under natural disturbance regimes  -  . This idea is further supported by the significant increase in annual and monocot species presence in RDX plots as well as the increase in annuals CWM in TNT plots. Such shifts are common on disturbed sites and thought to have increased tolerance for contaminant presence  -  . There were shifts towards graminoid dominance in RDX plots. Graminoid species are highly resilient to RDX presence  , however, graminoids are sensitive to TNT and other explosive mixtures   . TNT plots showed a novel shift in functional dominance, via CWM values, with the significant increase in C4 plant presence. This dominance of C4 species in TNT plots is particularly interesting as TNT induces significant impacts to photosynthetic mechanisms yet C4 species appear to be more tolerant overall to presence of explosive compounds, and may have led to this shift  .
Simpson’s index of similarity showed that species composition was ~60% or less similar for all treatment groups when all species were considered. Given the proximity of the plots (5 - 50 m) and the homogeneity of both soil and terrain the logical cause of this difference was the presence of contaminants. Observed similarities in community composition appeared to be largely driven by herbaceous species as there was little change in index values when woody species were removed from the analysis. In contrast woody species composition exhibited high variance in similarity across groups with values being as low as 14% and as high as 43%. Reduced species similarity and alterations to dominant functional groups among plots in conjunction with a lack of change for species richness and diversity suggest that compositionally plots were also different across treatment groups. It appeared that the field plots exhibited an altered successional trajectory emphasizing the need to investigate broad scale effects of anthropogenic activities on natural systems.
Community composition produced the clearest differentiation of plots with three distinct groups 1) reference, 2) single contaminant, and 3) multiple contaminant plots. Functionally differences were observed for TNT and RDX in terms of representative species count (Figure 4) and CWM values, yet there were no significant impacts for Comp B plots. Similar response patterns have been observed at the individual scale for both morphological and physiological characteristics  . TNT appears to influence belowground structures and photosynthetic operation while RDX impacts aboveground structures as well as photosynthetic operation  -  . Unlike the two constituent compounds, Comp B tends to produce reduced impacts, requiring larger concentrations to induce the same level of impairment at the species level. Differences were present for CWM values denoting shifts in dominant traits between treatment groups  and potentially altering ecosystem processes   .
Small scale impacts of disturbance can have significant direct and indirect influences on larger scale processes  . At the individual level explosives have an array of morphological and physiological impacts on vegetation  which can influence success and fecundity of colonizing species  . By limiting intolerant species and pro-
Figure 4. Graphical representation of community shifts across treatment area.
viding niches for tolerant ones contaminants act the same way that environmental filters do  . Considerable effort has been put into understanding short term impacts of anthropogenic contaminants on vegetative health, but more work is needed investigating long term impacts of exotoxins.
Species diversity and composition cannot fully explain ecosystem level impacts of a disturbance. By quantifying functional traits and composition of a community a much more complete understanding of disturbance effects and influences on recolonization trends can be gained   -  . From a functional trait perspective there was less distinction among treatment groups, with most shifts occurring in RDX contaminated areas. RDX plots had far more annuals, monocots, graminoids and vines, as well as fewer animal dispersed species. Increases in annuals and graminoids have been observed in the presence of contaminants    . Compositional changes away from zoochorous species has been connected with harsh habitats; particularly young or disturbed ones   . Functional richness and divergence have often been linked to community assembly processes  -  or ecosystem functioning  -  .
Our results suggest that long-term impacts of explosives contamination affect community composition and functional traits from chronic leaching of point source toxic compounds over nearly two decades. The dummy mines were designed and placed in the field to closely mimic anti-personnel landmines allowing for direct comparisons to locations with UXO presence. Given the global distribution of UXOs and the continued use of landmines  understanding the ecological impacts of these devices is of great importance. Concentration of UXOs at the field site is relatively small suggesting that vegetative community condition on contaminated systems (military training grounds, past and current warzones, industrial landscapes, etc.) may exhibit more significant responses. While UXOs have devastating impacts to human life and health they also pose a serious ecological risk both from physical damage via explosive potential and from toxic compounds that may enter the soil.
Legacy community level impacts of explosive compounds are capable of being identified using various species and functional trait metrics. Presence of explosive compounds had impacts on community species and functional composition yet elicited no response in standard species diversity metrics. RDX and TNT contaminated plots had shifts in dominant functional traits, suggesting an influx of more tolerant species. RDX possessed the most unique combination of functional traits relative to reference and other treatment groups. Comp B had no significant impact on any metric which agreed with findings in recent literature. More long term broad scale investigation of contaminants and disturbance are needed to fully understand anthropogenic impacts on natural systems.
The authors would like to thank Paul Manley for assisting with data collection and species identification.
 Hansen, W.D., Romme, W.H., Ba, A. and Turner, M.G. (2016) Shifting Ecological Filters Mediate Postfire Expansion of Seedling Aspen (Populus tremuloides) in Yellowstone. Forest Ecology and Management, 362, 218-230.
 Via, S.M. and Zinnert, J.C. (2016) Impacts of Explosive Compounds on Vegetation: A Need for Community Scale Investigations. Environmental Pollution, 208, 495-505.
 Newman, M.C. and Clements, W. (2008) Biotic and Abiotic Factors That Regulate Communities. In: Newman, M.C. and Clements, W., Eds., Ecotoxicology. A Comprehensive Treatment, Taylor and Francis/CRC Press, Boca Raton, 379-407.
 Kim, T.N., Spiesman, B.J., Buchanan, A.L., Hakes, A.S., Halpern, S.L., Inouye, B.D. and Kilanowski, A.L. (2015) Selective Manipulation of a Non-Dominant Plant and Its Herbivores Affects an Old-Field Plant Community. Plant Ecology, 216, 1029-1045.
 Prach, K., Rehounková, K., Lencová, K., Jirová, A., Konvalinková, P., Mudrák, O., Student, V., Vaněcek, Z., Tichy, L., Petrík, P. and Smilauer, P. (2014) Vegetation Succession in Restoration of Disturbed Sites in Central Europe: The Direction of Succession and Species Richness across 19 Seres. Applied Vegetation Science, 17, 193-200.
 Brady, C.J. and Noske, R.A. (2010) Succession in Bird and Plant Communities over a 24-Year Chronosequence of Mine Rehabilitation in the Australian Monsoon Tropics. Restoration Ecology, 18, 855-864.
 Zhang, H. and Chu, L.M. (2011) Plant Community Structure, Soil Properties and Microbial Characteristics in Revegetated Quarries. Ecological Engineering, 37, 1104-1111.
 Pandey, B., Agrawal, M. and Singh, S. (2014) Coal Mining Activities Change Plant Community Structure Due to Air Pollution and Soil Degradation. Ecotoxicology, 23, 1474-1483.
 Woodwell, G.M. and Rebuck, A.L. (1967) Effects of Chronic Gamma Radiation on the Structure and Diversity of an Oak-Pine Forest. Ecological Monographs, 37, 53-69.
 Cardinale, B.J., Srivastava, D.S., Duffy, J.E., Wright, J.P., Downing, A.L., Sankaran, M. and Jouseau, C. (2006) Effects of Biodiversity on the Functioning of Trophic Groups and Ecosystems. Nature, 443, 989-992.
 Cadotte, M.W., Cavender-Bares, J., Tilman, D. and Oakley, T.H. (2009) Using Phylogenetic, Functional and Trait Diversity to Understand Patterns of Plant Community Productivity. PLoS ONE, 4, e5695.
 Flynn, D.F.B., Mirotchnick, N., Jain, M., Palmer, M.I. and Naeem, S. (2011) Functional and Phylogenetic Diversity as Predictors of Biodiversity-Ecosystem Function Relationships. Ecology, 92, 1573-1581.
 Feng, G., Svenning, J.C., Mi, X., Jia, Q., Rao, M., Ren, H., Bebber, D.P. and Ma, K. (2014) Anthropogenic Disturbance Shapes Phylogenetic and Functional Tree Community Structure in a Subtropical Forest. Forest Ecology and Management, 313, 188-198.
 Klein, W. and Scheunert, I. (1982) Bound Pesticide Residues in Soil, Plants and Food with Particular Emphasis on the Application of Nuclear Techniques. International Symposium on Agrochemicals: Fate in Food and the Environment Using Isotope Techniques, Rome, 7-11 June 1982, 177-205.
 Rylott, E.L. and Bruce, N.C. (2009) Plants Disarm Soil: Engineering Plants for the Phytoremediation of Explosives. Trends in Biotechnology, 27, 73-81.
 Ansari, A.A., Gill, S.S., Gill, R., Lanza, G.R. and Newman, L. (Eds.) (2015) Phytoremediation: Management of Environmental Contaminants. Vol. 1, Springer International Publishing, Gewerbestrasse.
 Travis, E.R., Bruce, N.C. and Rosser, S.J. (2008) Microbial and Plant Ecology of a Long-Term TNT-Contaminated Site. Environmental Pollution (Barking, Essex: 1987), 153, 119-126.
 Huopalainen, M., Tuittila, E., Vanha-Majamaa, I., Nousiainen, J. and Vasander, H. (2000) The Potential of Soil Seed Banks for Revegetation of Bogs in SW Finland after Long-Term Aerial Pollution. Annales Botanici Fennici, 37, 1-9.
 Huopalainen, M., Tuittila, E.S., Vanha-Majamaa, I., Nousiainen, H., Laine, J. and Vasander, H. (2001) Effects of Long-Term Aerial Pollution on Soil Seed Banks in Drained Pine Mires in Southern Finland. Water, Air, and Soil Pollution, 125, 69-79.
 Walsh, M.R., Walsh, M.E., Ramsey, C.A., Thiboutot, S. and Ampleman, G. (2015) Energetics Residues Deposition from Training with Large Caliber Weapon Systems. European Conference of Defence and the Environment Conference Proceedings, 100-109.
 Pilon-Smits, E.A. and Freeman, J.L. (2006) Environmental Cleanup Using Plants: Biotechnological Advances and Ecological Considerations. Frontiers in Ecology and the Environment, 4, 203-210.
 Zhuang, J., Yu, H.Q., Henry, T.B. and Sayler, G.S. (2015) Fate and Toxic Effects of Environmental Stressors: Environmental Control. Ecotoxicology, 24, 2043-2048.
 Flynn, D.F.B., Gogol-Prokurat, M., Nogeire, T., Molinari, N., Richers, B.T., Lin, B.B., Simpson, N., Mayfield, M.M. and DeClerck, F. (2009) Loss of Functional Diversity under Land Use Intensification across Multiple Taxa. Ecology Letters, 12, 22-33.
 Vandewalle, M., Purschke, O., de Bello, F., Reitalu, T., Prentice, H.C., Lavorel, S., Johansson, L.J. and Sykes, M.T. (2014) Functional Responses of Plant Communities to Management, Landscape and Historical Factors in Semi-Natural Grasslands. Journal of Vegetation Science, 25, 750-759.
 Kusumoto, B., Shiono, T., Miyoshi, M., Maeshiro, R., Fujii, S.J., Kuuluvainen, T. and Kubota, Y. (2015) Functional Response of Plant Communities to Clearcutting: Management Impacts Differ between Forest Vegetation Zones. Journal of Applied Ecology, 52, 171-180.
 Donaldson, J.R., Stevens, M.T., Barnhill, H.R. and Lindroth, R.L. (2006) Age-Related Shifts in Leaf Chemistry of Clonal Aspen (Populus tremuloides). Journal of Chemical Ecology, 32, 1415-1429.
 Juvany, M., Müller, M. and Munné-Bosch, S. (2013) Plant Age-Related Changes in Cytokinins, Leaf Growth and Pigment Accumulation in Juvenile Mastic Trees. Environmental and Experimental Botany, 87, 10-18.
 Verkleij, J.A.C., Golan-Goldhirsh, A., Antosiewisz, D.M., Schwitzguébel, J.-P. and Schroder, P. (2009) Dualities in Plant Tolerance to Pollutants and Their Uptake and Translocation to the Upper Plant Parts. Environmental and Experimental Botany, 67, 10-22.
 Hawari, J., Halasz, A., Sheremata, T., Beaudet, S., Groom, C., Paquet, L., Rhofir, C., Ampleman, G. and Thiboutot, S. (2000) Characterization of Metabolites during Biodegradation of Hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) with Municipal Anaerobic Sludge. Applied and Environmental Microbiology, 66, 2652-2657.
 Khatisashvili, G., Gordeziani, M., Adamia, G., Kvesitadze, E., Sadunishvili, T. and Kvesitadze, G. (2009) Higher Plants Ability to Assimilate Explosives. World Academy of Science, Engineering and Technology, 33, 256-270.
 Pichtel, J. (2012) Distribution and Fate of Military Explosives and Propellants in Soil: A Review. Applied and Environmental Soil Science, 2012, Article ID: 617236.
 Oxley, J.C., Smith, J.L., Resende, E. and Pearce, E. (2003) Quantification and Aging of the Post-Blast Residue of TNT Landmines. Journal of Forensic Science, 48, 742-753.
 Morris, E.K., Caruso, T., Buscot, F., Fischer, M., Hancock, C., Maier, T.S., Meiners, T., Müller, C., Obermaier, E., Prati, D. and Socher, S.A. (2014) Choosing and Using Diversity Indices: Insights for Ecological Applications from the German Biodiversity Exploratories. Ecology and Evolution, 4, 3514-3524.
 Garnier, E., Cortez, J., Billes, G., Navas, M.L., Roumet, C., Debussche, M., Laurent, G., Blanchard, A., Aubry, D., Bellmann, A., Neill, C. and Toussaint, J.P. (2004) Plant Functional Markers Capture Ecosystem Properties during Secondary Succession. Ecology, 85, 2630-2637.
 Barbaro, L., Dutoit, T., Anthelmea, F. and Corcket, E. (2004) Respective Influence of Habitat Conditions and Management Regimes on Pre-Alpine Calcareous Grasslands. Journal of Environmental Management, 72, 261-275.
 Stevens, C.J., Thompson, K., Grime, J.P., Long, C.J. and Gowing, D.J.G. (2010) Contribution of Acidification and Eutrophication to Declines in Species Richness of Calcifuge Grasslands along a Gradient of Atmospheric Nitrogen Deposition. Functional Ecology, 24, 478-484.
 Kurek, P., Kapusta, P. and Holeksa, J. (2014) Burrowing by Badgers (Meles meles) and Foxes (Vulpes vulpes) Changes Soil Conditions and Vegetation in a European Temperate Forest. Ecological Research, 29, 1-11.
 Schnoor, J.L., Light, L.A., McCutcheon, S.C., Wolfe, N.L. and Carreia, L.H. (1995) Phytoremediation of Organic and Nutrient Contaminants. Environmental Science & Technology, 29, 318A-323A.
 Quist, M.C., Fay, P.A., Guy, C.S., Knapp, A.K. and Rubenstein, B.N. (2003) Military Training Effects on Terrestrial and Aquatic Communities on a Grassland Military Installation. Ecological Applications, 13, 432-442.
 Winfield, L.E., Rodger, J.H. and D’surney, S.J. (2004) The Responses of Selected Terrestrial Plants to Short (< 12 Days) and Long Term (2, 4 and 6 Weeks) Hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) Exposure. Part I: Growth and Developmental Effects. Ecotoxicology, 13, 335-347.
 Panz, K., Miksch, K. and Sójka, T. (2013) Synergetic Toxic Effect of an Explosive Material Mixture in Soil. Bulletin of Environmental Contamination and Toxicology, 91, 555-559.
 Zhang, H. and Chu, L.M. (2013) Changes in Soil Seed Bank Composition during Early Succession of Rehabilitated Quarries. Ecological Engineering, 55, 43-50.
 Best, E.P., Smith, T., Hagen, F.L., Dawson, J. and Torrey, A.J. (2007) Candidate Herbaceous Plants for Phytoremediation of Energetics on Ranges.
 Peterson, M.M., Horst, G.L., Shea, P.J. and Comfort, S.D. (1998) Germination and Seedling Development of Switchgrass and Smooth Bromegrass Exposed to 2,4,6-Trinitrotoluene. Environmental Pollution, 99, 53-59.
 Via, S.M., Zinnert, J.C., Butler, A.D. and Young, D.R. (2014) Comparative Physiological Responses of Morella cerifera to RDX, TNT, and Composition B Contaminated Soils. Environmental and Experimental Botany, 99, 67-74.
 Via, S.M., Zinnert, J.C. and Young, D.R. (2015) Differential Effects of Two Explosive Compounds on Seed Germination and Seedling Morphology of a Woody Shrub, Morella cerifera. Ecotoxicology (London, England), 24, 194-201.
 Zinnert, J.C., Via, S.M. and Young, D.R. (2013) Distinguishing Natural from Anthropogenic Stress in Plants: Physiology, Fluorescence and Hyperspectral Reflectance. Plant and Soil, 366, 133-141.
 Villéger, S., Mason, N.W.H. and Mouillot, D. (2008) New Multidimensional Functional Diversity Indices for a Multifaceted Framework in Functional Ecology. Ecology, 89, 2290-2301.
 Vile, D., Shipley, B. and Garnier, E. (2006) Ecosystem Productivity Can Be Predicted from Potential Relative Growth Rate and Species Abundance. Ecology Letters, 9, 1061-1067.
 Mokany, K., Ash, J. and Roxburgh, S. (2008) Functional Identity Is More Important than Diversity in Influencing Ecosystem Processes in a Temperate Native Grassland. Journal of Ecology, 96, 884-893.
 Heffernan, J.B., Soranno, P.A., Angilletta, M.J., Buckley, L.B., Gruner, D.S., Keitt, T.H., Kellner, J.R., Kominoski, J.S., Rocha, A.V., Xiao, J. and Harms, T.K. (2014) Macrosystems Ecology: Understanding Ecological Patterns and Processes at Continental Scales. Frontiers in Ecology and the Environment, 12, 5-14.
 Ding, Y., Zang, R.G., Letcher, S.G., Liu, S.R. and He, F.L. (2012) Disturbance Regime Changes the Trait Distribution, Phylogenetic Structure and Community Assembly of Tropical Rain Forests. Oikos, 121, 1263-1270.
 Letcher, S.G. (2010) Phylogenetic Structure of Angiosperm Communities during Tropical Forest Succession. Proceedings of the Royal Society of London: Biological Sciences, 277, 97-104.
 Letcher, S.G., Chazdon, R.L., Andrade, A.C.S., Bongers, F., Breugel, M., Finegan, B., Laurance, S.G., Mesquita, R.C.G., Martínez-Ramos, M. and Williamson, G.B. (2012) Phylogenetic Community Structure during Succession: Evidence from Three Neotropical Forest Sites. Perspectives in Plant Ecology, Evolution and Systematics, 14, 79-87.
 Norden, N., Letcher, S.G., Boukili, V., Swenson, N.G. and Chazdon, R. (2012) Demographic Drivers of Successional Changes in Phylogenetic Structure across Life-History Stages in Plant Communities. Ecology, 93, S70-S82.
 Rehounková, K. and Prach, K. (2010) Life-History Traits and Habitat Preferences of Colonizing Plant Species in Long-Term Spontaneous Succession in Abandoned Gravel-Sand Pits. Basic and Applied Ecology, 11, 45-53.
 Mouchet, M.A., Villéger, S., Mason, N.W.H. and Mouillot, D. (2010) Functional Diversity Measures: An Overview of Their Redundancy and Their Ability to Discriminate Community Assembly Rules. Functional Ecology, 24, 867-876.
 Mason, N.W.H., Richardson, S.J., Peltzer, D.A., de Bello, F., Wardle, D.A. and Allen, R.B. (2012) Changes in Coexistence Mechanisms along a Long-Term Soil Chronosequence Revealed by Functional Trait Diversity. Journal of Ecology, 100, 678-689.
 Spasojevic, M.J. and Suding, K.N. (2012) Inferring Community Assembly Mechanisms from Functional Diversity Patterns: The Importance of Multiple Assembly Processes. Journal of Ecology, 100, 652-661.
 Mason, N.W.H., De Bello, F., Mouillot, D., Pavoine, S. and Dray, S. (2013) A Guide for Using Functional Diversity Indices to Reveal Changes in Assembly Processes along Ecological Gradients. Journal of Vegetation Science, 24, 794-806.
 Mouillot, D., Graham, N.A.J., Villéger, S., Mason, N.W.H. and Bellwood, D.R. (2013) A Functional Approach Reveals Community Responses to Disturbances. Trends in Ecology and Evolution, 28, 167-177.