Influence of grass soil cover on water runoff and soil detachment under rainfall simulation in a sub‐humid South African degraded rangeland

In most regions of the world overgrazing plays a major role in land degradation and thus creates a major threat to natural ecosystems. Several feedbacks exist between overgrazing, vegetation, soil infiltration by water and soil erosion that need to be better understood. In this study of a sub‐humid overgrazed rangeland in South Africa, the main objective was to evaluate the impact of grass cover on soil infiltration by water and soil detachment. Artificial rains of 30 and 60 mm h^?1 were applied for 30 min on 1 m² micro‐plots showing similar sandy‐loam Acrisols with different proportions of soil surface coverage by grass (Class A: 75–100%; B: 75–50%; C: 50–25%; D: 25–5%; E: 5–0% with an outcropping A horizon; F: 0% with an outcropping B horizon) to evaluate pre‐runoff rainfall (Pr), steady state water infiltration (I), sediment concentration (SC) and soil losses (SL). Whatever the class of vegetal cover and the rainfall intensity, with the exception of two plots probably affected by biological activity, I decreased regularly to a steady rate <2 mm h^?1 after 15 min rain. There was no significant correlation between I and Pr with vegetal cover. The average SC computed from the two rains increased from 0·16 g L^?1 (class A) to 48·5 g L^?1 (class F) while SL was varied between 4 g m^?2 h^?1 for A and 1883 g m^?2 h^?1 for F. SL increased significantly with decreasing vegetal cover with an exponential increase while the removal of the A horizon increased SC and SL by a factor of 4. The results support the belief that soil vegetation cover and overgrazing plays a major role in soil infiltration by water but also suggest that the interrill erosion process is self‐increasing. Abandoned cultivated lands and animal preferred pathways are more vulnerable to erosive processes than simply overgrazed rangelands. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling of the effect of flow depth on sediment discharged by rain‐impacted flows from sheet and interrill erosion areas: a review

Sediment, nutrients and pollutants discharged from sheet and interrill erosion areas by rain‐impacted flows may influence water quality in streams and rivers. The depth of water on the soil surface influences the capacity of raindrop impacts to detach soil material underlying rain‐impacted flows, and a number of so‐called process‐based and mechanistic models erroneously use equations on the basis of the effect of water depth on splash erosion to account for this effect. Also, a number of these models require complex mathematical solutions to make them operate and can only predict sediment composition and discharges well if many of their parameters are calibrated specifically to the situations where they are being applied. Experiments with rain‐impacted flows, where flow depth and velocity over eroding surfaces have been controlled, have been reported in the literature and provide more appropriate equations to account for the drop size – flow depth interactions that affect detachment and transport of particles in rain‐impacted flows. There is a need to develop modeling approaches that rely on relevant data obtained under well‐controlled flow conditions where flow depths and velocities are known. Copyright © 2012 John Wiley & Sons, Ltd.     

Laboratory soil piping and internal erosion experiments: evaluation of a soil piping model for low‐compacted soils

Mechanistic models have been proposed for soil piping and internal erosion on well‐compacted levees and dams, but limited research has evaluated these models in less compacted (more erodible) soils typical of hillslopes and streambanks. This study utilized a soil box (50 cm long, 50 cm wide and 20 cm tall) to conduct constant‐head, soil pipe and internal erosion experiments for two soils (clay loam from Dry Creek and sandy loam from Cow Creek streambanks) packed at uniform bulk densities. Initial gravimetric moisture contents prior to packing were 10, 12 and 14% for Dry Creek soil and 8, 12, and 14% for Cow Creek soil. A 1‐cm diameter rod was placed horizontally along the length of the soil bed during packing and carefully removed after packing to create a continuous soil pipe. A constant head was maintained at the inflow end. Flow rates and sediment concentrations were measured from the pipe outlet. Replicate submerged jet erosion tests (JETs) were conducted to derive erodibility parameters for repacked samples at the same moisture contents. Flow rates from the box experiments were used to calibrate the mechanistic model. The influence of the initial moisture content was apparent, with some pipes (8% moisture content) expanding so fast that limited data was collected. The mechanistic model was able to estimate equivalent flow rates to those observed in the experiments, but had difficulty matching observed sediment concentrations when the pipes rapidly expanded. The JETs predicted similar erodibility coefficients compared to the mechanistic model for the more erodible cases but not for the less erodible cases (14% moisture content). Improved models are needed that better define the changing soil pipe cross‐section during supply‐ and transport‐limited internal erosion, especially for piping through lower compacted (more erodible) soils as opposed to more well‐compacted soils resulting from constructing levees and dams. Copyright © 2013 John Wiley & Sons, Ltd.     

