Disturbance Alters the Relative Importance of Topographic and Biogeochemical Controls on Microbial Activity in Temperate Montane Forests

Fire and pathogen-induced tree mortality are the two dominant forms of disturbance in Western U.S. montane forests. We investigated the consequences of both disturbance types on the controls of microbial activity in soils from 56 plots across a topographic gradient one year after the 2012 High Park wildfire in Colorado. Topsoil biogeochemistry, soil CO2 efflux, potential exoenzyme activities, and microbial biomass were quantified in plots that experienced fire disturbance, beetle disturbance, or both fire and beetle disturbance, and in plots where there was no recent evidence of disturbance. Soil CO2 efflux, N-, and P-degrading exoenzyme activities in undisturbed plots were positively correlated with soil moisture, estimated from a topographic wetness index; coefficient of determinations ranged from 0.5 to 0.65. Conversely, the same estimates of microbial activities from fire-disturbed and beetle-disturbed soils showed little correspondence to topographically inferred wetness, but demonstrated mostly negative relationships with soil pH (fire only) and mostly positive relationships with DOC/TDN (dissolved organic carbon/total dissolved nitrogen) ratios for both disturbance types. The coefficient of determination for regressions of microbial activity with soil pH and DOC/TDN reached 0.8 and 0.63 in fireand beetle-disturbed forests, respectively. Drivers of soil microbial activity change as a function of disturbance type, suggesting simple mathematical models are insufficient in capturing the impact of disturbance in forests.


Introduction
Forests account for half of Earth's terrestrial organic carbon (C) stocks by storing carbon as plant biomass, forest litter, and soil organic matter (SOM) [1,2].Wildfire, insect and disease outbreaks, droughts, and other natural disasters are an integral component of forest dynamics, but changes to the frequency and severity of these disturbances world-wide present a global challenge for monitoring and managing these forests and the ecosystem services they provide [3,4].Forests in the Western U.S. are experiencing larger wildfire disturbances attributed to warmer temperatures and earlier spring snowmelts [5], with a 60% increase in large fires from the mid-1980s in the Northern Rocky Mountains and significantly greater extents of severe fires in the Southern Rocky Mountains [6].Furthermore, unprecedented mountain pine beetle (MPB; Dendroctonous ponderosae) outbreaks in Western North America killed >10 Mha of trees over the last 2-3 decades [7,8], with major impacts on water resources [9,10], nutrient cycling [11], and forest carbon storage [12].These disturbance-caused changes require consideration when assessing global change risks to the stability of forest carbon [13,14].Landscape history, disturbance severity, and disturbance type represent interactive factors that reshape the composition, structure, and carbon budget of forests following disturbance events [15].Topography also plays an integral role in the distribution of water, nutrients, and soil carbon in forested systems [16]; it is therefore likely that hillslope-scale controls would influence soil biogeochemistry and microbial activity in disturbed systems.
Wildfire disturbance in forested systems alters vegetation communities and soil biogeochemical cycling as a function of fire duration, frequency, and intensity [17]-factors that are directly or indirectly shaped by topography [18][19][20].Hillslope-scale topography exerts important control on fire severity and burn frequency in Western U.S. landscapes: more frequent, low intensity wildfires occur on drier south-facing slopes in moderate drought years compared to densely vegetated north-facing slopes that experience fewer yet more severe fires during extreme drought years [20].Organic and mineral surface soils often display similar post-fire changes in geochemistry.These include a reduction in soil organic carbon (SOC) and total nitrogen (N) through the partial to complete removal of overlying organic matter [17,21,22]; the production of inorganic nitrogen including ammonium (NH 4 + ), a product of combustion, and nitrate (NO 3 − ), a subsequent product of nitrification, that occurs in the days to months following a wildfire [23,24]; and the concentration of non-volatized base cations (e.g., Mg 2+ , Ca 2+ , K + ) into ash lost via post-fire leaching and surface runoff [25,26].The extent of post-fire changes in soil geochemistry reflects ecosystem type, burn severity, and fire frequency [27][28][29], all of which are influenced by topography [30].Forest diebacks driven by MPB outbreaks alter C and N cycling pathways to varying degrees, depending on percent tree mortality [8] and phase of infestation [9,10,31,32], among other factors.The changes in soil biogeochemistry following MPB infestation include greater soil moisture availability likely due to reduced evapotranspiration associated with high tree mortality [33]; a reduction in phosphate (PO 4 3 ) concentrations [34]; lower C/N ratios in the leaf litter compared to little to no change in soil C/N ratios [34]; greater aromaticity of dissolved organic carbon (DOC) [8] and elevated concentrations of NH 4 + in the soil environment given the cessation of nitrogen uptake and more litter available following tree mortality [11,35,36].Nitrate concentrations following MPB outbreaks vary, with the strongest NO 3 − responses in forests that receive substantial N deposition (e.g., Europe, Asia) compared to N-limited system in Western North America [37][38][39].Studies of terrain-driven changes in soil biogeochemistry following MPB outbreaks are more limited, yet terrain attributes, such as slope, elevation, and aspect, serve as important factors for modeling hillslope-scale water availability following MPB outbreaks [9,10].Microbial communities in forested systems demonstrate strong responses to disturbance-driven changes in soil geochemistry and substrate inputs in the months to years following wildfire [40,41] and MPB outbreaks [12,42].Wildfires heat surface soils from 60 to 400 • C [43,44] and change the physical, chemical, and biological attributes of soils [45].Fire directly impacts microbial function by reducing microbial biomass (via mortality) and changes microbial communities by filtering microbial taxa [40,46].Fire indirectly impacts microbial function through altered substrate stoichiometry and edaphic conditions [47,48].Over time, fire changes soil texture, bulk density, and hydrophobicity, leading to increased soil erosion and further depletion of soil SOM [23,49,50].In contrast, microbial responses to MPB disturbance may be more driven by changes to root and organic matter inputs into soil.Mycorrhizal fungi exhibit up to a 12-fold reduction in biomass within a year of MPB-associated forest dieback coupled with an increase in the prevalence of saprotrophic fungal communities [34].Higher bacterial/fungal biomass ratios occur in the absence of fungal-tree symbioses after tree mortality [34], which suggest that fungi and bacteria respond differently; bacterial biomass inputs remain potentially unchanged or increase in soils following MPB disturbance [12].
Forest disturbances alter rates of organic matter decomposition by microbes that introduce uncertainty to the fate of soil carbon and aboveground carbon stocks.Quantifying soil microbial response to disturbance will reduce this uncertainty in dynamic, topographically complex ecosystems that are expected to experience greater degrees of disturbance given drought, warmer temperatures, and the lower snowpack associated with climate change [51].Soil C, N, and phosphorus (P) biogeochemistry vary based on the time since disturbance and are contingent on moisture availability [32].The dynamic responses of microbial enzyme activities to linked and compounded disturbance interactions requires more consideration [52], particularly in forested mountain systems where complex terrain contributes to variability in microbial responses to wildfire and MPB outbreak events.
Here, we examine how topography and two disturbance agents, wildfire and mountain pine beetle, interactively regulate soil biogeochemical responses and potential exoenzyme activities in montane forests of Colorado.Our goal is to assess how microbial functions, as indicated by soil CO 2 efflux and potential exoenzyme activities, and the resulting soil biogeochemistry in Western forested landscapes, are changed when forests are impacted by insect outbreaks, fire, and both insects and fire.We sampled 56 plots subjected to beetle disturbance from 2007 to 2009 (>90% tree mortality) and differing wildfire severity in pine-dominated forests of Colorado.Topsoil biogeochemistry, potential exoenzyme activity, and microbial biomass were quantified within each plot.We examined terrain-induced variability in moisture availability using the SAGA (System for Automated Geoscientific Analyses) wetness index generated from LiDAR-derived digital terrain models (DTMs), aspect, slope, catchment area, and solar radiation.Our sampling design incorporated combinations of beetle disturbance and fire, which allowed tests of whether microbial responses to disturbance interactions might be additive, synergistic, or negative [52].We hypothesized that microbial activities would correlate with topography and available moisture in soils of undisturbed forests, and questioned whether this relationship would be robust to fire, insect outbreak, and the combination of both disturbances.

