Developing Management Guidelines that Balance Cattle and Timber Production with Ecological Interests in the Black Hills of South Dakota

Forested lands contribute to the United States (US) economy by providing livestock and timber production. Livestock grazing of forested lands has been widespread throughout the western US since the settlement era, and currently occurs on 51.4 million hectares (ha) representing 16% of all US grazing land and 22% of all US forested land (Nickerson et al. 2011). While livestock grazing and timber harvest are occurring on a substantial amount of forested land, relationships between management practices, tree stocking, timber production, forage production, livestock grazing, wildlife, aesthetics, and ecological integrity are not well documented. Whether considering timber or cattle, finding a balance between production and resource conservation is a fundamental challenge to agricultural producers, and is often a tradeoff between short term gains and long term sustainability. This dissertation aims to identify livestock and timber management practices that optimize production and are ecologically conservative. Specifically, I focused on three objectives. First, I reviewed the published literature and summarized what is known about best-practices for concurrent management of livestock and timber production in pine forests in the US. I found most studies came from the southeastern and western US where timber and livestock production on the same land unit are common. The relationship between pine cover and forage seemed fairly consistent across the US, and production was optimized when cattle grazed open canopy forests with basal areas between 5 and 14 m2 ha-1 (15-35% tree canopy cover). Second, I developed forest cover maps to estimate forage production in the Black Hills, South Dakota (SD) for the period from 1999 to 2015. I developed a regression model based on Landsat and Ikonos satellite imagery and was able to detect large changes in forest cover over time. I then used these maps in combination with maps of soil type and Palmer Drought Severity Index (PDSI) to update forage production estimates for the region. These changes in forest cover have large implications for forage production in the Black Hills. Over the 15 year period, mean tree cover decreased in 181 pastures in the Mystic Ranger District by 17.6 ± 0.6%, and there was a corresponding 15.5 ± 0.6% increase in mean forage production. Third, I conducted a 2 -year field experiment in the Black Hills, SD to study the relationships between management practices such as livestock stocking rates, grazing pressure, and timber harvest history, and aspects of resource condition such as tree regeneration, forage production, and plant community composition. From 2014-2015, I visited 44 pastures across a spectrum of management practices and measured seedling regeneration (590 plots), plant species richness (393 plots), primary production (246 plots), and visual obstruction (120 transects). I found that cattle grazing did not affect ponderosa pine regeneration. Grazing did affect plant diversity, and I found the highest plant diversity in areas of moderate grazing pressure. This work suggests that moderate stocking rates should have no effect on the timber industry but could positively affect native plant diversity. In the conclusion, I summarize what I learned from the literature review, mapping exercise, and field study and provide some management recommendations based on this work. Overall, I found that updated forage production estimates based on satellite imagery, and using grazing pressure index (GPI) to identify optimal stocking rates are tools that can facilitate management of livestock and timber production in the Black Hills, SD

Process-level controls on CO2 fluxes from a seasonally snow-covered subalpine meadow soil, Niwot Ridge, Colorado

Fluxes of CO2 during the snow-covered season contribute to annual carbon budgets, but our understanding of the mechanisms controlling the seasonal pattern and magnitude of carbon emissions in seasonally snow-covered areas is still developing. In a subalpine meadow on Niwot Ridge, Colorado, soil CO2 fluxes were quantified with the gradient method through the snowpack in winter 2006 and 2007 and with chamber measurements during summer 2007. The CO2 fluxes of 0.71 μmol m−2 s−1 in 2006 and 0.86 μmol m−2 s−1 in 2007 are among the highest reported for snow-covered ecosystems in the literature. These fluxes resulted in 156 and 189 g C m−2 emitted over the winter, ~30% of the annual soil CO2 efflux at this site. In general, the CO2 flux increased during the winter as soil moisture increased. A conceptual model was developed with distinct snow cover zones to describe this as well as the three other reported temporal patterns in CO2 flux from seasonally snow-covered soils. As snow depth and duration increase, the factor controlling the CO2 flux shifts from freeze–thaw cycles (zone I) to soil temperature (zone II) to soil moisture (zone III) to carbon availability (zone IV). The temporal pattern in CO2 flux in each zone changes from periodic pulses of CO2 during thaw events (zone I), to CO2 fluxes reaching a minimum when soil temperatures are lowest in mid-winter (zone II), to CO2 fluxes increasing gradually as soil moisture increases (zone III), to CO2 fluxes decreasing as available carbon is consumed. This model predicts that interannual variability in snow cover or directional shifts in climate may result in dramatically different seasonal patterns of CO2 flux from seasonally snow-covered soils

