Impacts of Material Engineering Properties on Slope Wash and Stability in Fine-Grained Bedrock Slopes at Fossil-Bearing Sites , Badlands National Park , South Dakota , USA

Engineering properties of bedrock materials at Badlands National Park were used to develop models for Park managers to assess slope erosion and stability for fossil resource protection. Six fully instrumented sites were used to document slope conditions. Bedrock consisted of Oligocene White River Group rocks. Bulk erosion rates correlated to grain size with silty-sandy materials producing higher mass erosion rates as a function of the silt-to-clay ratio and plastic index. Data indicated that as grain size decreased, plastic index increased leading to a decrease in erodibility. These parameters were used to construct a grain-size proxy, ψ, that was substituted for grain size, D, in Bagnold’s entrainment equation and provided significant improvement in calculation of critical entrainment velocities for fine-grained materials. Hydraulic analyses of slope and pediment surface processes indicated surface roughness was a controlling factor and materials washed from rough steep slopes were effectively transported across smooth low-angle pediments with slope-to-pediment angle ratios of nearly 6:1. Slope stability modeling of ten slopes produced high factors of safety for all slopes, even under saturated conditions and was attributable to clay cohesion. All results were used to construct models that predicted years until net slope erosion equaled 2.5 cm (1 inch). Using these results, Park managers were advised to visit erosion-prone sites on a 1to 6-year schedule, based on site geology and slope aspect, to adequately protect critical fossil resources from destruction.


Introduction
Badlands formations are a worldwide phenomenon viewed as a group of systems with location-specific climatic and lithologic drivers [1].Studies of badlands erosional processes in Europe since ~2000 have utilized standard field methods, erosion pins and climatic data, to quantify climatic-landform interactions, that is, linking precipitation intensity, duration and their variations to erosion rates [2][3][4][5][6][7].Resulting erosion rates were determined to primarily be climatically driven and were based on trends, or changing trends, in precipitation or fluvial activity.These and additional studies [8][9][10] have indicated varying erosion rates ranging from negligible to ~0.5 cm/year with higher rates up 1.0 cm/year [11] and >3.0 cm/year in specific areas [12].
In the humid eastern U.S., that is, Perth Amboy, NJ, small badlands-like morphologies resulting from erosion of clay-sand fill created badland forms that compared favorably to morphologies at Badlands National Park (BNP) in South Dakota [13].In South Dakota, the Tertiary White River Badlands have played a leading role in the development of vertebrate paleontology in North America, beginning with the discovery and description of a titanothere mandible [14].This and similar work during the next 80 years led to the establishment of Badlands National Monument in 1929 and Geosciences 2018, 8, 267; doi:10.3390/geosciences8070267www.mdpi.com/journal/geosciencessubsequent renaming to Badlands National Park in 1978, with the mission of preserving the scenic and scientific value of a portion of the White Rive Badlands, while remaining open and accessible to the public.Subsequent to Prout's initial description, thousands of fossils have been collected, described and catalogued and have served to define the Tertiary White River geologic deposits, which contains one of the world's most abundant vertebrate fossil records from the late Eocene and early Oligocene [15].Scientific research focused on erosion at BNP has been limited.Schumm [16] studied the retreat rates of miniature pediments over an ~8-year period from the mid 1950s to early 1960s.Slopes were measured by reinforcing rods driven into the ground at the base of a slope and monitoring exposure rates over the duration of the study.He used Manning's equation to determine that flow velocities generated across both slope and pediment surfaces were nearly identical.He concluded that due to the lower angle of the pediment surface and the equal flow velocity, more erosive energy was directed toward the pediment and these surfaces effectively transported all debris away from the slope.Slope retreat was replaced by development of a new pediment surface between 6 cm and 12 cm wide at the base of the slopes.Uniformly applying this rate over the study period, slope retreat was determined to be 0.75 cm/year to 1.5 cm/year.Later, Stoffer [17] estimated an average erosion rate in BNP of ~2.5 cm/year.These minimal efforts expended at BNP led to a Park-wide study to determine erosion rates at fossil-bearing sites.Initial results indicated slope retreat rates, using erosion pins, were ~2 mm during the initial 2.5 months of the project [18].Final project results [19] showed site-wide average erosion rates were 0.41 cm/year and ranged from 0.25 cm to 1.5 cm as a function of geologic formation, material engineering properties, slope aspect and precipitation.

Research Goal
The BNP mission of resource preservation led to establishing the goal of constructing a fossil-site monitoring schedule.The schedule was based on statistical analysis of measured erosion rates as functions of the physical and mechanical properties of the fine-grained bedrock units in the Park.These properties were identified from analysis of collected data from six fully instrumented erosion sites located in the North Unit of BNP.Three sites were located in the east and three sites in the west (Figure 1).Bedrock and storm-based slope wash samples were collected and analyzed to define engineering and sedimentological properties.Resulting slope wash and stability models were used to assess the effects individual properties had on erosion processes at Badlands.This paper will describe field and analytical techniques utilized to determine erosion rates and development of a new grain size proxy term used to calculate fluid entrainment velocities for fine-grained bedrock materials.
America, beginning with the discovery and description of a titanothere mandible [14].This and similar work during the next 80 years led to the establishment of Badlands National Monument in 1929 and subsequent renaming to Badlands National Park in 1978, with the mission of preserving the scenic and scientific value of a portion of the White Rive Badlands, while remaining open and accessible to the public.Subsequent to Prout's initial description, thousands of fossils have been collected, described and catalogued and have served to define the Tertiary White River geologic deposits, which contains one of the world's most abundant vertebrate fossil records from the late Eocene and early Oligocene [15].
Scientific research focused on erosion at BNP has been limited.Schumm [16] studied the retreat rates of miniature pediments over an ~8-year period from the mid 1950s to early 1960s.Slopes were measured by reinforcing rods driven into the ground at the base of a slope and monitoring exposure rates over the duration of the study.He used Manning's equation to determine that flow velocities generated across both slope and pediment surfaces were nearly identical.He concluded that due to the lower angle of the pediment surface and the equal flow velocity, more erosive energy was directed toward the pediment and these surfaces effectively transported all debris away from the slope.Slope retreat was replaced by development of a new pediment surface between 6 cm and 12 cm wide at the base of the slopes.Uniformly applying this rate over the study period, slope retreat was determined to be 0.75 cm/year to 1.5 cm/year.Later, Stoffer [17] estimated an average erosion rate in BNP of ~2.5 cm/year.These minimal efforts expended at BNP led to a Park-wide study to determine erosion rates at fossil-bearing sites.Initial results indicated slope retreat rates, using erosion pins, were ~2 mm during the initial 2.5 months of the project [18].Final project results [19] showed site-wide average erosion rates were 0.41 cm/year and ranged from 0.25 cm to 1.5 cm as a function of geologic formation, material engineering properties, slope aspect and precipitation.

