Search USGSSearch

SEARCH · Search USGS

Results for “Photogrammetric Record”

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.

34 records · Page 2Linked to original sources

Computation with physical values from Landsat digital data

Landsat digital images are commonly analyzed by using the digital numbers for each pixel recorded on a computer-compatible magnetic tape. Although this procedure may be satisfactory when only a single, internally consistent image is used, the procedure may produce incorrect results if more than one image is used for analysis as in mosaics or temporal overlays. The digital numbers for each pixel should be converted to their dimensioned equivalents such as radiance, as measured at the satellite, in milliwatts per square centimetre per steradian, or reflectance.-from Author

Photogrammetric Engineering and Remote Sensing

The relationship of acquisition systems to automated stereo correlation.

Today a concerted effort is being made to expedite the mapping process through automated correlation of stereo data. Stereo correlation involves the comparison of radiance (brightness) signals or patterns recorded by sensors. Conventionally, two-dimensional area correlation is utilized but this is a rather slow and cumbersome procedure. Digital correlation can be performed in only one dimension where suitable signal patterns exist, and the one-dimensional mode is much faster. Electro-optical (EO) systems, suitable for space use, also have much greater flexibility than film systems. Thus, an EO space system can be designed which will optimize one-dimensional stereo correlation and lead toward the automation of topographic mapping.-from Author

Photogrammetric Engineering and Remote Sensing

Twenty years of explosive-effusive activity at El Reventador volcano (Ecuador) recorded in its geomorphology

Shifts in activity at long-active, open-vent volcanoes are difficult to forecast because precursory signals are enigmatic and can be lost in and amongst daily activity. Here, we propose that crater and vent morphologies, along with summit height, can help us bring some insights into future activity at one of Ecuador’s most active volcanoes El Reventador. On 3 November 2002, El Reventador volcano experienced the largest eruption in Ecuador in the last 140 years and has been continuously active ever since with transitions between and coexistence of explosive and effusive activity, characterized by Strombolian and Vulcanian behavior. Based on the analysis of a large dataset of thermal and visual images, we determined that in the last 20 years of activity, the volcano faced three destructive events: A. Destruction of the upper part of the summit leaving a north-south breached crater (3 November 2002), B. NE border crater collapse (2017), and C. NW flank collapse (2018), with two periods of reconstruction of the edifice: Period 1. Refill of the crater (2002-early 2018) and Period 2. Refill of the 2018 scar (April 2018–December 2022). Through photogrammetric analysis of visual and thermal images acquired in 11 overflights of the volcano, we created a time-series of digital elevation models (DEMs) to determine the maximum height of the volcano at each date, quantify the volume changes between successive dates, and characterize the morphological changes in the summit region. We estimate that approximately 34.1x10 6 m 3 of volcanic material was removed from the volcano due to destructive events, whereas 64.1x10 6 m 3 was added by constructive processes. The pre-2002 summit height was 3,560 m and due to the 2002 eruption it decreased to 3,527 m; it regained its previous height between 2014 and 2015 and the summit crater was completely filled by early April 2018. Event A resulted from an intrusion of magma that erupted violently; we proposed that Events B and C could be a result of an intrusion as well but may also be due to a lack of stability of the volcano summit which occurs when it reaches its maximum height of approximately 3,590 and 3,600 m.

El Reventador

Resource understanding: a challenge to aerial methods

Aerial survey methods are speeding acquisition of survey data needed to provide and manage the nation's resources. These methods have been applied to topographic mapping for a number of years and the record clearly shows their advantages in terms of cost and speed in contrast to the ground methods that have been historically employed. Limited use is now being made of aerial methods to assist cadastral surveys, in location, acquisition and development of National Parks, in mapping the geology of the nation, in locating and developing water resources, and in surveys of the oceans. It is the purpose of this paper to call attention to these uses and to encourage the scientific community to further refine aerial methods so that their use may be increased and the veracity of data improved.

Photogrammetric Engineering