Process-level controls on CO2 fluxes from a seasonally snow-covered subalpine meadow soil, Niwot Ridge, Colorado

Fluxes of CO2 during the snow-covered season contribute to annual carbon budgets, but our understanding of the mechanisms controlling the seasonal pattern and magnitude of carbon emissions in seasonally snow-covered areas is still developing. In a subalpine meadow on Niwot Ridge, Colorado, soil CO2 fluxes were quantified with the gradient method through the snowpack in winter 2006 and 2007 and with chamber measurements during summer 2007. The CO2 fluxes of 0.71 μmol m−2 s−1 in 2006 and 0.86 μmol m−2 s−1 in 2007 are among the highest reported for snow-covered ecosystems in the literature. These fluxes resulted in 156 and 189 g C m−2 emitted over the winter, ~30% of the annual soil CO2 efflux at this site. In general, the CO2 flux increased during the winter as soil moisture increased. A conceptual model was developed with distinct snow cover zones to describe this as well as the three other reported temporal patterns in CO2 flux from seasonally snow-covered soils. As snow depth and duration increase, the factor controlling the CO2 flux shifts from freeze–thaw cycles (zone I) to soil temperature (zone II) to soil moisture (zone III) to carbon availability (zone IV). The temporal pattern in CO2 flux in each zone changes from periodic pulses of CO2 during thaw events (zone I), to CO2 fluxes reaching a minimum when soil temperatures are lowest in mid-winter (zone II), to CO2 fluxes increasing gradually as soil moisture increases (zone III), to CO2 fluxes decreasing as available carbon is consumed. This model predicts that interannual variability in snow cover or directional shifts in climate may result in dramatically different seasonal patterns of CO2 flux from seasonally snow-covered soils

Winter and summer nitrous oxide and nitrogen oxides fluxes from a seasonally snow-covered subalpine meadow at Niwot Ridge, Colorado

The soil emission rates (fluxes) of nitrous oxide (N2O) and nitrogen oxides (NO + NO2 = NO x ) through a seasonal snowpack were determined by a flux gradient method from near-continuous 2-year measurements using an automated system for sampling interstitial air at various heights within the snowpack from a subalpine site at Niwot Ridge, Colorado. The winter seasonal-averaged N2O fluxes of 0.047–0.069 nmol m−2 s−1 were ~15 times higher than observed NO x fluxes of 0.0030–0.0067 nmol m−2 s−1. During spring N2O emissions first peaked and then dropped sharply as the soil water content increased from the release of snowpack meltwater, while other gases, including NO x and CO2 did not show this behavior. To compare and contrast the winter fluxes with snow-free conditions, N2O fluxes were also measured at the same site in the summers of 2006 and 2007 using a closed soil chamber method. Summer N2O fluxes followed a decreasing trend during the dry-out period after snowmelt, interrupted by higher values related to precipitation events. These peaks were up to 2–3 times higher than the background summer levels. The integrated N2O-N loss over the summer period was calculated to be 1.1–2.4 kg N ha−1, compared to ~0.24–0.34 kg N ha−1 for the winter season. These wintertime N2O fluxes from subniveal soil are generally higher than the few previously published data. These results are of the same order of magnitude as data from more productive ecosystems such as fertilized grasslands and high-N-cycling forests, most likely because of a combination of the relatively well-developed soils and the fact that subnivean biogeochemical processes are promoted by the deep, insulating snowpack. Hence, microbially mediated oxidized nitrogen emissions occurring during the winter can be a significant part of the N-cycle in seasonally snow-covered subalpine ecosystems