Search USGS⌕ Search

Geology topics

Skye C. Corbett

Publications and source records attributed to Skye C. Corbett.

17 recordsLinked to original sources

Widespread landslide activity in an extreme wet season and implications for regional sediment management, eastern San Francisco Bay area, California

Watershed sediment production is expected to increase in a warmer future with more extreme rain, with cascading effects throughout drainage and sediment-transport networks. This study investigated landscape-scale sediment movement in the eastern San Francisco Bay area, California, USA, during the extreme 2016–2017 wet season that brought major rainfall, landslides, and flooding. Mapping 8,928 landslides across a 1,050-km 2 study area revealed new sediment yield of 510–956 t/km 2 , equivalent to denudation of 193–361 mm/ky. These results correspond closely to long-term denudation rates in the northern and central California Coast Ranges, indicating that mass wasting in very wet years dominates long-term sediment mobilization. However, due to long residence times in drainage networks, the 2017 landslides contributed at most ∼1%–2% of the estimated locally derived fluvial sediment transport to San Francisco Bay. Although the amount of sediment mobilized did not threaten municipal water supplies, small rangeland impoundments in this mixed-use landscape lost storage capacity to new sedimentation. Considering regional sediment supply and demand, even the exceptionally large sediment delivery in an extreme wet year cannot meet the need for sediment to accrete tidal wetlands in the bay. To keep pace with rising sea levels, this abnormally high terrestrial sediment input would need to occur in 50 of the next 75 years, an unlikely occurrence due to the prevalence of recent drought years. Shoreline protection and restoration in the bay would need additional sources of sediment, such as through management of dredged sediment through beneficial-reuse programs.

California↗

Assessment of extreme subsurface hydrologic conditions captured during atmospheric river storms in the San Francisco Bay area (California, USA) with applications to shallow landslide early warning

An increase in soil pore water pressure is the typical trigger for the majority of landslides caused by rainfall. Atmospheric river storms, common to the west coast of North America during the winter season, can deliver landslide triggering rainfall resulting in severe impacts to coastal communities. Using a network of hydrologic monitoring stations situated within landslide-prone terrain in the San Francisco Bay area of California (USA), we assess the meteorologic conditions and resulting hydrologic and landslide response resulting from eight consecutive storm events that caused thousands of shallow landslides during the winter of 2022–2023. We find disparate hydrological responses and resultant degrees of observed landsliding ranging from < 1 landslide/km2 to 18 landslides/km2 at the monitoring sites that reflect the interplay and differences between rainfall delivery, subsurface hydrological characteristics, and geotechnical properties at each site. Antecedent soil moisture from both early season rainfall and the first storm in the sequence played a critical role in setting up some hillslopes for failure. Subsequent storms then generated elevated pore water pressures for several hours with associated landsliding. However, we find that the occurrence of widespread landsliding required not only sufficient pore water pressure magnitude in susceptible hillslopes, but also full and prolonged effective soil saturation throughout hillslope profiles. Landslides may still occur at lower values and durations of effective saturation but are likely to be less extensive regionally. We present these findings within the context of research directions and improvements to landslide early warning systems first suggested by researchers 40 years ago.

California↗

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗

An open-source workflow for scaling burn severity metrics from drone to satellite to support post-fire watershed management

Wildfires are increasing in size and severity across much of the western United States, exposing vulnerable wildland-urban interfaces to post-fire hazards. The Mediterranean chaparral region of Northern California contains many high sloping watersheds prone to hazardous post-fire flood events and identifying watersheds at high risk of soil loss and debris flows is a priority for post-fire response and management. Uncrewed Aerial Systems (UAS; aka drones) offer post-fire management teams the ability to quickly mobilize and survey burned areas with very high-resolution imagery (∼1 cm), facilitating emergency management and post-fire hazard assessment. However, adoption of this technology by hazard response teams may be hindered by complicated workflows for UAS data acquisition, image processing and analysis. We present an open-source workflow using mature Geographic Information Systems (GIS) software and Python packages in a Jupyter Notebook environment that guides users through classification of true-color UAS imagery to generate high resolution burn severity maps which can then be scaled across larger watersheds using Sentinel-2 normalized burn ratio (NBR) images. Soil burn severity classifications using a weighted brightness (WB) image and Char Index (CI) generated from UAS imagery were validated with in-situ data and random stratified points, resulting in the CI having the highest overall accuracy of 87.5%. CI also displayed a marginally stronger relationship over the WB with the post-fire Sentinel-2 NBR, R 2 = 0.79 and R 2 = 0.78 respectively. Our methods offer the unique opportunity to standardize GIS workflows, promoting replication through transparency, while improving the user's understanding of scientific GIS functionality.