Geologic map of the Cochiti Dam quadrangle, Sandoval County, New Mexico

The Cochiti Dam quadrangle is located in the southern part of the Española Basin and contains sedimentary and volcanic deposits that record alluvial, colluvial, eolian, tectonic and volcanic processes over the past seventeen million years. The geology was mapped from 1997 to 1999 and modified in 2004 to 2008. The primary mapping responsibilities were as follows: Dethier mapped the surficial deposits, basin-fill sedimentary deposits, Miocene to Quaternary volcanic deposits of the Jemez volcanic field, and a preliminary version of fault distribution. Thompson and Hudson mapped the Pliocene and Quaternary volcanic deposits of the Cerros del Rio volcanic field. Thompson, Minor, and Hudson mapped surface exposures of faults and Hudson conducted paleomagnetic studies for stratigraphic correlations. Thompson prepared the digital compilation of the geologic map. The mapped distribution of units is based primarily on interpretation of 1:16,000-scale, color aerial photographs taken in 1992, and 1:40,000-scale, black-and-white, aerial photographs taken in 1996. Most of the contacts on the map were transferred from the aerial photographs using a photogrammetric stereo-plotter and subsequently field checked for accuracy and revised based on field determination of allostratigraphic and lithostratigraphic units. Determination of lithostratigraphic units in volcanic deposits was aided by geochemical data, 40 Ar/ 39 Ar geochronology, aeromagnetic and paleomagnetic data. Supplemental revision of mapped contacts was based on interpretation of USGS 1-meter orthoimagery.

New Mexico

The relationship of channel planform and point bar architecture on a reach of the Wabash River near Grayville, Illinois

The erosional and depositional characteristics of meandering rivers lead to the formation and maintenance of point bars along the inner banks of meander bends. Point bars are composed of sediment layers in patterns resulting from the rate and style of channel migration, hydrodynamics, and sediment transport and deposition within the river system (e.g. Jackson, 1976; Dietrich and Smith, 1984; Dietrich, 1987, Abad and Garcia, 2009). The distribution of the sediments preserved in the internal architecture of a river point bar provides a record of channel planform evolution. Geophysical methods are used to gain a large-scale visualization of the subsurface and aid in the interpretation of historic channel patterns (Best et al., 2003; Woodward et al., 2003; Sambrook Smith et al., 2006). Comparing known surficial extents of the point bar to features identified in the subsurface can also enhance the understanding of historic channel planform. This study investigates two point bars along bends with different styles of migration, Maier and TB3, in a well-documented reach of the Wabash River near Grayville, IL. Evidence from historic aerial photography, modern lidar, photogrammetry, and geophysical surveys were used to determine the relationship between the point bar architecture and channel planform. Airborne lidar was flown in 2011 and is used to create the 2011 point bar surface. In 2017, a terrestrial lidar and Real Time Kinematic-Global Navigation Satellite System (RTK-GNSS) topographic survey were combined to create the surface for TB3. In 2018, a photogrammetric survey collected with a small unmanned aerial system (sUAS) was used to create a structure-frommotion (SfM) derived surface for Maier bend. The 2018 survey-based point bar surfaces were differenced from the 2011 point bar surfaces to get a Digital Elevation Model (DEM) of difference (DoD) to visualize areas of erosion and deposition. In addition, geophysical surveys using ground penetrating radar (GPR) were conducted in transverse and streamwise lines across the point bars in 2018. Elevation profiles from the 2011 point bar surfaces are extracted and overlain onto the 2018 GPR images to determine how the point bar is preserving the structure of sediments previously deposited. Results from this study provide an update to current models of point bar architecture.

Illinois

Characterizing precariously balanced rocks (PBRs) in the eastern United States for estimating maximum past earthquake ground motions

