Next Article in Journal
Development of Somatic Embryo Maturation and Growing Techniques of Norway Spruce Emblings towards Large-Scale Field Testing
Next Article in Special Issue
Modest Effects of Host on the Cold Hardiness of Emerald Ash Borer
Previous Article in Journal
Effects of Vegetation Management on Wood Properties and Plant Water Relations of Four Conifer Species in the Pacific Northwest of the USA
Previous Article in Special Issue
Biotic and Abiotic Drivers of Sap Flux in Mature Green Ash Trees (Fraxinus pennsylvanica) Experiencing Varying Levels of Emerald Ash Borer (Agrilus planipennis) Infestation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Response of Black Ash Wetland Gaseous Soil Carbon Fluxes to a Simulated Emerald Ash Borer Infestation

1
Department of Earth, Environmental and Geographical Sciences, Northern Michigan University, Marquette, MI 49855, USA
2
School of Forest Resources & Environmental Science, Michigan Technological University, Houghton, MI 49931, USA
3
Hawkesbury Institute for the Environment, Western Sydney University, Penrith, NSW 2753, Australia
4
USDA Forest Service, Pacific Southwest Research Station, Arcata, CA 95521, USA
5
USDA Forest Service, Northern Research Station, Grand Rapids, MN 55744, USA
6
Department of Natural Resource Sciences, Thompson Rivers University, Kamloops, BC V2C 0C8, Canada
*
Author to whom correspondence should be addressed.
Forests 2018, 9(6), 324; https://doi.org/10.3390/f9060324
Submission received: 31 March 2018 / Revised: 9 May 2018 / Accepted: 23 May 2018 / Published: 4 June 2018
(This article belongs to the Special Issue Understanding and Managing Emerald Ash Borer Impacts on Ash Forests)

Abstract

:
The rapid and extensive expansion of emerald ash borer (EAB) in North America since 2002 may eliminate most existing ash stands, likely affecting critical ecosystem services associated with water and carbon cycling. To our knowledge, no studies have evaluated the coupled response of black ash (Fraxinus nigra Marsh.) wetland water tables, soil temperatures, and soil gas fluxes to an EAB infestation. Water table position, soil temperature, and soil CO2 and CH4 fluxes were monitored in nine depressional headwater black ash wetlands in northern Michigan. An EAB disturbance was simulated by girdling (girdle) or felling (ash-cut) all black ash trees with diameters greater than 2.5 cm within treated wetlands (n = 3 per treatment). Soil gas fluxes were sensitive to water table position, temperature, and disturbance. Soil CO2 fluxes were significantly higher, and high soil CH4 fluxes occurred more frequently in disturbed sites. Soil CH4 fluxes in ash-cut were marginally significantly higher than girdle during post-treatment, yet both were similar to control sites. The strong connection between depressional black ash wetland study sites and groundwater likely buffered the magnitude of disturbance-related impact on water tables and carbon cycling.

1. Introduction

Disturbance events are known to impact the natural function of wetland ecosystems by altering nutrient, carbon, energy, and hydrologic fluxes [1,2]. Recently, emerald ash borer (EAB) (Coleoptera: Buprestidae, Agrilus planipennis Fairmaire) was introduced as an invasive pest to ash (Fraxinus spp.) forests in North America [3]. Hundreds of millions of ash trees were estimated to have been killed in North America within the decade following the initial detection of EAB in southeastern Michigan in July 2002 [4,5]. The expansion of EAB results in significant economic costs [6,7] and causes significant perturbation to forest ecosystems [5,8,9], and will likely continue to decimate existing ash stands throughout North America.
Forested wetlands sequester large quantities of carbon from the atmosphere [10], and their carbon storage capacity is largely a result of prolonged periods of elevated water table position that limits decomposition [11,12]. Forest and soil gaseous-carbon fluxes are sensitive to disturbance of overstory trees [13,14], and fluxes from wetland soils are also regulated by water table position and soil temperature [15,16]. Disturbance to the dominant living woody component of a forested wetland has been shown to affect wetland water table positions [17,18] and may affect soil temperatures due to the decreased amount of transpiration and increased amount of incoming solar radiation. Ecohydrological responses to a simulated EAB infestation were evaluated concurrently with gaseous soil carbon fluxes in depressional black ash wetland study sites examined in this study. Significantly lower site transpiration, smaller rates of water table drawdown during the growing season, and elevated summertime water tables were detected following a simulated EAB infestation treatment in depressional black ash wetland study sites [17]. Similarly, mineral soil black ash wetlands in northern Minnesota treated with a simulated EAB infestation had significantly elevated water tables when compared to control sites [18]. Ultimately, wetland soil carbon pools and fluxes will be concomitantly influenced by overstory disturbance and by the response of water table position and soil temperature.
Most research on wetland carbon storage and fluxes has focused on boreal or sub-boreal peatland ecosystems dominated by Sphagnum species and sparse coniferous tree species [19,20], yet large stores of carbon also exist in mixed deciduous forested wetlands with histic, histic-mineral, and mineral soils [21,22]. Histic-mineral and mineral soil wetlands are five times more abundant on earth’s surface than histosol wetlands, and comprise an estimated 22% of the global wetland carbon pool [10,23]. Freshwater mineral soil wetlands may contain significant organic soil deposits (20 to 40 cm depth), and have been inconsistently classified as either peatlands or fresh water mineral soil wetlands in regional- or global-scale wetland carbon cycle studies [24]. To our knowledge, no studies have specifically examined carbon cycling processes in deciduous forested peatlands.
Biogeochemophysical conditions such as temperature, pH, water table position, soil moisture, redox potential, and soil microbial communities affect the production and fluxes of gaseous soil carbon constituents (CO2, CH4) [19]. There are complex relationships between soil microbial community composition, physical and chemical conditions, the quality of available organic matter, and rates of soil decomposition [25]. As a result, it is often difficult to separate the influences of individual biogeochemophysical factors from one another. Peatland decomposition processes have consistently been shown to be sensitive to both soil temperature and redox conditions [26,27].
Carbon dioxide fluxes from wetlands generally increase following episodic short-term water table fluctuations; however, the magnitude of the response can be highly variable [26]. Specifically, CO2 production and flux may increase following episodic water table fluctuations due to redox reaction–induced chemical breakdown and activation of exo-enzymes following short-term exposure to oxygen [11,28,29]. In contrast, the quality of organic substrate and presence of certain microbial communities can inhibit decomposition rates following episodic oxidation-state modifying events [19]. It is unclear how EAB-induced alteration of water tables may affect oxidation–reduction processes and ultimately decomposition rates in black ash wetlands.
Temperature controls on CO2 production from wetlands have been well documented in both field and laboratory settings [15,30]. In laboratory studies, a 10 °C temperature increase resulted in a positive change in CO2 production (Q10) of two- to threefold [15]. However, in situ Q10 values have been inconsistent across vertical and lateral gradients within peatlands [30,31]. The uncertainty associated with the response of CO2 production rates to experimental temperature manipulations in peatlands has likely resulted from the complexity of in situ biogeochemical conditions. Additionally, the need to characterize soil warming impacts on peatland ecosystems remains critical for the development of accurate global carbon cycle models [32]. Removal of the overstory following an EAB infestation may increase insolation, increase wetland soil temperatures, and ultimately alter gaseous carbon fluxes.
Similar to CO2 production and fluxes, there are high degrees of spatial and temporal variability associated with CH4 production and fluxes from peatland soils [33]. Because both fermentation and CO2 reduction processes occur under anoxic conditions, the location of production and fate of CH4 once produced are highly sensitive to water table position [34]. Specific physical and chemical conditions need to be present for methane production to occur. Among the most important of these is the combined depletion of oxygen and absence of other, more efficient electron acceptors (e.g., N O 3 , F e 3 + , M n 4 + , and S O 4 2 ) that must be present before methanogenesis and CO2 reduction can proceed [28,35]. The initiation or rate of methane production has also been shown to be sensitive to pH [36], temperature [15], and substrate quality and availability [19,37].
Due to the co-occurrence of methanogens and methanotrophs in peatland soils, water table position is the single most influential physical parameter that regulates the quantity of CH4 fluxes from the soil surface. Several studies have evaluated the response of CH4 fluxes in peatlands to temperature and water table conditions [33,38,39,40]. Yet, to our knowledge, no study has examined the relationship between water table, soil temperature, and soil CO2 and CH4 fluxes in deciduous forested peatlands, nor has any study evaluated the response of gaseous soil carbon fluxes to a simulated invasive pest disturbance in deciduous forested peatlands. Temperature-dependent incubation experiments of CH4 production in cores extracted from hummock, lawn, and hollow locations in a peatland in southwestern Scotland revealed that CH4 production was similar among sites, yet in situ CH4 fluxes were highly variable and closely associated with the position of the water table [33]. The response of water tables and soil temperatures to an EAB-induced disturbance in black ash (Fraxinus nigra Marsh.) wetlands will likely alter CH4 production processes, and ultimately the amount of CH4 emitted from the soil surface.
This study examines the relationship between, and the coupled response of, water table, soil temperature, and gaseous soil carbon fluxes within the context of the previously established altered hydrologic regime in depressional black ash wetlands. The specific objectives of this study were to (1) characterize the relationships between water table position, soil temperature, and soil CO2 and CH4 fluxes; and (2) evaluate the responses of soil temperature and gaseous soil carbon fluxes to disturbances in black ash wetlands. Ultimately, this study will augment the currently limited knowledge of carbon cycling processes in deciduous forested peatlands and examine the ecosystem-level implications of a looming emerald ash borer infestation. We hypothesize that (1) water table position and soil temperature will be major influences on soil CO2 and CH4 fluxes from these soils; (2) soil temperature will increase because of increased insolation in disturbed sites; (3) soil CO2 fluxes will increase with the increased soil temperature in disturbed sites; and (4) soil CH4 fluxes will be higher because of the expected increased temperatures and elevated water tables detected in response to a simulated EAB infestation.