California↗

Investigating the atmospheric conditions associated with impactful shallow landslides in California (USA)

Shallow landslides are often triggered during rainfall events, which can increase subsurface soil water pressure and destabilize hillslopes. The likelihood of regional shallow landslide initiation is often assessed through a comparison of rainfall intensity and duration to pre-established thresholds. While informative for landslide warning, this exclusive focus on rainfall exceeding thresholds does not consider the meteorological conditions producing the rainfall. Here, we ask the question, are there common meteorological characteristics that lead to landslide-triggering precipitation? We develop a catalog of 18 post-1995 widespread, impactful shallow landslide events occurring within 13 storms across California, USA, where initiation time could be constrained to a ≤6-h window. We examine storm characteristics during the landslide initiation window using atmospheric reanalysis products, radar observations, and quantitative precipitation estimates. We find that, while there are some common atmospheric characteristics across landslide events, they can occur under a range of atmospheric conditions. For example, all Northern California landslide events assessed are associated with moderate to strong atmospheric rivers (ARs), while Southern California landslides feature non-AR to strong AR conditions. The storm events evaluated herein share many characteristics of hydrologically important storms in California that did not necessarily result in landslides; thus, atmospheric characteristics alone may not be sufficient to determine whether landslides will occur. However, documenting the characteristics of landslide-triggering storms defines the conditions under which landslides tend to occur, provides analog events that can be useful in forecast applications, helps define future research directions relating to atmospheric conditions and landslides, and supports interdisciplinary research efforts.

California↗

Landslides triggered by the 2016–2017 storm season, eastern San Francisco Bay region, California

Discussion The winter rainy season of 2016–2017 brought abundant rainfall to the State of California and to the San Francisco Bay region. In January and February of 2017, intense rainfall from strong winter storms saturated soils in the region and triggered thousands of shallow landslides. The highest concentration of these landslides was in the eastern part of the bay region, where landslides in the hills east of the Cities of Richmond, Berkeley, Oakland, Hayward, and Fremont (see main map for locations) damaged homes, displaced a major electrical transmission-line tower, and blocked several heavily traveled roadways. This map shows 8,928 landslides manually mapped from rectified high-resolution (0.25-meter [m]) satellite imagery (March 11, 2017, from Google Earth) in a three-dimensional geographic information system (GIS) framework. The map area encompasses approximately 1,050 square kilometers (km 2 ), bounded by the Carquinez Strait and San Francisco Bay to the north and west, respectively, and extending to the Interstate Highway 680 corridor to the south and east. Individual landslides were mapped as polygons, but for ease of display, they are shown here as points denoting the highest elevation of each landslide headscarp. The greatest calculated landslide concentration (measured as the total number of landslides per unit area) exceeded 80 landslides per 0.25 km 2 in the hills east of the City of Berkeley. This zone is illustrated on the map in red, a color which denotes areas of greater than or equal to 30 landslides per 0.25 km 2 . The source area and deposit polygon datasets are available as a data release (Corbett and Collins, 2023). Complementary field investigations at more than 150 landslides of those displayed here indicate that most mapped landslides are shallow (less than 1 m deep) debris slides and debris flows. Although most landslides occurred in undeveloped areas (for example, parks and open space), about 1,400 landslides, or 16 percent of the total number, affected structures or infrastructure (including paved and dirt roads) in some way. The inset figures (figs. 1–5) show photographs taken in January and February 2017 that depict some areas of damage that occurred during and after the storms; the locations of where these photographs were taken are shown on the map. In the San Francisco Bay region, landslides usually take place each year during the winter rainy season, but widespread, intense rainfall events such as those whose effects are depicted here typically occur as a result of considerably above-average precipitation conditions. Although the landslides of early 2017 are not the most damaging to affect the region historically, they are still a potent reminder of the potential for landslide hazards present in the San Francisco Bay region.