Fragile geologic features (FGFs) can be used to estimate an upper bound to the strength of earthquake ground shaking in the past thousands to tens of thousands of years. FGFs could be especially useful in the eastern United States, where few active faults are available for paleoseismic studies and where low seismic attenuation means that an earthquake could damage FGFs at large distances. This article describes a pilot study of one form of FGF, precariously balanced rocks (PBRs), in the Blue Ridge of Virginia, ∼160 km from the central Virginia seismic zone (CVSZ) and ∼100 km from Virginia’s Giles County seismic zone (GCSZ). The main goal is to develop and refine methods for studying PBRs in the eastern United States. Photogrammetric 3D models of the rocks were used to locate their centers of mass and delineate their basal contacts, with the latter corroborated by wrapping a wire around the bases to verify the circumferences. Relations between hill geometry and amplification were used to estimate topographic amplification. Seismometers on the PBRs provided the frequencies, amplitudes, and directions of rocking after giving the rocks gentle pushes. The levels of ground shaking likely to topple the rocks were compared with the median hazard curves from the U.S. Geological Survey’s National Seismic Hazard Model and were used to map magnitudes of earthquakes in the nearby region likely to topple the rocks. Results indicate that an M w ≥ 7.5 earthquake is needed in the CVSZ or near the 2020 Sparta, North Carolina, earthquake to topple the two most fragile PBRs and that an earthquake of ~ M w 7–7.5 in the GCSZ is needed. The fragilities and estimated ages of the PBRs (20–40 ka) are thus consistent with the absence of such events in the historical record in the nearby region.

eastern United States

Accuracy assessment for the U.S. Geological Survey Regional Land-Cover Mapping Program: New York and New Jersey Region

The U.S. Geological Survey, in cooperation with other government and private organizations, is producing a conterminous U.S. land-cover map using Landsat Thematic Mapper 30-meter data for the Federal regions designated by the U.S. Environmental Protection Agency. Accuracy assessment is to be conducted for each Federal region to estimate overall and class-specific accuracies. In Region 2, consisting of New York and New Jersey, the accuracy assessment was completed for 15 land-cover and land-use classes, using interpreted 1:40,000-scale aerial photographs as reference data. The methodology used for Region 2 features a two-stage, geographically stratified approach, with a general sample of all classes (1,033 sample sites), and a separate sample for rare classes (294 sample sites). A confidence index was recorded for each land-cover interpretation on the 1:40,000-scale aerial photography The estimated overall accuracy for Region 2 was 63 percent (standard error 1.4 percent) using all sample sites, and 75.2 percent (standard error 1.5 percent) using only reference sites with a high-confidence index. User's and producer's accuracies for the general sample and user's accuracy for the sample of rare classes, as well as variance for the estimated accuracy parameters, were also reported. Narrowly defined land-use classes and heterogeneous conditions of land cover are the major causes of misclassification errors. Recommendations for modifying the accuracy assessment methodology for use in the other nine Federal regions are provided.

Photogrammetric Engineering and Remote Sensing

CLICK: The new USGS center for LIDAR information coordination and knowledge

Elevation data is rapidly becoming an important tool for the visualization and analysis of geographic information. The creation and display of three-dimensional models representing bare earth, vegetation, and structures have become major requirements for geographic research in the past few years. Light Detection and Ranging (lidar) has been increasingly accepted as an effective and accurate technology for acquiring high-resolution elevation data for bare earth, vegetation, and structures. Lidar is an active remote sensing system that records the distance, or range, of a laser fi red from an airborne or space borne platform such as an airplane, helicopter or satellite to objects or features on the Earth’s surface. By converting lidar data into bare ground topography and vegetation or structural morphologic information, extremely accurate, high-resolution elevation models can be derived to visualize and quantitatively represent scenes in three dimensions. In addition to high-resolution digital elevation models (Evans et al., 2001), other lidar-derived products include quantitative estimates of vegetative features such as canopy height, canopy closure, and biomass (Lefsky et al., 2002), and models of urban areas such as building footprints and three-dimensional city models (Maas, 2001).

Photogrammetric Engineering and Remote Sensing