Materials and Methods
We selected three field areas in Larimer County, Colorado that spanned the High Park wildfire and previous beetle kill disturbance gradients (Figure 1): Stove Prairie Ranch (SPR), Buckhorn United Methodist Camp (UMC) and, Lory State Park (LSP).The sites are located within 15 km of each other and are situated at elevations of 2090-2435 m in the Ponderosa pine (P.ponderosae Douglas ex C.Lawson) dominated montane forest that composes about three quarters of the tree cover in the region.Douglas fir (P.menziesii (Mayr) Franco), quaking aspen (P.tremuloides Michx.) and lodgepole pine (P.contorta Douglas) are also present and make up an additional quarter of tree cover.The majority of precipitation falls as snow; the area receives on average 425 cm of snow with significant snowfall from October through May, peaking in April.The average annual rainfall total is 56 cm.The average air temperature is 7.8 • C (min: −7.8 • C; max: 22.2 • C) with peak monthly averages in August and September during our sampling period.All three field areas overlay the Wetmore-Boyle-Rock Outcrop Complex and consist of shallow, well to excessively drained, upland Alfisol and Mollisol soils.Soils in this complex are classified as frigid Lithic Haplustalfs and frigid shallow Aridic Argiustolls [53].Underlying bedrock units were determined using 1:24,000 geologic quadrangle maps from the USGS (United States Geological Survey) National Map Geologic Database.All bedrock types were Precambrian (Early Proterozoic) in age, ranging from metasedimentary rocks composed of knotted mica schist at SPR to knotted mica schist and quartzofeldspathic mica schist at UMC [54].The LSP site is underlain by Boulder Creek granodiorite and metavolcanic hornblende-plagioclase gneiss or schist interlayered with massive amphibolite [54].

Verification of Cause of Mortality
The High Park wildfire burned a total of 353 km 2 in June 2012, including all three field areas in our study.Despite the extensive burn perimeter, we identified areas that were unburned and burned at each study site.Prior to the wildfire, all three sites experienced patchy MPB infestation in 2007-2009.We verified the cause of mortality as follows.
To distinguish trees killed by fire from those killed by MPB, we used multi-proxy methods because, in some cases, the wildfire burned away external evidence of MPB infestation.Briefly, to locate sites where MPB and wildfire both occurred, we first created a disturbance progress map using United States Forest Service (USFS) MPB outbreak data and overlaid the High Park wildfire boundary [55].To ascertain the location of MPB-killed trees that were affected by the High Park wildfire, we pinpointed clusters of red phase pines before the wildfire using historical aerial images on Google Earth to compare with images after the wildfire to see if they were in the wildfire zone.To confirm MPB infestation in the field, we peeled back the tree bark to reveal beetle galleries and blue stain fungus in the wood.We cored all dead trees in the candidate plots measuring >10 cm in diameter with a Haglof 5.1 mm diameter increment borer (Haglof Company Group, Langsele, Sweden).To ensure consistency and to minimize the chance of capturing growth caused by anything other than climate, such as slope, two cores per tree were taken parallel to the slope at a constant height of 1.3 m [56].The outer ring was required to determine the year-of-death (YOD) of the tree, so care was taken to ensure that the outer ring was captured intact.Cores were labeled and stored in paper straws until analysis.