California↗

Multi-stage soil-hydraulic recovery and limited ravel accumulations following the 2017 Nuns and Tubbs wildfires in Northern California

Wildfire can impact soil-hydraulic properties by reducing saturated hydraulic conductivity and sorptivity, making recently burned landscapes prone to debris flows and flash floods. The post-fire hazard window can range from years to decades. In Northern California, where wildfire frequency is steadily increasing, the impact and soil-hydraulic recovery from wildfires is unknown. Following the October 2017 Nuns and Tubbs fires in the Northern Bay Area of California, we established 41 monitoring sites for repeat tension-disc infiltrometer measurements of field-saturated hydraulic conductivity ( K fs ) over 3.5 years. Our site arrays, which encompass grasslands, chaparral, and oak and conifer forests across a range in lithology, show a marked decrease in K fs following the wildfires and a swift partial recovery following the initial post-fire rainy season. Our time series reveals a complex path to soil-hydraulic recovery marked by distinct seasonal stages. Analysis of changing K fs , sorptivity, and infiltration model residuals collectively suggests that these stages are related to transitions between soil-hydraulic processes like structural soil sealing from rainsplash, thermal cracking of bare soil, and vegetation regrowth. While soil infiltration rates were strongly impacted by the 2017 fires, dry ravel estimates are an order of magnitude less for similar slopes than the 2009 Station fire in the San Gabriel mountains of Southern California, suggesting that limited ravel flux may insufficiently load channels for debris flows that initiate from within-channel failure. Our analysis suggests that burned landscapes in the Northern Bay Area of California may experience rapid soil-hydraulic recovery and limited pathways toward post-fire debris flow initiation.

California↗

Rockfall kinematics from massive rock cliffs: Outlier boulders and flyrock from Whitney Portal, California, rockfalls

Geologic conditions and topographic setting are among the most critical factors for assessing rockfall hazards. However, other subtle features of rockfall motion may also govern the runout of rockfall debris, particularly for those sourced from massive cliffs where debris can have substantial momentum during transport. Rocks may undergo collisions with trees and talus boulders, with the latter potentially generating flyrock—launched rock pieces resulting from boulder collisions that follow distinctively different paths than the majority of debris. Collectively, these intricacies of rockfall kinematics may substantially govern the hazards expected from rockfall to both persons and infrastructure located beneath steep cliffs. Here, we investigate the kinematics, including outlier boulder and flyrock trajectories, of seismically triggered rockfalls on 24 June 2020 that damaged campground facilities near Whitney Portal, CA, a heavily used outdoor recreation gateway to the Sierra Nevada mountains. Our results, obtained in part by rockfall runout model simulations, indicate that outlier boulder trajectories resulted from opportunities provided by less steep terrain beyond the talus edge. The influence of trees, initially thought to have served a protective capacity in attenuating rockfall energy, appears to have been negligible for the large boulder volumes (>50 m 3 ) mobilized, although they did potentially deflect the trajectory of flyrock debris. Rockfall outlier boulders from the event were comparable in volume and runout distance to prehistoric boulders located beyond the talus slope, thereby providing some level of confidence in the use of a single rockfall shadow angle for estimating future rockfall hazards at the site.

California↗

Linking mesoscale meteorology with extreme landscape response: Effects of narrow cold frontal rainbands (NCFR)

Landscapes evolve in response to prolonged and/or intense precipitation resulting from atmospheric processes at various spatial and temporal scales. Whereas synoptic (large‐scale) features (e.g., atmospheric rivers and hurricanes) govern regional‐scale hydrologic hazards such as widespread flooding, mesoscale features such as thunderstorms or squall lines are more likely to trigger localized geomorphic hazards such as landslides. Thus, to better understand relations between hydrometeorological drivers and landscape response, a knowledge of mesoscale meteorology and its impacts is needed. Here we investigate the extreme geomorphic response associated with one type of mesoscale meteorological feature, the narrow cold frontal rainband (NCFR). Resulting from low‐level convergence and shallow convection along a cold front, NCFRs are narrow bands of high‐intensity rainfall that occur in midlatitude areas of the world. Our study examines an NCFR impacting the Sierra Nevada foothills (California, USA) that initiated over 500 landslides, mobilized ~360,000 metric tons of sediment to the fluvial system (as much as 16 times the local annual sediment yield), and severely damaged local infrastructure and regional water transport facilities. Coupling geomorphological field investigations with meteorological analyses, we demonstrate that precipitation associated with the NCFR was both intense (maximum 15 min intensity of 70 mm/hr) and localized, resulting in a highly concentrated band of shallow landsliding. This meteorological phenomenon likely plays an important role in landscape evolution and hazard initiation. Other types of mesoscale meteorological features also occur globally and offer new avenues for understanding the effects of storms on landscapes.

