Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research, Earth Surface”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3Linked to original sources
USGS2023

Time-lapse seafloor surveys reveal how turbidity currents and internal tides in Monterey Canyon interact with the seabed at centimeter-scale

Here we show how ultra-high resolution seabed mapping using new technology can help to understand processes that sculpt submarine canyons. Time-lapse seafloor surveys were conducted in the axis of Monterey Canyon, ∼50 km from the canyon head (∼1,840 m water depth) over an 18-month period. These surveys comprised 5-cm resolution multibeam bathymetry, 1-cm resolution lidar bathymetry, and 2-mm resolution stereophotographic imagery. Bathymetry data reveal centimeter-scale textures that would be undetectable by more traditional survey methods. Upward-looking Acoustic Doppler Current Profilers at the site recorded the flow character of internal tides and the passage of three turbidity currents, while sediment cores collected from the site record flow deposits. Combined with flow and core data, the bathymetry shows how turbidity currents and internal tides modify the seabed. The turbidity currents drape sediment across the site, infilling bedform troughs and smoothing erosional features carved by the internal tides (e.g., rippled scours). Turbidity currents with speeds of 0.9–3.3 m/s failed to cause notable bedform movement, which is surprising given that flows with similar speeds produced rapid bedform migration elsewhere, including the upper Monterey Canyon. The lack of migration may be related to the character of the underlying substrate or indicate that turbidity currents at the site lack dense, near-bed layers. The scale of scours produced by the internal tides (≤0.7 m/s) approaches the scale of features recorded in the ancient rock record. Thus, these results illustrate how the scale gap between seabed mapping technology and the rock record may eventually be bridged.

California↗
USGS2023

Mapping landslide susceptibility over large regions with limited data

Landslide susceptibility maps indicate the spatial distribution of landslide likelihood. Modeling susceptibility over large or diverse terrains remains a challenge due to the sparsity of landslide data (mapped extent of known landslides) and the variability in triggering conditions. Several different data sampling strategies of landslide locations used to train a susceptibility model are used to mitigate this challenge. However, to our knowledge, no study has systematically evaluated how different sampling strategies alter a model's predictor effects (i.e., how a predictor value influences the susceptibility output) critical to explaining differences in model outputs. Here, we introduce a statistical framework that examines the variation in predictor effects and the model accuracy (measured using receiver operator characteristics) to highlight why certain sampling strategies are more effective than others. Specifically, we apply our framework to an array of logistic regression models trained on landslide inventories collected at sub-regional scales over four terrains across the United States. Results show significant variations in predictor effects depending on the inventory used to train the models. The inconsistent predictor effects cause low accuracies when testing models on inventories outside the domain of the training data. Grouping test and training sets according to physiographic and ecological characteristics, which are thought to share similar triggering mechanisms, does not improve model accuracy. We also show that using limited landslide data distributed uniformly over the entire modeling domain is better than using dense but spatially isolated data to train a model for applications over large regions.

California, Kentucky, North Carolina, West Virgini↗
USGS2022

The influence of drying on the aeolian transport of river-sourced sand

Transgression and regression of water levels (stages) have impacted the evolution of aeolian landforms and sedimentary deposits throughout geologic history. We studied this phenomenon over a five-day period of reduced flow on the Colorado River in Grand Canyon National Park, AZ, USA, in March 2021. These transient low flows exposed river-channel sand deposits to the air, causing progressive desiccation (drying) and thereby making these deposits susceptible to aeolian transport. We measured aeolian threshold friction velocities ( u *t ) for sand saltation and PM10 dust emissions, as well as other characteristics, on a subaerially exposed sandbar and downwind aeolian dunefield during each day of the low river flow. The sandbar transitioned from supply-limited to transport-limited aeolian sediment transport conditions during the regression in river water stage. A possible tipping point between the two transport conditions occurred approximately 48 hours after the drop in river flow. The empirically measured u *t decreased as the sandbar sediment dried with increased subaerial exposure time. Theoretical estimates and empirical measurements of u *t corresponded closely on the aeolian dunefield and on the sandbar when it was drier during the third and fourth day of the experiment. Eighty-seven percent of the variability in u *t was explained by empirical models that provide practical estimates of aeolian transport potential of subaerial river sediment deposits using monitoring data that are commonly available in this and other river systems. The work provides theoretical insight into the response of aeolian processes to sediment supply changes driven by periods of anthropogenic activity, drought, and climate change.