2. Materials and Methods

2.1. Study Area

The Ottawa National Forest (ONF) is located in the western portion of the upper peninsula of Michigan, USA (Figure 1). The ONF contains approximately 525 km2 of wetlands within 4000 km2 of mixed-hardwood forested area. Based upon US Department of Agriculture Forest Service Forest Inventory Analysis data, it is estimated that 70% of the 180 km2 deciduous forested wetlands within the ONF are dominated by black ash [41].
Black ash wetlands in the ONF commonly occur as landform depressions surrounded by mixed-hardwood upland tree species such as sugar maple (Acer saccharum L.) and white pine (Pinus strobus L.). Other tree species, such as red maple (Acer rubrum L.), yellow birch (Betula alleghaniensis Britt.), balsam fir (Abies balsamea L. (Mill)), and white cedar (Thuja occidentalis L.), are common, though typically at low abundance within these black ash–dominated wetlands [42]. Shrub, fern, and sedge species are abundant but unevenly distributed in the understory.
Black ash wetland soils in this study were composed of woody peats. The study site soils were, on average, 75% organic matter and 47% organic carbon [17]. Peat thicknesses in the study sites ranged from 40 cm to greater than 690 cm, with an average peat depth of 140 cm. A clay lens or a poorly sorted clay loam was commonly detected at the bottom of the peat layer. Black ash wetlands in the ONF are also characterized by seasonally inundated conditions and commonly have discernible surface drainage outlet channels. Most soils in the ONF originate from glacial landforms and may exist as thin horizons overlying near-surface bedrock [43].

2.2. Experimental Study Design

Nine study sites within the Ottawa National Forest, ranging in size from 0.25 to 1.25 ha, were established in 2011 in isolated landform depressions within first-order watersheds [17] (Figure 1). Black ash basal area ranged from 26% to 85% of total basal area among study sites (mean = 66%), and black ash was always the most abundant overstory tree [42]. Treatments were assigned randomly to 3 sites in 3 geographically distinct locations (n = 3 per treatment). Each distinct location included 1 untreated control, 1 girdle treatment, and 1 ash-cut treatment. A 1-year baseline data-collection period (2012) preceded the 2-year posttreatment monitoring period (2013–2014) of the study. The girdle treatment was used to simulate an active EAB infestation, wherein near-complete stand mortality could take 3–4 years [45]. In the girdled sites, all black ash trees greater than 2.5 cm in diameter at a height of 1.37 m had the bark, cambium, and phloem removed in a 15–30 cm tall circumferential band. The ash-cut treatment was used to simulate post-disturbance ecosystem conditions that might occur following complete loss of black ash from the overstory. In the ash-cut sites, all black ash trees greater than 2.5 cm at a height of 1.37 m were cut with chainsaws within 1.2 m of the soil surface during February 2013 and left on site. Treatments were applied between November 2012 and February 2013.

2.3. Instrumentation and Monitoring

Three random locations within each wetland site were determined with ArcGIS (Esri, Redlands, CA, USA), and were located at least 8 m from the wetland boundary. Two 13 cm tall gas collars constructed of 25.4 cm diameter PVC pipe were co-located 6 m apart and driven into the soil to a depth of 5 cm at each random location (n = 6 per site). Portable static dark chambers were used to measure both CO2 and CH4 fluxes [46,47], and were constructed of 25.4 cm ID PVC slip caps. A rim of closed cell foam was adhered to a PVC ring mounted within the chamber to ensure that there was an airtight seal between the gas collar and the chamber. The dark static chambers used to measure CH4 fluxes were equipped with a small fan to mix the chamber headspace [40]. When water table elevations exceeded 12 cm above the ground surface (6.3% of samples), high-density polyethylene extensions 50 cm long and 28 cm in diameter were used to connect the soil gas collar with the dark chamber.
CO2 and CH4 fluxes were measured at approximately monthly intervals during snow-free seasons (June–October) from 2012–2014. All measurements were collected on calm days between 11:00 and 17:30 h (EDT) to reduce insolation variability. Measurements were not collected if more than 6.4 mm of precipitation fell within the previous 8 h. All 18 gas collars within 1 geographically distinct location were measured on the same day, and all sites were measured within a 4-day period each month.
An infrared gas analyzer (IRGA) (PP Systems EGM-4, PP Systems International Inc., Amesbury, MA, USA) was used with a static dark chamber, equipped with a small fan to mix the chamber headspace, to measure CO2 fluxes [39,40]. CO2 concentrations were logged by the IRGA at 2–3 s intervals for a total of 124 s. The slope of CO2 concentration versus time (t, s) was estimated by linear regression [48]. Portions of nonlinear responses that occurred within the first 80 s of the measurement period were manually removed prior to calculating the slope. CO2 fluxes were then estimated according to Equation (1):
CO 2 F l u x = d CO 2 d t × PV ART
where d CO 2 d t is the slope from the regression, P is the atmospheric pressure (Pascal), V is the volume of the headspace within the chamber (m3), A is the area of the soil encircled by the gas collar (m2), R is the gas constant (J K−1 mol−1), and T is the air temperature (K). All dark chamber headspace volume estimates used in Equation (1) were corrected for water table position and extension headspace volume when necessary.
Gas-tight syringes were used to remove five 20 mL gas samples at consecutive 10 min intervals from each gas chamber [46]. Gas samples were immediately injected into pre-evacuated 6 mL EXETAINER® glass vials (Labco Ltd., Lampeter, Wales, UK), and the septa were coated with 100% silicone for added leakage protection. Gas samples were analyzed within 30 days using gas chromatography (Varian 3800, Agilent Technologies Inc., Santa Clara, CA, USA) at Michigan Technological University’s Hydrology and Wetlands Laboratory. As a quality-control measure, concentration standards were injected into EXETAINER® vials in the laboratory, transported to study sites during sampling excursions, and analyzed alongside field samples. The slope of the CH4 concentration versus time response was calculated using least-squares regression as described for CO2 above, and the slopes of the CH4 responses ( d CH d t ) were used in place of d CO 2 d t in Equation (1) to estimate CH4 fluxes [15]. Gaseous soil carbon fluxes were reported as mg m−2 day−1 of CO2 and CH4 molecules. Approximately 8% of samples with nonlinear slopes were flagged as possible human-induced ebullition events and were not included in the analysis.
Water levels were manually measured during monthly collection intervals and were recorded at 15 min intervals (Levelogger Junior M5, Solinst Canada Ltd., Georgetown, ON, Canada) in a steel monitoring well in each wetland [17]. The water level relative to the wetland soil surface ( W L m w , ( cm ) ) was calculated for each gas collar ( W L c o l , ( cm ) ) using total station surveyed elevations [17] at gas collars and metal monitoring wells, respectively. Soil temperature was measured within the gas collar 5 cm below the ground surface ( S T 5 c m ,   ° C ) immediately following each CO2 flux measurement. Soil temperature probes (107, Campbell Scientific, Logan, UT, USA) connected to on-site data loggers (CR800, Campbell Scientific) were installed to record every 15 min in all sites in 2011, but repeat damage by animals resulted in large amounts of missing data during 2012 and 2013. Temperature loggers (iButtons, Embedded Data Systems LLC, Lawrenceburg, KY, USA) set to record at 2 h intervals were subsequently installed 5 cm beneath the soil surface and within PVC solar shields located 1 m above the soil surface near the wetland monitoring well from June through November of 2014.