Verification of Cause of Mortality
The High Park wildfire burned a total of 353 km 2 in June 2012, including all three field areas in our study.Despite the extensive burn perimeter, we identified areas that were unburned and burned at each study site.Prior to the wildfire, all three sites experienced patchy MPB infestation in 2007-2009.We verified the cause of mortality as follows.
To distinguish trees killed by fire from those killed by MPB, we used multi-proxy methods because, in some cases, the wildfire burned away external evidence of MPB infestation.Briefly, to locate sites where MPB and wildfire both occurred, we first created a disturbance progress map using United States Forest Service (USFS) MPB outbreak data and overlaid the High Park wildfire boundary [55].To ascertain the location of MPB-killed trees that were affected by the High Park wildfire, we pinpointed clusters of red phase pines before the wildfire using historical aerial images on Google Earth to compare with images after the wildfire to see if they were in the wildfire zone.To confirm MPB infestation in the field, we peeled back the tree bark to reveal beetle galleries and blue stain fungus in the wood.We cored all dead trees in the candidate plots measuring >10 cm in diameter with a Haglof 5.1 mm diameter increment borer (Haglof Company Group, Langsele, Sweden).To ensure consistency and to minimize the chance of capturing growth caused by anything other than climate, such as slope, two cores per tree were taken parallel to the slope at a constant height of 1.3 m [56].The outer ring was required to determine the year-of-death (YOD) of the tree, so care was taken to ensure that the outer ring was captured intact.Cores were labeled and stored in paper straws until analysis.
The year of death of trees killed by the MPB before the High Park wildfire in 2012 was verified using standard dendrochronological procedures [56,57].We measured the width of each growth ring using the LINTAB tree-ring measurement station (LINTAB TM 5, Rinntech Inc., Heidelberg, Germany).We compiled the live tree data into a site master with known outside ring dates and crossdated the dead trees.We performed Pearson's correlation coefficient analysis (p < 0.001, CI = 0.3281) of ring widths in the samples and the master chronology using COFECHA 2.111 (The Laboratory of Tree Ring Research, Tucson, AZ, USA [58]).Crossdating was verified against the standardized ring-width series of the master chronology TSAP (TSAP-Win [59]; Rinntech Inc., Heidelberg, Germany [60,61]).

Establishment of Study Plots
During summer 2013, we established a total of fifty-six sampling plots across the three field sites: twenty at SPR, twenty at LSP, and sixteen at UMC (Figure 1).At each of the three field sites, we located sampling plots in four disturbance types: Burned, Beetle Present (BP, n = 20); Burned, Beetle Absent (BA, n = 20); Unburned, Beetle Present (UP; n = 8), and Unburned, Beetle Absent (UA, n = 8).In any ecological sampling scheme, spatial autocorrelation is a concern.To reduce the possibility that the sampling strategy could bias our conclusions, we established a categorical study by site by sampling disturbance plots across each of the three locations.We also indicate location (SPR, LSP, UMC) in the figures reporting the results so that any potential bias can be assessed visually.Each plot was established by placing a flag in the center of a cluster of 3-5 pines that were affected by the same disturbance type and which varied in size between 10 and 20 cm in diameter.For each plot we identified a minimum of three mature pine focal trees of each disturbance type in isolation (the edge of the plot was at least 1 m distance from any live trees or other significant vegetation or disturbance type in all directions).Working with the heterogeneity of the landscape and disturbance patterns, no plots were closer than 2 m apart and most were more than 10 m apart.Only trees above 10 cm in diameter were considered.We designated the center of the plot as the point equidistant from the focal trees.GPS coordinates, elevation for the center of the plot were recorded.
Each plot was approximately 2 m 2 in size.Locations of trees, shrubs, stumps, logs, and PVC collars (see below) were mapped within each plot to facilitate subsampling.Briefly, we designated North-South and East-West transect lines and marked these points at the plot boundary with labeled flags to create quadrants.For each tree within the plots, we measured the diameter at breast height (dbh; Forestry Suppliers Metric Fabric Diameter Tape Model 283D/5M) and heights (Suunto PM-5/360 PC Clinometer).Trees were labelled with a unique designation and disturbance status.Shrub cover was mapped but not identified to species.Finally, percent cover was estimated by two people and averaged for each quadrant.

Soil Efflux Measurements
To establish the overall metabolic activity of soils that was coincident in time with soil biogeochemical and enzyme measurements, we measured the soil CO 2 efflux from four locations within each 2 m 2 plot within 24 h of collecting samples for laboratory analysis (described below).Measurements were made prior to soil sampling for laboratory analysis between the window of 10 a.m. and 3 p.m.This is a period of stable soil efflux measurements for similar soils, with vegetation, climate and site histories in neighboring counties in CO where we have previously published studies [31,32,62].
Four PVC soil collars were installed per plot, selecting a random location in each quadrant.The location of the soil collar was at least 30 cm within the plot boundary, on relatively level ground that was free of vegetation.After 48 h, soil efflux was measured in the quadrant of each plot, using LI-COR LI-6400 infa-red gas analyzers (IRGA) equipped with 6400-09 soil CO 2 flux chambers (LI-COR Incorporated, Lincoln, NE, USA).The IRGA uses the characteristic infrared absorption of gases to measure the concentration of CO 2 building up in the chamber over time and calculates a soil efflux rate.Soil temperature in each quadrant was measured concurrently with soil efflux using the LI-COR LI-6400 thermocouple attachment, which measures temperature integrated over a depth of 0-15 cm.
IRGAs are calibrated at LI-COR every 2 two years, with the last factory calibrations in 2012.Prior to the field season IRGA zeros and spans were checked in the laboratory using tank standards.
Prior to each set of measurements, the IRGA was allowed to warm up for at least 10 min.Standard checks were made to verify pressure, water and CO 2 readings were within normal bounds.Prior to each measurement, the chamber cap was removed and the chamber was laid adjacent to the soil collar to allow the user to set the average (ambient) [CO 2 ] for the instrument.The soil thermocouple probe was inserted fully into the soil adjacent to the soil collar prior to the initiation of the measurement cycle.The collar depth was measured and the instrument was set such that the insertion depth was zero.In cases when it was not possible for the chamber to be flush with the surface, an insertion depth was calculated and entered into the IRGA so that the chamber volume was correctly estimated.Three efflux measurements were made at each collar using automated sampling mode for this instrument.Briefly, air was pumped through a soda lime scrub to bring the chamber CO 2 concentration to 10 ppm below the ambient [CO 2 ] and the rate of CO 2 build up was used to calculate the soil efflux up to 10 ppm above ambient [CO 2 ].In cases where variation in efflux at a single collar exceeded 10% during a measurement cycle, data were discarded and the cycle was repeated.Soil temperature measurements were recorded simultaneously by the LI6400.