Arizona↗
USGS2023

Impacts of spontaneous waterfall development on bedrock river longitudinal profile morphology

River profiles are shaped by climatic and tectonic history, lithology, and internal feedbacks between flow hydraulics, sediment transport and erosion. In steep channels, waterfalls may self-form without changes in external forcing (i.e., autogenic formation) and erode at rates faster or slower than an equivalent channel without waterfalls. We use a 1-D numerical model to investigate how self-formed waterfalls alter the morphology of bedrock river longitudinal profiles. We modify the standard stream power model to include a slope threshold above which waterfalls spontaneously form and a rate constant allowing waterfalls to erode faster or slower than other fluvial processes. Using this model, we explore how waterfall formation alters both steady state and transient longitudinal profile forms. Our model predicts that fast waterfalls create km-scale reaches in a dynamic equilibrium with channel slope held approximately constant at the threshold slope for waterfall formation, while slow waterfalls can create local channel slope maxima at the location of slow waterfall development. Furthermore, slow waterfall profiles integrate past base level histories, leading to multiple possible profile forms, even at steady-state. Consistency between our model predictions and field observations of waterfall-rich rivers in the Kings and Kaweah drainages in the southern Sierra Nevada, California, supports the hypothesis that waterfall formation can modulate river profiles in nature. Our findings may help identify how bedrock channels are influenced by waterfall erosion and aid in distinguishing between signatures of external and internal perturbations, thereby strengthening our ability to interpret past climate and tectonic changes from river longitudinal profiles.

Journal of Geophysical Research - Earth Surface↗
USGS2023

Kinematic evolution of a large paraglacial landslide in the Barry Arm fjord of Alaska

Our warming climate is adversely affecting cryospheric landscapes via glacial retreat, permafrost degradation, and associated slope destabilization. In Prince William Sound, Alaska, the rapid retreat of Barry Glacier has destabilized the slopes flanking the glacier, resulting in numerous landslides. The largest of these landslides (∼500 Mm 3 in volume) is more than 2 km wide and has the potential to generate a tsunami that could affect nearby recreationists, marine traffic, infrastructure, natural and cultural resources, and the community of Whittier, located 60 km from the landslide. Here, we combine landslide structural and kinematic element mapping with data acquired from bi-yearly airborne lidar, multi-week satellite-based synthetic aperture radar (SAR), sub-hourly ground-based SAR, and seismic monitoring from 2020 to 2022 to characterize this landslide and examine its evolution. While some methods serve as a snapshot in time that is a culmination of events, others emphasize the ever-evolving nature of the landslide and associated hazards. Four major kinematic elements define the overall structure of the landslide, which vary in deformation type and rate, from creep (5 mm per day over several months) to episodic movement (2 m in 30 days) and landslide-wide to localized events. In some areas of the landslide, short-term deformation deviates from structures formed by cumulative movement, implying structural and kinematic evolution associated with glacier retreat. These insights are important for assessing landslide hazards and hazard evolution for large, slow-moving bedrock landslides in actively deglaciating environments.

Alaska↗
USGS2023

Tree-ring derived avalanche frequency and climate associations in a high-latitude, maritime climate

Snow avalanches are a natural hazard in mountainous areas worldwide with severe impacts that include fatalities, damage to infrastructure, disruption to commerce, and landscape disturbance. Understanding long-term avalanche frequency patterns, and associated climate and weather influences, improves our understanding of how climate change may affect avalanche activity. We used dendrochronological techniques to evaluate the historical frequency of large magnitude avalanches (LMAs) in the high-latitude climate of southeast Alaska, United States. We collected 434 cross sections throughout six avalanche paths near Juneau, Alaska. This resulted in 2706 identified avalanche growth disturbances between 1720 and 2018, which allowed us to reconstruct 82 years with LMA activity across three sub-regions. By combining this tree-ring-derived avalanche data set with a suite of climate and atmospheric variables and applying a generalized linear model to fit a binomial regression, we found that February and March precipitation and the Oceanic Niño Index (ONI) were significant predictors of LMA activity in the study area. Specifically, LMA activity occurred during winters with substantial February and March precipitation and neutral or negative (cold) ONI values, while years not characterized by LMAs occur more frequently during warm winters (positive ONI values). Our examination of the climate-avalanche relationship in southeast Alaska sheds light on important climate variables and physical processes associated with LMA years. These results can be used to inform long-term infrastructure planning and avalanche mitigation operations in an urban area, such as Juneau, where critical infrastructure is subject to substantial avalanche hazard.