Raindrop‐impact‐induced erosion processes and prediction: a review

Raindrop‐impact‐induced erosion is initiated when detachment of soil particles from the surface of the soil results from an expenditure of raindrop energy. Once detachment by raindrop impact has taken place, particles are transported away from the site of the impact by one or more of the following transport processes: drop splash, raindrop‐induced flow transport, or transport by flow without stimulation by drop impact. These transport processes exhibit varying efficiencies. Particles that fall back to the surface as a result of gravity produce a layer of pre‐detached particles that provides a degree of protection against the detachment of particles from the underlying soil. This, in turn, influences the erodibility of the eroding surface. Good understanding of rainfall erosion processes is necessary if the results of erosion experiments are to be properly interpreted. Current process‐based erosion prediction models do not deal with the issue of temporal variations in erodibility during a rainfall event or variabilities in erodibility associated with spatial changes in dominance of the transport processes that follow detachment by drop impact. Although more complex erosion models may deal with issues like this, their complexity and high data requirement may make them unsuitable for use as general prediction tools. Copyright © 2005 John Wiley & Sons, Ltd.     

Development of a rain‐on‐snow detection algorithm using passive microwave radiometry

Currently observed climate warming in the Arctic has numerous consequences. Of particular relevance, the precipitation regime is modified where mixed and liquid precipitation can occur during the winter season leading to rain‐on‐snow (ROS) events. This phenomenon is responsible for ice crust formation, which has a significant impact on ecosystems (such as biological, hydrological, ecological and physical processes). The spatially and temporally sporadic nature of ROS events makes the phenomenon difficult to monitor using meteorological observations. This paper focuses on the detection of ROS events using passive microwave (PMW) data from a modified brightness temperature (T_B) gradient approach at 19 and 37 GHz. The approach presented here was developed empirically for observed ROS events with coincident ground‐based PMW measurements in Sherbrooke, Quebec, Canada. It was then tested in Nunavik, Quebec, with the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR‐E). We obtained a detection accuracy of 57, 71 and 89% for ROS detection for three AMSR‐E grid cells with a maximum error of 7% when considering all omissions and commissions with regard to the total number of AMSR‐E passes throughout the winter period. Copyright © 2016 John Wiley & Sons, Ltd.     

Role of biological soil crust cover in bioweathering and protection of sandstones in a semi‐arid landscape (Torrollones de Gabarda,Huesca, Spain)

Sandstone structural landscapes in the semi‐arid Torrollones de Gabarda area (Province of Huesca, NE Spain) are often covered by a well developed biological soil crust of lichens, mosses and cyanobacteria and black coatings on vertical surfaces. By using scanning electron microscopy with backscattered detector imaging, the biological soil crust studied evidenced high activity in the sandstone–crust interface. Processes such as physical disintegration, etching and dwelling as well as biomineralization by calcium oxalate and ?xation of mineral particles by extracellular polymeric substances were observed. On the horizontal sandstone surfaces these processes may cause the occurrence of gnammas and the development of a protective coating that favours intense ?aking when the crust is disturbed. On the sandstone cliffs, columnar and tafoni weathering development is clearly guided by the protective action of the biological soil crust. These qualitative observations are important to develop methodologies to address their quantitative importance in geomorphological processes in semi‐arid landscapes. Copyright © 2004 John Wiley & Sons, Ltd.     

The Daily Erosion Project – daily estimates of water runoff,soil detachment,and erosion