Research Goal
The BNP mission of resource preservation led to establishing the goal of constructing a fossilsite monitoring schedule.The schedule was based on statistical analysis of measured erosion rates as functions of the physical and mechanical properties of the fine-grained bedrock units in the Park.These properties were identified from analysis of collected data from six fully instrumented erosion sites located in the North Unit of BNP.Three sites were located in the east and three sites in the west (Figure 1).Bedrock and storm-based slope wash samples were collected and analyzed to define engineering and sedimentological properties.Resulting slope wash and stability models were used to assess the effects individual properties had on erosion processes at Badlands.This paper will describe field and analytical techniques utilized to determine erosion rates and development of a new grain size proxy term used to calculate fluid entrainment velocities for fine-grained bedrock materials.

Site Geology
Rocks exposed within BNP consisted primarily of the White River Group and surrounding sediments that span from late Cretaceous to mid Tertiary (early Oligocene) time [15].Cretaceous-Paleocene rocks span a timeframe from 74 to 55.5 Ma.Cretaceous rocks have been divided into several facies members [17] of Campanian and Maestricthian age.Sediments consisted of marine facies deposited in the Cretaceous Western Interior Seaway and were ~110 m thick in the Park.These rocks were primarily gray shale with abundant bentonite beds, limestone concretions and chalk and marl beds hosting abundant marine fauna.After the sea drained away, the upper Cretaceous Fox Hills Formation marine rocks were subjected to intensive tropical weathering conditions and pedogenesis.This resulted in development of the Paleocene Yellow Mounds paleosol, a series of goldcolored alternating fine-grained sandstone and shale beds ubiquitous throughout the Park (Figure 2).
Above the marine Cretaceous rocks, White River Group sediments were mainly fluvial in origin and consisted of poorly consolidated fluvial clay and silt deposits that included volcanoclastic and

Site Geology
Rocks exposed within BNP consisted primarily of the White River Group and surrounding sediments that span from late Cretaceous to mid Tertiary (early Oligocene) time [15].Cretaceous-Paleocene rocks span a timeframe from 74 to 55.5 Ma.Cretaceous rocks have been divided into several facies members [17] of Campanian and Maestricthian age.Sediments consisted of marine facies deposited in the Cretaceous Western Interior Seaway and were ~110 m thick in the Park.These rocks were primarily gray shale with abundant bentonite beds, limestone concretions and chalk and marl beds hosting abundant marine fauna.After the sea drained away, the upper Cretaceous Fox Hills Formation marine rocks were subjected to intensive tropical weathering conditions and pedogenesis.This resulted in development of the Paleocene Yellow Mounds paleosol, a series of gold-colored alternating fine-grained sandstone and shale beds ubiquitous throughout the Park (Figure 2).
Above the marine Cretaceous rocks, White River Group sediments were mainly fluvial in origin and consisted of poorly consolidated fluvial clay and silt deposits that included volcanoclastic and bentonite clay layers.Illite and montmorillonite have been identified as the predominant minerals [21] in these fine-grained deposits.White River Group rocks consisted of the Chamberlain Pass, Chadron and Brule Formations [15].The caprock for the entire sequence consists of the overlying late Oligocene Sharps Formation.
Eocene rocks represent deposition from 55.5 to 36 Ma The lowermost Chamberlain Pass Formation consisted of overbank sediments, conglomerates and channel sands and was ~4 m thick.Overbank fines have been weathered into the bright red Interior paleosol immediately overlying and providing color contrast to the Yellow Mounds paleosol.Above the paleosols, Chadron Formation rocks ranged up to 55 m in thickness and were deposited within the pre-existing valley and hill topography.These rocks were mainly gray fine-grained claystone beds containing abundant bentonite units that swell when wet and becomes very slippery.When dry, it forms a hard surface marked by desiccation features resembling small aggregated popcorn-shaped mounds (Figure 3).Weathering of Chadron rocks produced generally low relief convexo-concave hills and were quite distinct in occurrence and ubiquitous across the Park.Overlying Oligocene rocks were deposited from 36 to 24 Ma and consisted of the Brule Formation, a mainly fluvial deposit of clay, silt and sand up to 140 m thick.The formation has been divided into the Scenic (Orellan) and Poleslide (Whitneyan) members.The lower Scenic consisted of yellowish gray to buff arenaceous claystones, paleosols and limestone beds near the top.The upper Poleslide consisted of gray to creamy-colored clay, dominated by light gray massive cliff forming sandstone that contained carbonate nodules near the top [17,22].
Geosciences 2018, 8, x FOR PEER REVIEW 4 of 17 bentonite clay layers.Illite and montmorillonite have been identified as the predominant minerals [21] in these fine-grained deposits.White River Group rocks consisted of the Chamberlain Pass, Chadron and Brule Formations [15].The caprock for the entire sequence consists of the overlying late Oligocene Sharps Formation.Eocene rocks represent deposition from 55.5 to 36 Ma The lowermost Chamberlain Pass Formation consisted of overbank sediments, conglomerates and channel sands and was ~4 m thick.Overbank fines have been weathered into the bright red Interior paleosol immediately overlying and providing color contrast to the Yellow Mounds paleosol.Above the paleosols, Chadron Formation rocks ranged up to 55 m in thickness and were deposited within the pre-existing valley and hill topography.These rocks were mainly gray fine-grained claystone beds containing abundant bentonite units that swell when wet and becomes very slippery.When dry, it forms a hard surface marked by desiccation features resembling small aggregated popcorn-shaped mounds (Figure 3).Weathering of Chadron rocks produced generally low relief convexo-concave hills and were quite distinct in occurrence and ubiquitous across the Park.Overlying Oligocene rocks were deposited from 36 to 24 Ma and consisted of the Brule Formation, a mainly fluvial deposit of clay, silt and sand up to 140 m thick.The formation has been divided into the Scenic (Orellan) and Poleslide (Whitneyan) members.The lower Scenic consisted of yellowish gray to buff arenaceous claystones, paleosols and limestone beds near the top.The upper Poleslide consisted of gray to creamy-colored clay, dominated by light gray massive cliff forming sandstone that contained carbonate nodules near the top [17,22].Where erosion has removed the softer overlying clay, sandstone channels have weathered out producing distinctive columns of underlying clay capped by the protective sands, characterized by a step-and-riser topography with steep sided cliffs and sharp pinnacles (Figure 2).Texturally, the Brule rocks ranged from lower sections that were ultra fine-grained, forming low-angle slopes [22,23], to an upper section that was more silty, thus able to hold near vertical slopes that approach heights of 100 m.It is this textural variation that leads to the formation of the many spires and pinnacles that Badlands National Park is known for.
Miocene Arikaree Group rock unconformably overlies this entire sequence and forms the cap rock.Locally, it consisted of the Rockyford Ash, derived from eolian ash deposits that originated from Great Basin volcanism and the Sharps Formation, consisting mainly of light gray tuffaceous Where erosion has removed the softer overlying clay, sandstone channels have weathered out producing distinctive columns of underlying clay capped by the protective sands, characterized by a step-and-riser topography with steep sided cliffs and sharp pinnacles (Figure 2).Texturally, the Brule rocks ranged from lower sections that were ultra fine-grained, forming low-angle slopes [22,23], to an upper section that was more silty, thus able to hold near vertical slopes that approach heights of 100 m.It is this textural variation that leads to the formation of the many spires and pinnacles that Badlands National Park is known for.
Miocene Arikaree Group rock unconformably overlies this entire sequence and forms the cap rock.Locally, it consisted of the Rockyford Ash, derived from eolian ash deposits that originated from Great Basin volcanism and the Sharps Formation, consisting mainly of light gray tuffaceous sandstone and mudrock deposited in a relict topography of channels and flats.
Regolith was usually very thin to absent throughout the Park and ranged from unweathered bedrock to quasi developed soils ~30 cm deep [24].Generally, Brule outcrops had thin and compact 1.0 cm to 5.0 cm of regolith and Chadron outcrops had a loose but thicker regolith that could be >10 cm.Regolith accumulation was related to slope angle where steep slope angles hold less weathered material.