2.4. Statistical Analysis

Linear mixed-effects models (LMMs) were fit with the restricted maximum likelihood procedure in R [48,49] to assess the impact of water level, temperature, and treatment on soil CO2 and CH4 fluxes. Zero-concentration data (n = 400) from gas chromatograph calibration curves were used to determine the 0.36 ppm instrument detection limit for CH4 [50]. Based upon the 40 min collection period, the lower detectable soil CH4 flux limit was determined to be 1.25 (mg m−2 day−1). Approximately 20% of soil CH4 flux estimates were below the lower detection limit. Distribution parameters of the natural-log transformed CH4 flux data were calculated using maximum likelihood estimation to adjust for the censored values. Estimates for censored values were then imputed using random quantiles below the detection limit drawn from the fitted distribution as replacement values [51,52]. Following natural log-transformation of estimates of CO2 and CH4 from Equation (1), zero-mean assumptions of normality and homogeneity of variance were determined to be acceptable based on visual examination of residual plots for all linear mixed-effects models [48,49]. LMM methods were used to determine significance [53], the Tukey method was used to adjust p-values for all pairwise comparisons [54], and degrees of freedom were calculated according to the Kenward–Roger approximation [55]. Significance was reported where p < 0.05, and marginal significance was reported where p < 0.1. Following the transformation of CO2 and CH4 fluxes, the back-transformed regression equations took the form of Equation (2):
C F l u x = β ° e β 1 T s o i l
where β ° is the back-transformed model intercept coefficient, and β 1 is the model slope coefficient.
Logged temperature and water level data from 2014 were used to test the response of maximum daily soil temperatures ( S T m a x ) to maximum daily air temperatures ( A T m a x ) and W L m w . The influence of W L m w on the relationship between S T m a x and A T m a x was evaluated by comparing a two-level LMM ( A T m a x · W L m w ) with a single-level LMM ( A T m a x ), where A T m a x and W L m w were used as fixed effects, and site and day were used as crossed random effects. Additionally, the response of S T m a x to A T m a x , W L m w , and treatment during the posttreatment study period was tested with a three-level LMM where A T m a x , W L m w , and treatment were used as fixed effects, and site and day were used as crossed random effects.
LMM selection criteria, including the conditional variance ( r c 2 ) and marginal variance ( r m 2 ), the Akaike information criterion (AIC), the residual proportion change in variance (PCV), and residual standard error (RSE) were used to evaluate model performance [53,56]. Multivariate LMM response surfaces were predicted [57] by generating slopes for four water level positions: 10 cm above, 10 cm below, 30 cm below, and 50 cm below the wetland soil surface. Observations were binned with equidistant spacing between the four water level positions and used for visual interpretation only.
The influence of W L c o l on the relationship between soil gas flux and S T 5 c m was evaluated by comparing a two-level LMM ( S T 5 c m · W L c o l ) with a single-level LMM ( S T 5 c m ) for CO2 and CH4, respectively (Table 1). Both S T 5 c m and W L c o l were used as fixed effects, and correlated random intercepts and slopes were generated for each collar nested within site. Posttreatment CO2 and CH4 fluxes and responses to S T 5 c m and W L c o l were modeled with a three-level LMM using S T 5 c m , W L c o l , and treatment as interacting fixed effects with a correlated random intercept and slope for each collar nested within site. Annual differences in CO2 and CH4 soil fluxes by treatment were examined using a four-level LMM with S T 5 c m , W L c o l , treatment, and year as interacting fixed-effects terms, with a correlated random intercept and slope for each collar nested within site. A four-level LMM was constructed to compare CO2 fluxes between treatments and study periods using S T 5 c m , W L c o l , treatment, and study period (pre- or posttreatment) as interacting fixed effects, with a correlated random intercept and slope for each collar nested within site. Finally, mean untransformed soil gas fluxes were calculated for each site during each sampling event, and Kruskal–Wallis test was used to compare between-treatment means within each year.

3. Results

3.1. Regional Climate and Soil Temperature

Precipitation data for Ironwood, Michigan [44], indicated that the 83 cm of precipitation received in the ONF during 2012 was below the 30-year average of 89 cm, whereas the 127 cm and 111 cm received during 2013 and 2014, respectively, were above the 30-year average. The ONF received approximately 150% more snow-water equivalent (SWE) precipitation during winter 2012–2013 and 100% more SWE during winter 2013–2014 compared to winter 2011–2012. In addition, the ONF received approximately 25% more rain during 2013 and 15% more rain during 2014 compared to 2012.
A comparison of model selection criteria indicated that the full LMM ( A T m a x · W L m w ) for S T m a x performed better than the null model ( A T m a x only) by improving the relative quality (AIC), increasing the residual proportion change in variance (PCV), and decreasing the residual standard error (RSE). The full LMM explained 0.98 of the conditional variance ( r c 2 ) and 0.73 of the marginal variance ( r m 2 ) (Table 1). The 0.88 (°C °C−1) slope predicted for the S T m a x response to A T m a x when W L m w was farthest beneath the soil surface was larger than the 0.39 (°C °C−1) slope predicted when W L m w was above the soil surface (Figure 2). The mean daily soil temperatures of 13.2 and 13.6 °C in the girdle and ash-cut sites, respectively, were both greater than the 12.7 °C soil temperature detected in the controls, although these differences were not significant (Table 2).

3.2. Soil CO2 and CH4 Fluxes

As with soil temperature, a comparison of model statistics indicated that the full LMM ( S T 5 c m · W L c o l ) performed better than the null LMM ( S T 5 c m   o n l y ) for CO2 and CH4, respectively, by increasing the r c 2 and r m 2 , increasing the PCV, and decreasing the RSE (Table 1). The predicted slope of the CO2 flux response to soil temperature was 0.18 (mg m−2 day−1 °C−1) when W L c o l was farthest beneath the soil surface, and this value was larger than the 0.06 (mg m−2 day−1 °C−1) response slope predicted when W L c o l was above the soil surface (Figure 2) [57]. Conversely, the predicted slope of the CH4 flux response to soil temperature was 0.05 (mg m−2 day−1 °C−1) when W L c o l was above the soil surface, and this value was larger than the 0.006 (mg m−2 day−1 °C−1) response slope predicted when W L c o l was farthest below the soil surface (Figure 2).
Soil temperature ( S T 5 c m ) and water level ( W L c o l ) were used to predict soil CO2 and CH4 fluxes for each treatment during the posttreatment study period (Figure 3). The CO2 flux LMM had a 0.52 r c 2 and 0.35 r m 2 and a residual standard error of 2.1 mg m−2 day−1, whereas the CH4 flux LMM had a 0.37 r c 2 and 0.24 r m 2 and a residual standard error of 0.9 mg m−2 day−1. Due to the limited number of CO2 flux observations collected when soil temperatures were less than 10 °C during the pretreatment study period and the absence of pretreatment CH4 flux observations, the pretreatment response slope estimates were either not reliable or nonexistent and therefore were excluded from the three-level LMM. The girdle, ash-cut, and control response slopes for soil CO2 flux (0.10, 0.08, and 0.10 mg m−2 day−1 °C−1) and soil CH4 flux (0.02, 0.04, and 0.06 mg m−2 day−1 °C−1) were not significantly different from one another during the posttreatment study period. The mean CO2 flux detected in ash-cut sites was marginally significantly higher than the controls (p = 0.094), and the mean CH4 fluxes detected in the girdle and ash-cut sites were not significantly different (p > 0.1) than the control sites during the posttreatment study period (Figure 3).
There were no significant differences in mean soil CO2 flux detected among treatments during 2012 (pretreatment), whereas the 6599 mg m−2 day−1 ash-cut mean soil CO2 flux was marginally significantly higher (p = 0.076) than the 3193 mg m−2 day−1 control mean soil CO2 flux during 2013, and the 6722 and 6835 mean soil CO2 fluxes in the girdle (p = 0.042) and ash-cut (p = 0.042), respectively, were significantly higher than the 2979 mg m−2 day−1 control mean soil CO2 flux during 2014 (Figure 4). Neither the girdle nor the ash-cut soil CH4 fluxes were significantly different than the control during 2013 and 2014 (Figure 4). When the entire pre- and posttreatment periods were considered, the 5625 and 6545 mg m−2 day−1 mean soil CO2 fluxes detected in the girdle (p = 0.096) and ash-cut sites (p = 0.046), respectively, were significantly higher than the 3038 mg m−2 day−1 detected in the control sites during the posttreatment study period (Table 2). The nonparametric Kruskal–Wallis test shows that higher CO2 fluxes occurred more frequently in the ash-cut than the controls (p = 0.085) (Table 2). The 11.9 mg m−2 day−1 ash-cut mean soil CH4 flux was almost twice as high as the 6.6 mg m−2 day−1 control mean soil CH4 flux, and marginally significantly (p = 0.079) higher than the 2.4 mg m−2 day−1 girdle mean soil CH4 flux during the posttreatment study period (Table 2). While the mean CH4 fluxes detected in girdle and ash-cut sites were similar to controls during the posttreatment period, the nonparametric Kruskal–Wallis test shows that higher CH4 fluxes occurred more frequently in the ash-cut sites than the controls (p = 0.026) (Table 2). Within-treatment soil CO2 flux comparisons were not considered due to the limited number of flux observations collected when soil temperatures were less than 10 °C during the pretreatment study period, and no pretreatment CH4 flux data were available.