Water runoff and sediment transport from agricultural uplands are substantial threats to water quality and sustained crop production. To improve soil and water resources, farmers, conservationists, and policy‐makers must understand how landforms, soil types, farming practices, and rainfall interact with water runoff and soil erosion processes. To that end, the Iowa Daily Erosion Project (IDEP) was designed and implemented in 2003 to inventory these factors across Iowa in the United States. IDEP utilized the Water Erosion Prediction Project (WEPP) soil erosion model along with radar‐derived precipitation data and government‐provided slope, soil, and management information to produce daily estimates of soil erosion and runoff at the township scale (93 km² [36 mi²]). Improved national databases and evolving remote sensing technology now permit the derivation of slope, soil, and field‐level management inputs for WEPP. These remotely sensed parameters, along with more detailed meteorological data, now drive daily WEPP hillslope soil erosion and water runoff estimates at the small watershed scale, approximately 90 km² (35 mi²), across sections of multiple Midwest states. The revisions constitute a substantial improvement as more realistic field conditions are reflected, more detailed weather data are utilized, hill slope sampling density is an order of magnitude greater, and results are aggregated based on surface hydrology enabling further watershed research and analysis. Considering these improvements and the expansion of the project beyond Iowa it was renamed the Daily Erosion Project (DEP). Statistical and comparative evaluations of soil erosion simulations indicate that the sampling density is adequate and the results are defendable. The modeling framework developed is readily adaptable to other regions given suitable inputs. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Soil detachment and transport on field‐ and laboratory‐scale interrill areas: erosion processes and the size‐selectivity of eroded sediment

Field‐ and laboratory‐scale rainfall simulation experiments were carried out in an investigation of the temporal variability of erosion processes on interrill areas, and the effects of such variation upon sediment size characteristics. Poorly aggregated sandy soils from the semi‐arid environment of Senegal, West Africa, were used on both a 40 m² field plot and a 0·25 m² laboratory plot; rainfall intensity for all experiments was 70 mm h^?1 with a duration of 1 to 2 hours. Time‐series measurements were made of the quantity and the size distribution of eroded material: these permitted an estimate of the changing temporal balance between the main erosion processes (splash and wash). Results from both spatial scales showed a similar temporal pattern of runoff generation and sediment concentration. For both spatial scales, the dominant erosional process was detachment by raindrops; this resulted in a dynamic evolution of the soil surface under raindrop impact, with the rapid formation of a sieving crust followed by an erosion crust. However, a clear difference was observed between the two scales regarding the size of particles detached by both splash and wash. While all measured values were lower than the mean weight diameter (MWD) value of the original soil (mean 0·32 mm), demonstrating the size‐selective nature of wash and splash processes, the MWD values of washed and splashed particles at the field scale ranged from 0·08 to 0·16 mm and from 0·12 to 0·30 mm respectively, whereas the MWD values of washed and splashed particles at the laboratory scale ranged from 0·13 to 0·29 mm and from 0·21 to 0·32 mm respectively. Thus only at the field scale were the soil particles detached by splash notably coarser than those transported by wash. This suggests a transport‐limited erosion process at the field scale. Differences were also observed between the dynamics of the soil loss by wash at the two scales, since results showed wider scatter in the field compared to the laboratory experiments. This scatter is probably related to the change in soil surface characteristics due to the size‐selectivity of the erosion processes at this spatial scale. Copyright © 2006 John Wiley & Sons, Ltd.     

SCALES: a large‐scale assessment model of soil erosion hazard in Basse‐Normandie (northern‐western France)

The cartography of erosion risk is mainly based on the development of models, which evaluate in a qualitative and quantitative manner the physical reproduction of the erosion processes (CORINE, EHU, INRA). These models are mainly semi‐quantitative but can be physically based and spatially distributed (the Pan‐European Soil Erosion Risk Assessment, PESERA). They are characterized by their simplicity and their applicability potential at large temporal and spatial scales. In developing our model SCALES (Spatialisation d'éChelle fine de l'ALéa Erosion des Sols/large‐scale assessment and mapping model of soil erosion hazard), we had in mind several objectives: (1) to map soil erosion at a regional scale with the guarantee of a large accuracy on the local level, (2) to envisage an applicability of the model in European oceanic areas, (3) to focus the erosion hazard estimation on the level of source areas (on‐site erosion), which are the agricultural parcels, (4) to take into account the weight of the temporality of agricultural practices (land‐use concept). Because of these objectives, the nature of variables, which characterize the erosion factors and because of its structure, SCALES differs from other models. Tested in Basse‐Normandie (Calvados 5500 km²) SCALES reveals a strong predisposition of the study area to the soil erosion which should require to be expressed in a wet year. Apart from an internal validation, we tried an intermediate one by comparing our results with those from INRA and PESERA. It appeared that these models under estimate medium erosion levels and differ in the spatial localization of areas with the highest erosion risks. SCALES underlines here the limitations in the use of pedo‐transfer functions and the interpolation of input data with a low resolution. One must not forget however that these models are mainly focused on an interregional comparative approach. Therefore the comparison of SCALES data with those of the INRA and PESERA models cannot result on a convincing validation of our model. For the moment the validation is based on the opinion of local experts, who agree with the qualitative indications delivered by our cartography. An external validation of SCALES is foreseen, which will be based on a thorough inventory of erosion signals in areas with different hazard levels. Copyright © 2010 John Wiley & Sons, Ltd.     

