Patrick R. Kormos

Hydrologic Vulnerability of Sagebrush Steppe Following Pinyon and Juniper Encroachment

Abstract Woodland encroachment on United States rangelands has altered the structure and function of shrub steppe ecosystems. The potential community structure is one where trees dominate, shrub and herbaceous species decline, and rock cover and bare soil area increase and become more interconnected. Research from the Desert Southwest United States has demonstrated areas under tree canopies effectively store water and soil resources, whereas areas between canopies (intercanopy) generate significantly more runoff and erosion. We investigated these relationships and the impacts of tree encroachment on runoff and erosion processes at two woodland sites in the Intermountain West, USA. Rainfall simulation and concentrated flow methodologies were employed to measure infiltration, runoff, and erosion from intercanopy and canopy areas at small-plot (0.5 m2) and large-plot (13 m2) scales. Soil water repellency and vegetative and ground cover factors that influence runoff and erosion were quantified. Runoff and erosion from rainsplash, sheet flow, and concentrated flow processes were significantly greater from intercanopy than canopy areas across small- and large-plot scales, and site-specific erodibility differences were observed. Runoff and erosion were primarily dictated by the type and quantity of ground cover. Litter offered protection from rainsplash effects, provided rainfall storage, mitigated soil water repellency impacts on infiltration, and contributed to aggregate stability. Runoff and erosion increased exponentially (r2  =  0.75 and 0.64) where bare soil and rock cover exceeded 50%. Sediment yield was strongly correlated (r2  =  0.87) with runoff and increased linearly where runoff exceeded 20 mm·h−1. Measured runoff and erosion rates suggest tree canopies represent areas of hydrologic stability, whereas intercanopy areas are vulnerable to runoff and erosion. Results indicate the overall hydrologic vulnerability of sagebrush steppe following woodland encroachment depends on the potential influence of tree dominance on bare intercanopy expanse and connectivity and the potential erodibility of intercanopy areas.

Hydrologic and Erosion Responses of Sagebrush Steppe Following Juniper Encroachment, Wildfire, and Tree Cutting

Abstract Extensive woodland expansion in the Great Basin has generated concern regarding ecological impacts of tree encroachment on sagebrush rangelands and strategies for restoring sagebrush steppe. This study used rainfall (0.5 m2 and 13 m2 scales) and concentrated flow simulations and measures of vegetation, ground cover, and soils to investigate hydrologic and erosion impacts of western juniper (Juniperus occidentalis Hook.) encroachment into sagebrush steppe and to evaluate short-term effects of burning and tree cutting on runoff and erosion responses. The overall effects of tree encroachment were a reduction in understory vegetation and formation of highly erodible, bare intercanopy between trees. Runoff and erosion from high-intensity rainfall (102 mm · h−1, 13 m2 plots) were generally low from unburned areas underneath tree canopies (13 mm and 48 g · m−2) and were higher from the unburned intercanopy (43 mm and 272 g · m−2). Intercanopy erosion increased linearly with runoff and exponentially where bare ground exceeded 60%. Erosion from simulated concentrated flow was 15- to 25-fold greater from the unburned intercanopy than unburned tree canopy areas. Severe burning amplified erosion from tree canopy plots by a factor of 20 but had a favorable effect on concentrated flow erosion from the intercanopy. Two years postfire, erosion remained 20-fold greater on burned than unburned tree plots, but concentrated flow erosion from the intercanopy (76% of study area) was reduced by herbaceous recruitment. The results indicate burning may amplify runoff and erosion immediately postfire. However, we infer burning that sustains residual understory cover and stimulates vegetation productivity may provide long-term reduction of soil loss relative to woodland persistence. Simply placing cut-downed trees into the unburned intercanopy had minimal immediate impact on infiltration and soil loss. Results suggest cut-tree treatments should focus on establishing tree debris contact with the soil surface if treatments are expected to reduce short-term soil loss during the postcut understory recruitment period.

Hydrologic Response to Mechanical Shredding in a Juniper Woodland