Alaska↗
USGS2024

Predicting characteristic length scales of barrier island segmentation in microtidal environments

Segmented barrier islands can be found in regions with small tidal ranges. In contrast to tidally dominated barriers, where inlet dynamics are thought to control island length scales, the controls on barrier island length scales in wave-dominated environments have not been quantified. These microtidal barriers typically have a curved shoreline, suggesting the influence of wave-driven alongshore sediment transport. Microtidal barriers are also typically hydrodynamically isolated from one another, as weak tidal flows limit interactions between adjoining barriers. To better understand the controls on and scales of barrier segmentation in the relative absence of tides, here we develop a theoretical framework to estimate the alongshore length scales at which a barrier will either breach or heal following a disturbance in the barrier morphology. The non-dimensional framework compares the timescales of overwash (advective) and alongshore sediment transport (diffusive) processes along barrier island chains. We then apply this framework to modern barrier islands in the microtidal Gulf of Mexico using wave hindcast data and the lengths, widths, heights, and lagoon depths measured from remotely sensed geospatial data and topobathymetric data. We find that most of these barriers are currently longer than their critical length scale, often as a result of coastal restoration efforts. Our critical length scale analysis suggests that most of the Gulf of Mexico barriers are vulnerable to segmentation despite coastal restoration efforts intended to protect fisheries and the mainland coasts.

Florida, Louisiana, Mississippi↗
USGS2024

Vortex trapping of suspended sand grains over ripples