Induced biological soil crust controls on wind erodibility and dust (PM10) emissions

Inducing biological soil crust (biocrust) development is an appealing approach for dust mitigation in drylands due to the resistance biocrusts can provide against erosion. Using a portable device, we evaluated dust emissions from surfaces either inoculated with biocrust, amended with a plant-based soil stabilizer, or both at varying wind friction velocities. Four months after application, emissions from all treatments were either indistinguishable from or greater than controls, despite evidence of biocrust establishment. All treatments had greater surface roughness and showed more evidence of entrapment of windblown sediment than controls, factors which may have been partially responsible for elevated emissions. There was a synergistic effect of inoculation and stabilizer addition, resulting in a nearly two-fold reduction in estimated emissions compared to either treatment alone. Stepwise regression analysis indicated that variables associated with surface crust strength (aggregate stability, penetration resistance) were negatively associated with emissions and variables associated with sediment supply (sand content, loose sediment cover) were positively associated with emissions. With more time to develop, the soil-trapping activity and surface integrity of biocrust inoculum and soil stabilizer mixtures is expected to increase with the accumulation of surface biomass and enhancement of roughness through freeze–thaw cycles. © 2019 John Wiley & Sons, Ltd.     

Simulations demonstrating interaction between coarse and fine sediment loads in rain‐impacted flow

Rain‐impacted flows dominate sheet and interrill erosion and are important in eroding soil rich in nutrients and other chemicals which may have deleterious effects on water quality. Erosion in rain‐impacted flow is associated with raindrop detachment followed by transport either by the combination of flow velocity and raindrop impact (raindrop‐induced flow transport, RIFT) or the inherent capacity of the flow to transport detached material. Coarse particles tend to be transported by RIFT, while fine particles tend to be transported without any assistance from raindrop impact. Because the transport process associated with coarse particles is not 100 per cent efficient, it generates a layer of loose particles on the soil surface and this layer protects the underlying soil from detachment. Simulations were performed by modelling the uplift and downstream movement of both fine and coarse particles detached from the soil surface by individual raindrop impacts starting with a surface where no loose material was present. The simulations produced a flush of fine material followed by a decline in the discharge of fine material as the amount of loose material built up on the bed. The decline in the discharge of fine material was accompanied by an increase in the discharge of coarse material. The relative amounts of coarse and fine material discharged in the flow varied with flow velocity and cohesion in the surface of the soil matrix. The results indicate that the discharge of various sized sediments is highly dependent on local soil, rain and flow conditions and that extrapolating the results from one situation to another may not be appropriate. Copyright © 2006 John Wiley & Sons, Ltd.     

Spatio‐temporal controls of snowmelt and runoff generation during rain‐on‐snow events in a mid‐latitude mountain catchment

A network of 30 standalone snow monitoring stations was used to investigate the snow cover distribution, snowmelt dynamics, and runoff generation during two rain‐on‐snow (ROS) events in a 40 km² montane catchment in the Black Forest region of southwestern Germany. A multiple linear regression analysis using elevation, aspect, and land cover as predictors for the snow water equivalent (SWE) distribution within the catchment was applied on an hourly basis for two significant ROS flood events that occurred in December 2012. The available snowmelt water, liquid precipitation, as well as the total retention storage of the snow cover were considered in order to estimate the amount of water potentially available for the runoff generation. The study provides a spatially and temporally distributed picture of how the two observed ROS floods developed in the catchment. It became evident that the retention capacity of the snow cover is a crucial mechanism during ROS. It took several hours before water was released from the snowpack during the first ROS event, while retention storage was exceeded within 1 h from the start of the second event. Elevation was the most important terrain feature. South‐facing terrain contributed more water for runoff than north‐facing slopes, and only slightly more runoff was generated at open compared to forested areas. The results highlight the importance of snowmelt together with liquid precipitation for the generation of flood runoff during ROS and the large temporal and spatial variability of the relevant processes. Copyright © 2015 John Wiley & Sons, Ltd.     