Eruptive stratigraphy of the Tatara-San Pedro complex, 36°S, sourthern volcanic zone, Chilean Andes: Reconstruction method and implications for magma evolution at long-lived arc volcanic centers

The Quaternary Tatara–San Pedro volcanic complex (36°S, Chilean Andes) comprises eight or more unconformity-bound volcanic sequences, representing variably preserved erosional remnants of volcanic centers generated during ∼930 ky of activity. The internal eruptive histories of several dominantly mafic to intermediate sequences have been reconstructed, on the basis of correlations of whole-rock major and trace element chemistry of flows between multiple sampled sections, but with critical contributions from photogrammetric, geochronologic, and paleomagnetic data. Many groups of flows representing discrete eruptive events define internal variation trends that reflect extrusion of heterogeneous or rapidly evolving magma batches from conduit–reservoir systems in which open-system processes typically played a large role. Long-term progressive evolution trends are extremely rare and the magma compositions of successive eruptive events rarely lie on precisely the same differentiation trend, even where they have evolved from similar parent magmas by similar processes. These observations are not consistent with magma differentiation in large long-lived reservoirs, but they may be accommodated by diverse interactions between newly arrived magma inputs and multiple resident pockets of evolved magma and/or crystal mush residing in conduit-dominated subvolcanic reservoirs. Without constraints provided by the reconstructed stratigraphic relations, the framework for petrologic modeling would be far different. A well-established eruptive stratigraphy may provide independent constraints on the petrologic processes involved in magma evolution—simply on the basis of the specific order in which diverse, broadly cogenetic magmas have been erupted. The Tatara–San Pedro complex includes lavas ranging from primitive basalt to high-SiO 2 rhyolite, and although the dominant erupted magma type was basaltic andesite (∼52–55 wt % SiO 2 ) each sequence is characterized by unique proportions of mafic, intermediate, and silicic eruptive products. Intermediate lava compositions also record different evolution paths, both within and between sequences. No systematic long-term pattern is evident from comparisons at the level of sequences. The considerable diversity of mafic and evolved magmas of the Tatara–San Pedro complex bears on interpretations of regional geochemical trends. The variable role of open-system processes in shaping the compositions of evolved Tatara–San Pedro complex magmas, and even some basaltic magmas, leads to the conclusion that addressing problems such as arc magma genesis and elemental fluxes through subduction zones on the basis of averaged or regressed reconnaissance geochemical datasets is a tenuous exercise. Such compositional indices are highly instructive for identifying broad regional trends and first-order problems, but they should be used with extreme caution in attempts to quantify processes and magma sources, including crustal components, implicated in these trends.

Andes Mountains, Tatara-San Pedro Volcanic Complex

Validation of the ASTER instrument level 1A scene geometry

An independent assessment of the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) instrument geometry was undertaken by the U.S. ASTER Team, to confirm the geometric correction parameters developed and applied to Level 1A (radiometrically and geometrically raw with correction parameters appended) ASTER data. The goal was to evaluate the geometric quality of the ASTER system and the stability of the Terra spacecraft. ASTER is a 15-band system containing optical instruments with resolutions from 15- to 90-meters; all geometrically registered products are ultimately tied to the 15-meter Visible and Near Infrared (VNIR) sub-system. Our evaluation process first involved establishing a large database of Ground Control Points (GCP) in the mid-western United States; an area with features of an appropriate size for spacecraft instrument resolutions. We used standard U.S. Geological Survey (USGS) Digital Orthophoto Quads (DOQS) of areas in the mid-west to locate accurate GCPs by systematically identifying road intersections and recording their coordinates. Elevations for these points were derived from USGS Digital Elevation Models (DEMS). Road intersections in a swath of nine contiguous ASTER scenes were then matched to the GCPs, including terrain correction. We found no significant distortion in the images; after a simple image offset to absolute position, the RMS residual of about 200 points per scene was less than one-half a VNIR pixel. Absolute locations were within 80 meters, with a slow drift of about 10 meters over the entire 530-kilometer swath. Using strictly simultaneous observations of scenes 370 kilometers apart, we determined a stereo angle correction of 0.00134 degree with an accuracy of one microradian. The mid-west GCP field and the techniques used here should be widely applicable in assessing other spacecraft instruments having resolutions from 5 to 50-meters. ?? 2008 American Society for Photogrammetry and Remote Sensing.