Abstract We investigated soil compaction and hydrologic responses from mechanically shredding Utah juniper (Juniperus ostesperma [Torr.] Little) to control fuels in a sagebrush/bunchgrass plant community (Artemisia nova A. Nelson, Artemisia tridentata Nutt. subsp. wyomingensis Beetle & Young/Pseudoroegneria spicata [Pursh] A. Löve, Poa secunda J. Presl) on a gravelly loam soil with a 15% slope in the Onaqui Mountains of Utah. Rain simulations were applied on 0.5-m2 runoff plots at 64 mm · h−1 (dry run: soil initially dry) and 102 mm · h−1 (wet run: soil initially wet). Runoff and sediment were collected from runoff plots placed in five blocks, each containing four microsites (juniper mound, shrub mound, vegetation-free or bare interspace, and grass interspace) with undisturbed or tracked treatments for each microsite type and a residue-covered treatment for grass and bare interspace microsites. Soil penetration resistance was measured at the hill slope scale, and canopy and ground cover were measured at the hill slope and runoff plot scale. Although shredding trees at a density of 453 trees · ha−1 reduced perennial foliar cover by 20.5%, shredded tree residue covered 40% of the ground surface and reduced non–foliar-covered bare ground and rock by 17%. Tire tracks from the shredding operation covered 15% of the hill slope and increased penetration resistance. For the wet run, infiltration rates of grass interspaces were significantly decreased (39.8 vs. 66.1 mm · h−1) by tire tracks, but infiltration rates on juniper mounds and bare interspaces were unchanged. Bare interspace plots covered with residue had significantly higher infiltration rates (81.9 vs. 26.7 mm · h−1 ) and lower sediment yields (38.6 vs. 313 g · m−2 ) than those without residue. Because hydrologic responses to treatments are site- and scale-dependent, determination of shredding effects on other sites and at hill slope or larger scales will best guide management actions.

Ecosystem Water Availability in Juniper versus Sagebrush Snow-Dominated Rangelands

ABSTRACT Western Juniper (Juniperus occidentalis Hook.) has greatly expanded in the past 150+years and now dominates over 3.6 million ha of rangeland in the Intermountain Western United States. The impacts of juniper encroachment on critical ecohydrological relationships among snowdistribution, water budgets, plant community transitions, and habitat requirements for wildlife, such as the greater sage grouse (Centrocercus urophasianus), remain poorly understood. The goal of this study is to better understand how juniper encroachment affects water availability for ecohydrologic processes and associatedwildlife habitat in snow-dominated sagebrush (Artemisia spp.) steppe ecosystems. A 6-yr combined measurement and modeling study is conducted to explore differences in snow distribution, water availability, and annual water balances between juniper-dominated and sagebrushdominated catchments. Although there is large interannual variability in both measured weather data and modeled hydrologic fluxes during the study, results indicate that juniper-dominated catchments have greater peak accumulations of snow water equivalent, earlier snow melt, and less streamflow relative to sagebrushdominated catchments. Water delivery is delayed by an average of 9 days in the sagebrush-dominated scenario comparedwith the juniper-dominated scenario as a result of increasedwater storage in snow drifts. The delayed water input to sagebrush-dominated ecosystems in typical water years has wide-ranging implications for available surface water, soil water, and vegetation dynamics associated with wildlife habitat for sagebrush obligates such as sage grouse. Results from this study imply that the retention of high-elevation, sagebrush-dominated landscapes may become crucial for sage grouse habitat management if mid- and low-elevation precipitation continues to transition from snow to rain dominated.

Short-Term Effects of Tree Removal on Infiltration, Runoff, and Erosion in Woodland-Encroached Sagebrush Steppe

Abstract Land owners and managers across the western United States are increasingly searching for methods to evaluate and mitigate the effects of woodland encroachment on sagebrush steppe ecosystems. We used small-plot scale (0.5 m2) rainfall simulations and measures of vegetation, ground cover, and soils to investigate woodland response to tree removal (prescribed fire and mastication) at two late-succession woodlands. We also evaluated the effects of burning on soil water repellency and effectiveness of aggregate stability indices to detect changes in erosion potential. Plots were located in interspaces between tree and shrub canopies and on undercanopy tree and shrub microsites. Erosion from untreated interspaces in the two woodlands differed more than 6-fold, and erosion responses to prescribed burning differed by woodland site. High-intensity rainfall (102 mm · h−1) on the less erodible woodland generated amplified runoff and erosion from tree microsites postfire, but erosion (45–75 g · m−2) was minor relative to the 3–13-fold fire-induced increase in erosion on tree microsites at the highly erodible site (240–295 g · m−2). Burning the highly erodible woodland also generated a 7-fold increase in erosion from shrub microsites (220–230 g · m−2) and 280–350 g · m−2 erosion from interspaces. High levels of runoff (40–45 mm) and soil erosion (230–275 g · m−2) on unburned interspaces at the more erodible site were reduced 4–5-fold (10 mm and 50 g · m−2) by masticated tree material. The results demonstrate that similarly degraded conditions at woodland-encroached sites may elicit differing hydrologic and erosion responses to treatment and that treatment decisions should consider inherent site-specific erodibility when evaluating tree-removal alternatives. Strong soil water repellency was detected from 0 cm to 3 cm soil depth underneath unburned tree canopies at both woodlands and its strength was not altered by burning. However, fire removal of litter exacerbated repellency effects on infiltration, runoff generation, and erosion. The aggregate stability index method detected differences in relative soil stability between areas underneath trees and in the intercanopy at both sites, but failed to provide any indication of between-site differences in erodibility or the effects of burning on soil erosion potential.