Climate Setting
Climate data were available for the Park from 1998 to present utilizing the Remote Automated Weather Stations (RAWS).Badlands climate is typical intercontinental characterized by hot dry summers and typically cool wetter winters.Precipitation averaged 430 mm/year ranging from a low of 249 mm/year to a high of 712 mm/year.Yearly distribution typically contained late-spring, mid-summer precipitation peaks where single events greater than 20 mm randomly occurred.Predominate spring and summer precipitation patterns originated from the north to northwest.The average annual temperature was ~9.5 • C with an average annual high temperature of 17.4 • C and average annual low of 2.8 • C. Annual average wind speed was ~18 kph with maximum average gusts of 48 kph.Typical spring and summer precipitation is derived from predominantly W-NW directions.
Both spring and fall seasons contained multiple weeks of diurnal freeze-thaw of the surface materials.The effects were most pronounced in fine-grained clay rocks such as the Chadron Formation.Ice crystal formation effectively splits the top few centimeters of depth into small spear-like slivers, disintegrating the rock into distinct small pieces that become subjective to downslope creep processes.These splits also allow rapid disaggregation to occur during subsequent precipitation wetting and slope wash processes.The typical popcorn texture (Figure 3) of the Chadron is in part due to frost action.

Materials and Methods
Several different measurement methods were used and included stainless steel erosion pins measured using digital calipers, stainless steel rules, aluminum engineering scales, or some combination thereof, that were set at multiple locations along a slope.Full descriptions of site-layouts have been provided elsewhere [18].Both rules and scales were imaged with a 14 MP Sony digital camera (Tokyo, Japan).One-meter-long sediment trays placed flush with the toe of monitored slopes were used to collect event-based sediment runoff from a known upslope contributing area.Precipitation data were collected at each site using a digital tipping-bucket rain gauge.Initial slope geometry was measured using a Brunton compass (Riverton, WY, USA) and tape and subsequently using 3D photogrammetry.Slope data included documented slope profile and surface changes and surface changes around monitored fossils.Twenty-four bedrock samples were collected and analyzed for water content, grain size distribution (wet-sieve and laser particle sizing), Atterberg limits and undrained cohesion.Slope modeling at each site, using Slide 6.0, a 2D limit equilibrium slope stability model (rocscience, Toronto, ON, Canada, www.rocscience.com/products/8/Slide),resulted in factor-of-safety (FS) determination for natural and saturated conditions.

Materials Testing Procedures
Slope wash sediment collected from trays at slope toes were weighed, oven dried at 105 • C for 24 h and re-weighed to obtain water content and eroded mass.All bedrock samples were prepared by dispersing the material using sodium hexametaphosphate [25] and wet-sieving.The silt and finer fraction was determined using the Modified Wentworth system screen sizes 230 (63 µm) to 850 (10 µm) under a low constant water flow for an eight-minute period.All material exiting the 10 µm screen were subsequently evaluated using a laser-based liquid particle size analyzer.Additional laboratory testing of whole bedrock included determining specific gravity, Atterberg limits and unconfined compressive strength.Atterberg limits testing resulted in determining liquid and plastic limits (LL and PL) and the plasticity index (PI = LL − PL).

Slope Stability Modeling Procedures
Physical slope dimensions were measured directly and by using photogrammetry and were used to construct model profiles for stability analysis.Each slope consisted of a single bedrock layer that was weathered on the slope surface into a varying thickness mantle of loose, disintegrated material that appeared much like a soil regolith.All slopes measured consisted of this weathered surface (thickness dependent upon slope angle) and was considered less cohesive than solid bedrock.Thus, analyses were performed using properties of normal regolith instead of bedrock and were assumed to simulate the worst weathered case where bedrock had reduced to a loose soil-like material on the surface, which was visual confirmed by field observation and testing.
Three limit equilibrium methods were utilized to solve for the FS, Bishop's simplified method, Janbu's simplified method and GLE/Morgenstern-Price method [26,27].A sensitivity analysis was performed by varying the ranges of cohesion and friction angle values used in the models to determine the influence these parameters had on slope stability.
Rounded convexo-concave slopes in clay-rich badlands formations (primarily Chadron rocks) have been attributed to be primarily a function of creep processes [13,28,29].Formations having clays subjective to shrink-swell experience large volume changes with reduction/addition of water which tends to expand the surface as weathering progresses.The resulting top layer of loose, soil-like material has a lower bulk strength than the underlying bedrock and moves downslope by sag and creep [24].As a result, slopes are at, or close to, their maximum angle of stability.Such conditions can be modeled by the infinite slope model for saturated regolith [30][31][32].The general infinite slope model [24] was used to determine the FS for homogeneous soil acting under undrained shear strength for the two slope cases specified below.
Case 1: when where d w = depth of water above the bedrock surface (depends on soil permeability), d r = depth of regolith (weathered) layer, θ = slope angle determined for each slope profile, γ sat = saturated unit weight, γ = total unit weight, UCS = unconfined compressive strength, all in consistent units.The depth of water above the bedrock surface typically changes with precipitation and soil permeability which varied from 10 −3 to 10 −7 cm s −1 for fine sand, silt and mixtures of sand, silt and clay to <10 −7 cm s −1 for homogeneous clay [33].The depth of regolith layer was estimated based on Howard [24] and visual estimations from the field that were consistent with these values.The slope angle was a measured value and remained invariant for each individual slope.γ sat and γ were obtained from laboratory work.Using Equations ( 1) and ( 2), FS values were calculated for slopes at Badlands.