Photogrammetric Engineering and Remote Sensing

Historical record of Landsat global coverage

The long-term, 34+ year record of global Landsat remote sensing data is a critical resource to study the Earth system and human impacts on this system. The National Satellite Land Remote Sensing Data Archive (NSLRSDA) is charged by public law to: “maintain a permanent, comprehensive Government archive of global Landsat and other land remote sensing data for long-term monitoring and study of the changing global environment” (U.S. Congress, 1992). The advisory committee for NSLRSDA requested a detailed analysis of observation coverage within the U.S. Landsat holdings, as well as that acquired and held by International Cooperator (IC) stations. Our analyses, to date, have found gaps of varying magnitude in U.S. holdings of Landsat global coverage data, which appear to reflect technical or administrative variations in mission operations. In many cases it may be possible to partially fill these gaps in U.S. holdings through observations that were acquired and are now being held at International Cooperator stations.

Photogrammetric Engineering and Remote Sensing

Biomass modeling of four water intensiveleading world crops using hyperspectral narrowbands in support of HyspIRI Mission

New satellite missions are expected to record high spectral resolution information globally and consistently for the first time, so it is important to identify modeling techniques that take advantage of these new data. In this paper, we estimate biomass for four major crops using ground-based hyperspectral narrowbands. The spectra and their derivatives are evaluated using three modeling techniques: two-band hyperspectral vegetation indices (HVIs), multiple band-HVIs (MB-HVIs) developed from Sequential Search Methods (SSM), and MB-HVIs developed from Principal Component Regression. Overall, the two-band HVIs and MB-HVIs developed from SSMs using first derivative transformed spectra in the visible blue and green and NIR explained more biomass variability and had lower error than the other approaches or transformations; however a better search criterion needs to be developed in order to reflect the true ability of the two-band HVI approach. Short-Wave Infrared 1 (1000 to 1700 nm) proved less effective, but still important in the final models.

Photogrammetric Engineering and Remote Sensing

Unmanned aircraft systems (UAS) activities at the Department of the Interior

The U.S. Department of the Interior (DOI) is responsible for protecting and managing the natural resources and heritage on almost 20% of the land in the United States. The DOI’s mission requires access to remotely sensed data over vast lands, including areas that are remote and potentially dangerous to access. Unmanned Aircraft Systems (UAS) technology has the potential to enable the DOI to be a better steward of the land by: (1) Improving natural hazard forecasting and the analysis of the impacts. (2) Improving the understanding of climate change to better plan for likely impacts. (3) Developing precipitation and evaporation forecasting to better manage water resources. (4) Monitoring Arctic ice change and its impacts on ecosystems, coasts, and transportation. (5) Increasing safety and effectiveness of wildland fire management. (6) Enhancing search and rescue capabilities. (7) Broadening the abilities to monitor environmental or landscape conditions and changes. (8) Better understanding and protecting the Nation’s ecosystems. The initial operational testing and evaluations performed by the DOI have proven that UAS technology can be used to support many of the Department’s activities. UAS technology provides scientists a way to look longer, closer and more frequently at some of Earth’s most remote areas—places that were previously too dangerous or expensive to monitor in detail. The flexibility of operations and relative low cost to purchase and operate Small Unmanned Aerial System (sUAS) enhances the ability to track long-term landscape and environmental change. The initial testing indicates the operational costs are approximately 10% of traditional manned aircraft. In addition, users can quickly assess landscape-altering events such as wildland fires, floods and volcanoes. UAS technology will allow the DOI to do more with less and in the process enhance the Department’s ability to provide unbiased scientific information to help stakeholders make informed decisions. It will also provide a digital baseline record that can be archived and used when monitoring future events or conditions. One possible future scenario has scientists carrying sUAS into the field allowing quick deployment and operation to observe the environment or for emergency response. This scenario could also include a persistent monitoring capability provided by a UAS that can stay airborne over a small geographic area for days or weeks, or possibly longer. While the DOI focus is on sUAS, the Department recognizes that larger UAS systems will also play a role in meeting its mission. The Department anticipates meeting long-duration or specialized acquisition commitments, such as state or national aerial photography, by collaboration with other agencies or through commercial contracts. Even though the DOI continues to evaluate UAS and sensor technology to meet the Department’s mission, some of its bureaus are already moving towards an operational capability. The authors fully anticipate that by 2020 UAS will emerge as one of the primary platforms for DOI remote sensing applications.