Trends and sensitivities of low streamflow extremes to discharge timing and magnitude in Pacific Northwest mountain streams

Path analyses of historical streamflow data from the Pacific Northwest indicate that the precipitation amount has been the dominant control on the magnitude of low streamflow extremes compared to the air temperature-affected timing of snowmelt runoff. The relative sensitivities of low streamflow to precipitation and temperature changes have important implications for adaptation planning because global circulation models produce relatively robust estimates of air temperature changes but have large uncertainties in projected precipitation amounts in the Pacific Northwest U.S. Quantile regression analyses indicate that low streamflow extremes from the majority of catchments in this study have declined from 1948 to 2013, which may significantly affect terrestrial and aquatic ecosystems, and water resource management. Trends in the 25th percentile of mean annual streamflow have declined and the center of timing has occurred earlier. We quantify the relative influences of total precipitation and air temperature on the annual low streamflow extremes from 42 stream gauges using mean annual streamflow as a proxy for precipitation amount effects and streamflow center of timing as a proxy for temperature effects on low flow metrics, including 7q10 summer (the minimum 7 day flow during summer with a 10 year return period), mean August, mean September, mean summer, 7q10 winter, and mean winter flow metrics. These methods have the benefit of using only readily available streamflow data, which makes our results robust against systematic errors in high elevation distributed precipitation data. Winter low flow metrics are weakly tied to both mean annual streamflow and center of timing.

RISK ASSESSMENT OF EROSION FROM CONCENTRATED FLOW ON RANGELANDS USING OVERLAND FLOW DISTRIBUTION AND SHEAR STRESS PARTITIONING

Abstract. Erosion rates of overland flow on rangelands tend to be relatively low, but under certain conditions where flow is concentrated, soil loss can be significant. Therefore, a rangeland site can be highly vulnerable to soil erosion where overland flow is likely to concentrate and exert high shear stress on soil grains. This concept is commonly applied in cropland and wildland soil erosion modeling using predictions of flow effective shear stress (shear stress applied on soil grains). However, historical approaches to partition shear stress in erosion models are computationally complex and require extensive parameterization. Furthermore, most models are not capable of predicting the conditions in which concentrated flow occurs on rangelands. In this study, we investigated the rangelands conditions at which overland flow is more likely to become concentrated and developed equations for partitioning the shear stress of concentrated flow on rangelands. A logistic equation was developed to estimate the probability of overland flow to become concentrated. Total shear stress of rangeland overland flow was partitioned into components exerted on soil, vegetation, and rock cover using field experimental data. In addition, we investigated the vegetation cover limit at which the effective shear stress component is substantially reduced, limiting the erosion rate. The results from the partitioning equations show that shear stress exerted on soil grains was relatively small in sheet flow. Shear stress exerted on soil grains in concentrated flow was significantly higher when bare soil exceeded 60% of the total surface area but decreased significantly when the bare soil area was less than 25% or when the plant base cover exceeded 20%. These percentages could be used as relative measures of hydrologic recovery for disturbed rangelands or as triggers that indicate a site is crossing a threshold beyond which soil erosion might accelerate due to the high effective shear stress.

Modeling the isotopic evolution of snowpack and snowmelt : Testing a spatially distributed parsimonious approach

Abstract Use of stable water isotopes has become increasingly popular in quantifying water flow paths and travel times in hydrological systems using tracer‐aided modeling. In snow‐influenced catchments, snowmelt produces a traceable isotopic signal, which differs from original snowfall isotopic composition because of isotopic fractionation in the snowpack. These fractionation processes in snow are relatively well understood, but representing their spatiotemporal variability in tracer‐aided studies remains a challenge. We present a novel, parsimonious modeling method to account for the snowpack isotope fractionation and estimate isotope ratios in snowmelt water in a fully spatially distributed manner. Our model introduces two calibration parameters that alone account for the isotopic fractionation caused by sublimation from interception and ground snow storage, and snowmelt fractionation progressively enriching the snowmelt runoff. The isotope routines are linked to a generic process‐based snow interception‐accumulation‐melt model facilitating simulation of spatially distributed snowmelt runoff. We use a synthetic modeling experiment to demonstrate the functionality of the model algorithms in different landscape locations and under different canopy characteristics. We also provide a proof‐of‐concept model test and successfully reproduce isotopic ratios in snowmelt runoff sampled with snowmelt lysimeters in two long‐term experimental catchment with contrasting winter conditions. To our knowledge, the method is the first such tool to allow estimation of the spatially distributed nature of isotopic fractionation in snowpacks and the resulting isotope ratios in snowmelt runoff. The method can thus provide a useful tool for tracer‐aided modeling to better understand the integrated nature of flow, mixing, and transport processes in snow‐influenced catchments.