Results
Table 1 contains the results of the materials testing for Badlands formations.These data were utilized to evaluate slope erosivity and stability, the latter by constructing stability models for assessing the effects that water content, cohesion and internal friction angle had on FS analysis.LL = liquid limit; PL = plastic limit; PI = plasticity index; SG = specific gravity.

Physical and Engineering Properties
The finest-grained material identified was the clay-rich gray-green mudstone from the Peanut Peak member of the Chadron Formation (Site BL-05) with 84% passing the 850 screen (10 µm).It had the smallest averaged mean diameter (8.09 µm), the highest Atterberg limits, highest initial water content and smallest unsaturated and saturated unit weight.Conversely, the coarsest material was the Scenic member of the Brule Formation at Site BL-02 which consisted of gray silty sandstone with 43% of the material passing the 850 screen.It had the largest averaged mean diameter (33 µm), the lowest Atterberg limits and largest unsaturated and saturated unit weight.Across the North Unit, the 3 sites in the east (BL-01, BL-02, BL-03) were coarser-grained (average 23.6 µm), had the lowest Atterberg limits (average PI = 27.9),largest unit weights and the smallest undrained cohesion.The 3 sites in the west (BL-04, BL-05, BL-06) were finer-grained (average 11.1 µm), had the greatest Atterberg limits (average PI = 69.3),smallest unit weights and the greatest undrained cohesion.One likely reason for the distinctions in parameters, even in the same member unit, was that the regional dip of the formations (~2-3 • E) exposed different portions of the formations from west to east.These data trends were consistent with the collected erosion data.
Specific gravity (SG) values were determined using the eroded materials collected from the sediment trays following ASTM methods [34].Scenic and Poleslide members were ~2.57g cm −3 and the Chadron was 4.7% higher at 2.68 g cm −3 .Specific gravity values had a standard deviation of 0.049.

Slope Erosivity and Stability
Slope profiles represent physical process responses between the forces generated on the slope driving downslope motion of bedrock and loose debris and the engineering properties of the slope materials resisting motion.Forces of erosion may be generated by either fluid or gravity.Howard [24] described erosional processes acting on Badlands slopes and divided them between wash (fluid) and mass wasting (gravity).Wash processes primarily occurred from rainstorms including runoff and splash erosion.Mass wasting resulted mainly from unstable regolith or bedrock moving downslope under gravity including creep, sliding and flows [17].Both situations were modeled during this study using engineering parameters for exposed bedrock at erosion sites.