Soil Sampling, Moisture, and Biogeochemical Analysis
Within each plot, several soil samples were collected to estimate a series of different properties; nitrogen, carbon and microbial biomass and more general characteristics of the plot (pH, soil moisture, bulk density).During soil collection, litter depth and the depth of the organic layer and burn layer were evaluated visually using a plastic ruler (mm).All soil sampling was conducted while wearing nitrile gloves.Soils were homogenized, stored on ice, and returned to the laboratory to be processed within six hours.Soil was passed through a 2-mm sieve to remove debris and coarse plant material.
To calculate bulk density, we cored two randomized locations within each plot using a standard soil corer of known radius (0.7106 cm).The core length was measured and returned to the laboratory, dried at 110 • C for 48 h, and weighed to determine mass.
Approximately 200 g of burned, organic, and mineral soil was collected from the top 5 cm of two randomized locations within the plot for pH and texture.To measure pH, we used an Oakton pH Testr 3 digital pH meter (Oakton instruments, Vernon Hills, IL, USA) on a 1:1 solution of 30 g soil to 30 mL deionized water.We weighed 5 g of homogenized soil, and dried the sample in a drying oven at 60 • C for 48 h before reweighing the dried soil to calculate percent soil moisture.
To determine soil biogeochemical metrics and microbial biomass, we collected approximately 10 g soil from the top 5 cm from four randomized locations.We extracted for dissolved organic carbon (DOC), dissolved organic nitrogen (DON), and dissolved inorganic nitrogen (DIN) following established procedures [31,63].Briefly, 25 mL of 0.5 M potassium sulfate (K 2 SO 4 ) was added to 5 g of homogenized sample, agitated on an orbital shaker for one hour at 120 rpm and then vacuum-filtered through glass fiber filters.Concurrently, an additional 5 g of each sample was fumigated with 2 mL of ethanol-free chloroform in Erlenmeyer flask for 24 h, vented for one hour and then extracted as above to determine microbial biomass carbon and nitrogen pools.All extracts were frozen at −20 • C until analysis.
DOC and total dissolved nitrogen (TDN) in the K 2 SO 4 extracts and the fumigated soil K 2 SO 4 extracts were quantified using the non-purgable-organic-C protocol on a Shimadzu TOC analyzer (TOC 5000) equipped with a total dissolved nitrogen module (Shimadzu Scientific Instruments, Inc., Columbia, MD, USA).Soil microbial biomass carbon (MBC) and nitrogen (MBN) were estimated by subtracting concentrations of DOC and DON in the un-fumigated samples as well as soil free controls from the fumigated samples [64,65].The efficiency factors for microbial biomass carbon (kEC = 0.45 [66]) and microbial biomass nitrogen (kEN = 0.54 [64]) were used to calculate the respective biomass as the difference between fumigated and non-fumigated samples.

Potential Exoenzyme Activity
We measured potential activities of seven hydrolytic enzymes: β-D-cellubiosidase (CB), α-Glucosidase (AG), β-Glucosidase (BG), and β-Xylosidase (XYL), which hydrolyze carbon-rich substrates; leucine aminopeptidase (LAP) and N-acetyl-β-Glucosaminidase (NAG), which hydrolyze nitrogen rich substrates, and Phosphatase (PHOS), which hydrolyzes phosphorus rich substrates.Potential exoenzyme activities were measured using a modified fluorimetric deep-well microplate technique [67].Soil slurries were prepared with 2.75 g of soil that was stored at 4 • C for no longer than 3 weeks and 91 mL of 50 mM Tris Buffer titrated to match soil pH.Assays were run with two internal standards: dilution series (0-100 µM) of 4-methylumbelliferone (MUB) or 7-amino-4methylcoumarin (MUC) for LAP, each mixed with soil homogenate [68].Standards for standard curves and assays were incubated at 25 • C for 1.25 and 1.5 h; 100 µL of 200 µM fluorimetric substrate was added to 900 µL of each soil slurry.Fluorescence was measured on Synergy™ 4 Multi-Mode microplate reader with an excitation wavelength of 365 nm and an emission wavelength of 450 nm.Total enzyme activities were calculated from fluorescence values as the rate of substrate converted in nmol h −1 g −1 soil [69].We calculated specific enzyme activities as activity per unit microbial biomass C.