4. Discussion

4.1. Soil Temperature and Water Level

Wetland water level ( W L m w ) significantly altered the response of soil temperature ( S T m a x ) to air temperature ( A T m a x ; Figure 2). When W L m w was lower and the soil presumably drier, S T m a x more closely paralleled A T m a x . In contrast, when W L m w was higher (near or above surface), S T m a x was less correlated to A T m a x . The combined influence of the higher specific heat capacity of water and the consistent groundwater discharge into black ash wetlands [17] likely decreased the sensitivity of S T m a x to A T m a x when soils were saturated. Furthermore, soil temperature response to treatment was likely mitigated by the consistent influx of groundwater because of the expected increase in magnitude of latent heat exchange and reduction in magnitude of sensible heat exchange between soil and atmosphere [59].
The wetter conditions (Table 2) and higher wetland water levels during the posttreatment period [17] may have buffered the magnitude of S T m a x response to disturbance. Because of the wetter-than-average posttreatment water years and the resultant elevated water tables, the average soil temperatures in girdle and ash-cut sites were not significantly greater than the controls, despite the presumed increase in canopy openness and insolation. We theorized that the increased amount of insolation likely received by soils in the disturbed study sites would have a greater influence on the response of S T m a x to A T m a x during drier periods. It is likely that average S T m a x differences would have been larger during the posttreatment study period under drier conditions.

4.2. Responses of Soil CO2 and CH4 Fluxes to Disturbance

The water level position relative to the wetland soil surface ( W L c o l ) regulated the response of soil CO2 and CH4 fluxes to soil temperature ( S T 5 c m ; Figure 2), and co-varied soil CO2 and CH4 flux responses to S T 5 c m and W L c o l have been observed in numerous wetland types [15,19,60]. Soil CO2 fluxes tended to be higher and increase more rapidly in response to increased temperature when W L c o l was below the wetland soil surface, whereas soil CH4 fluxes tended to be higher and increase more rapidly in response to increased temperature when W L c o l was at or above the wetland soil surface. The magnitude of observed soil CO2 and CH4 fluxes from depressional black ash wetlands was similar to fluxes observed in wetlands under field and laboratory conditions [15,61,62,63].
Given that both the magnitude of soil CO2 and CH4 fluxes and the correlation of S T m a x to A T m a x were controlled by wetland water levels, it is likely that interannual climatic conditions also affect the relative magnitude of mean annual soil CO2 and CH4 fluxes. Specifically, mean annual soil CO2 fluxes may be higher during drier periods and lower during wetter periods. Conversely, soil CH4 fluxes would be expected to be higher during wetter periods and lower during drier periods. Therefore, it is possible that the soil CO2 flux response to disturbance was buffered and the soil CH4 flux response was enhanced because of the wetter-than-average conditions during the posttreatment study period [17]. Black ash wetland water levels tended to be higher during the two wetter water years [17], and soil temperature tended to be less responsive to air temperature when water levels were higher (Figure 2). Moreover, wetland soil CO2 fluxes are typically lower when water tables are above the surface compared to fluxes measured during lower water table positions [19,62] Therefore, soil CO2 fluxes from disturbed sites may have been higher and possibly more dissimilar than control sites if drier climatic conditions had occurred during the two posttreatment monitoring years.
The magnitude of soil CO2 fluxes was significantly higher in disturbed sites, and high soil CH4 fluxes occurred more frequently in the ash-cut than the control during the posttreatment study period (Figure 3 and Figure 4; Table 2). As evidenced by the similarity among slopes, the sensitivity of soil CO2 fluxes to environmental variables ( S T 5 c m and W L c o l ) did not change following disturbance (Figure 3). Therefore, the higher soil CO2 fluxes from disturbed sites were not caused by an altered sensitivity to environmental variables. Similarly, soil CH4 response slopes for the three treatments were not significantly different from one another during the posttreatment study period (Figure 3), and the significantly more frequent occurrence of high CH4 fluxes in ash-cut sites may not have been caused by an altered sensitivity to environmental variables.
It is possible that the increased rate of heterotrophic respiration from disturbed soils stemmed from an alteration in the organic chemical composition of peatland soil [64]. Carbon mineralization rates were positively correlated with labile carbon compounds and negatively correlated with recalcitrant carbon compounds in a northern Minnesota peatland [65]. Increased rates of CH4 production were detected following the addition of glucose to lignin-rich peat in northern North American peatlands matter [66], and CO2 and CH4 production was significantly affected by pH and C:N ratios in Dutch peat grasslands [67]. Fine root respiration may also account for a substantial proportion of soil CO2 fluxes, and increased soil CO2 fluxes have been attributed to turnover and decomposition of fine roots [68,69,70]. Black ash trees have a shallow, fibrous root system [71], and the majority of this fine root system was likely available for decomposition following disturbance. Unlike the rapid infilling from co-occurring tree species that occurred following EAB infestation in green ash (Fraxinus pennsylvanica Marsh.) stands in northwest Ohio [72], neither overstory or understory stems of co-occurring tree species within these black ash wetland study sites responded positively following treatment [42]. Therefore, the increased magnitude of CO2 fluxes (Figure 4) from disturbed black ash wetland soils may be at least partially attributed to the transient large influx of decomposable fine roots and an associated response in microbial activity and biomass [63,73], and not likely the result of increased autotrophic respiration from the release of co-occurring tree species. Inversely, the loss of autotrophic root respiration from mature individuals may lead to a decrease in soil CO2 flux following disturbance [74,75]. A lack of response of inorganic soil nitrogen availability following treatment in these sites during this study period [76] supports the explanation that increased microbial activity drove the observed CO2 response.
The study sites were consistently connected to mineral-rich groundwater during the study period [17], and this connection likely influenced carbon mineralization dynamics within depressional black ash wetland soils. The relative standard error of soil CH4 fluxes (18.0%) was considerably larger than the relative standard error of soil CO2 fluxes (<0.1%; Table 2). Because the highest rates of methane oxidation are expected near oxic/anoxic interfaces [28], the large variation of soil CH4 fluxes that occurred when wetlands were saturated and/or inundated (Figure 3) was not likely caused by increased rates of methanotrophy within the peat profile. This further suggests that methanogenesis may have been constrained by the availability of electron acceptors [77], inhibited by chemical compounds [28] or low pH [78], or transported to the atmosphere through vegetation with arenchyma [79,80].
When compared to control sites, and because of the warmer and wetter conditions observed in the disturbed sites after treatment, we expected to observe higher posttreatment mean soil CH4 fluxes in ash-cut and girdle sites. However, the persistent connection with mineral-rich groundwater during the wetter posttreatment study period may, in fact, have buffered the relative magnitude of CH4 flux response to disturbance by inhibiting methanogenesis. Specifically, it is likely that study sites were consistently being replenished with more efficient electron-accepting redox species such as F e 2 + , M n 4 + , or S O 4 2 during the wetter posttreatment period, which effectively inhibited methanogenesis. It is also possible that methane production within the soil increased following disturbance, but increased CH4 fluxes were not detected from the soil surface due to plant-mediated transport.
The large increase of herbaceous cover following disturbance in these sites may have limited natural woody regeneration recruitment to larger seedling size classes or saplings [42], and canopy removal has been shown to convert black ash–dominated wetlands to herbaceous or scrub-shrub–dominated wetlands [81]. As a result, canopy replacement through natural regeneration may take several decades [82], where planting mitigation efforts will influence the tree species composition, likely decrease the duration of canopy recovery, and help maintain ecosystem function following disturbance [83,84,85]. Though the hydrogeomorphic setting is distinctly different from the black ash wetlands considered here, these findings may have important implications for potential effects in other ash-dominated wetlands such as those dominated by green and pumpkin ash (Fraxinus profunda Bush). The contribution of herbaceous vegetation to soil respiration following disturbance may be substantial [86], and the contribution of the herbaceous community to autotrophic soil respiration may have increased, given the increased amount of herbaceous cover detected in these sites post treatment [42]. In addition, the increased magnitude of heterotrophic respiration associated with the elevated soil temperatures and transient large influx of decomposable fine roots likely contributed to the significantly higher girdle and ash-cut mean soil CO2 fluxes during the posttreatment period when compared to controls. Warmer and wetter conditions are expected in the upper Great Lakes region under elevated atmospheric CO2 conditions, which may also affect soil CH4 fluxes in EAB-affected wetlands [87]. Therefore, it is likely that successful planting mitigation efforts will effectively decrease the duration of EAB infestation–induced increases of soil CO2 fluxes and buffer the losses of soil C storage within depressional black ash wetlands.