California↗

Terrestrial lidar data of the February 14, 2019 Sausalito Boulevard Landslide, Sausalito, California

On February 14, 2019, just before 2:56 am local time (Pacific Standard Time), a landslide initiated from the hillslopes above the Hurricane Gulch section of the City of Sausalito, Marin County, California. The landslide, specifically classified as a debris flow, overran a road (Sausalito Boulevard) immediately below the landslide source area and impacted a residential structure that subsequently toppled downslope and collided with another residential structure. The second structure then crossed a lower road (Crescent Avenue) that runs along the base of the slope before the mixture of soil and structural debris came to rest in and near the valley axis that drains the lower area of Hurricane Gulch. The U.S. Geological Survey responded to this event within hours of the landslide and provided situational awareness of possible secondary landslide hazards associated with the unstable slope. The USGS also rapidly mobilized its topographic surveying capabilities (specifically, GPS and terrestrial lidar devices) and collected a three-dimensional point cloud model of the landslide source area and surrounding terrain to capture the as-failed condition of the slope for use in potential future studies. This report summarizes the methods and available data collected during this response.

California↗

Thermal influences on spontaneous rock dome exfoliation

Rock domes, with their onion-skin layers of exfoliation sheets, are among the most captivating landforms on Earth. Long recognized as integral in shaping domes, the exact mechanism(s) by which exfoliation occurs remains enigmatic, mainly due to the lack of direct observations of natural events. In August 2014, during the hottest days of summer, a granitic dome in California, USA, spontaneously exfoliated; witnesses observed extensive cracking, including a ~8000 kg sheet popping into the air. Subsequent exfoliation episodes during the following two summers were recorded by instrumentation that captured—for the first time—exfoliation deformation and stress conditions. Here we show that thermal cycling and cumulative dome surface heating can induce subcritical cracking that culminates in seemingly spontaneous exfoliation. Our results indicate that thermal stresses—largely discounted in dome formation literature—can play a key role in triggering exfoliation and therefore may be an important control for shaping domes worldwide.

Nature Communications↗

Seismic velocity site characterization of 10 Arizona strong-motion recording stations by spectral analysis of surface wave dispersion

Vertical one-dimensional shear wave velocity ( V S ) profiles are presented for strong-motion sites in Arizona for a suite of stations surrounding the Palo Verde Nuclear Generating Station. The purpose of the study is to determine the detailed site velocity profile, the average velocity in the upper 30 meters of the profile ( V S30 ), the average velocity for the entire profile ( V SZ ), and the National Earthquake Hazards Reduction Program (NEHRP) site classification. The V S profiles are estimated using a non-invasive continuous-sine-wave method for gathering the dispersion characteristics of surface waves. Shear wave velocity profiles were inverted from the averaged dispersion curves using three independent methods for comparison, and the root-mean-square combined coefficient of variation ( COV ) of the dispersion and inversion calculations are estimated for each site.

Arizona↗

Conditions and processes affecting sand resources at archeological sites in the Colorado River corridor below Glen Canyon Dam, Arizona

