Abstract
Restoration of polluted sites via in situ bioremediation relies heavily on the indigenous microbes and their activities. Spatial heterogeneity of microbial populations, contaminants and soil chemical parameters on such sites is a major hurdle in optimizing and implementing an appropriate bioremediation regime. We performed a grid-based sampling of an aged creosote-contaminated site followed by geostatistical modelling to illustrate the spatial patterns of microbial diversity and activity and to relate these patterns to the distribution of pollutants. Spatial distribution of bacterial groups unveiled patterns of niche differentiation regulated by patchy distribution of pollutants and an east-to-west pH gradient at the studied site. Proteobacteria clearly dominated in the hot spots of creosote pollution, whereas the abundance of Actinobacteria, TM7 and Planctomycetes was considerably reduced from the hot spots. The pH preferences of proteobacterial groups dominating in pollution could be recognized by examining the order and family-level responses. Acidobacterial classes came across as generalists in hydrocarbon pollution whose spatial distribution seemed to be regulated solely by the pH gradient. Although the community evenness decreased in the heavily polluted zones, basal respiration and fluorescein diacetate hydrolysis rates were higher, indicating the adaptation of specific indigenous microbial populations to hydrocarbon pollution. Combining the information from the kriged maps of microbial and soil chemistry data provided a comprehensive understanding of the long-term impacts of creosote pollution on the subsurface microbial communities. This study also highlighted the prospect of interpreting taxa-specific spatial patterns and applying them as indicators or proxies for monitoring polluted sites.