5. Conclusions

Control, girdle, and ash-cut treatments were applied to three sites and used to simulate an EAB disturbance. Water level, soil temperature, and soil CO2 and CH4 fluxes were monitored during pre- and posttreatment study periods. Mean soil temperatures were higher in girdle and ash-cut sites compared to the controls. The relative water table position influenced the degree of soil temperature sensitivity to air temperature, where soil temperature more closely paralleled air temperature when water levels were below the wetland soil surface. Consequently, the wetter-than-average posttreatment climatic conditions [17] likely inhibited the magnitude of soil temperature response to disturbance where average or drier-than-average climatic conditions may have yielded even larger differences between treated and untreated sites.
Soil CO2 fluxes were significantly higher and higher soil CH4 fluxes occurred more frequently in disturbed sites, but mean soil CH4 fluxes were similar in disturbed and undisturbed sites. Soil temperature and wetland water levels controlled soil CO2 and CH4 fluxes, where warmer and drier conditions yielded higher soil CO2 fluxes, and warmer and wetter conditions yielded higher soil CH4 fluxes. The relationship between soil CO2 and CH4 fluxes and soil temperature and water level was not influenced by a simulated EAB disturbance. The combination of higher mean soil temperatures, increased herbaceous cover, and the transient large influx of decomposable fine roots following canopy suppression may have contributed to the higher soil CO2 fluxes detected in disturbed sites. Conversely, the consistent connection with mineral-rich groundwater flow likely constrained methanogenesis, and consequently reduced the relative magnitude of soil CH4 flux response to a simulated EAB disturbance.
Black ash wetland soil carbon mineralization processes were sensitive, but also showed resilience to EAB-induced impacts. Both CO2 and CH4 fluxes were regulated by soil temperature and water level position. The consistent connection with groundwater flow systems and wetter-than-average weather conditions during the posttreatment period likely buffered the relative magnitude of soil CH4 response to disturbance.

Author Contributions

M.V.G., J.S., J.D., R.K., and T.P. conceived and designed the experiments; M.V.G., J.S., N.B., and J.D. performed the experiments; M.V.G., J.S., and N.B. analyzed the data; and all authors contributed to writing the paper.

Acknowledgments