Terrain Analysis
Terrain attributes were calculated for each plot in the study areas using a 1-m LiDAR-derived digital surface model produced by the National Ecological Observatory Network.The 1-m digital surface model was resampled using the cubic convolution technique to eliminate microtopographic effects on the landscape and preprocessed to fill sinks [70].Raster datasets were generated in SAGA GIS (v.2.1.4,2014, System for Automated Geoscientific Analyses [71]) for aspect, local slope, solar radiation, catchment area, topographic wetness index, and stream power index.Aspect and local slope were produced using the Slope, Aspect, and Curvature Module whereas total, direct, and diffuse solar radiation data sets were determined from the Potential Incoming Solar Radiation Module (v.2.1.4,2014, System for Automated Geoscientific Analyses [71]), Cary, NC, USA.Catchment area was calculated with the triangular multiple flow direction algorithm which supports multidirectional downslope flow [72].The topographic wetness index represents putatively long-term soil moisture conditions across the study sites as determined by the likelihood that water will accumulate in or flow out of a raster cell.We calculated the TWI in Raster Calculator (ArcMap, v.10.2, 2014, Environmental Systems Research Institute, Redlands, CA, USA) using the SAGA-derived catchment area and slope raster imagery [73,74]: where A is the specific catchment area and β is the local slope.High TWIs occur in concave, water-gathering landscapes compared to low TWIs which are more common for steep, convex, water-shedding hillslopes.
The stream power index (SPI) represents potential flow erosion at a given point on the landscape as a function of catchment area and slope [75]: Higher stream power indices suggest a greater potential for erosion due to greater contributing catchment areas, higher slope gradients, and the resultant increase in water flow velocity [76].

Statistical Analyses
All statistical tests were performed on natural-log transformed data in JMP Pro (Version 13.0.0,SAS Institute, Cary, NC, USA).We employed oneway ANOVA with Tukey-Kramer (HSD) posthoc test to compare means of terrain and geochemical attributes among the 4 disturbance groupings at a 95% confidence interval.Linear regressions were used to compare potential exoenzyme activities with dependent variables of terrain and geochemistry (e.g., SAGA wetness index; pH; geochemical ratios).We compared responses of exoenzyme activities and specific exoenzyme activities (activity per unit biomass C).There were no differences in responses and therefore only exoenzyme activities are discussed.

Biogeochemical and Topographic Relationships by Disturbance Type
Wildfire disturbance had the largest impact on soil pH (Figure 2e), DOC/TDN (Figure 2f), N geochemistry (Figure A1), and potential exoenzyme activities, with less identifiable differences attributed to the presence or absence of MPB (Table 1).Mean soil pH was significantly higher in burned plots compared to unburned plots regardless of beetle presence.Greater NH 4 + , TDN, and NO 3 − concentrations occurred in the burned sites versus the unburned plots irrespective of MPB (Figure A1a-c), whereas biomass N was highest in unburned locations (Figure A1d).The unburned, beetle present plots exhibited the lowest NH 4 + concentrations of the four sample groupings (Figure A1a).Soil C biogeochemistry varied little among the disturbance types, except for significantly higher DOC concentrations in burned, beetle present plots compared to beetle absent samples in both burned and unburned landscapes (Figure A2).Phosphatase (PHOS) was the only one of seven exoenzymes to exhibit significant differences by disturbance type, with lower potential activities in soils from burned, beetle absent plots compared to unburned beetle absent and beetle present plots.
Forests 2018, 9, x FOR PEER REVIEW 8 of 23 per unit biomass C).There were no differences in responses and therefore only exoenzyme activities are discussed.

Biogeochemical and Topographic Relationships by Disturbance Type
Wildfire disturbance had the largest impact on soil pH (Figure 2e), DOC/TDN (Figure 2f), N geochemistry (Figure A1), and potential exoenzyme activities, with less identifiable differences attributed to the presence or absence of MPB (Table 1).Mean soil pH was significantly higher in burned plots compared to unburned plots regardless of beetle presence.Greater NH4 + , TDN, and NO3 − concentrations occurred in the burned sites versus the unburned plots irrespective of MPB (Figure A1a-c), whereas biomass N was highest in unburned locations (Figure A1d).The unburned, beetle present plots exhibited the lowest NH4 + concentrations of the four sample groupings (Figure A1a).Soil C biogeochemistry varied little among the disturbance types, except for significantly higher DOC concentrations in burned, beetle present plots compared to beetle absent samples in both burned and unburned landscapes (Figure A2).Phosphatase (PHOS) was the only one of seven exoenzymes to exhibit significant differences by disturbance type, with lower potential activities in soils from burned, beetle absent plots compared to unburned beetle absent and beetle present plots.Aspect, or the occurrence of sites on north or south-facing hillslopes, played a role in the distribution of burned and unburned sites in our study, whereas localized surface soil moisture varied, in part, with the presence or absence of MPB.We observed that the burned plots occurred consistently on the north to northwest-facing hillslopes that receive lower influxes of average total solar radiation versus unburned landscapes with southern aspects and higher solar radiation inputs (Figure 2a,b and Figure A2a).Soil temperature did not significantly vary among the disturbance types (Figure 2b), but a one-time measure of soil moisture documented more moisture in surface soils from beetle present plots compared to beetle absent sites (Figure 2d).Specifically, plots with beetle-kill contained significantly higher surface soil moisture percentages in burned and unburned plots compared to plots where beetles were absent (Figure 2d).The SAGA wetness index did not differ among the four disturbance types as observed with aspect (Figure 2a) and total solar radiation (Figure A2b), leading to its selection as a primary terrain attribute for examining topographic control on soil CO 2 efflux, geochemistry, and potential exoenzyme activity in our study.

Topographic Controls on Soil CO 2 Efflux and Potential Exoenzyme Activity
The SAGA wetness index was identified as a significant hillslope-level control in unburned, beetle absent soils that accounted for more than 50% of variability in soil CO 2 efflux and potential exoenzyme activities compared to under 10% of variability explained by this terrain attribute in burned, beetle absent soils.The SAGA wetness index for unburned, beetle absent soils exhibited strong positive relationships with soil CO 2 efflux (Figure 3a), PHOS activity (Figure 3b), NAG activity (Figure 3c), and LAP activity (Figure 3d).The SAGA wetness index did not significantly correlate with soil efflux or potential exoenzyme activities in the burned, beetle present plots, with the exception of LAP (Figure 3d).Soil C mineralizing enzyme activities and microbial biomass C did not coincide with wetness indices in either plot type (Figure 3e,f).Wetness indices from unburned, beetle present plots showed no correspondence with soil efflux or exoenzyme activities, aside for LAP, despite strong relationships among SAGA wetness index, pH, and biomass C (Table 1).