Erosivity from Wash
Slope wash produces eroded sediment as a function of mass and depth of flowing water, slope angle and length and the properties (size, shape, mechanical attributes) of particles resting on the slope.
For non-cohesive particles, downslope transport occurs when the critical shear stress, τ cr , exceeds material strength and is dependent upon grain size, D, that is, τ cr α D [35].τ cr is a difficult parameter to quantify for cohesive slit and clay (<16 µm diameter) [36], in that erosion resistance is increased above that for granular material (0.1-10 mm) and depends on cohesive attraction (valency charge forces) generated between adjacent particles rather than on the physical parameters of the particles themselves.Thus, a critical shear stress does not exist in the same meaning as for non-cohesive materials.
At Badlands, slopes <35-40 • consisted uniformly of weathered bedrock materials (regolith) that mantled slopes to depths of up to ~30 cm but were typically thinner.The material appeared loose, fluffy and airy with significant cohesion, particularly in finer-grained units.When wetted, slope wash physically removed and transported some of this material toward the slope toe depositing it as a type of fan, or pediment (also noted by others) [13,16,22,24] that was typically smooth and compacted.The same runoff was likely to carry some portion of the eroded material away from the slope resulting in a net lateral slope retreat over time.
Physically, the erosional mechanism for regolith proceeds as a grain-fluid interaction.Wetting of the regolith effectively disaggregated the non-lithified fine-grained bedrock particles and slope wash consisted of individual particles, not aggregated clumps.This process produced the smooth and compacted pediment surfaces noted by Schumm, Smith and Howard [16,22,24] and were observed during the present study.Additionally, eroded materials collected in sampling trays consisted of disaggregated sediment with occasional aggregated clumps (Figure 3).These clumps were always attributable to non-precipitation event gravity falls and were easily distinguished from slope wash as they were always resting on top of the sediment in the tray and not incorporated into it.Immediately after precipitation events and until slopes desiccated, no collected eroded materials were aggregated.Thus, for these bedrock materials, slope wash occurred mainly by entrainment and transport of individual silt and clay particles.
Utley and Wynn [37] evaluated wash processes for cohesive materials using a method that related critical shear stress, τ cr (Pa), to several soil properties including PI. Critical shear stress was then used to determine a soil erodibility coefficient, k d (cm 3 /N-s): τ cr = 0.16(PI) 0.84  (3) where k d is an indicator of erosion resistance, where the smaller the value, the greater the erosional resistance.Mechanically, as PI increased, cohesive forces between adjacent particles increased (Figure 4A) requiring higher shear stresses to mobilize the grain which resulted in an increased resistance to erosion.The plasticity index is also reflective of the silt to clay ratio, S:M, based on the sieve analysis and defined as the ratio between the percentage of grains at 60% and 10% finer, in that smaller ratios indicate greater clay content and an increased PI value.On Figure 4A, the coarse-grained material (closed symbols) had a S:M of 0.19 and the fine-grained materials (open symbols) had a S:M of 0.05.These data indicate that higher clay content resulted in a lower S:M ratio, an increased PI and a decrease in k d , resulting in increased erosion resistance.
Resistance to erosion was determined by directly comparing calculated k d values to the mass flux, ṁ (g m −2 ), of eroded sediment from monitored slopes (Figure 4B).Mass flux was determined by dividing the total mass (g) of materials collected in a sediment tray by the contributing slope area (m 2 ).These results suggest that erosivity was a function of grain size, that is, finer-grained materials were less erosive than coarser-grained materials and these relations represent engineering properties of the materials.
As PI decreased, k d increased logarithmically in response to a decrease in erosion resistance.Increasing S:M ratios resulted in greater amounts of slope wash, however, mass flux decreased logarithmically through the process.This suggests the existence of an upper limit to PI effects for erosion resistance and lies most likely in the LL, itself a function of the S:M ratio.These results suggested that erodibility assessment for cohesive materials are possible if engineering properties of the material are known, particularly PI.In most cases, attaining engineering properties for representative materials is a rapid process and has been used effective for mitigation efforts rather than determining site-specific erodibility conditions.At the least, stability results based on materials engineering properties would enhance site selection where long-term monitoring may be required.Although these results appear to be robust for calculating erosivity potential, they do not provide an actual value for τ cr .Therefore, a method for determining τ cr for cohesive materials is needed.
Hjulström [38] and Shields [39] demonstrated that erosion (entrainment) resistance increased for cohesive materials above those of granular non-cohesive materials.Shields used dimensionless parameters for grain size and shear stress and determined regions of motion and no motion.Hjulström's diagram (Figure 5) segregated process-based regions from erosional to depositional conditions for grains as small as 0.001 mm (1 µm).Velocities required to entrain grains ranged widely between unconsolidated and consolidated materials and was defined as the region between the dashed lines.For example, an unconsolidated particle of size 0.01 mm would move with a velocity of ~30 cm s −1 rising to >100 cm s −1 if it were consolidated.Although equations exist to determine entrainment velocities for granular materials >0.1 mm, few equations exist for clay particles.
Using Figure 5, average-sized particles for Badlands formations require fluid entrainment velocities between ~70-300 cm s −1 .The standard equation used to calculate critical entrainment velocity, u * c , for granular materials was developed by Bagnold [40] and has been widely used: where k is the fluid coefficient (~0.2 for water), σ is the particle density, ρ is the fluid density, g is the gravitational constant and D is the mean particle diameter.This equation is valid for particles >0.1 mm (100 µm) and, thus, cannot be used directly for the fine-grained materials at BNP.For smaller particles, entrainment velocities increase as a function of decreasing grain size attributable to chemical force attraction between particles, or cohesion and fundamental particle shape changes.The complexity of the particle force attractions versus motive forces makes this problem non-trivial.That k d has been shown to correlate to material engineering properties suggests these same properties may be useful to construct a grain size proxy for use in Equation (5).For cohesive soils, PI has been shown to be a function of grain size and to control erosion resistance through k d (Figure 4).Thus, a proxy grain size parameter, ψ, was developed based on the engineering parameters of cohesive materials.The relative strength of these parameters relies on the physical distribution of particles represented by the S:M ratio.The grain size proxy was determined as where C is the cohesion.The plasticity index and cohesion hold the material together by generating bulk strength properties and are primarily related to the small particle sizes for clay-rich material.Materials defined by small S:M ratios provide conditions where grain platelets arrange in a more stacked configuration with decreasing size where increased molecular force interactions on grain edges and surfaces dominate.Wetting and hydrating the particles will further increase attractive forces [41].The combined effect is to generate material strength that is reflected by increases in PI and cohesion.Changes in the grain-size distribution alters this function.Materials having a larger S:M ratio leads to decreased grain-to-grain contact and reduced molecular charge attraction.The more granular the mixture becomes, the less bulk strength the material has and is reflected by lower PI and cohesion.The effect of the S:M ratio provides the basic physical mechanism that defines erodibility of fine-grained materials.
Substituting ψ as a proxy for the mean grain diameter, D, in ( 5) and plotted on the Hjulström diagram (Figure 5) indicated that Badlands materials behave primarily as consolidated, that is, cohesive, material.Decreases in the S:M ratio increased bulk strength increasing erosive velocity in fine-grained clay-dominant materials.Materials that contained larger S:M ratios (higher silt content) behaved oppositely.The range of critical velocities required to entrain these materials were calculated to be between 75 and 300 cm s −1 .Not only does u * c plot correctly in terms of material type but is it significantly correlated to k d , the resistance parameter (Figure 6).
For the range of particle-size distributions for Badlands slopes, entrainment velocities were calculated using: where k d is based on PI.This equation remains valid over the range of PI values for the materials that were tested.Schumm [16] used Manning's equation for open channel flow to evaluate velocity on both slopes and pediments at BNP.In his analysis, the flow depth was assumed to be the same on slopes and pediments thereby eliminating it as a parameter in the equation.The result was a relative velocity measure, V r : V r = 1.49n S 0.5 (8) where n is the roughness coefficient and S is the slope inclination.Roughness values were selected from Chow [42] and were 0.06 and 0.02 for the slope and pediment, respectively (slopes were rougher).
Geosciences 2018, 8, x FOR PEER REVIEW 10 of 17 where n is the roughness coefficient and S is the slope inclination.Roughness values were selected from Chow [42] and were 0.06 and 0.02 for the slope and pediment, respectively (slopes were rougher).
(A) (B) All measurements were averaged and Vr was 23 and 24 for slopes and pediments, respectively.That Vr was the same on both surfaces was a result of decreasing roughness on the pediment being compensated by a corresponding decrease in the slope angle away from the slope.All measurements were averaged and V r was 23 and 24 for slopes and pediments, respectively.That V r was the same on both surfaces was a result of decreasing roughness on the pediment being compensated by a corresponding decrease in the slope angle away from the slope.Hjulström's diagram relates processes to velocities.Particles <0.1 mm require higher entrainment velocities as functions of grain size, shape and particle attractive forces.Badlands materials are plotted that were determined using the grain-size proxy.Open symbols had a smaller S:M ratio and were more clay-rich than closed symbol materials that had a larger S:M ratio.
Using measurements of slopes and pediments obtained from 3D photogrammetry and the same roughness values used by Schumm [16], Vr comparisons were calculated.Average relative velocity values were the same on slopes and pediment surfaces and were 26.Although the angle ratio between slope and pediment was 5.8, the decrease in roughness on the pediment produced a consistent velocity and energy condition to effectively transport particles across the lower angle surfaces.Sediment delivered to the toe of the slope, was effectively transported by efficiency of flow energy across the pediment and in time resulted in a general lateral slope retreat.This process was facilitated by the disaggregation of the bedrock when wetted and Atterberg limits were above the LL and particle strength was reduced.Sediment transport processes were activated during runoff events and the more precipitation received on a slope, the higher the mass flux of sediment (Figure 7).These data show net slope retreat for 2 adjacent slopes at site BL-01 as functions of precipitation.Slope retreat values were measured as mm of erosion (surface lowering) between successive measurement periods.Slope retreat was accelerated during summer 2011 during a relatively large precipitation event.Slope movement was mostly subdued until the next precipitation event in spring 2012.However, retreat was less in response to a lower intensity event.These data also show effects of slope aspect, or direction of the slope to the direction of the precipitation.Slope 05R was south-facing and was shielded from the primary intensity of the precipitation from the northwest.The north-facing slope, 06R, received the full impact of the precipitation and produced greater sediment mass flux and larger net slope retreat.These slopes were separated by a ~2 m wide pediment and were in the same units of the same formation.These processes acting through time are responsible for the generation of the low-angle smooth pediment surfaces leading away from slopes that are ubiquitous features across the Badlands.Hjulström's diagram relates processes to velocities.Particles <0.1 mm require higher entrainment velocities as functions of grain size, shape and particle attractive forces.Badlands materials are plotted that were determined using the grain-size proxy.Open symbols had a smaller S:M ratio and were more clay-rich than closed symbol materials that had a larger S:M ratio.
Using measurements of slopes and pediments obtained from 3D photogrammetry and the same roughness values used by Schumm [16], V r comparisons were calculated.Average relative velocity values were the same on slopes and pediment surfaces and were 26.Although the angle ratio between slope and pediment was 5.8, the decrease in roughness on the pediment produced a consistent velocity and energy condition to effectively transport particles across the lower angle surfaces.Sediment delivered to the toe of the slope, was effectively transported by efficiency of flow energy across the pediment and in time resulted in a general lateral slope retreat.This process was facilitated by the disaggregation of the bedrock when wetted and Atterberg limits were above the LL and particle strength was reduced.Sediment transport processes were activated during runoff events and the more precipitation received on a slope, the higher the mass flux of sediment (Figure 7).These data show net slope retreat for 2 adjacent slopes at site BL-01 as functions of precipitation.Slope retreat values were measured as mm of erosion (surface lowering) between successive measurement periods.Slope retreat was accelerated during summer 2011 during a relatively large precipitation event.Slope movement was mostly subdued until the next precipitation event in spring 2012.However, retreat was less in response to a lower intensity event.These data also show effects of slope aspect, or direction of the slope to the direction of the precipitation.Slope 05R was south-facing and was shielded from the primary intensity of the precipitation from the northwest.The north-facing slope, 06R, received the full impact of the precipitation and produced greater sediment mass flux and larger net slope retreat.These slopes were separated by a ~2 m wide pediment and were in the same units of the same formation.These processes acting through time are responsible for the generation of the low-angle smooth pediment surfaces leading away from slopes that are ubiquitous features across the Badlands.