We appreciate the support of the US Forest Service, in particular Mark Fedora, and the rest of the Ottawa National Forest. Funding for this project was provided by the USDA, the US-EPA Great Lakes Restoration Initiative (DW-12-92429101-0), Michigan Tech’s Center for Water and Society, and Michigan Tech’s Ecosystem Science Center. Special thanks to Stephen Sebestyen, Alex Mayer, Robert Froese, Nicholas Schreiner, Jarrod Nelson, Leah Harrison, Erica Jones, Jon Bontrager, Dan Hutchinson, John Hrjiblan, Lynette Potvin, and Timothy Veverica for their technical and logistical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McLaughlin, J.W.; Gale, M.R.; Jurgensen, M.F.; Trettin, C.C. Soil organic matter and nitrogen cycling in response to harvesting, mechanical site preparation, and fertilization in a wetland with a mineral substrate. For. Ecol. Manag. 2000, 129, 7–23. [Google Scholar] [CrossRef]
  2. Sun, G.; McNulty, S.G.; Shepard, J.P.; Amatya, D.M.; Riekerk, H.; Comerford, N.B.; Skaggs, W.; Swift, L. Effects of timber management on the hydrology of wetland forests in the southern united states. For. Ecol. Manag. 2001, 143, 227–236. [Google Scholar] [CrossRef]
  3. MacFarlane, D.W.; Meyer, S.P. Characteristics and distribution of potential ash tree hosts for emerald ash borer. For. Ecol. Manag. 2005, 213, 15–24. [Google Scholar] [CrossRef]
  4. Kashian, D.M.; Witter, J.A. Assessing the potential for ash canopy tree replacement via current regeneration following emerald ash borer-caused mortality on southeastern Michigan landscapes. For. Ecol. Manag. 2011, 261, 480–488. [Google Scholar] [CrossRef]
  5. Flower, C.E.; Knight, K.S.; Gonzalez-Meler, M.A. Impacts of the emerald ash borer (Agrilus planipennis fairmaire) induced ash (Fraxinus spp.) mortality on forest carbon cycling and successional dynamics in the eastern united states. Biol. Invasions 2013, 15, 931–944. [Google Scholar] [CrossRef]
  6. Sydnor, T.D.; Bumgardner, M.; Todd, A. The potential economic impacts of emerald ash borer (Agrilus planipennis) on Ohio, US, Communities. Arboricult. Urban For. 2007, 33, 48–54. [Google Scholar]
  7. Kovacs, K.F.; Haight, R.G.; McCullough, D.G.; Mercader, R.J.; Siegert, N.W.; Liebhold, A.M. Cost of potential emerald ash borer damage in U.S. Communities, 2009–2019. Ecol. Econ. 2010, 69, 569–578. [Google Scholar] [CrossRef]
  8. Gandhi, K.K.; Herms, D. Direct and indirect effects of alien insect herbivores on ecological processes and interactions in forests of eastern north America. Biol. Invasions 2010, 12, 389–405. [Google Scholar] [CrossRef]
  9. Lovett, G.M.; Canham, C.D.; Arthur, M.A.; Weathers, K.C.; Fitzhugh, R.D. Forest ecosystem responses to exotic pests and pathogens in eastern north America. BioScience 2006, 56, 395–405. [Google Scholar] [CrossRef]
  10. Kimble, J.M. The Potential of U.S. Forest Soils to Sequester Carbon and Mitigate the Greenhouse Effect; CRC Press: Boca Raton, FL, USA, 2003; p. 429. [Google Scholar]
  11. Freeman, C.; Ostle, N.; Kang, H. An enzymic ‘latch’ on a global carbon store–A shortage of oxygen locks up carbon in peatlands by restraining a single enzyme. Nature 2001, 409, 149. [Google Scholar] [CrossRef] [PubMed]
  12. Gorham, E. Northern peatlands: Role in the carbon cycle and probable responses to climatic warming. Ecol. Appl. 1991, 1, 182–195. [Google Scholar] [CrossRef] [PubMed]
  13. Hogberg, P.; Nordgren, A.; Buchmann, N.; Taylor, A.F.S.; Ekblad, A.; Hogberg, M.N.; Nyberg, G.; Ottosson-Lofvenius, M.; Read, D.J. Large-scale forest girdling shows that current photosynthesis drives soil respiration. Nature 2001, 411, 789–792. [Google Scholar] [CrossRef] [PubMed]
  14. Zerva, A.; Mencuccini, M. Short-term effects of clearfelling on soil CO2, CH4, and N2O fluxes in a sitka spruce plantation. Soil Biol. Biochem. 2005, 37, 2025–2036. [Google Scholar] [CrossRef]
  15. Moore, T.R.; Dalva, M. The influence of temperature and water-table position on carbon-dioxide and methane emissions from laboratory columns of peatland soils. J. Soil Sci. 1993, 44, 651–664. [Google Scholar] [CrossRef]
  16. Macdonald, J.A.; Fowler, D.; Hargreaves, K.J.; Skiba, U.; Leith, I.D.; Murray, M.B. Methane emission rates from a northern wetland; response to temperature, water table and transport. Atmos. Environ. 1998, 32, 3219–3227. [Google Scholar] [CrossRef]
  17. Van Grinsven, M.J.; Shannon, J.P.; Davis, J.C.; Bolton, N.W.; Wagenbrenner, J.W.; Kolka, R.K.; Grant Pypker, T. Source water contributions and hydrologic responses to simulated emerald ash borer infestations in depressional black ash wetlands. Ecohydrology 2017, 10. [Google Scholar] [CrossRef]
  18. Slesak, R.A.; Lenhart, C.F.; Brooks, K.N.; D’Amato, A.W.; Palik, B.J. Water table response to harvesting and simulated emerald ash borer mortality in black ash wetlands in Minnesota, USA. Can. J. For. Res. 2014, 44, 961–968. [Google Scholar] [CrossRef]
  19. Blodau, C. Carbon cycling in peatlands—A review of processes and controls. Environ. Rev. 2002, 10, 111–134. [Google Scholar] [CrossRef]
  20. Nancy, D.; Narasinha, S.; Peter, W.; Shashi, V.; Elon, V.; Eville, G.; Patrick, C.; Robert, H.; Cheryl, K.; Joseph, Y.; et al. Carbon emissions from peatlands. In Peatland Biogeochemistry and Watershed Hydrology at the Marcell Experimental Forest; CRC Press: Boca Raton, FL, USA, 2011; pp. 297–347. [Google Scholar]
  21. McFee, W.W.; Kelly, J.M. Soil carbon in northern forested wetlands: Impacts of silvicultural practices. Carbon 1995, 1, 437–461. [Google Scholar]
  22. Trettin, C.C.; Jurgensen, M.F.; Grigal, D.F.; Gale, M.R.; Jeglum, J.R. Northern Forested Wetlands Ecology and Management; CRC Press: Boca Raton, FL, USA, 1996. [Google Scholar]
  23. Eswaran, H.; Van Den Berg, E.; Reich, P. Organic carbon in soils of the world. Soil Sci. Soc. Am. J. 1993, 57, 192–194. [Google Scholar] [CrossRef]
  24. Bridgham, S.; Megonigal, J.P.; Keller, J.; Bliss, N.; Trettin, C. The carbon balance of north American wetlands. Wetlands 2006, 26, 889–916. [Google Scholar] [CrossRef]
  25. Freeman, C.; Kim, S.Y.; Lee, S.H.; Kang, H. Effects of elevated atmospheric CO2 concentrations on soil microorganisms. J. Microbiol. 2004, 42, 267–277. [Google Scholar] [PubMed]
  26. Blodau, C.; Basiliko, N.; Moore, T.R. Carbon turnover in peatland mesocosms exposed to different water table levels. Biogeochemistry 2004, 67, 331–351. [Google Scholar] [CrossRef] [Green Version]
  27. Dalva, M.; Moore, T.R.; Arp, P.; Clair, T.A. Methane and soil and plant community respiration from wetlands, Kejimkujik National Park, Nova Scotia: Measurements, predictions, and climatic change. J. Geophys. Res. Atmos. 2001, 106, 2955–2962. [Google Scholar] [CrossRef] [Green Version]
  28. Segers, R. Methane production and methane consumption: A review of processes underlying wetland methane fluxes. Biogeochemistry 1998, 41, 23–51. [Google Scholar] [CrossRef]
  29. Aller, R.C. Bioturbation and remineralization of sedimentary organic-matter—Effects of redox oscillation. Chem. Geol. 1994, 114, 331–345. [Google Scholar] [CrossRef]
  30. van Hulzen, J.B.; Segers, R.; van Bodegom, P.M.; Leffelaar, P.A. Temperature effects on soil methane production: An explanation for observed variability. Soil Biol. Biochem. 1999, 31, 1919–1929. [Google Scholar] [CrossRef]
  31. Fissore, C.; Giardina, C.P.; Kolka, R.K.; Trettin, C.C. Soil organic carbon quality in forested mineral wetlands at different mean annual temperature. Soil Biol. Biochem. 2009, 41, 458–466. [Google Scholar] [CrossRef]
  32. Hagerty, S.B.; Van Groenigen, K.J.; Allison, S.D.; Hungate, B.A.; Schwartz, E.; Koch, G.W.; Kolka, R.K.; Dijkstra, P. Accelerated microbial turnover but constant growth efficiency with warming in soil. Nat. Clim. Chang. 2014, 4, 903–906. [Google Scholar] [CrossRef]
  33. Daulat, W.E.; Clymo, R.S. Effects of temperature and watertable on the efflux of methane from peatland surface cores. Atmos. Environ. 1998, 32, 3207–3218. [Google Scholar] [CrossRef]
  34. Kettunen, A.; Kaitala, V.; Lehtinen, A.; Lohila, A.; Alm, J.; Silvola, J.; Martikainen, P.J. Methane production and oxidation potentials in relation to water table fluctuations in two boreal mires. Soil Biol. Biochem. 1999, 31, 1741–1749. [Google Scholar] [CrossRef]
  35. Peters, V.; Conrad, R. Sequential reduction processes and initiation of CH4 production upon flooding of oxic upland soils. Soil Biol. Biochem. 1996, 28, 371–382. [Google Scholar] [CrossRef]
  36. Moore, T.R.; Heyes, A.; Roulet, N.T. Methane emissions from wetlands, southern hudson-bay lowland. J Geophys. Res. Atmos. 1994, 99, 1455–1467. [Google Scholar] [CrossRef]
  37. Amaral, J.A.; Knowles, R. Methane metabolism in a temperate swamp. Appl. Environ. Microbiol. 1994, 60, 3945–3951. [Google Scholar] [PubMed]
  38. Shannon, R.D.; White, J.R. A three-year study of controls on methane emissions from two Michigan peatlands. Biogeochemistry 1994, 27, 35–60. [Google Scholar] [CrossRef]
  39. Ballantyne, D.; Hribljan, J.; Pypker, T.; Chimner, R. Long-term water table manipulations alter peatland gaseous carbon fluxes in northern Michigan. Wetl. Ecol. Manag. 2014, 22, 35–47. [Google Scholar] [CrossRef]
  40. Johnson, C.P.; Pypker, T.G.; Hribljan, J.A.; Chimner, R.A. Open top chambers and infrared lamps: A comparison of heating efficacy and CO2/CH4 dynamics in a northern Michigan peatland. Ecosystems 2013, 16, 736–748. [Google Scholar] [CrossRef]
  41. Prasad, A.M.; Iverson, L.R. Little’s Range and FIA Importance Value Database for 135 Eastern US Tree Species; Northeastern Research Station: Delaware, OH, USA, 2003. [Google Scholar]
  42. Davis, J.C.; Shannon, J.P.; Bolton, N.W.; Kolka, R.K.; Pypker, T. Vegetation responses to simulated emerald ash borer infestation in Fraxinus nigra dominated wetlands of Upper Michigan, USA. Can. J. For. Res. 2017, 47, 319–330. [Google Scholar] [CrossRef]
  43. Soil Survey Staff. Web Soil Survey. Natural Resources Conservation Service, United States Department of Agriculture. Available online: http://websoilsurvey.Nrcs.Usda.Gov/ (accessed on 14 November 2014).
  44. NOAA-NCEI. National Oceanic and Atomospheric Administration. National Centers for Environmental Information: Asheville, NC, USA. Available online: https://www.ncei.noaa.gov/ (accessed on 15 November 2017).
  45. McCullough, D.G.; Schneeberger, N.F.; Katovich, S.A.; Siegert, N.W. Pest Alert: Emerald Ash Borer. USDA For. Serv., Northeastern Area, State & Private Forestry: Newtown Square, PA, USA, 2015. NA-PR-02-04. Available online: http://www.na. fs.fed.us/spfo/pubs/pest_al/eab/eab.pdf (accessed on 5 December 2017).
  46. Tuittila, E.-S.; Komulainen, V.-M.; Vasander, H.; Nykänen, H.; Martikainen, P.J.; Laine, J. Methane dynamics of a restored cut-away peatland. Glob. Chang. Biol. 2000, 6, 569–581. [Google Scholar] [CrossRef]
  47. Carroll, P.; Crill, P. Carbon balance of a temperate poor fen. Glob. Biogeochem. Cycles 1997, 11, 349–356. [Google Scholar] [CrossRef] [Green Version]
  48. R Core Team. R: A Language and Environment for Statistical Computing [Computer Software]; R: Foundation for Statistical Computing: Vienna, Austria, 2016. [Google Scholar]
  49. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting linear mixed-effects models using LME4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  50. Skoog, D.A.; Holler, J.; Crouch, S.R. Liquid chromatography. In Principles of Instrumental Analysis, 6th ed.; Thomson Cengage Learning: Boston, MA, USA, 2006; pp. 816–855. [Google Scholar]
  51. Hewett, P.; Ganser, G.H. A comparison of several methods for analyzing censored data. Ann. Occup. Hyg. 2007, 51, 611–632. [Google Scholar] [PubMed]
  52. Millard, S.P. Envstats, an R Package for Environmental Statistics; Springer Science & Business Media: Berlin, Germany, 2013. [Google Scholar]
  53. Nakagawa, S.; Schielzeth, H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol. Evol. 2013, 4, 133–142. [Google Scholar] [CrossRef]
  54. Lenth, R.V. Lsmeans: Least-Squares Means. R Package Version 2013, Volume 2. Available online: https://www.jstatsoft.org/article/view/v069i01/0 (accessed on 18 October 2016).
  55. Halekoh, U.; Højsgaard, S.; Højsgaard, M.S.; Matrix, I. Package ‘pbkrtest’. 2014. Available online: https://cran.r-project.org/web/packages/pbkrtest/index.html (accessed on 23 September 2016).
  56. Gałecki, A.; Burzykowski, T. Linear Mixed-Effects Models Using R; Springer: New York, NY, USA, 2013. [Google Scholar]
  57. Fox, J. Effect displays in R for generalised linear models. J. Stat. Softw. 2003, 8, 1–27. [Google Scholar] [CrossRef]
  58. Wickham, H.; Chang, W. Ggplot2: An Implementation of the Grammar of Graphics. 2015. Available online: http://ftp.auckland.ac.nz/software/CRAN/src/contrib/Descriptions/ggplot.html (accessed on 11 August 2016).
  59. Walko, R.L.; Band, L.E.; Baron, J.; Kittel, T.G.; Lammers, R.; Lee, T.J.; Ojima, D.; Pielke, R.A., Sr.; Taylor, C.; Tague, C. Coupled atmosphere–biophysics–hydrology models for environmental modeling. J. Appl. Meteorol. 2000, 39, 931–944. [Google Scholar] [CrossRef]
  60. Kolka, R.K.; Sebestyen, S.D.; Verry, E.S.; Brooks, K.N. Peatland Biogeochemistry and Watershed Hydrology at the Marcell Experimental Forest; CRC Press Taylor & Francis Group: Boca Raton, FL, USA, 2011. [Google Scholar]
  61. Noh, N.J.; Shannon, J.; Bolton, N.; Davis, J.; Van Grinsven, M.; Pypker, T.; Kolka, R.; Wagenbrenner, J. Carbon dioxide fluxes from coarse dead wood in black ash wetlands. Forests 2018, in press. [Google Scholar]
  62. Dinsmore, K.; Skiba, U.; Billett, M.; Rees, R. Effect of water table on greenhouse gas emissions from peatland mesocosms. Plant Soil 2009, 318, 229–242. [Google Scholar] [CrossRef] [Green Version]
  63. Chimner, R.A.; Cooper, D.J. Influence of water table levels on CO2 emissions in a colorado subalpine fen: An in situ microcosm study. Soil Biol. Biochem. 2003, 35, 345–351. [Google Scholar] [CrossRef]
  64. Leifeld, J.; Steffens, M.; Galego-Sala, A. Sensitivity of peatland carbon loss to organic matter quality. Geophys. Res. Lett. 2012, 39. [Google Scholar] [CrossRef] [Green Version]
  65. Updegraff, K.; Pastor, J.; Bridgham, S.D.; Johnston, C.A. Environmental and substrate controls over carbon and nitrogen mineralization in northern wetlands. Ecol. Appl. 1995, 5, 151–163. [Google Scholar] [CrossRef]
  66. Yavitt, J.B.; Williams, C.J.; Wieder, R.K. Production of methane and carbon dioxide in peatland ecosystems across north America: Effects of temperature, aeration, and organic chemistry of pent. Geomicrobiol. J. 1997, 14, 299–316. [Google Scholar] [CrossRef]
  67. Best, E.P.H.; Jacobs, F.H.H. The influence of raised water table levels on carbon dioxide and methane production in ditch-dissected peat grasslands in the netherlands. Ecol. Eng. 1997, 8, 129–144. [Google Scholar] [CrossRef]
  68. Ewel, K.C.; Cropper, W.P., Jr.; Gholz, H.L. Soil CO2 evolution in florida slash pine plantations. Ii. Importance of root respiration. Can. J. For. Res. 1987, 17, 330–333. [Google Scholar] [CrossRef]
  69. Hanson, P.J.; Edwards, N.T.; Garten, C.T.; Andrews, J.A. Separating root and soil microbial contributions to soil respiration: A review of methods and observations. Biogeochemistry 2000, 48, 115–146. [Google Scholar] [CrossRef]
  70. Edwards, N.T.; Harris, W.F. Carbon cycling in a mixed deciduous forest floor. Ecology 1977, 58, 431–437. [Google Scholar] [CrossRef]
  71. Burns, R.M.; Honkala, B.H. Silvics of North America; United States Department of Agriculture: Washington, DC, USA, 1990; Volume 2. [Google Scholar]
  72. Costilow, K.; Knight, K.; Flower, C. Disturbance severity and canopy position control the radial growth response of maple trees (Acer spp.) in forests of northwest Ohio impacted by emerald ash borer (Agrilus planipennis). Ann. For. Sci. 2017, 74, 10. [Google Scholar] [CrossRef]
  73. Thomas, K.L.; Benstead, J.; Davies, K.L.; Lloyd, D. Role of wetland plants in the diurnal control of CH4 and CO2 fluxes in peat. Soil Biol. Biochem. 1996, 28, 17–23. [Google Scholar] [CrossRef]
  74. Nave, L.; Gough, C.; Maurer, K.; Bohrer, G.; Hardiman, B.; Le Moine, J.; Munoz, A.; Nadelhoffer, K.; Sparks, J.; Strahm, B. Disturbance and the resilience of coupled carbon and nitrogen cycling in a north temperate forest. J. Geophys. Res. Biogeosci. 2011, 116. [Google Scholar] [CrossRef] [Green Version]
  75. Nuckolls, A.E.; Wurzburger, N.; Ford, C.R.; Hendrick, R.L.; Vose, J.M.; Kloeppel, B.D. Hemlock declines rapidly with hemlock woolly adelgid infestation: Impacts on the carbon cycle of southern Appalachian forests. Ecosystems 2009, 12, 179–190. [Google Scholar] [CrossRef]
  76. Davis, J.C. Vegetation Dynamics and Nitrogen Cycling Responses to Simulated Emerald Ash Borer Infestation in Fraxinus nigra-Dominated Wetlands of Upper Michigan, USA; Michigan Technological University: Houghton, MI, USA, 2016. [Google Scholar]
  77. Zehnder, A.; Stumm, W. Geochemistry and biogeochemistry of anaerobic habitats. In Biology of Anaerobic Microorganisms; John Wiley & Sons, Inc.: New York, NY, USA, 1988; pp. 1–38. [Google Scholar]
  78. Dunfield, P.; Knowles, R.; Dumont, R.; Moore, T.R. Methane production and consumption in temperate and subarctic peat soils: Response to temperature and PH. Soil Biol. Biochem. 1993, 25, 321–326. [Google Scholar] [CrossRef]
  79. Moore, T.R.; Dalva, M. Methane and carbon dioxide exchange potentials of peat soils in Aerobic and Anaerobic laboratory incubations. Soil Biol. Biochem. 1997, 29, 1157–1164. [Google Scholar] [CrossRef]
  80. Gauci, V.; Gowing, D.J.G.; Hornibrook, E.R.C.; Davis, J.M.; Dise, N.B. Woody stem methane emission in mature wetland alder trees. Atmos. Environ. 2010, 44, 2157–2160. [Google Scholar] [CrossRef]
  81. Erdmann, G.G.; Crow, T.R.; Ralph, M., Jr.; Wilson, C.D. Managing Black Ash in the Lake States; General Technical Report NC-115; US Dept. of Agriculture, Forest Service, North Central Forest Experiment Station: St. Paul, MN, USA, 1987; Volume 115.
  82. Flower, C.E.; Gonzalez-Meler, M.A. Responses of temperate forest productivity to insect and pathogen disturbances. Annu. Rev. Plant Biol. 2015, 66, 547–569. [Google Scholar] [CrossRef] [PubMed]
  83. Looney, C.E.; D’Amato, A.W.; Palik, B.J.; Slesak, R.A. Canopy treatment influences growth of replacement tree species in Fraxinus nigra forests threatened by the emerald ash borer in Minnesota, USA. Can. J. For. Res. 2016, 47, 183–192. [Google Scholar] [CrossRef]
  84. Bolton, N.; Shannon, J.; Davis, J.; Grinsven, M.V.; Noh, N.J.; Schooler, S.; Kolka, R.; Pypker, T.; Wagenbrenner, J. Methods to improve survival and growth of planted alternative species seedlings in black ash ecosystems threatened by emerald ash borer. Forests 2018, 9, 146. [Google Scholar] [CrossRef]
  85. Iverson, L.; Knight, K.S.; Prasad, A.; Herms, D.A.; Matthews, S.; Peters, M.; Smith, A.; Hartzler, D.M.; Long, R.; Almendinger, J. Potential species replacements for black ash (Fraxinus nigra) at the confluence of two threats: Emerald ash borer and a changing climate. Ecosystems 2016, 19, 248–270. [Google Scholar] [CrossRef]
  86. Zehetgruber, B.; Kobler, J.; Dirnböck, T.; Jandl, R.; Seidl, R.; Schindlbacher, A. Intensive ground vegetation growth mitigates the carbon loss after forest disturbance. Plant Soil 2017, 420, 239–252. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Hayhoe, K.; VanDorn, J.; Croley, T., II; Schlegal, N.; Wuebbles, D. Regional climate change projections for Chicago and the US great lakes. J. Great Lakes Res. 2010, 36, 7–21. [Google Scholar] [CrossRef]