Potential Exoenzyme Activities in Fire-Disturbed Landscapes
Soil pH and DOC/TDN emerged as significant variables in burned, beetle absent plots that accounted for 30 to 80% of variability in C, N, and P potential exoenzyme activities (Figures 4a-d and 5a-d), suggesting that soil pH and geochemistry strongly influence exoenzyme activities in burned systems.Soils from burned, beetle absent plots showed strong, negative associations between pH and Sum of C activities (Figure 4a), PHOS activity (Figure 4b), and NAG activity (Figure 4c).Conversely, positive relationships occurred between soil pH and LAP activity in burned and unburned beetle absent sites (Figure 4d).DOC/TDN ratios presented positive strong relationships with Sum of C (Figure 5a), PHOS (Figure 5b), and NAG activities (Figure 5c) in burned, beetle absent plots compared to no observed relationships between DOC/TDN and exoenzyme activities in unburned, beetle absent sites (Figure 5).LAP activity did not significantly correspond to DOC/TDN in either disturbance type (Figure 5d).

Potential Exoenzyme Activities in Fire-Disturbed Landscapes
Soil pH and DOC/TDN emerged as significant variables in burned, beetle absent plots that accounted for 30 to 80% of variability in C, N, and P potential exoenzyme activities (Figures 4a-d and  5a-d), suggesting that soil pH and geochemistry strongly influence exoenzyme activities in burned systems.Soils from burned, beetle absent plots showed strong, negative associations between pH and Sum of C activities (Figure 4a), PHOS activity (Figure 4b), and NAG activity (Figure 4c).Conversely, positive relationships occurred between soil pH and LAP activity in burned and unburned beetle absent sites (Figure 4d).DOC/TDN ratios presented positive strong relationships with Sum of C (Figure 5a), PHOS (Figure 5b), and NAG activities (Figure 5c) in burned, beetle absent plots compared to no observed relationships between DOC/TDN and exoenzyme activities in unburned, beetle absent sites (Figure 5).LAP activity did not significantly correspond to DOC/TDN in either disturbance type (Figure 5d).

