Rates of erosion and landscape change along the Blue Ridge escarpment,southern Appalachian Mountains,estimated from in situ cosmogenic 10Be

The Blue Ridge escarpment, located within the southern Appalachian Mountains of Virginia and North Carolina, forms a distinct, steep boundary between the lower‐elevation Piedmont and higher‐elevation Blue Ridge physiographic provinces. To understand better the rate at which this landform and the adjacent landscape are changing, we measured cosmogenic beryllium‐10 (¹⁰Be) in quartz separated from sediment samples (n = 50) collected in 32 streams and from three exposed bedrock outcrops along four transects normal to the escarpment, allowing us to calculate erosion rates integrated over 10⁴–10⁵ years. These basin‐averaged erosion rates (5.4–49 m Myr^?1) are consistent with those measured elsewhere in the southern Appalachain Mountains and show a positive relationship between erosion rate and average basin slope. Erosion rates show no relationship with basin size or relative position of the Brevard fault zone, a fundamental structural element of the region. The cosmogenic isotopic data, when considered along with the distribution of average basin slopes in each physiographic province, suggest that the escarpment is eroding on average more rapidly than the Blue Ridge uplands, which are eroding more rapidly than the Piedmont lowlands. This difference in erosion rates by geomorphic setting suggests that the elevation difference between the uplands and lowlands adjacent to the escarpment is being reduced but at extremely slow rates. Copyright © 2016 John Wiley & Sons, Ltd.     

Controls on deep critical zone architecture: a historical review and four testable hypotheses

The base of Earth's critical zone (CZ) is commonly shielded from study by many meters of overlying rock and regolith. Though deep CZ processes may seem far removed from the surface, they are vital in shaping it, preparing rock for infusion into the biosphere and breaking Earth materials down for transport across landscapes. This special issue highlights outstanding challenges and recent advances of deep CZ research in a series of articles that we introduce here in the context of relevant literature dating back to the 1500s. Building on several contributions to the special issue, we highlight four exciting new hypotheses about factors that drive deep CZ weathering and thus influence the evolution of life‐sustaining CZ architecture. These hypotheses have emerged from recently developed process‐based models of subsurface phenomena including: fracturing related to subsurface stress fields; weathering related to drainage of bedrock under hydraulic head gradients; rock damage from frost cracking due to subsurface temperature gradients; and mineral reactions with reactive fluids in subsurface chemical potential gradients. The models predict distinct patterns of subsurface weathering and CZ thickness that can be compared with observations from drilling, sampling and geophysical imaging. We synthesize the four hypotheses into an overarching conceptual model of fracturing and weathering that occurs as Earth materials are exhumed to the surface across subsurface gradients in stress, hydraulic head, temperature, and chemical potential. We conclude with a call for a coordinated measurement campaign designed to comprehensively test the four hypotheses across a range of climatic, tectonic and geologic conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

Morphodynamics of a bedrock‐alluvial meander bend that incises as it migrates outward: approximate solution of permanent form

In meandering rivers cut into bedrock, erosion across a channel cross‐section can be strongly asymmetric. At a meander apex, deep undercutting of the outer bank can result in the formation of a hanging cliff (which may drive hillslope failure), whereas the inner bank adjoins a slip‐off slope that connects to the hillslope itself. Here we propose a physically‐based model for predicting channel planform migration and incision, point bar and slip‐off slope formation, bedrock abrasion, the spatial distribution of alluvial cover, and adaptation of channel width in a mixed bedrock‐alluvial channel. We simplify the analysis by considering a numerical model of steady, uniform bend flow satisfying cyclic boundary conditions. Thus in our analysis, ‘sediment supply’, i.e. the total volume of alluvium in the system, is conserved. In our numerical simulations, the migration rate of the outer bank is a specified parameter. Our simulations demonstrate the existence of an approximate state of dynamic equilibrium corresponding to a near‐solution of permanent form in which a bend of constant curvature, width, cross‐sectional shape and alluvial cover distribution migrates diagonally downward at constant speed, leaving a bedrock equivalent of a point bar on the inside of the bend. Channel width is set internally by the processes of migration and incision. We find that equilibrium width increases with increasing sediment supply, but is insensitive to outer bank migration rate. The slope of the bedrock point bar varies inversely with both outer bank migration rate and sediment supply. Although the migration rate of the outer bank is externally imposed here, we discuss a model modification that would allow lateral side‐wall abrasion to be treated in a manner similar to the process of bedrock incision. Copyright © 2016 John Wiley & Sons, Ltd.     

Late Pleistocene outburst floods from Issyk Kul,Kyrgyzstan?