More rain,less soil: long‐term changes in rainfall intensity with climate change

This commentary discusses the role of long‐term climate change in driving increases in soil erosion. Assuming that land use and management remain effectively constant, we discuss changes in the ability of rainfall to cause erosion (erosivity), using long daily rainfall data sets from southeast England. An upward trend in mean rainfall per rain day is detected at the century‐plus timescale. Implications for soil erosion and sediment delivery are discussed and evidence from other regions reviewed. We conclude that rates of soil erosion may well increase in a warmer, wetter world. Copyright © 2015 John Wiley & Sons, Ltd.     

Spatial and temporal variations in soil erosion and deposition due to land‐levelling in a semi‐arid area of Basilicata (Southern Italy)

This paper examines the changes from 1955 to 2002 in soil erosion and deposition due to changes in land‐use patterns in the semi‐arid territory of Craco, which is characterized by landsliding and badland erosion. The area underwent continuous degradation during the last century due not only to its lithological vulnerability but also to the anthropic pressure favoured by the introduction of Common Agricultural Policy (CAP) measures, which has led to the reclamation of scrub lands and badlands for durum wheat cultivation. Our analysis integrates the Unit Stream Power Erosion Deposition (USPED) model with a geographic information system (GIS) to quantify erosion risk and predict deposition patterns. Soil data, land use inventory, digital elevation data and climatic atlases were used as resource data sets to generate USPED factor values. The obtained results correlate well with field measured erosion data by other researchers. In the investigated 47 years, stable areas decreased by about 280 ha (3·8% of the total surface area), largely attributable to the increase of the low and moderate erosion intensity without significant change in sedimentation. Results from this study have implications related to understanding the geomorphic response of sites that were abandoned following remodelling due to the application of the F measure of Regulation CEE 2078/92. The average annual erosion rates estimated for abandoned and remodelled sites are respectively 15·99 and 10·64 t ha^?1, meaning that the total amount of erosion in 20 years could be estimated at around 100 t ha^?1. Copyright © 2007 John Wiley & Sons, Ltd.     

Magnetic enhancement in wildfire‐affected soil and its potential for sediment‐source ascription

Intense rainfall following wildfire can cause substantial soil and sediment redistribution. With concern for the increasing magnitude and frequency of wildfire events, research needs to focus on hydrogeomorphological impacts of fire, particularly downstream fluxes of sediment and nutrients. Here, we investigate variation in magnetic enhancement of soil by fire in burnt eucalypt forest slopes to explore its potential as a post‐fire sediment tracer. Low‐frequency magnetic susceptibility values (χ_lf) of <10 µm material sourced from burnt slopes (c. 8·0–10·4 × 10^?6 m³ kg^?1) are an order of magnitude greater than those of <10 µm material derived from long‐unburnt areas (0·8 × 10^?6 m³ kg^?1). Susceptibility of anhysteretic remanent magnetization (χARM) and saturation isothermal remanent magnetization (SIRM) values are similarly enhanced. Signatures are strongly influenced by soil and sediment particle size and storage of previously burnt material in footslope areas. Whilst observations indicate that signatures based on magnetic enhancement show promise for post‐fire sediment tracing, problems arise with the lack of dimensionality in such data. Magnetic grain size indicators χ_fd%, χARM/SIRM and χ_fd/χARM offer further discrimination of source material but cannot be included in numerical unmixing models owing to non‐linear additivity. This leads to complications in quantitatively ascribing downstream sediment to source areas of contrasting burn severity since sources represent numerical multiples of each other, indicating the need to involve additional indicators, such as geochemical evidence, to allow a more robust discrimination. Copyright © 2005 John Wiley & Sons, Ltd.     

Soil formation and weathering in a permafrost environment of the Swiss Alps: a multi‐parameter and non‐steady‐state approach