Discussion
Soil moisture is an important control of microbial activity in undisturbed forests.Soil CO 2 efflux along with N-and P-degrading enzymes in undisturbed forests indicated greater microbial activities in sites with higher SAGA wetness indices (Figure 3).Importantly, topography and soil moisture correlate with NPP, soil mineral nutrient, and organic matter content [77,78], with aspect playing an important role.Our unburned study sites are situated on south-facing hillslopes in low to moderate elevation forests (2090-2435 m) where moisture would be of importance.Earlier work supports this interpretation.In moderate elevation (1840-2421 m) lodgepole-dominated forests of the northern Rocky Mountains, Montana, greater seasonal CO 2 efflux in upland forest sites corresponded to soils that exhibited higher upslope accumulated areas, a proxy for lateral flow of soil water [79].Greater fine root densities and lower C/N ratios also occurred in these wetter landscape positions that the authors attributed to greater soil CO 2 efflux and higher litter decomposition rates, respectively [79].Soil moisture was also identified as a limiting factor for soil respiration in low to moderate elevation pine forests of Boulder Creek, Colorado, where soil CO 2 efflux corresponded to topographic wetness indices along a low to high forested elevation gradient [80].Thus, the SAGA wetness index serves as a useful geospatial proxy for microbial activity in undisturbed forests.
In our study, conducted one year after a fire, we observed lower microbial biomass C and N in certain burned plots, and greater variability compared to undisturbed plots (Figures 2f and A1d).The plots with the lowest microbial biomass were likely where the fire was most severe.This interpretation is supported by the corresponding changes in biomass, soil pH, soil moisture, and DOC/TDN that we observed (Figure 2).Both the dominant species and the quantity of ash added to the plot after the fire would influence soil pH [81][82][83], variation in pH here is a reasonable proxy for fire intensity (Figure A3).At higher pH, organic compounds become more soluble and therefore more prone to leaching after wildfires in western forests [84].As expected [23,40,49,50] pH increased (Figure 2e) and DOC declined relative to TDN (Figure 2f) in burned plots.
Fire influences microbial activity by altering soil microclimate conditions (soil temperature, soil moisture, solar insolation); soil chemistry and physical structure (pH, nutrient availability, soil structure, soil texture); and increases hydrophobicity in soils [45].Overall, plots showing the greatest physical and biogeochemical impacts of fire displayed different microbial function from undisturbed plots (Figures 2-4).We found that, in contrast to undisturbed forest, soil CO 2 efflux and exoenzyme activities in burned plots did not correlate with the SAGA wetness index (Table 1; Figure 3b-d).The burned sites occurred on wetter, north-facing hillslopes, suggesting that microbial communities may not have experienced the same moisture limitations that result on drier south-facing sites in the unburned environments (Figure 2a).Alterations to the soil microclimate, primarily the fire-induced loss of vegetation and reduced thickness of organic soil horizons, may also explain the variability in soil temperature seen in burned sites given the potential for an increased degree of solar insolation (Figure 2b).Soil carbon and nitrogen pools change dramatically post-fire both in terms of availability and chemical structure, and microbial enzyme activity has been shown to respond to these changes to soil C/N availability both rapidly and persistently [41].The lack of correspondence between soil moisture and microbial activities in the burned sites of our study could therefore indicate stronger control exerted by shifting nutrient limitation and microbial responses to biogeochemical properties in the burned sites (Table 1) or the influence of aspect on soil moisture conditions.
In our study, soil pH varies with fire intensity and emerges as an important explanatory variable for all exoenzyme activities, except LAP, in burned sites (Table 1; Figure 4).The pH of unburned soils was more acidic due to organic matter inputs from pine trees, whereas soils from burned forests were more alkaline but also had higher range in pH likely due to differences in burn severity (Figures 2e and 4).Soil pH represents a well-known descriptor of microbial community composition at landscape and biome scales [85][86][87].Field observations of ash depth, a proxy for burn intensity was positively correlated with pH (r = 0.7, p < 0.05) and patterns of litter depth were also consistent with this interpretation (Figure A3).Associated with microbial N acquisition, LAP was positively correlated with pH in unburned sites and remained positively correlated in burned sites (Figure 4d).A study comparing effects of fire in Ponderosa pine and mixed conifer forests found that, regardless of initial soil conditions, soil bacteria converged on similar community composition after severe fire [46].Changes in substrate availability likely play a role in this shift; SOC is more soluble in basic soils [88,89].In our study, greater fire intensity resulted in lower DOC/TDN ratios (Figure 5), which is consistent with reduced organic matter.All exoenzyme activities except LAP were positively correlated with DOC/TDN ratios (Table 1; Figure 5).Severe fires cause significant losses of SOM [90] and suggest potential greater C limitation for microbial communities in severe burn soils relative to greater N-limitation for microbial communities in less intense burns [40].
The effects of tree mortality from MPB on soil microbial function depend on the timing of beetle induced tree mortality.In contrast to the burned plots, there is no change in pH in MPB plots (Figure A3), which is consistent with several previous studies [32,36,91,92].Unlike fire, which increases soil pH and potential for rapid leaching of soluble organic material, beetle-induced tree mortality reduces C and nutrient inputs to the soil from root exudates and litterfall in the short term, which is followed by increased C inputs from leaf, branch and root litter necromass [31,32].Mountain Pine Beetles kill trees in association with fungi that colonize the tree when beetles bore into the phloem which disrupts transport of sugars belowground, effectively girdling the tree [90].After tree mortality, soil extractable dissolved C and N, inorganic P and microbial biomass (C and N) decline initially but tend to increase between 4 and 5 years after tree mortality [12,32].This loss of exudates causes increased N mineralization [32,63] and lower C/N ratios and SOC [12,31].In this study, we see few detectable changes in the biogeochemical quantities measured (Figures 2 and A1).This is consistent with the variability we would expect, given that beetle kill occurred from 2007 to 2009 (6 to 4 years before the study); spanning the period where previous studies indicate that losses of organic inputs from live trees is transitioning to gains in organic inputs from litter decomposition [12,32].Nevertheless, beetle disturbance still influences microbial functional controls; while PHOS, LAP and NAG are all correlated with moisture in the undisturbed plots, only LAP is statistically related to moisture in the beetle killed plots (Table 1).This might be explained by the tendency to have higher soil moisture in plots impacted by beetles (Figure 2).The regression analysis also indicates that NAG, involved with N-degradation, is correlated with the DOC/TDN ratio in beetle killed plots which contrasts with undisturbed soils.This could be because of the tendency to find higher TDN and a lower NO 3 − in the beetle killed plots compared to the undisturbed forest, indicating greater substrate availability (chitin and organic N sources) in the undisturbed forest.
Prior beetle disturbance may have modified the effect of fire on the controls of microbial activity.Soil moisture was higher in plots where beetle kill preceded fire compared to plots that were recently burned (Figure 2d).This interaction between beetle kill and fire resulted in intermediate responses of pH and DOC/TDN ratios compared to no disturbance or beetle-only plots and the plots that were burned (Figure 2e,f).We have no evidence that burn severity was lower in these plots; it is more likely that that the higher soil moisture and thicker needle layer [31,32] provided some buffering from the effects of fire (Figure A3).Exoenzyme activities were not strongly regulated by moisture in these plots, and the relationships between potential exoenzyme activities and DOC/TDN were similar to plots that experienced only fire, although soil CO 2 efflux did show a stronger response to DOC/TDN in burned sites with prior beetle disturbance compared to burned, beetle absent sites (Table 1).Also, due to the narrower range in pH, the strong relationships between exoenzyme activities and pH were absent in beetle-killed plots (Table 1).
Comparing microbial function after the fire and MPB disturbances independently and jointly allows us to evaluate whether these disturbances act in a compound manner; that the response of the second consecutive disturbance is modified by the first [52].We have tentative support for the compound disturbance hypothesis.Both disturbances contribute enough variation to the environmental drivers of microbial function to reduce the importance of topography and its correlates as determinants of the rates of CO 2 efflux and potential enzyme activities (Table 1).While beetle mortality alone may be complicated by the timing of when that mortality occurred [31,93,94], the environmental correlates of microbial function after fire alone (pH and DOC/TDN) are not the same as the correlates after compounded beetles and fire disturbance (DOC/TDN).In our previous work [32] and this current study, we have found limited effects of beetle mortality on soil pH itself, though we find some evidence that beetle kill apparently buffers soil pH (Figures 2e and A3).How these compound effects impact later recruitment and longer-term ecosystem function [95,96] requires longer monitoring.However, it is clear that fire does not completely erase the effects of beetle kill on microbial function and could result in longer-term changes to these ecosystems.Our analysis suggests that these disturbances are compound, rather than independent.
Montane forests recovering from historical logging and managed for some degree of fire suppression are responsible for most of the carbon uptake and storage in the Western U.S. [97,98].Starting in the early 1990s, catastrophic levels of tree mortality have been observed across millions of hectares of bark beetle-infested forest in North America [99][100][101].Forest disturbance is a fundamental driver of terrestrial carbon cycle dynamics [102], but the long-term effects of different disturbance types are poorly understood and characterized in today's land surface models [103].Attempts to scale the impacts of different types of tree mortality remain challenging [104].Heterotrophic function in land surface models is broadly controlled by bulk live and dead biomass, temperature, moisture, and in some cases N, and is strongly dependent on plant productivity [93,105].Our study suggests that the environmental factors limiting microbial function shift depending on whether fire or pathogeninduced mortality is the disturbance agent.Our results are consistent with studies limited to single disturbances [31,32,40,50,92,94].Soil biogeochemical processes can be represented through a variety of modeling approaches including a simple bulk substrate, moisture and Q10 response [105], a more complex treatment of soil substrate suitability [62,106], or even trait-based approaches [107,108].If the type of disturbance influences biogeochemical cycling by changing availability of organic substrates and the structure or function of the microbial community, then a more complete representation of microbial processes is required.