This study examined links among fluvial, aeolian, and hillslope geomorphic processes that affect archeological sites and surrounding landscapes in the Colorado River corridor downstream from Glen Canyon Dam, Arizona. We assessed the potential for Colorado River sediment to enhance the preservation of river-corridor archeological resources through aeolian sand deposition or mitigation of gully erosion. By identifying locally prevailing wind directions, locations of modern sandbars, and likely aeolian-transport barriers, we determined that relatively few archeological sites are now ideally situated to receive aeolian sand supply from sandbars deposited by recent controlled floods. Whereas three-fourths of the 358 river-corridor archeological sites we examined include Colorado River sediment as an integral component of their geomorphic context, only 32 sites currently appear to have a high degree of connectivity (coupled interactions) between modern fluvial sandbars and sand-dominated landscapes downwind. This represents a substantial decrease from past decades, as determined by aerial-photograph analysis. Thus, we infer that recent controlled floods have had a limited, and declining, influence on archeological-site preservation. Within the study area, overland-flow (gully) erosion is less severe in sand landscapes with active aeolian sand than in landscapes that lack aeolian transport; gullies terminate more commonly in active sand (sand that is mobile by wind rather than stabilized by biologic soil crust). We infer that these characteristics largely result from aeolian sand transport being an effective gully-limiting and gully-annealing mechanism. Aeolian sand activity in the river corridor varies substantially as a function of reach morphology and dominant wind direction relative to the river-corridor orientation, factors that control accommodation space for river-derived sand and the modern sand supply to aeolian dunes. These attributes, together with an inverse correlation between aeolian sand activity and gully occurrence, define varying degrees of net long-term gully-erosion risk for sediment deposits and associated archeological sites in different regions of the river corridor. Over most of the river corridor, including some of the archeologically richest regions, sand is too inactive with respect to aeolian transport to anneal gullies effectively. At eight selected archeological sites that we studied with high-resolution terrestrial lidar scans for more than a year, sand loss by overland flow (gully erosion) and aeolian deflation generally exceeded deposition, such that erosion dominated over most monitoring intervals&mdash;even at four sites with strong connectivity to modern sand supply. The Glen Canyon reach of the river corridor appears especially vulnerable to gully erosion. Among the sites that we monitored in detail, erosion generally dominated over deposition to a greater degree at four Glen Canyon sites with no modern sand supply than at four Marble&ndash;Grand Canyon sites with aeolian sand supply from controlled-flood sandbars. Although gross annual-scale erosion rates were similar among the Glen Canyon sites and among the Marble&ndash;Grand Canyon sites, a relative lack of depositional processes led to greater net erosion at the Glen Canyon sites. Having found no differences in weather patterns to suggest greater erosive forcing in Glen Canyon, and no conclusively influential differences in the slope or watershed area contributing to gully formation, we attribute the greater erosion at the Glen Canyon sites to a combination of inherent geomorphic context (high terraces that do not receive modern sediment supply) and pronounced effects of postdam sediment-supply limitation. We conclude that most of the river-corridor archeological sites are at elevated risk of net erosion under present dam operations. In the present flow regime, controlled floods do not simulate the magnitude or frequency of natural floods, and are not large enough to deposit sand at elevations that were flooded at annual to decadal intervals in predam time. For archeological sites that depend upon river-derived sand, we infer elevated erosion risk owing to a combination of reduced sand supply (both fluvial and aeolian) through (1) the lower-than-natural flood magnitude, frequency, and sediment supply of the controlled-flooding protocol; (2) reduction of open, dry sand area available for wind redistribution under current normal (nonflood) dam operations, which do not include flows as low as natural seasonal low flows and do include substantial daily flow fluctuations; and (3) impeded aeolian sand entrainment and transport owing to increased riparian vegetation growth in the absence of larger, more-frequent floods. If dam operations were to increase the supply of sand available for windblown transport&mdash;for example, through larger floods, sediment augmentation, or increased fluvial sandbar exposure by low flows&mdash;and also decrease riparian vegetation, the prevalence of active aeolian sand could increase over time, and the propensity for unmitigated gully erosion could decrease. Although the evolution of river-corridor landscapes and archeological sites has been altered fundamentally by the lack of large, sediment-rich floods (flows on the order of 5,000 m 3 /s), some combination of sediment-rich flows above 1,270 m 3 /s, seasonal flows below 226 m 3 /s, and riparian-vegetation removal might increase the preservation potential for sand-dependent archeological resources in the Colorado River corridor.

Arizona↗

Relations between rainfall–runoff-induced erosion and aeolian deposition at archaeological sites in a semi-arid dam-controlled river corridor