Slope Stability Modeling
Landslides represent normal geomorphic processes indicating instability in slope profiles or force balance.Acting through time, these processes produced equilibrium slope profiles as a function of slope material, slope height and climate intensity [22].At Badlands, landslides have been mapped in the Cedar Pass, Sage Creek and Norbeck Pass areas [23] that are contained in the red circles on Figure 1.Howard [20] showed that mass wasting rates rapidly increased when slope gradients approached the limiting slope angle of ~34-36°.Slope angles greater than the limiting value typically occur in Brule formation rocks which forms the Badlands Wall (Figure 2).Brule rocks contain sufficient silt and fine sand that with the clay needed to provide cohesive strength has the inherent ability to weather into the characteristic tall and steep sloped forms [22].Expectedly, the greatest number of landslides in the Park have been identified in the Brule.Instability leading to failure has occurred due to tension failure from both desiccation and/or undercutting.

Slope Stability Modeling
Landslides represent normal geomorphic processes indicating instability in slope profiles or force balance.Acting through time, these processes produced equilibrium slope profiles as a function of slope material, slope height and climate intensity [22].At Badlands, landslides have been mapped in the Cedar Pass, Sage Creek and Norbeck Pass areas [23] that are contained in the red circles on Figure 1.Howard [20] showed that mass wasting rates rapidly increased when slope gradients approached the limiting slope angle of ~34-36°.Slope angles greater than the limiting value typically occur in Brule formation rocks which forms the Badlands Wall (Figure 2).Brule rocks contain sufficient silt and fine sand that with the clay needed to provide cohesive strength has the inherent ability to weather into the characteristic tall and steep sloped forms [22].Expectedly, the greatest number of landslides in the Park have been identified in the Brule.Instability leading to failure has occurred due to tension failure from both desiccation and/or undercutting.