Figure 1. Map of Ottawa National Forest and black ash wetland study site locations within the Great Lakes region of North America (lower-left inset). Three geographically distinct locations (squares) were each subdivided into three treatment sites (diamonds). Regional precipitation data were obtained from the National Oceanic and Atmospheric Administration’s National Centers for Environmental Information (NOAA-NCEI) [44] station in Ironwood, Michigan (star).
Figure 1. Map of Ottawa National Forest and black ash wetland study site locations within the Great Lakes region of North America (lower-left inset). Three geographically distinct locations (squares) were each subdivided into three treatment sites (diamonds). Regional precipitation data were obtained from the National Oceanic and Atmospheric Administration’s National Centers for Environmental Information (NOAA-NCEI) [44] station in Ironwood, Michigan (star).
Forests 09 00324 g001
Figure 2. (A) Responses of daily maximum soil temperatures to daily maximum air temperatures and water levels; (B) Responses of soil CO2 fluxes to soil temperatures and water levels; (C) Responses of soil CH4 fluxes to soil temperatures and water levels. The lines in each panel represent the predicted [57] responses for four water level positions relative to the wetland soil surface. Observations (points) were color coordinated using equidistant binning for the four water level positions, and are used for visual interpretation only.
Figure 2. (A) Responses of daily maximum soil temperatures to daily maximum air temperatures and water levels; (B) Responses of soil CO2 fluxes to soil temperatures and water levels; (C) Responses of soil CH4 fluxes to soil temperatures and water levels. The lines in each panel represent the predicted [57] responses for four water level positions relative to the wetland soil surface. Observations (points) were color coordinated using equidistant binning for the four water level positions, and are used for visual interpretation only.
Forests 09 00324 g002
Figure 3. (A) Predicted soil CO2 flux response to soil temperature and water level ( W L c o l ) for each treatment, with shaded 95% confidence intervals [57]; (B) Predicted soil CH4 flux response to soil temperature and water level ( W L c o l ) for each treatment, with shaded 95% confidence intervals. W L c o l was embedded, using a correlated random intercept and slope, within each LMM as an interacting fixed effect. Significant differences between treatment means are indicated by (a) and (b) at p < 0.05, where (′) indicates detection of marginally significant between-treatment means differences at p < 0.1.
Figure 3. (A) Predicted soil CO2 flux response to soil temperature and water level ( W L c o l ) for each treatment, with shaded 95% confidence intervals [57]; (B) Predicted soil CH4 flux response to soil temperature and water level ( W L c o l ) for each treatment, with shaded 95% confidence intervals. W L c o l was embedded, using a correlated random intercept and slope, within each LMM as an interacting fixed effect. Significant differences between treatment means are indicated by (a) and (b) at p < 0.05, where (′) indicates detection of marginally significant between-treatment means differences at p < 0.1.
Forests 09 00324 g003
Figure 4. (A) Soil CO2 flux for each treatment during 2012, 2013, and 2014; (B) Soil CH4 flux for each treatment during 2013 and 2014. Significant between-treatment differences in mean flux at p > 0.05 are indicated by (a) and (b), where (′) indicates detection of marginally significant between-treatment mean differences at p < 0.1. Boxplot lines represent median, 25%, and 75% quantiles, while whiskers represent quantiles ± 1.5 ∙ interquantile range, respectively [58].
Figure 4. (A) Soil CO2 flux for each treatment during 2012, 2013, and 2014; (B) Soil CH4 flux for each treatment during 2013 and 2014. Significant between-treatment differences in mean flux at p > 0.05 are indicated by (a) and (b), where (′) indicates detection of marginally significant between-treatment mean differences at p < 0.1. Boxplot lines represent median, 25%, and 75% quantiles, while whiskers represent quantiles ± 1.5 ∙ interquantile range, respectively [58].
Forests 09 00324 g004
Table 1. Model evaluation statistics for soil temperature, CO2 flux, and CH4 flux linear mixed models (LMMs; Figure 2). Residual standard error (RSE) was reported for soil temperature for soil CO2 and CH4 flux LMMs. AIC, Akaike information criterion; PCV, proportion change in variance.
Table 1. Model evaluation statistics for soil temperature, CO2 flux, and CH4 flux linear mixed models (LMMs; Figure 2). Residual standard error (RSE) was reported for soil temperature for soil CO2 and CH4 flux LMMs. AIC, Akaike information criterion; PCV, proportion change in variance.
LMMPredictors r c 2 r m 2 AICPCVRSE
Soil Temp.Null: x = A T m a x 0.970.7317630.62
Full: x = A T m a x · W L m w 0.980.7312250.6%0.39
CO2 FluxNull: x = S T 5 c m 0.340.1415642.48
Full: x = S T 5 c m · W L c o l 0.510.3114890.2%2.27
CH4 FluxNull: x = S T 5 c m 0.490.056210.86
Full: x = S T 5 c m · W L c o l 0.500.106530.1%0.85
Table 2. Summary statistics, including degree of freedom (df), standard error (SE), Tukey honest significant difference (Mean Test), and Kruskal–Wallis rank means comparisons (KW Test), for each soil gas (CO2 and CH4), study period (pre- and posttreatment), and treatment.
Table 2. Summary statistics, including degree of freedom (df), standard error (SE), Tukey honest significant difference (Mean Test), and Kruskal–Wallis rank means comparisons (KW Test), for each soil gas (CO2 and CH4), study period (pre- and posttreatment), and treatment.
GasStudy PeriodTreatmentdfMean Flux (mg m−2 d−1)SE (mg m−2 d−1)Mean TestnRank MeanKW Test
CO2PreControl18.93.0 × 1031.26a1013.5a
Girdle16.54.0 × 1031.25a1016.5a
Ash-Cut30.95.1 × 1031.32a1016.5a
PostControl6.23.1 × 1031.18a2430.8a
Girdle7.55.7 × 1031.19b’2437.4ab
Ash-Cut8.76.7 × 1031.24b2441.3b’
CH4PostControl5.36.61.10ab2433.3a
Girdle5.52.41.11a2426.6a
Ash-Cut5.911.91.11b’2447.6b
Significant between-treatment mean differences are indicated by (a) and (b) at p < 0.05, where (′) indicates detection of marginally significant between-treatment mean differences at p < 0.1.

Share and Cite

MDPI and ACS Style

Van Grinsven, M.; Shannon, J.; Bolton, N.; Davis, J.; Noh, N.J.; Wagenbrenner, J.; Kolka, R.; Pypker, T. Response of Black Ash Wetland Gaseous Soil Carbon Fluxes to a Simulated Emerald Ash Borer Infestation. Forests 2018, 9, 324. https://doi.org/10.3390/f9060324

AMA Style

Van Grinsven M, Shannon J, Bolton N, Davis J, Noh NJ, Wagenbrenner J, Kolka R, Pypker T. Response of Black Ash Wetland Gaseous Soil Carbon Fluxes to a Simulated Emerald Ash Borer Infestation. Forests. 2018; 9(6):324. https://doi.org/10.3390/f9060324

Chicago/Turabian Style

Van Grinsven, Matthew, Joseph Shannon, Nicholas Bolton, Joshua Davis, Nam Jin Noh, Joseph Wagenbrenner, Randall Kolka, and Thomas Pypker. 2018. "Response of Black Ash Wetland Gaseous Soil Carbon Fluxes to a Simulated Emerald Ash Borer Infestation" Forests 9, no. 6: 324. https://doi.org/10.3390/f9060324

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop