Most available biogeochemical models focus within a soil profile and cannot adequately resolve contributions of the lighter size fractions of organic rich soils for enrichment ratio (ER) estimates, thereby causing unintended errors in soil organic carbon (SOC) storage predictions. These models set ER as constant, usually equal to unity. The goal of this study is to provide spatiotemporal predictions of SOC stocks at the hillslope scale that account for the selective entrainment and deposition of lighter size fractions. It is hypothesized herein that ER values may vary depending on hillslope location, Land Use/Land Cover (LULC) conditions, and magnitude of the hydrologic event. An ER module interlinked with two established models, CENTURY and Watershed Erosion Prediction Project, is developed that considers the effects of changing runoff coefficients, bare soil coverage, tillage depth, fertilization, and soil roughness on SOC redistribution and storage. In this study, a representative hillslope is partitioned into two control volumes (CVs): a net erosional upslope zone and a net depositional downslope zone. We first estimate ER values for both CVs I and II for different hydrologic and LULC conditions. Second, using the improved ER estimates for the two CVs, we evaluate the effects that management practices have on SOC redistribution during different crop rotations. Overall, LULC promoting less runoff generally yielded higher ER values, which ranged between 0.97 and 3.25. Eroded soils in the upland CV were up to 4% more enriched in SOC than eroded soils in the downslope CV due to larger interrill contributions, which were found to be of equal importance to rill contributions. The chronosequence in SOC storage for the erosional zone revealed that conservation tillage and enhanced crop yields begun in the 1980s reversed the downward trend in SOC losses, causing nearly 26% of the lost SOC to be regained.

We introduce a new theoretical framework for knickpoint migration in unventilated flows prevalent in channelized semicohesive first‐ and second‐order streams characterized by a sequence of constrictions and expansions in the channel cross section. A key feature of the framework is the inclusion of shear‐induced fluvial erosion due to flow suction at the nappe as the key mechanism driving migration. Channel surveys, water stage, time‐lapse photography, and laser scans from Mud Creek, Iowa, USA, confirm this unventilated erosive mechanism and reveal a dual advective‐diffusive mode of retreat. We treat the fluvial bed shear stress as being similar to the flow on the lee side of a submerged obstacle and derive a governing equation which is a generalized Burgers' equation. The equation is solved to successfully simulate knickpoints monitored in two sites (Iowa and Mississippi), and a modified Peclet number is used to analyze the advective‐diffusive nature of the knickpoint migration.

An improved modeling framework for capturing the effects of space and time‐variant resistance to overland flow is developed for intensively managed landscapes. The framework builds on the WEPP model but it removes the limitations of the “equivalent” plane and time‐invariant roughness assumption. The enhanced model therefore accounts for spatiotemporal changes in flow resistance along a hillslope due to changes in roughness, in profile curvature, and downslope variability. The model is used to quantify the degree of influence—from individual soil grains to aggregates, “isolated roughness elements,” and vegetation—on overland flow characteristics under different storm magnitudes, downslope gradients, and profile curvatures. It was found that the net effects of land use change from vegetation to a bare surface resulted in hydrograph peaks that were up to 133% larger. Changes in hillslope profile curvature instead resulted in peak runoff rate changes that were only up to 16%. The stream power concept is utilized to develop a taxonomy that relates the influence of grains, isolated roughness elements, and vegetation, on overland flow under different storm magnitudes and hillslope gradients. Critical storm magnitudes and hillslope gradients were found beyond which the effects of these landscape attributes on the peak stream power were negligible. The results also highlight weaknesses of the space/time‐invariant flow resistance assumption and demonstrate that assumptions on landscape terrain characteristics exert a strong control both on the shape and magnitude of hydrographs, with deviations reaching 65% in the peak runoff when space/time‐variant resistance effects are ignored in some cases.

The role of tillage practices on soil aggregate properties has been mainly addressed at the pedon scale (i.e., soilscape scale) by treating landscape elements as disconnected. However, there is observed heterogeneity in aggregate properties along flowpaths, suggesting that landscape scale hydraulic processes are also important. This study examines this supposition using field, laboratory and modeling analysis to assess aggregate size and stability along flowpaths under different management conditions: (1) tillage-induced abrasion effects on aggregate size were evaluated with the dry mean weight diameter (DMWD); (2) raindrop impact effects were evaluated with small macroaggregate stability (SMAGG_STAB) using rainfall simulators; and (3) these aggregate proxies were studied in the context of connectivity through the excess bed shear stress (δ), quantified using a physically-based landscape model. DMWD and SMAGG_STAB decreased along the flowpaths for all managements, and a negative correspondence between the proxies and δ was observed. δ captured roughness effects on connectivity along the flowpaths: highest connectivity was noted for parallel-ridge-till flowpaths, where δ ranged from 0–8.2 Pa, and lowest connectivity for contour-ridge-till flowpaths, where δ ranged from 0–1.1 Pa. High tillage intensity likely led to an increase in aggregate susceptibility to hydraulic forcing, reflected in the higher gradients of aggregate size and stability trendlines with respect to δ. Finally, a linear relationship between DMWD and SMAGG_STAB was established.

This study aimed to better understand how tillage row orientation with respect to dominant flow‐pathway along hillslope impacts runoff and the transport of different sediment size fractions. Experimental plots were constructed in contour ridge till (CRT) and parallel ridge till (PRT) sites to monitor runoff and sediment fluxes. Particle size fractions of the sediment, along with organic C and N contents were measured to quantify physical and chemical enrichment ratios for the two tillage orientations. In the CRT plot, tillage produced large oriented roughness elements along the contours, which acted as little check dams, while in the PRT site, the roughness elements helped confine and concentrate the runoff. In the CRT, the “check dams” resulted in a runoff coefficient of .03, while in the PRT, the flow confinement between rows produced a runoff coefficient of .82. Moreover, the erosion rates at the CRT site were 97% less than those in the PRT plot. Large contours produced finer sediment fractions due to the selective sorting in ponded furrows. More aggregate sediment fractions were present in the PRT site, which was dominated by rill erosion. Physical enrichment revealed selective entrainment of finer sediment particles during erosion. Finer‐grain particles with higher specific surface areas, attached more organic C resulting in chemical enrichment. Physical and chemical enrichment methods were found to be in good agreement. These findings suggest transport models that can simulate size fraction updates to the soil active layer can be used to estimate SOC redistribution and hence more accurate C budgets.