Elevated shorelines and lake sediments surrounding Issyk Kul, the world's second largest mountain lake, record fluctuating lake levels during Quaternary times. Together with bathymetric and geochemical data, these markers document alternating phases of lake closure and external drainage. The uppermost level of lake sediments requires a former damming of the lake's western outlet through the Boam gorge. We test previous hypothesised ice or landslide dam failures by exploring possible links between late Quaternary lake levels and outbursts. We review and recompile the chronology of reported changes in lake site, and offer new ages of abandoned shorelines using ¹⁴C in bivalve and gastropod shells, and plant detritus, as well as sand lenses in delta and river sediments using Infrared Stimulated Luminescence. Our dates are consistent with elevated lake levels between ~45 ka and 22 ka. Cosmogenic ¹⁰Be and ²⁶Al exposure ages of fan terraces containing erratic boulders (>3 m) downstream of the gorge constrain the timing of floods to 20.5–18.5 ka, postdating a highstand of Issyk Kul. A flow‐competence analysis gives a peak discharge of >10⁴ m³ s^–1 for entraining and transporting these boulders. Palaeoflood modelling, however, shows that naturally dammed lakes unconnected to Issyk Kul could have produced such high discharges upon sudden emptying. Hence, although our data are consistent with hypotheses of catastrophic outburst floods, average lake‐level changes of up to 90 mm yr^–1 in the past 150 years were highly variable without any outbursts, so that linking lake‐level drops to catastrophic dam breaks remains ambiguous using sedimentary archives alone. This constraint may readily apply to other Quaternary lakes of that size elsewhere. Nonetheless, our reconstructed Pleistocene floods are among the largest reported worldwide, and motivate further research into the palaeoflood hydrology of Central Asia. Copyright © 2017 John Wiley & Sons, Ltd.     

Sediment transport in high‐speed flows over a fixed bed: 2. Particle impacts and abrasion prediction

Single bed load particle impacts were experimentally investigated in supercritical open channel flow over a fixed planar bed of low relative roughness height simulating high‐gradient non‐alluvial mountain streams as well as hydraulic structures. Particle impact characteristics (impact velocity, impact angle, Stokes number, restitution and dynamic friction coefficients) were determined for a wide range of hydraulic parameters and particle properties. Particle impact velocity scaled with the particle velocity, and the vertical particle impact velocity increased with excess transport stage. Particle impact and rebound angles were low and decreased with transport stage. Analysis of the particle impacts with the bed revealed almost no viscous damping effects with high normal restitution coefficients exceeding unity. The normal and resultant Stokes numbers were high and above critical thresholds for viscous damping. These results are attributed to the coherent turbulent structures near the wall region, i.e. bursting motion with ejection and sweep events responsible for turbulence generation and particle transport. The tangential restitution coefficients were slightly below unity and the dynamic friction coefficients were lower than for alluvial bed data, revealing that only a small amount of horizontal energy was transferred to the bed. The abrasion prediction model formed by Sklar and Dietrich in 2004 was revised based on the new equations on vertical impact velocity and hop length covering various bed configurations. The abrasion coefficient k_v was found to be vary around k_v ~ 10⁵ for hard materials (tensile strength f_t > 1 MPa), one order of magnitude lower than the value assumed so far for Sklar and Dietrich's model. Copyright © 2017 John Wiley & Sons, Ltd.     

Potential wind erosion rate response to climate and land‐use changes in the watershed of the Ningxia–Inner Mongolia reach of the Yellow River,China, 1986–2013

Climate and land‐use changes could strongly affect wind erosion and in turn cause a series of environmental problems. Thus, the objective of this study was to assess potential wind erosion rate (PWER) response to climate and land‐use changes in the watershed of the Ningxia–Inner Mongolia Reach of the Yellow River (NIMRYR), China. The watershed of NIMRYR suffers from serious wind erosion hazards, and over recent decades, wind erosion intensity and distribution has changed, following climate and land‐use changes. To understand these processes in the NIMRYR watershed, the Integrated Wind Erosion Modelling System (IWEMS) and the Revised Wind Erosion Equation (RWEQ) were used to calculate the PWER under different climate conditions and land‐use scenarios, and to assess the influences of climate and land‐use changes on the PWER. The results show the PWER in the whole watershed had a significant declining trend from 1986 to 2013. The results of the relationship among PWER, climate change, and land‐use changes showed that climate change was the dominant control on the PWER change in this watershed. Compared to the period 1986–1995, the average PWER decreased 23.32% and 64.98% as a result of climate change in the periods 1996–2005 and 2006–2013, respectively. In contrast with climate change, the effects of land‐use changes on the average PWER were much lower, and represented a change in PWER of less than 3.3% across the whole watershed. The study method we used could provide some valuable reference for wind erosion modelling, and the research results should help climate and land‐use researchers to develop strategies to reduce wind erosion. Copyright © 2017 John Wiley & Sons, Ltd.     