Spatially discontinuous permafrost conditions frequently occur in the European Alps. How soils under such conditions have evolved and how they may react to climate warming is largely unknown. This study focuses on the comparison of nearby soils that are characterised by the presence or absence of permafrost (active‐layer thickness: 2–3 m) in the alpine (tundra) and subalpine (forest) range of the Eastern Swiss Alps using a multi‐method (geochemical and mineralogical) approach. Moreover, a new non‐steady‐state concept was applied to determine rates of chemical weathering, soil erosion, soil formation, soil denudation, and soil production. Long‐term chemical weathering rates, soil formation and erosion rates were assessed by using immobile elements, fine‐earth stocks and meteoric ¹⁰Be. In addition, the weathering index (K + Ca)/Ti, the amount of Fe‐ and Al‐oxyhydroxides and clay minerals characteristics were considered. All methods indicated that the differences between permafrost‐affected and non‐permafrost‐affected soils were small. Furthermore, the soils did not uniformly differ in their weathering behaviour. A tendency towards less intense weathering in soils that were affected by permafrost was noted: at most sites, weathering rates, the proportion of oxyhydroxides and the weathering stage of clay minerals were lower in permafrost soils. In part, erosion rates were higher at the permafrost sites and accounted for 79–97% of the denudation rates. In general, soil formation rates (8.8–86.7 t/km²/yr) were in the expected range for Alpine soils. Independent of permafrost conditions, it seems that the local microenvironment (particularly vegetation and subsequently soil organic matter) has strongly influenced denudation rates. As the climate has varied since the beginning of soil evolution, the conditions for soil formation and weathering were not stable over time. Soil evolution in high Alpine settings is complex owing to, among others, spatio‐temporal variations of permafrost conditions and thus climate. This makes predictions of future behaviour very difficult. Copyright © 2016 John Wiley & Sons, Ltd.     

The increase of rainfall erosivity and initial soil erosion processes due to rainfall acidification

The drastic growth of population in highly industrialized urban areas, as well as fossil fuel use, is increasing levels of airborne pollutants and enhancing acid rain. In rapidly developing countries such as Iran, the occurrence of acid rain has also increased. Acid rain is a driving factor of erosion due to the destructive effects on biota and aggregate stability; however, little is known about its impact on specific rates of erosion at the pedon scale. Thus, the present study aimed to investigate the effect of acid rain at pH levels of 5.25, 4.25, and 3.75 for rainfall intensities of 40, 60, and 80 mm h^?1 on initial soil erosion processes under dry and saturated soil conditions using rainfall simulations. The results were compared using a two‐way ANOVA and Duncan tests and showed that initial soil erosion rates with acidic rain and non‐acidic rain under dry soil conditions were significantly different. The highest levels of soil particle loss due to splash effects in all rainfall intensities were observed with the most acidic rain (pH = 3.75), reaching maximum values of 16 g m^?2 min^?1. The lowest levels of particle losses were observed in the control plot where non‐acidic rain was used, with values ranging from 3.8 to 8.1 g m^?2 min^?1. Similarly, under saturated soil conditions, the lowest level of soil particle loss was observed in the control plot, and the highest peaks of soil loss were observed for the most acidic rains (pH = 3.75 and pH = 4.25), reaching maximum average values of 40 g m^?2 min^?1. However, for saturated soils with acidic water but with non‐acidic rain, the highest soil particle loss was observed for the control plot for all the rainfall intensities. In conclusion, acidic rain has a negative impact on soils, which can be more intense with a concomitant increase in rainfall intensity. Rapid solutions, therefore, need to be found to reduce the emission of pollutants into the air, otherwise, rainfall erosivity may drastically increase.     

Diagnosing snow accumulation errors in a rain‐snow transitional environment with snow board observations

Diagnosing the source of errors in snow models requires intensive observations, a flexible model framework to test competing hypotheses, and a methodology to systematically test the dominant snow processes. We present a novel process‐based approach to diagnose model errors through an example that focuses on snow accumulation processes (precipitation partitioning, new snow density, and snow compaction). Twelve years of meteorological and snow board measurements were used to identify the main source of model error on each snow accumulation day. Results show that modeled values of new snow density were outside observational uncertainties in 52% of days available for evaluation, while precipitation partitioning and compaction were in error 45% and 16% of the time, respectively. Precipitation partitioning errors mattered more for total winter accumulation during the anomalously warm winter of 2014–2015, when a higher fraction of precipitation fell within the temperature range where partition methods had the largest error. These results demonstrate how isolating individual model processes can identify the primary source(s) of model error, which helps prioritize future research.