Conclusions
Terrain attributes, such as the SAGA wetness index, function as useful predictors for soil CO 2 efflux and microbial activities when upscaling point analyses to watershed scales of study; however, microclimate conditions, seasonality, precipitation type, and disturbance history complicate topographic relationships with microbial activities in complex terrain.In this study, we show that the SAGA wetness index predicts exoenzyme activities in undisturbed landscapes, but that terrain-driven relationships with soil CO 2 efflux and exoenzyme activities are altered by disturbance events (e.g., fire, beetle kill) that are predicted to intensify with climate change [5,8].Soil moisture was not a strong predictor of exoenzyme activities after both disturbance types studied here; only potential LAP activity, involved in N-degradation, was significantly correlated with the SAGA wetness index following beetle infestation (Figure 3, Table 1).Rather, pH and indicators of microbial substrate availability emerged as strong correlates of activity (Table 1; Figures 3-5).Our research demonstrates the effectiveness of using an integrated approach that combines measures of potential exoenzyme activities, soil biogeochemistry, and terrain attributes for examining the dynamic nature of heterotrophic activity and soil respiration in forest ecosystems following disturbance events.The suite of enzyme and biogeochemical assays used here enabled us to see that the breakdown of environmental drivers for soil CO 2 efflux was most likely caused by changing terrain attributes and physical and biogeochemical factors that limit the activity of the microbial community that would not be provided by soil efflux alone.Addressing soil decomposition processes in complex terrain requires that we assess the topographic, biological and biogeochemical variables that drive microbial activities, particularly in western U.S. forests where the compounded effects of wildfire and insect outbreaks are expected to increase in intensity and frequency because of global climate change.

Figure 1 .
Figure 1.A map of the study area in Larimer County, CO, USA.The area burned by the 2012 wildfire is shaded red.Fifty-six study plots were established across the three study locations Stove Prairie Ranch (SPR), Buckhorn United Methodist Camp (UMC) and, Lory State Park (LSP) indicated with black circles.Inset map: the open black square contains the area of Colorado depicted by the main map.

Figure 1 .
Figure 1.A map of the study area in Larimer County, CO, USA.The area burned by the 2012 wildfire is shaded red.Fifty-six study plots were established across the three study locations Stove Prairie Ranch (SPR), Buckhorn United Methodist Camp (UMC) and, Lory State Park (LSP) indicated with black circles.Inset map: the open black square contains the area of Colorado depicted by the main map.

Figure 4 .
Figure 4. Burned, Beetle Absent (BA) and Unburned, Beetle Absent (UA) plots for (a) Sum of C activities; (b) PHOS activity; (c) NAG activity, and (d) LAP activity as a function of pH.Symbology in the figure represents data points for the following field sites: UMC (triangles), LSP (circles), and SPR (squares).

Figure 5 .
Figure 5. Burned, Beetle Absent and Unburned, Beetle Absent plots for (a) Sum of C EEA (exoenzyme activity); (b) PHOS activity; (c) NAG activity; and (d) LAP activity as a function of DOC/TDN ratios.Symbology in the figure represents data points for the following field sites: UMC (triangles), LSP (circles), and SPR (squares).

Figure 4 .
Figure 4. Burned, Beetle Absent (BA) and Unburned, Beetle Absent (UA) plots for (a) Sum of C activities; (b) PHOS activity; (c) NAG activity, and (d) LAP activity as a function of pH.Symbology in the figure represents data points for the following field sites: UMC (triangles), LSP (circles), and SPR (squares).

Figure 4 .
Figure 4. Burned, Beetle Absent (BA) and Unburned, Beetle Absent (UA) plots for (a) Sum of C activities; (b) PHOS activity; (c) NAG activity, and (d) LAP activity as a function of pH.Symbology in the figure represents data points for the following field sites: UMC (triangles), LSP (circles), and SPR (squares).

Figure 5 .
Figure 5. Burned, Beetle Absent and Unburned, Beetle Absent plots for (a) Sum of C EEA (exoenzyme activity); (b) PHOS activity; (c) NAG activity; and (d) LAP activity as a function of DOC/TDN ratios.Symbology in the figure represents data points for the following field sites: UMC (triangles), LSP (circles), and SPR (squares).

Figure 5 .
Figure 5. Burned, Beetle Absent and Unburned, Beetle Absent plots for (a) Sum of C EEA (exoenzyme activity); (b) PHOS activity; (c) NAG activity; and (d) LAP activity as a function of DOC/TDN ratios.Symbology in the figure represents data points for the following field sites: UMC (triangles), LSP (circles), and SPR (squares).

Figure A3 .
Figure A3.The relationship between volumetric soil moisture (%), litter depth (cm), soil pH and ash depth (cm) in plots with and without beetle mortality prior to the High Park fire in 2012.Symbology in the figure represents data points for the following field sites: UMC (triangles), LSP (circles), and SPR (squares).

Figure A3 .
Figure A3.The relationship between volumetric soil moisture (%), litter depth (cm), soil pH and ash depth (cm) in plots with and without beetle mortality prior to the High Park fire in 2012.Symbology in the figure represents data points for the following field sites: UMC (triangles), LSP (circles), and SPR (squares).

Table 1 .
Reported r 2 and p-values for environmental control and microbial function relationships from undisturbed and disturbed sites in Larimer County, Colorado.