Photogrammetric Engineering and Remote Sensing

Shore erosion as a sediment source to the tidal Potomac River, Maryland and Virginia

The shoreline of the tidal Potomac River attained its present form as a result of the Holocene episode of sea-level rise; the drowned margins of the system are modified by wave activity in the shore zone and by slope processes on banks steepened by basal-wave erosion. Shore erosion leaves residual sand and gravel in shallow water and transports silt and clay offshore to form a measurable component of the suspended-sediment load of the tidal Potomac River. Erosion rates were measured by comparing digitized historical shoreline maps and modern maps, and by comparing stereopairs of aerial photographs taken at different points in time, with the aid of an interactive computer-graphics system and a digitizing stereoplotter. Cartographic comparisons encompassed 90 percent of the study reach and spanned periods of 38 to 109 years, with most measurements spanning at least 84 years. Photogrammetric comparisons encompassed 49 percent of the study reach and spanned 16 to 40 years. Field monitoring of erosion rates and processes at two sites, Swan Point Neck, Maryland, and Mason Neck, Virginia, spanned periods of 10 to 18 months. Estimated average recession rates of shoreline in the estuary, based on cartographic and photogrammetric measurements, were 0.42 to 0.52 meter per annum (Virginia shore) and 0.31 to 0.41 meter per annum (Maryland shore). Average recession rates of shoreline in the tidal river and transition zone were close to 0.15 meter per annum. Estimated average volume-erosion rates along the estuary were 1.20 to 1.87 cubic meters per meter of shoreline per annum (Virginia shore) and 0.56 to 0.73 cubic meter per meter of shoreline per annum (Maryland shore); estimated average volume-erosion rates along the shores of the tidal river and transition zone were 0.55 to 0.74 cubic meter per meter of shoreline per annum. Estimated total sediment contributed to the tidal Potomac River by shore erosion was 0.375 x 10 6 to 0.565 x 10 6 metric tons per annum; of this, the estimated amount of silt and clay ranged from 0.153x10 6 to 0.226x10 6 metric tons per annum. Between 49 and 60 percent of the sediment was derived from the Virginia shore of the estuary; 14 to 18 percent was derived from the Maryland shore of the estuary; and 23 to 36 percent was derived from the shores of the tidal river and transition zone. The adjusted modern estimate of sediment eroded from the shoreline of the estuary is about 55 percent of the historical estimate. Sediment eroded from the shoreline accounted for about 6 to 9 percent of the estimated total suspended load for the tidal Potomac River during water years 1979 through 1981 and for about 11 to 18 percent of the suspended load delivered to the estuary during the same period. Annual suspended-sediment loads derived from upland source areas fluctuated by about an order of magnitude during the 3 years of record (1979-81); shore erosion may have been a more important component of the sediment budget during periods of low flow than during periods of higher discharges. Prior to massive land clearance during the historical period of intensive agriculture in the 18th and 19th centuries, annual sediment loads from upland sources probably were smaller than they are at present; under these circumstances shore erosion would have been an important component of the sediment budget. At current rates of sediment supply, relative sea-level rise, and shoreline recession, the landward parts of the tidal Potomac River are rapidly being filled by sediment. If these rates were to remain constant over time, and no sediment were to escape into Chesapeake Bay, the tidal river and transition zone would be filled within 600 years, and the total system would be filled in less than 4,000 years. Given a slower rate of sediment supply, comparable to the measured rate during the low-flow 1981 water year, the volume of the tidal Potomac River might remain relatively stable or even increase over time. Changes in rates