Variation in reach‐scale thalweg‐migration intensity in a braided reach of the lower Yellow River in 1986–2015

Thalweg migration of an alluvial river plays a key role in channel evolution, which may influence the effect of existing river training works and biodiversity on floodplains, and cause losses in riparian land and property. The braided reach of the Lower Yellow River underwent continuous channel aggradation during the period from 1986 to 1999, and then remarkable channel degradation in 1999–2015 owing to the state of operation of the Xiaolangdi Reservoir in 1999. Here we quantify associated thalweg migration changes and identify the key influencing factor in the braided reach. Thalweg‐migration distances and intensities at section‐ and reach‐scales were calculated during the past 30 years from 1986 to 2015, in order to investigate the characteristics of thalweg migration in the reach. There was a 47% reduction in the reach‐scale thalweg‐migration distance and a 35% reduction in the corresponding migration intensity after the reservoir operation. It is also revealed that fluvial erosion intensity is a dominant factor in controlling the thalweg migration, based on the investigation into various influencing factors in the study reach. The thalweg‐migration intensity of the braided reach can be expressed as a power function of the previous four‐year average fluvial erosion intensity. The calculated thalweg‐migration intensities in 1986–2015 using the proposed relation generally agree with the observed data. Copyright © 2017 John Wiley & Sons, Ltd.     

Frequency analysis of nonstationary annual maximum flood series using the time‐varying two‐component mixture distributions

The most popular practice for analysing nonstationarity of flood series is to use a fixed single‐type probability distribution incorporated with the time‐varying moments. However, the type of probability distribution could be both complex because of distinct flood populations and time‐varying under changing environments. To allow the investigation of this complex nature, the time‐varying two‐component mixture distributions (TTMD) method is proposed in this study by considering the time variations of not only the moments of its component distributions but also the weighting coefficients. Having identified the existence of mixed flood populations based on circular statistics, the proposed TTMD was applied to model the annual maximum flood series of two stations in the Weihe River basin, with the model parameters calibrated by the meta‐heuristic maximum likelihood method. The performance of TTMD was evaluated by different diagnostic plots and indexes and compared with stationary single‐type distributions, stationary mixture distributions and time‐varying single‐type distributions. The results highlighted the advantages of TTMD with physically‐based covariates for both stations. Besides, the optimal TTMD models were considered to be capable of settling the issue of nonstationarity and capturing the mixed flood populations satisfactorily. Copyright © 2016 John Wiley & Sons, Ltd.     

Synthesis of soil‐hydraulic properties and infiltration timescales in wildfire‐affected soils

We collected soil‐hydraulic property data from the literature for wildfire‐affected soils, ash, and unburned soils. These data were used to calculate metrics and timescales of hydrologic response related to infiltration and surface runoff generation. Sorptivity (S) and wetting front potential (Ψ_f) were significantly different (lower) in burned soils compared with unburned soils, whereas field‐saturated hydraulic conductivity (K_fs) was not significantly different. The magnitude and duration of the influence of capillarity during infiltration was greatly reduced in burned soils, causing faster ponding times in response to rainfall. Ash had large values of S and K_fs but moderate values of Ψ_f, compared with unburned and burned soils, indicating ash has long ponding times in response to rainfall. The ratio of S²/K_fs was nearly constant (~100 mm) for unburned soils but more variable in burned soils, suggesting that unburned soils have a balance between gravity and capillarity contributions to infiltration that may depend on soil organic matter, whereas in burned soils the gravity contribution to infiltration is greater. Changes in S and K_fs in burned soils act synergistically to reduce infiltration and accelerate and amplify surface runoff generation. Synthesis of these findings identifies three key areas for future research. First, short timescales of capillary influences on infiltration indicate the need for better measurements of infiltration at times less than 1 min to accurately characterize S in burned soils. Second, using parameter values, such as Ψ_f, from unburned areas could produce substantial errors in hydrologic modeling when used without adjustment for wildfire effects, causing parameter compensation and resulting underestimation of K_fs. Third, more thorough measurement campaigns that capture soil‐structural changes, organic matter impacts, quantitative water repellency trends, and soil‐water content along with soil‐hydraulic properties could drive the development of better techniques for numerically simulating infiltration in burned areas.     

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.