Process dynamics in fluvial-based dryland environments are highly complex with fluvial, aeolian, and alluvial processes all contributing to landscape change. When anthropogenic activities such as dam-building affect fluvial processes, the complexity in local response can be further increased by flood- and sediment-limiting flows. Understanding these complexities is key to predicting landscape behavior in drylands and has important scientific and management implications, including for studies related to paleoclimatology, landscape ecology evolution, and archaeological site context and preservation. Here we use multi-temporal LiDAR surveys, local weather data, and geomorphological observations to identify trends in site change throughout the 446-km-long semi-arid Colorado River corridor in Grand Canyon, Arizona, USA, where archaeological site degradation related to the effects of upstream dam operation is a concern. Using several site case studies, we show the range of landscape responses that might be expected from concomitant occurrence of dam-controlled fluvial sand bar deposition, aeolian sand transport, and rainfall-induced erosion. Empirical rainfall-erosion threshold analyses coupled with a numerical rainfall&ndash;runoff&ndash;soil erosion model indicate that infiltration-excess overland flow and gullying govern large-scale (centimeter- to decimeter-scale) landscape changes, but that aeolian deposition can in some cases mitigate gully erosion. Whereas threshold analyses identify the normalized rainfall intensity (defined as the ratio of rainfall intensity to hydraulic conductivity) as the primary factor governing hydrologic-driven erosion, assessment of false positives and false negatives in the dataset highlight topographic slope as the next most important parameter governing site response. Analysis of 4+ years of high resolution (four-minute) weather data and 75+ years of low resolution (daily) climate records indicates that dryland erosion is dependent on short-term, storm-driven rainfall intensity rather than cumulative rainfall, and that erosion can occur outside of wet seasons and even wet years. These results can apply to other similar semi-arid landscapes where process complexity may not be fully understood.

Arizona↗

Shear Wave Velocity and Site Amplification Factors for 25 Strong-Motion Instrument Stations Affected by the M5.8 Mineral, Virginia, Earthquake of August 23, 2011

Vertical one-dimensional shear wave velocity (Vs) profiles are presented for 25 strong-motion instrument sites along the Mid-Atlantic eastern seaboard, Piedmont region, and Appalachian region, which surround the epicenter of the M5.8 Mineral, Virginia, Earthquake of August 23, 2011. Testing was performed at sites in Pennsylvania, Maryland, West Virginia, Virginia, the District of Columbia, North Carolina, and Tennessee. The purpose of the study is to determine the detailed site velocity profile, the average velocity in the upper 30 meters of the profile (VS,30), the average velocity for the entire profile (VS,Z), and the National Earthquake Hazards Reduction Program (NEHRP) site classification. The Vs profiles are estimated using a non-invasive continuous-sine-wave method for gathering the dispersion characteristics of surface waves. A large trailer-mounted active source was used to shake the ground during the testing and produce the surface waves. Shear wave velocity profiles were inverted from the averaged dispersion curves using three independent methods for comparison, and the root-mean square combined coefficient of variation (COV) of the dispersion and inversion calculations are estimated for each site.

North Carolina, Pennslyvannia, Tennesee, West Virg↗

High-resolution topography and geomorphology of select archeological sites in Glen Canyon National Recreation Area, Arizona