Water Supply Paper

Channel change and sediment transport in the Puyallup River watershed through 2022

The Puyallup River drains a 990 square mile watershed in western Washington, with headwaters on the glacier-covered flanks of Mount Rainier. Major tributaries include the White, Carbon, and Mowich Rivers. In the levee-confined reaches of the lower watershed, loss of flood conveyance due to sand and gravel deposition has been a chronic issue. Over much of the 20th century, flood conveyance was maintained through sediment removal, but this practice ended in the late 1990s. Flood hazard management activities since the 1990s have primarily involved levee removal or setback projects. Assessments of 1984-2009 repeat cross sections suggested that sediment deposition rates were particularly high in reaches with recent levee setbacks. However, there have been no assessments of recent deposition rates since the 2009 surveys. There are also concerns that intensifying flood hydrology or increased sediment delivery from Mount Rainier may exacerbate deposition. However, assessment of those risks has been hindered by limited understanding of watershed-scale sediment delivery and routing, particularly for coarse sand and gravel. The U.S. Geological Survey, in cooperation with Pierce County, initiated this study to improve understanding of sediment deposition in the lower Puyallup River watershed. This work is primarily based on differencing of multiple aerial lidar datasets collected during 2002–2022, supplemented by early 1990 photogrammetric elevation datasets, geomorphic assessments of streamgage data, historical topographic surveys from 1907, and previously collected sediment transport measurements. Analyses cover the Puyallup, Carbon, and Mowich Rivers, but do not include the White River. During 2004–2020, repeat aerial lidar indicates that 1.3 ± 0.3 million yd3 of sediment accumulated in the lower 20 valley miles (VMs) of the Puyallup River, averaging 80,000 ± 20,000 cubic yards per year (yd3/yr). Deposition was observed during both 2004–11 and 2011–20 lidar differencing intervals. This continued a long-term depositional trend that extends back to at least 1977. From 2004 to 2011, deposition rates along the Soldiers Home levee setback reach, the only setback project downstream of VM 20 completed prior to 2011, were approximately four times higher than in adjacent unmodified reaches. From 2011 to 2020, two additional setback projects were completed; volumetric deposition rates over all three setback reaches were similar to adjacent unmodified reaches, suggesting elevated setback deposition in the 2004–11 interval may have been influenced by an extreme flood in November 2006. These levee setback projects increased the local cross-sectional area of the floodway, used as a rough proxy for relative flood conveyance, by 50 to 200 percent above 2004 conditions. If deposition continued at recent rates, cross-sectional area over the levee setback reaches would be reduced back to 2004 values by 2050-90. Deposition also occurred over the lower six VMs of the Carbon River during 2004–20, though volumes (0.15 ± 0.09 million yd3) were an order of magnitude lower than along the Puyallup River. Relatively lower deposition rates in the Carbon River are most likely the combined result of modestly lower incoming sediment loads, modestly steeper channel slope, and the additional sediment transport capacity provided by two large non-glacial tributaries that enter the Carbon River near VM 5. Upstream of the depositional reaches described above, 2002–22 sediment storage trends along the Puyallup, Carbon, and Mowich Rivers were predominately negative (net erosion) up to the Mount Rainier National Park boundary. Net erosion was the result of bank and bluff erosion exceeding deposition across wetted channel and bare gravel areas, as opposed to uniform vertical downcutting. Net erosion along these river valleys delivered 3.4 ± 0.6 million yd3 to the river system, equivalent to 190,000 ± 35,000 yd3/yr. Most of that volume was supplied by erosion of relatively low (4–10 ft) surfaces along the Puyallup and Mowich Rivers and tall (300 ft) glacial bluffs along the lower Carbon River. Substantial aggradation from 1984 to 2009 reported by Czuba and others (2010) along reaches of the Puyallup River (VM 19–22) where levee confinement has recently been removed was most likely an artifact of methodologic bias. The Puyallup, Mowich, and Carbon Rivers drain five distinct glaciated watersheds on the flanks of Mount Rainier, four of which were assessed in this study. All four watersheds were impacted by an extreme November 2006 rainstorm. Between 2002 and 2008, debris flows occurred in all four headwater areas, collectively eroding at least 2.1 million yd3 of sediment. These debris flows formed distinct deposits one to two miles downstream of source areas, depositing 30-50 percent of the material eroded upstream. From 2008 to 2022, no headwater debris flows were observed and overall rates of geomorphic change in the headwaters were low. Rivers eroded into debris flow deposits emplaced over the 2002–08 interval, but re-deposited equivalent volumes of material within a half mile downstream. Stage-discharge relations at five streamgages on upland rivers draining Mount Rainier show either net channel incision or dynamic variability with no long-term trend over the past 60–100 years. Observations of pervasive river valley erosion and stable or incising trends at long-term streamgages in the upper watershed do not support prior claims of widespread and accelerating aggradation of upland rivers draining Mount Rainier. Erosion and deposition volumes estimated in this report were combined with sediment transport estimates from limited suspended sediment and bedload measurements, estimates of sub-glacial erosion rates, and sediment delivery from non-glacial tributaries to construct watershed-scale sediment budgets for the Puyallup River watershed. During 2004–20, the estimated sediment load entering the depositional lowlands was well balanced by estimated inputs from, in order of relative magnitude, subglacial erosion (33–60 percent of total sediment load), erosion along the major river valleys (25–45 percent), erosion in recently deglaciated headwater areas (7–17 percent) and non-glacial tributaries (3–9 percent). These results are specific to the study period and represent total sediment loads, most of which is fine material carried in suspension. The relative sourcing of sand and gravel may be different than implied by this sediment budget. Downstream of VM 12, comparison of 1907 and 2009 channel surveys show net lowering of the channel thalweg of 4–12 ft. A long-term gage near VM 22 shows lowering of 4–5 ft through the 1960s. Lowering at both locations was inferred to be a channel response to the substantial straightening, and so steepening, of the river during major phases of levee construction through the early and mid-20th century. Application of a simple empirical bedload-discharge power-law relation to an ensemble of model-estimated daily mean discharge records in the lower Puyallup River between 1977 and 2100 projects that annual bedload transport capacity in the lower Puyallup River will increase by 20–60 percent by the middle of the 21st century. Actual changes in bedload transport and deposition rates will depend on concurrent changes in sediment supply and local hydraulics governing deposition. This report presents several key conclusions. First, the persistence and spatial patterns of sand and gravel deposition along the lower Puyallup River support prior claims that deposition is fundamentally caused by decreases in channel slope moving downstream. Given this underlying cause and the abundance of sand and gravel available to be transported downstream, deposition is likely to continue for the foreseeable future. Second, despite continued sediment deposition, recent levee setback projects in the lower Puyallup River will likely provide several decades of flood conveyance benefits relative to a no-action alternative. Third, while the rivers linking Mount Rainier to the Puget Sound lowlands have often been discussed as conduits that either pass or accumulate sediment from Mount Rainier, observations from 2002–22 show these river valleys acting as substantial sediment sources, delivering three times more sediment than recently deglaciated headwater areas on Mount Rainier. While the persistence and underlying cause of recent river valley erosion remain unknown, sediment storage dynamics along these river valleys are likely to be a major control on sand and gravel delivery to the lower watershed.

EarthArXiv