Coastal hydrodynamics and morphodynamics integrate the effects of small-scale fluid-sediment interactions; yet, these small-scale processes are not well understood. To investigate sediment trapping by turbulent coherent structures or vortices, the transport of coarse sand over ripples was analyzed in a small-oscillatory flow tunnel with phase-separated Particle Image and Tracking Velocimetry. Results from one of the first direct measurements of vortex-trapped sand grains under oscillatory flows are presented. The vortices mobilized sand grains along the ripple slopes just prior to flow reversal and transported the suspended sediment grains. During several flow cycles, some sand grains were temporarily trapped in the vortex, prescribing semi-circular trajectories off-center from the vortex core in quadrants of the vortex that were closest to the ripple slope, as illustrated by Nielsen (1992, https://doi.org/10.1142/1269 ). Comparisons of the horizontal sediment grain velocity with the horizontal fluid velocity yielded a linear relationship with a slope of 0.87. The vertical grain velocities also varied linearly with the vertical fluid velocity with a slope of approximately 1 and an offset of −0.08 m s −1 . The offset is close to the still water settling velocity for coarse sand grains, as hypothesized during vortex trapping. Additionally, estimates of the off-center distance, between the centers of the semi-circular sediment paths and vortex cores, compared well with the ratio of the settling velocity to the radian frequency of the vortex yielding a linear regression slope of 0.99. Improved understanding of vortex trapping effects on sediment dynamics may decrease uncertainty in model predictions of large-scale coastal hydrodynamics and sediment transport.

Journal of Geophysical Research: Earth Surface↗
USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

Oregon↗
USGS2024

Postfire sediment mobilization and its downstream implications across California, 1984 – 2021

Fire facilitates erosion through changes in vegetation and soil, with major postfire erosion commonly occurring even with moderate rainfall. As climate warms, the western United States (U.S.) is experiencing an intensifying fire regime and increasing frequency of extreme rain. We evaluated whether these hydroclimatic changes are evident in patterns of postfire erosion by modeling hillslope erosion following all wildfires larger than 100 km 2 in California from 1984 to 2021. Our results show that annual statewide postfire hillslope erosion has increased significantly over time. To supplement the hillslope erosion modeling, we compiled modeled and measured postfire debris-flow volumes. We find that, in northern California, more than 50% of fires triggering the top 20 values of sediment mass and sediment yield occurred in the most recent decade (between 2011 and 2021). In southern California, the postfire sediment budget was dominated by debris flows, which showed no temporal trend. Our analysis reveals that 57% of postfire sediment erosion statewide occurred upstream of reservoirs, indicating potential impacts to reservoir storage capacity and thus increased risk to water-resource security with ongoing climate change.

California↗
USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗
USGS2025

Uncertainty reduction for subaerial landslide-tsunami hazards

Subaerial rock slopes may generate a tsunami by rapidly moving into the water. Large uncertainty in landslide characteristics propagates into large uncertainty in tsunami hazard, making hazard assessment more difficult for land and emergency managers. Once a potentially tsunamigenic landslide is identified, it may not be clear which landslide characteristics contribute most significantly to uncertainty in the tsunami hazard. Our aim is to document the relative worth of different landslide characteristics (e.g., size, material properties) for reducing uncertainty in landslide-tsunami hazard assessments. Isolating the relative importance of specific landslide characteristics may inform prioritization of data collection and improve efficiency in understanding hazard. To accomplish this, we generated a set of 288 landslide-tsunami simulations in which we systematically varied the size and material properties of possible failure extents at the Barry Arm landslide complex in northwestern Prince William Sound, Alaska, USA. We find that for landslides smaller than the receiving waterbody, the landslide volume has the strongest effect on resulting wave characteristics and thus the highest leverage on reducing uncertainty in tsunami hazard. In contrast, for landslides substantially larger than the waterbody, the duration of rapid movement of the landslide has the highest leverage. Based on our results, we propose a classification scheme for subaerial landslides based on the relative size of the landslide and waterbody. Additionally, our results support the generation of a tsunami height transfer function between existing tide gages and a nearby coastal city. These results have direct implications for the practice of operational early warning.

Alaska↗
USGS2026

Rapid seismic and infrasound assessment of large landslides: A case study from Denali National Park and Preserve (Alaska)

Large, rapid landslides are a global hazard that can occur in remote, mountainous areas. Eyewitness reports of landslides and satellite imagery can often be limited or delayed, particularly during inclement weather. However, landslide-generated seismic and infrasound (low-frequency atmospheric sound) waves can be remotely detected in near real-time. This information can significantly expedite characterization and possible landslide response activities. Here, we highlight these capabilities using a > 4 million m³ ice–rock avalanche in Denali National Park and Preserve (Alaska). This event was detected via a landslide-specific seismic location and volume estimation algorithm deployed in Alaska, and — notably — by standard earthquake monitoring systems. Following rapid detection of this event, we combined its seismic and infrasound dataset with optical, synthetic aperture radar, and oblique aerial imagery, multitemporal digital elevation models, and a numerical flow model to reconstruct its failure timeline and dynamics. We apply array processing to infrasound signals traveling > 250 km and find that two precursory events occurred minutes prior to the main failure. We use long-period seismic signals to infer the force exerted by the landslide on the Earth and constrain the rheological parameters of our numerical flow simulation with this result and deposit morphology. The main failure produced a steeply-dipping impulsive initial downward force and reached speeds exceeding 60 m/s. This impulsive force generated relatively strong seismic body waves, which contributed to the earthquake system detection. This large, remote Alaska landslide underscores the key value of seismic and infrasound analysis for rapid landslide assessment and motivates efforts to further operationalize these approaches.

Alaska↗
USGS2017

A framework for modeling connections between hydraulics, water surface roughness, and surface reflectance in open channel flows

This paper introduces a framework for examining connections between the flow field, the texture of the air-water interface, and the reflectance of the water surface and thus evaluating the potential to infer hydraulic information from remotely sensed observations of surface reflectance. We used a spatial correlation model describing water surface topography to illustrate the application of our framework. Nondimensional relations between model parameters and flow intensity were established based on a prior flume study. Expressing the model in the spatial frequency domain allowed us to use an efficient Fourier transform-based algorithm for simulating water surfaces. Realizations for both flume and field settings had water surface slope distributions positively correlated with velocity and water surface roughness. However, most surface facets were gently sloped and thus unlikely to yield strong specular reflections; the model exaggerated the extent of water surface features, leading to underestimation of facet slopes. A ray tracing algorithm indicated that reflectance was greatest when solar and view zenith angles were equal and the sensor scanned toward the Sun to capture specular reflections of the solar beam. Reflected energy was concentrated in a small portion of the sky, but rougher water surfaces reflected rays into a broader range of directions. Our framework facilitates flight planning to avoid surface-reflected radiance while mapping other river attributes, or to maximize this component to exploit relationships between hydraulics and surface reflectance. This initial analysis also highlighted the need for improved models of water surface topography in natural rivers.

Journal of Geophysical Research F: Earth Surface↗
USGS2002

A bibliography of terrain modeling (geomorphometry), the quantitative representation of topography: Supplement 4.0

Terrain modeling, the practice of ground-surface quantification, is an amalgam of Earth science, mathematics, engineering, and computer science. The discipline is known variously as geomorphometry (or simply morphometry), terrain analysis, and quantitative geomorphology. It continues to grow through myriad applications to hydrology, geohazards mapping, tectonics, sea-floor and planetary exploration, and other fields. Dating nominally to the co-founders of academic geography, Alexander von Humboldt (1808, 1817) and Carl Ritter (1826, 1828), the field was revolutionized late in the 20th Century by the computer manipulation of spatial arrays of terrain heights, or digital elevation models (DEMs), which can quantify and portray ground-surface form over large areas (Maune, 2001). Morphometric procedures are implemented routinely by commercial geographic information systems (GIS) as well as specialized software (Harvey and Eash, 1996; Köthe and others, 1996; ESRI, 1997; Drzewiecki et al., 1999; Dikau and Saurer, 1999; Djokic and Maidment, 2000; Wilson and Gallant, 2000; Breuer, 2001; Guth, 2001; Eastman, 2002). The new Earth Surface edition of the Journal of Geophysical Research, specializing in surficial processes, is the latest of many publication venues for terrain modeling. This is the fourth update of a bibliography and introduction to terrain modeling (Pike, 1993, 1995, 1996, 1999) designed to collect the diverse, scattered literature on surface measurement as a resource for the research community. The use of DEMs in science and technology continues to accelerate and diversify (Pike, 2000a). New work appears so frequently that a sampling must suffice to represent the vast literature. This report adds 1636 entries to the 4374 in the four earlier publications1. Forty-eight additional entries correct dead Internet links and other errors found in the prior listings. Chronicling the history of terrain modeling, many entries in this report predate the 1999 supplement. Coverage is representative from about 1800 through early–mid 2002. Papers increasingly are published exclusively or in duplicate on the Internet's World Wide Web; the dates given here for Web addresses (URLs) that lack a print publication indicate a Web site's last update or my last access of it. The bibliography is arranged alphabetically and thus is not readily summarized. This introduction cites about 500 entries, a third of them grouped under 24 morphometric topics, as a guide to the listing's contents. Continuing the practice of previous bibliographies in the series to provide more information on a few applications (see summary of past topics in Pike, 2000a), this report elaborates further on topographic data, putative new parameters, tectonic geomorphology/neo-orometry, biogeography, ice-cap morphometry, results from the Mars Global DEM, landslide-hazard mapping, terrain modeling as physics, Hack's law, and broad-scale computer visualization. The literature of some of these subjects is large, and none of the summaries is intended to more than introduce the topic and comment on some of the current contributions of terrain modeling. Closing the essay is a discussion of pre-1900 papers that trace the evolution of ridge-line and watercourse quantification by descriptive geometry, as well as comments on some new books and an on-line bulletin board.

Open-File Report↗
USGS2010

A rain splash transport equation assimilating field and laboratory measurements

Process-based models of hillslope evolution require transport equations relating sediment flux to its major controls. An equation for rain splash transport in the absence of overland flow was constructed by modifying an approach developed by Reeve (1982) and parameterizing it with measurements from single-drop laboratory experiments and simulated rainfall on a grassland in East Africa. The equation relates rain splash to hillslope gradient, the median raindrop diameter of a storm, and ground cover density; the effect of soil texture on detachability can be incorporated from other published results. The spatial and temporal applicability of such an equation for rain splash transport in the absence of overland flow on uncultivated hillslopes can be estimated from hydrological calculations. The predicted transport is lower than landscape-averaged geologic erosion rates from Kenya but is large enough to modify short, slowly eroding natural hillslopes as well as microtopographic interrill surfaces between which overland flow transports the mobilized sediment. Copyright 2010 by the American Geophysical Union. Copyright 2010 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2013

Evaluation of a new model of aeolian transport in the presence of vegetation

Aeolian transport is an important characteristic of many arid and semiarid regions worldwide that affects dust emission and ecosystem processes. The purpose of this paper is to evaluate a recent model of aeolian transport in the presence of vegetation. This approach differs from previous models by accounting for how vegetation affects the distribution of shear velocity on the surface rather than merely calculating the average effect of vegetation on surface shear velocity or simply using empirical relationships. Vegetation, soil, and meteorological data at 65 field sites with measurements of horizontal aeolian flux were collected from the Western United States. Measured fluxes were tested against modeled values to evaluate model performance, to obtain a set of optimum model parameters, and to estimate the uncertainty in these parameters. The same field data were used to model horizontal aeolian flux using three other schemes. Our results show that the model can predict horizontal aeolian flux with an approximate relative error of 2.1 and that further empirical corrections can reduce the approximate relative error to 1.0. The level of error is within what would be expected given uncertainties in threshold shear velocity and wind speed at our sites. The model outperforms the alternative schemes both in terms of approximate relative error and the number of sites at which threshold shear velocity was exceeded. These results lend support to an understanding of the physics of aeolian transport in which (1) vegetation's impact on transport is dependent upon the distribution of vegetation rather than merely its average lateral cover and (2) vegetation impacts surface shear stress locally by depressing it in the immediate lee of plants rather than by changing the bulk surface's threshold shear velocity. Our results also suggest that threshold shear velocity is exceeded more than might be estimated by single measurements of threshold shear stress and roughness length commonly associated with vegetated surfaces, highlighting the variation of threshold shear velocity with space and time in real landscapes.

Journal of Geophysical Research F: Earth Surface↗
USGS2015

Role of ground ice dynamics and ecological feedbacks in recent ice wedge degradation and stabilization

Ground ice is abundant in the upper permafrost throughout the Arctic and fundamentally affects terrain responses to climate warming. Ice wedges, which form near the surface and are the dominant type of massive ice in the Arctic, are particularly vulnerable to warming. Yet processes controlling ice wedge degradation and stabilization are poorly understood. Here we quantified ice wedge volume and degradation rates, compared ground ice characteristics and thermal regimes across a sequence of five degradation and stabilization stages and evaluated biophysical feedbacks controlling permafrost stability near Prudhoe Bay, Alaska. Mean ice wedge volume in the top 3 m of permafrost was 21%. Imagery from 1949 to 2012 showed thermokarst extent (area of water-filled troughs) was relatively small from 1949 (0.9%) to 1988 (1.5%), abruptly increased by 2004 (6.3%) and increased slightly by 2012 (7.5%). Mean annual surface temperatures varied by 4.9°C among degradation and stabilization stages and by 9.9°C from polygon center to deep lake bottom. Mean thicknesses of the active layer, ice-poor transient layer, ice-rich intermediate layer, thermokarst cave ice, and wedge ice varied substantially among stages. In early stages, thaw settlement caused water to impound in thermokarst troughs, creating positive feedbacks that increased net radiation, soil heat flux, and soil temperatures. Plant growth and organic matter accumulation in the degraded troughs provided negative feedbacks that allowed ground ice to aggrade and heave the surface, thus reducing surface water depth and soil temperatures in later stages. The ground ice dynamics and ecological feedbacks greatly complicate efforts to assess permafrost responses to climate change.

Journal of Geophysical Research F: Earth Surface↗