Along the Colorado River corridor between Glen Canyon Dam and Lees Ferry, Arizona, located some 25 km downstream from the dam, archaeological sites dating from 8,000 years before present through the modern era are located within and on top of fluvial and alluvial terraces of the prehistorically undammed river. These terraces are known to have undergone significant erosion and retreat since emplacement of Glen Canyon Dam in 1963. Land managers and policy makers associated with managing the flow of the Colorado River are interested in understanding how the operations of Glen Canyon Dam have affected the archeological sites associated with these terraces and how dam-controlled flows currently interact with other landscape-shaping processes. In 2012, the U.S. Geological Survey initiated a research project in Glen Canyon to study the types and causes of erosion of the terraces. This report provides the first step towards this understanding by presenting comparative analyses on several types of high-resolution topographic data (airborne lidar, terrestrial lidar, and airborne photogrammetry) that can be used in the future to document and analyze changes to terrace-based archaeological sites. Herein, we present topographic and geomorphologic data of four archaeological sites within a 14 km segment of Glen Canyon using each of the three data sources. In addition to comparing each method’s suitability for adequately representing the topography of the sites, we also analyze the data within each site’s context and describe the geomorphological processes responsible for erosion. Our results show that each method has its own strengths and weaknesses, and that terrestrial and airborne lidar are essentially interchangeable for many important topographic characterization and monitoring purposes. However, whereas terrestrial lidar provides enhanced capacity for feature recognition and gully morphology delineation, airborne methods (whether by way of laser or optical sensors) are better suited for reach- and regional-scale mapping. Our site-specific geomorphic analyses of the four archeological sites indicate that their current topographical conditions are a result of different and sometimes competing erosional agents, including bedrock- and terrace-based overland flow, fluvial-induced terrace bank collapse, and alluvial-fan-generated debris flows. Although the influences of anthropogenic-induced erosion from dam operations are not specifically analyzed in this report, we do identify geomorphic settings where dam operations are either more or less likely to affect archeological site stability. This information can be used to assist with future monitoring efforts of these sites and identification of similar conditions for other archeological sites along the Colorado River corridor in Glen Canyon.

Arizona↗

Topographic change detection at select archeological sites in Grand Canyon National Park, Arizona, 2007-2010

Human occupation in Grand Canyon, Arizona, dates from at least 11,000 years before present to the modern era. For most of this period, the only evidence of human occupation in this iconic landscape is provided by archeological sites. Because of the dynamic nature of this environment, many archeological sites are subject to relatively rapid topographic change. Quantifying the extent, magnitude, and cause of such change is important for monitoring and managing these archeological sites. Such quantification is necessary to help inform the continuing debate on whether and how controlled releases from Glen Canyon Dam, located immediately upstream of Grand Canyon National Park, are affecting site erosion rates, artifact transport, and archeological resource preservation along the Colorado River in Grand Canyon. Although long-term topographic change resulting from a variety of natural processes is inherent in the Grand Canyon region, continued erosion of archeological sites threatens both the archeological resources and our future ability to study evidence of past cultural habitation. Thus, this subject is of considerable interest to National Park Service managers and other stakeholders in the Glen Canyon Dam Adaptive Management Program. Understanding the causes and effects of archeological site erosion requires a knowledge of several factors, including the location, timing, and magnitude of the changes occurring in relation to archeological resources, the rates of change, and the relative contribution of potential causes. These potential causes include sediment depletion associated with managed flows from Glen Canyon Dam, site-specific weather and overland flow patterns, visitor impacts, and long-term regional climate change. To obtain this information, highly accurate, spatially specific data are needed from sites undergoing change. Using terrestrial lidar techniques, and building upon three previous surveys of archeological sites performed in 2006 and 2007, we collected two new datasets in April and September 2010 and processed and improved upon existing methods to generate high-accuracy (3 to 5 cm vertical change threshold) topographic change-detection maps for 10 survey areas encompassing 9 archeological sites along the Colorado River corridor. We also used terrestrial lidar techniques to investigate several other metrics for studying archeological site stability, including monitoring cultural structures and artifacts and remotely measuring cryptobiotic soil crust areas. Our topographic change results indicate that 9 of 10 survey areas showed signs of either erosion, deposition, or both during the 2007&ndash;2010 time interval and that these changes can be linked to a variety of geomorphic processes, primarily overland flow gullying and aeolian sand transport. In several cases, large (>50 cm) vertical change occurred, and in one case, more than 100 m 3 of sediment was eroded. Further, for all sites monitored throughout the river corridor during this time period, the overall signal was related to erosion rather than deposition. These results highlight the potential for rapid archeological site change in Grand Canyon. Whereas the topographic change results presented herein provide the highest level of change detection yet performed on entire archeological sites in Grand Canyon, additional work in combining these results with site-specific weather, hydrology, and geomorphology data is needed to provide a more thorough understanding of the causes of the documented topographic changes. Linking lidar-derived measurements of topographic changes with these other data sources should provide land managers with a scientific basis for making management decisions regarding archeological resources in Grand Canyon National Park and assist in answering open questions regarding the influence that sediment-depleted flows from Glen Canyon Dam have on archeological site stability.

Arizona↗