Short-Term Impacts of Tree Removal on Runoff and Erosion From Pinyon- and Juniper-Dominated Sagebrush Hillslopes

abstract Tree removal is often applied to woodland-encroached rangelands to restore vegetation and improve hydrologic function, but knowledge is limited regarding effects of tree removal on hydrologic response. This study used artificial rainfall and overland flow experiments (9–13 m2) and measures of vegetation and ground cover to investigate short-term (1–2 yr) responses to tree removal at two woodland-encroached sites. Plots were located under trees (tree zone) and in the intercanopy (shrub-interspace zone, 75% of area). Before tree removal, vegetation and ground cover were degraded and intercanopy runoff and erosion rates were high. Cutting and placing trees into the intercanopy did not significantly affect vegetation, ground cover, runoff, or erosion 1 yr posttreatment. Whole-tree mastication as applied in this study did not redistribute tree mulch within the intercanopy, but the treatment did result in enhanced herbaceous cover and hydrologic function in the intercanopy. Fire removal of litter and herbaceous cover increased tree-zone runoff and erosion under high-intensity rainfall by 4- and 30-fold at one site but had minimal impact at the other site. Site response differences were attributed to variability in burn conditions and site-specific erodibility. Burning had minimal impact on shrub-interspace runoff and erosion from applied high-intensity rainfall. However, 1 yr postfire, erosion from concentrated overland flow experiments was 2- to 13-fold greater on burned than unburned tree-zone and shrub-interspace plots and erosion for burned tree zones was 3-fold greater for the more erodible site. Two yr postfire, overland flow erosion remained higher for burned versus unburned tree zones, but enhanced intercanopy herbaceous cover reduced erosion from shrub-interspace zones. The net impact of burning included an initial increase in erosion risk, particularly for tree zones, followed by enhanced herbaceous cover and improved hydrologic function within the intercanopy. The overall results suggest that erosion from late-succession woodlands is reduced primarily through recruitment of intercanopy herbaceous vegetation and ground cover.

Shear Stress Partitioning of Overland Flow on Disturbed and Undisturbed Rangelands

Physically-based hillslope erosion models commonly estimate soil detachment and transport capacity based on overland flow shear stress applied to soil aggregates. However, vegetation and rock cover counteract the shear stress of overland flow where they occur. Accordingly, partitioning of total shear stress into components exerted on soil, vegetation, and rock cover is a key element for the erosion models. The objective of this study is to estimate the components of shear stress of overland flow on disturbed and undisturbed rangelands using field experimental data. In addition, this study investigates the vegetation cover limit at which the soil shear stress component is substantially reduced, limiting the erosion rate. The soil shear stress component was estimated based on the assumption that the ratio of soil shear stress to the total shear stress is equal to the ratio of hydraulic friction factor of soil to the friction factor of the composite surface. The total friction factor of the composite surface was estimated using empirical equations developed based on field experimental data over diverse rangeland landscapes within the Great Basin Region, United States. This equation logarithmically correlates the composite surface friction to the vegetation cover (plant base and plant litter) and rock cover components. Moreover, the hydraulic friction factor of each cover element was estimated based on its parameter in that equation. The soil hydraulic friction portion was assumed to be the logarithmic difference between the total friction and the friction of the cover elements. The result of this assumption was used to develop empirical equations that predict the ratio of soil shear stress to the total shear stress of concentrated flow and sheet flow in terms of bare soil fraction of total area. The predicting equation of total friction factor was improved by adding the slope and the flow discharge variables. The predicting equations of soil shear stress as a function of bare soil fraction did not change significantly when changing the assumption of a rectangular shape of cross section to a parabolic shape. The developed shear stress partitioning equations in this study are applicable across a wide span of ecological sites, soils, slopes, and vegetation and ground cover conditions and can be used by physically-based rangeland hydrology and erosion models. The results from the developed equations show that shear stress exerted on soil grains is significantly higher when bare soil exceeds 60% of the total surface area, while reduced significantly when bare soil area is less than 25% or when the plant base cover exceeds 20%. These percentages could be used as relative measures of hydrologic recovery for disturbed rangelands or triggers that indicate that a site is crossing a threshold where soil erosion might accelerate due to the high soil shear stress.