Slope Stability Modeling
Landslides represent normal geomorphic processes indicating instability in slope profiles or force balance.Acting through time, these processes produced equilibrium slope profiles as a function of slope material, slope height and climate intensity [22].At Badlands, landslides have been mapped in the Cedar Pass, Sage Creek and Norbeck Pass areas [23] that are contained in the red circles on Figure 1.Howard [20] showed that mass wasting rates rapidly increased when slope gradients approached the limiting slope angle of ~34-36 • .Slope angles greater than the limiting value typically occur in Brule formation rocks which forms the Badlands Wall (Figure 2).Brule rocks contain sufficient silt and fine sand that with the clay needed to provide cohesive strength has the inherent ability to weather into the characteristic tall and steep sloped forms [22].Expectedly, the greatest number of landslides in the Park have been identified in the Brule.Instability leading to failure has occurred due to tension failure from both desiccation and/or undercutting.
Smith [22] also suggested creep-like movement predominated on lower gradient slopes, such as in Chadron formation slopes that typically have maximum angles ~26 • (Figure 2).Grain size distributions indicated Brule rocks contained 50% of the grains >25 µm diameter and were the coarsest units identified and also contained the tallest and steepest (~52 • ) slopes that were measured during the study.Chadron rocks were the finest-grained units identified with 50% of the grains >4.6 µm and ~80% of the grains <10 µm.These rocks contained the greatest cohesion but due to a greater clay content had the least steep slopes and weathered into characteristic low-profile convexo-concave mounds.
Slope stability modeling (Equations ( 1) and ( 2)) requires a known material shear strength.Typically, saturated clays have no internal friction and in undrained conditions, undrained cohesion, c u , can be considered equal to shear strength.Thus, the strength test results provided c u and unconfined compressive strength, q u , was determined using: Unconfined compressive strength was subsequently used in the FS equations for slope stability assessment.Stability analysis were conducted using three limit equilibrium methods (Section 2.2).The physical setting was based on the measured slope profile.Model parameters were determined either through direct material testing or parameter modeling using test results.All FS results from all three methods were similar indicating no significant differences between methods.Variance between method results were also provided from the modeling and were zero for coarse-grained sites (BL-01 through BL-03) and ranged from zero to 9 for the fine-grained sites (BL-04 through BL-06).FS and variance were unusually high for slopes BL-06 since these were low angle slopes (<10 • ).Slopes BL-04 through BL-06 had higher FS values for 2 reasons: (1) grain-size distributions on these slopes were significantly skewed toward finer grains (50% < 10 µm) and (2) increased cohesion (high PI values) of the materials on the slope enhanced erosion resistance.Lowest FS was obtained in all cases under a surface water table and saturated conditions.Even under these saturation condition, no slopes contained a FS less than threshold conditions.Slopes closest to failure conditions were BL-01 and was located in an area of ubiquitous failures in Brule rocks, mainly from a combination of saturation and tension failure.Tallest modeled slopes were ~5 m high and slope failures began to appear in slope with heights above ~20 m.
Sensitivity analysis for the 10 modeled slopes indicated that slope stability was more sensitive to cohesion than angle of internal friction.Varying cohesion between 0 and 50% of the true value resulted in FS variations from 11% to 31%.Conversely, variations in internal angle of friction between 0 and 50% of the true value resulted in FS variations from 0 to 1.42%.
Material properties were mainly derived from the fine grain size, particularly the S:M ratio.Increasing the clay content increased cohesion thereby limiting variations in FS values.These FS variations were more pronounced for the coarse-grained materials at sites BL-01 through BL-03 where FS changes averaged 22%.For these coarse-grain dominate slopes, variations in the clay content changed cohesion enough to result in a high variability in FS.For fine-grained slopes (BL-04 through BL-06), variations in clay content had a more subdued effect on cohesion owing to the overall finer textured materials and the result was a subdued variation in FS values.These sites had an average FS change of 14%, a net decrease of 64% in average FS variation between the coarse-and fine-grained sites.

Site Monitoring Schedule
Results of net slope retreat from ten monitored slopes as functions of slope angle, aspect and precipitation were modeled to produce a first-cut suggestion for a site monitoring schedule.Multiple regression models were used to determine the most sensitive modeling variables that controlled erosion and could be used to predict future site monitoring.Slope aspect was the most critical factor for determining erosion potential and was followed by formation, or bedrock type (Table 2).An erosion ratio was calculated between the most erosive condition and all other slopes.For example, for slope aspect, the north-facing slopes had the highest erosion rate and south-facing slopes had the least.The erosion ratio for north-facing slopes was 1.0 and for south-facing slopes was 0.2.Thus, slopes having a north aspect were about five times more erosive.The same procedure was used for formation and location.These results indicated that the Poleslide was the most erosive formation and given the equal distribution of formations throughout the study region, the location erosion ratio was essentially the same.Accordingly, slope aspect was determined to be the most critical parameter to use in the site visitation schedule.Regression models were used to determine the erosion rate per year and the number of years for a slope to experience 2.54 cm of erosion.Years between site visits vary from ~1 year for north-facing slopes to ~6 years for south-facing slopes.Park officials will also be required to assess fossil types and protection standards when making final decisions for site visits.

Conclusions
Landforms at Badlands National Park have developed through long-term interactions of climate and geology.Material properties of the exposed formations governed erosional processes and were highly dependent on grain-size distributions and mineralogy within stratigraphic units.Variations in the silt to clay content were determining factors for erosion resistance of each formation.Increased clay content has been determined to increase PI, unit cohesion and increase erosion resistance.Increased silt content resulted in lower PI and decreased unit cohesion based on the S:M ratio.However, the added mechanical stability from increased silt enables near vertical cliffs to form, although these units have a decreased resistance to erosion and generate a greater mass flux than the clay-dominate units.
Analysis of material engineering properties of the bedrock has resulted in development of a grain-size proxy indicator, ψ, based the plasticity index, unit cohesion and the S:M ratio.This parameter was successfully substituted into Bagnold's entrainment equation to calculate critical threshold velocity that was plotted on a Hjulström diagram.These results indicated that valid entrainment velocities were able to be calculated based on the engineering properties of a fine-grained material.Plasticity index was observed to positively correlate to an erosivity resistance factor, k d and as PI increased, the erosion resistance of the material increased, reflected as a decrease in k d .These data and results suggested that slope erosion processes functioned as individual particles were disaggregated from bedrock during wetting and were transported downslope as sheetwash.Loose sediment material delivered to the toe of the slope produced a low-angle pediment leading away from the slope.Variations in surface roughness between slope and pediment was responsible for maintaining high velocity flows on the pediment that carried sediment away from the slope and prohibited sediment accumulation at the toe.The result was a general retreat of the slope through time producing an equilibrium slope profile based on geomechanical and mineralogical attributes of the units being weathered.Thus, equilibrium slope profiles at Badlands have developed as a direct response to the S:M ratio.
Geomechanical properties of the bedrock on all slopes produced stable profiles even under saturated conditions.Modeled slopes were sensitive to cohesion and contained up to 31% variance in FS under a 50% change in cohesion.Slope instability was visually observed where steep slopes exceeded ~20 m in height and become unstable from desiccation and tension fractures at the top of the slope and/or undercutting at the toe.None of the slopes modeled met either condition and were mechanically stable.
All results were used to develop a paleontological monitoring program for BNP managers based on the time (years) calculated for net slope erosion to equal 2.5 cm.The schedule was based on an erosion ratio that indicated erosion potential.North-facing slopes were ~5 times more erosive than south-facing slopes.This led to suggested site visits for north slopes every year and every 6 years for south slopes.Overall, coarser-grained rocks, such as Brule rocks, were 1/3 more erosive than fine-grained rocks.Thus, paleontological monitoring of fossil resources was determined to be highly variable but somewhat predictable.Erosion at BNP remains somewhat unknown given the limited scope of this study but was shown to be dependent upon the physical setting, bedrock properties and the climate.

Figure 1 .
Figure 1.Badlands National Park located in southwestern South Dakota consists of North and South administrative units.Erosion sites were located within the circles at known fossil sites in the North unit.

Figure 1 .
Figure 1.Badlands National Park located in southwestern South Dakota consists of North and South administrative units.Erosion sites were located within the circles at known fossil sites in the North unit.

1. 2 .
Regional SettingBadlands National Park, located in southwestern South Dakota (Figure1), is ~80 km east of the Black Hills extending eastward another ~100 km.The Park consists of 987 km 2 of sharply eroded buttes, pinnacles and spires set within and surrounded by the largest, protected mixed grass prairie in the United States [20].Physiographically, the Park consists of nearly flay-lying bedrock containing significant local relief due to variable erosion.The highest and lowest elevations in the Park are 1000 m on Sheep Mountain in the west and 750 m along Sage Creek, also in the west.The largest morphological feature in the Park is the Badlands 'Wall' (Figure 2), a 50 to 150 m high erosional scarp trending ~W-E for 100 km defined by steep gullied ridges with knife-like edges and pinnacles.Paleosols within the eroded sediment appear in various hues of red, yellow and purple.Geosciences 2018, 8, x FOR PEER REVIEW 3 of 17 1.2.Regional Setting Badlands National Park, located in southwestern South Dakota (Figure 1), is ~80 km east of the Black Hills extending eastward another ~100 km.The Park consists of 987 km 2 of sharply eroded buttes, pinnacles and spires set within and surrounded by the largest, protected mixed grass prairie in the United States [20].Physiographically, the Park consists of nearly flay-lying bedrock containing significant local relief due to variable erosion.The highest and lowest elevations in the Park are 1000 m on Sheep Mountain in the west and 750 m along Sage Creek, also in the west.The largest morphological feature in the Park is the Badlands 'Wall' (Figure 2), a 50 to 150 m high erosional scarp trending ~W-E for 100 km defined by steep gullied ridges with knife-like edges and pinnacles.Paleosols within the eroded sediment appear in various hues of red, yellow and purple.

Figure 2 .
Figure 2. A complete stratigraphic sequence for Badlands geology includes: the lower-most Yellow Mounds in the upper Cretaceous-Paleocene rocks (below orange line) and the immediately overlying red Interior paleosol.Tertiary Chadron formation that are gray shales weather into low-profile haystack-like mounds (between the orange and blue lines), the steep cliffs above Chadron rocks formed in the Brule formation.Overlying cap rocks consist of Miocene Arikaree Group rocks.Brule rock form the highly gullied and steep Badlands Wall.

Figure 2 .
Figure 2. A complete stratigraphic sequence for Badlands geology includes: the lower-most Yellow Mounds in the upper Cretaceous-Paleocene rocks (below orange line) and the immediately overlying red Interior paleosol.Tertiary Chadron formation that are gray shales weather into low-profile haystack-like mounds (between the orange and blue lines), the steep cliffs above Chadron rocks formed in the Brule formation.Overlying cap rocks consist of Miocene Arikaree Group rocks.Brule rock form the highly gullied and steep Badlands Wall.

Figure 3 .
Figure 3. Ice crystal expansion (below the scale) of the clay-rich Chadron Formation promotes disaggregation and formation of the aggregated popcorn-like texture (above the scale).

Figure 3 .
Figure 3. Ice crystal expansion (below the scale) of the clay-rich Chadron Formation promotes disaggregation and formation of the aggregated popcorn-like texture (above the scale).

Figure 4 .
Figure 4. (A) Erosion resistance correlates to plasticity index (PI) indicating fine-grained materials (open symbols) have greater erosion resistance (lower indicator of erosion resistance (kd)) than coarsegrained material (closed symbols).(B) Sediment mass flux from slopes support the use of the erodibility coefficient based on Atterberg Limits testing.Slopes in the western Badlands (open symbols) were finer-grained than those in the eastern Badlands (closed symbols) and yielded less sediment mass flux.

Figure 4 .
Figure 4. (A) Erosion resistance correlates to plasticity index (PI) indicating fine-grained materials (open symbols) have greater erosion resistance (lower indicator of erosion resistance (k d )) than coarse-grained material (closed symbols).(B) Sediment mass flux from slopes support the use of the erodibility coefficient based on Atterberg Limits testing.Slopes in the western Badlands (open symbols) were finer-grained than those in the eastern Badlands (closed symbols) and yielded less sediment mass flux.

Geosciences 2018, 8 , 17 Figure 5 .
Figure 5.Hjulström's diagram relates processes to velocities.Particles <0.1 mm require higher entrainment velocities as functions of grain size, shape and particle attractive forces.Badlands materials are plotted that were determined using the grain-size proxy.Open symbols had a smaller S:M ratio and were more clay-rich than closed symbol materials that had a larger S:M ratio.

Figure 5 .
Figure 5.Hjulström's diagram relates processes to velocities.Particles <0.1 mm require higher entrainment velocities as functions of grain size, shape and particle attractive forces.Badlands materials are plotted that were determined using the grain-size proxy.Open symbols had a smaller S:M ratio and were more clay-rich than closed symbol materials that had a larger S:M ratio.

Figure 6 .
Figure 6.Soil erodibility coefficient, kd and critical shear velocity, * , based on ψ, the grain-size proxy, hold physical significance for fine-grained materials and are both anchored by engineering properties of the materials.Open symbols have lower S:M ratios.

Figure 7 .
Figure 7. Slope retreat data for 2 slopes at site BL-01 for approximately 2 years.Slope 05R utilized erosion pins and slope 06R utilized steel rules.Precipitation is shown as lines.

Figure 6 . 17 Figure 6 .
Figure 6.Soil erodibility coefficient, k d and critical shear velocity, u * c , based on ψ, the grain-size proxy, hold physical significance for fine-grained materials and are both anchored by engineering properties of the materials.Open symbols have lower S:M ratios.

Figure 7 .
Figure 7. Slope retreat data for 2 slopes at site BL-01 for approximately 2 years.Slope 05R utilized erosion pins and slope 06R utilized steel rules.Precipitation is shown as lines.

Figure 7 .
Figure 7. Slope retreat data for 2 slopes at site BL-01 for approximately 2 years.Slope 05R utilized erosion pins and slope 06R utilized steel rules.Precipitation is shown as lines.

Funding:
This research was funded by a research grant from the US Department of Interior, National Park Service utilizing a Cooperative Ecosystems Studies Unit (CESU), Task Agreement J1300110051.Continued field access has been granted by Badlands National Park through 2018.

Table 1 .
Average geomechanical properties of slope wash and bedrock samples obtained from laboratory tests from each of the six erosion sites.

Table 2 .
Erosion modelling results for slope aspect, bedrock formation and location in the Park.The erosion ratio is based upon the most erosive parameter.Suggested years between visits were normalized on ~years for 2.54 cm of erosion to occur.