USGS Science⌕ Search

SEARCH · USGS Science

Results for “Landslides”

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 1,369 records · Page 76Linked to original sources
USGS2008

Recent bright gully deposits on Mars: Wet or dry flow?

Bright gully sediments attributed to liquid water flow have been deposited on Mars within the past several years. To test the liquid water flow hypothesis, we constructed a high-resolution (1 m/pixel) photogrammetric digital elevation model of a crater in the Centauri Montes region, where a bright gully deposit formed between 2001 and 2005. We conducted one-dimensional (1-D) and 2-D numerical flow modeling to test whether the deposit morphology is most consistent with liquid water or dry granular How. Liquid water flow models that incorporate freezing can match the runout distance of the flow for certain freezing rates but fail to reconstruct the distributary lobe morphology of the distal end of the deposit. Dry granular flow models can match both the observed runout distance and the distal morphology. Wet debris flows with high sediment concentrations are also consistent with the observed morphology because their rheologies are often similar to that of dry granular flows. As such, the presence of liquid water in this flow event cannot be ruled out, but the available evidence is consistent with dry landsliding.

Geology↗
USGS2025

Forest-floor burial in 1507 by the largest Mount Rainier lahar of the past millennium

New dating of lahar-killed trees underscores volcano hazards in the Puget Sound metropolitan area. Beginning as a landslide from the west flank of Mount Rainier, Washington, USA, the Electron Mudflow, which was the largest lahar of the last millennium, swept more than 60 km down the Puyallup River drainage into areas now densely populated. Wiggle matching of seven radiocarbon ages from buried, bark-bearing Douglas-fir ( Pseudotsuga menziesii ) trees brackets the mudflow's age between 1477 and 1522 CE with 99.7% certainty. To narrow this date, we applied dendrochronology crossdating on samples collected from 21 trees killed by the lahar, measuring 86 time series for statistical verification. The four bark-bearing trees died the same year while the final rings in all other trees had decayed, exposing rings formed in earlier years. When averaged together, the crossdated measurements form a 475 yr master chronology that was correlated against absolutely dated tree-ring chronologies in the region. The Electron chronology best matched with chronologies from low-elevation sites, especially a Douglas-fir chronology from Vancouver Island, Canada, to show that the Electron trees died in 1507 CE. Latewood in the final ring was beginning to form, indicating the mudflow likely occurred in the late-summer months. What caused the Electron Mudflow is unknown, but this precise date will help to assess possible relationships with other events, assist in interpreting Indigenous narratives about the mudflow, and increase awareness of potential lahar hazards.

Washington↗
USGS2009

Saddle Mountain fault deformation zone, Olympic Peninsula, Washington: Western boundary of the Seattle uplift

The Saddle Mountain fault, first recognized in the early 1970s, is now well mapped in the Hoodsport area, southeastern Olympic Peninsula (northwestern United States), on the basis of light detection and ranging (LIDAR) surveys, aerial photography, and trench excavations. Drowned trees and trench excavations demonstrate that the Saddle Mountain fault produced a M W 6.5–7.0 earthquake 1000–1300 yr ago, likely contemporaneous with the M W 7.5 Seattle fault earthquake 1100 yr ago and with a variety of other fault and landslide activity over a wide region of the Olympic Peninsula and Puget Lowland. This near synchroneity suggests that the Saddle Mountain and Seattle fault may be kinematically linked. Aeromagnetic anomalies and LIDAR topographic scarps define an en echelon sequence of faults along the southeastern Olympic Peninsula of Washington, all active in Holocene time. A detailed analysis of aeromagnetic data suggests that the Saddle Mountain fault extends at least 35 km, from 6 km southwest of Lake Cushman northward to the latitude of the Seattle fault. A magnetic survey over Price Lake using a nonmagnetic canoe illuminated two east-dipping reverse faults with 20 m of vertical offset at 30 m depth associated with 2–4 m of vertical displacement at the topographic surface. Analysis of regional aeromagnetic data indicates that the Seattle fault may extend westward across Hood Canal and into the Olympic Mountains, where it terminates near the northward terminus of the Saddle Mountain fault. The en echelon alignment of the Saddle Mountain and nearby Frigid Creek and Canyon River faults, all active in late Holocene time, reflects a >45-km-long zone of deformation that may accommodate the northward shortening of Puget Lowland crust inboard of the Olympic massif. In this view, the Seattle fault and Saddle Mountain deformation zone form the boundaries of the northward-advancing Seattle uplift.

Washington↗
USGS2012

Eruptive history of Mount Katmai, Alaska

Mount Katmai has long been recognized for its caldera collapse during the great pyroclastic eruption of 1912 (which vented 10 km away at Novarupta in the Valley of Ten Thousand Smokes), but little has previously been reported about the geology of the remote ice-clad stratovolcano itself. Over several seasons, we reconnoitered all parts of the edifice and sampled most of the lava flows exposed on its flanks and caldera rim. The precipitous inner walls of the 1912 caldera remain too unstable for systematic sampling; so we provide instead a photographic and interpretive record of the wall sequences exposed. In contrast to the several andesite-dacite stratovolcanoes nearby, products of Mount Katmai range from basalt to rhyolite. Before collapse in 1912, there were two overlapping cones with separate vent complexes and craters; their products are here divided into eight sequences of lava flows, agglutinates, and phreatomagmatic ejecta. Latest Pleistocene and Holocene eruptive units include rhyodacite and rhyolite lava flows along the south rim; a major 22.8-ka rhyolitic plinian fall and ignimbrite deposit; a dacite-andesite zoned scoria fall; a thick sheet of dacite agglutinate that filled a paleocrater and draped the west side of the edifice; unglaciated leveed dacite lava flows on the southeast slope; and the Horseshoe Island dacite dome that extruded on the caldera floor after collapse. Pre-collapse volume of the glaciated Katmai edifice was ∼30 km 3 , and eruptive volume is estimated to have been 57±13 km 3 . The latter figure includes ∼40±6 km 3 for the edifice, 5±2 km 3 for off-edifice dacite pyroclastic deposits, and 12±5 km 3 for the 22.8-ka rhyolitic pyroclastic deposits. To these can be added 13.5 km 3 of magma that erupted at Novarupta in 1912, all or much of which is inferred to have been withdrawn from beneath Mount Katmai. The oldest part of the edifice exposed is a basaltic cone, which gave a 40 Ar/ 39 Ar plateau age of 89 ± 25 ka. The seismic record of caldera collapse includes 14 earthquakes of magnitude 6.0–7.0. By combining the times of earthquakes, the hours of downwind plinian-fall episodes from Novarupta, and the stratigraphic record of hydrothermal explosion breccia and phreatic mud layers ejected around the caldera rim and intercalated within the Novarupta pumice-fall sequence, it can be inferred that collapse began in the 11th hour of the 60-h-long eruption and continued fitfully for 3.5 days. Several big landslides and pumiceous debris flows shaken loose by the collapse-related seismicity are bracketed in time by their levels of intercalation within the Novarupta pumice-fall sequence. An intracaldera lake was ∼10 m deep by 1916, drained away in 1923, and has since deepened progressively to ∼250 m today. Compositionally, products of Mount Katmai represent an ordinary medium-K arc array, both tholeiitic and calcalkaline, that extends from 51.6% to 72.3% SiO 2 . Values of 87 Sr/ 86 Sr range from 0.70335 to 0.70372, correlating loosely with fractionation indices. The 5–6 km 3 of continuously zoned andesite-dacite magma (58%–68% SiO 2 ) that erupted at Novarupta in 1912 was withdrawn from beneath Mount Katmai and bears close compositional affinity with products of that edifice, not with pre-1912 products of the adjacent Trident cluster. Evidence is presented that the 7–8 km 3 of high-silica rhyolite (77% SiO 2 ) released in 1912 is unlikely to have been stored under Novarupta or Trident. Pre-eruptive contiguity with the andesite-dacite reservoir is suggested by (1) eruption of ∼3 km 3 of rhyolite magma first, followed by mutual mingling in fluctuating proportions; (2) thermal and redox continuity of the whole zoned sequence despite the wide compositional gap; (3) Nd, Sr, O isotopic, and rare earth element (REE) affinities of the whole array; (4) compositional continuity of the nearly aphyric rhyolite with the glass (melt) phase of the phenocryst-rich dacite; and (5) phase-equilibrium experiments that indicate similar shallow pre-eruptive storage depths (3–6 km) for rhyolite, dacite, and andesite.

Alaska↗
USGS2015

Tsunami recurrence in the eastern Alaska-Aleutian arc: A Holocene stratigraphic record from Chirikof Island, Alaska

Despite the role of the Alaska-Aleutian megathrust as the source of some of the largest earthquakes and tsunamis, the history of its pre–twentieth century tsunamis is largely unknown west of the rupture zone of the great (magnitude, M 9.2) 1964 earthquake. Stratigraphy in core transects at two boggy lowland sites on Chirikof Island’s southwest coast preserves tsunami deposits dating from the postglacial to the twentieth century. In a 500-m-long basin 13–15 m above sea level and 400 m from the sea, 4 of 10 sandy to silty beds in a 3–5-m-thick sequence of freshwater peat were probably deposited by tsunamis. The freshwater peat sequence beneath a gently sloping alluvial fan 2 km to the east, 5–15 m above sea level and 550 m from the sea, contains 20 sandy to silty beds deposited since 3.5 ka; at least 13 were probably deposited by tsunamis. Although most of the sandy beds have consistent thicknesses (over distances of 10–265 m), sharp lower contacts, good sorting, and/or upward fining typical of tsunami deposits, the beds contain abundant freshwater diatoms, very few brackish-water diatoms, and no marine diatoms. Apparently, tsunamis traveling inland over low dunes and boggy lowland entrained largely freshwater diatoms. Abundant fragmented diatoms, and lake species in some sandy beds not found in host peat, were probably transported by tsunamis to elevations of >10 m at the eastern site. Single-aliquot regeneration optically stimulated luminescence dating of the third youngest bed is consistent with its having been deposited by the tsunami recorded at Russian hunting outposts in 1788, and with the second youngest bed being deposited by a tsunami during an upper plate earthquake in 1880. We infer from stratigraphy, 14 C-dated peat deposition rates, and unpublished analyses of the island’s history that the 1938 tsunami may locally have reached an elevation of >10 m. As this is the first record of Aleutian tsunamis extending throughout the Holocene, we cannot estimate source earthquake locations or magnitudes for most tsunami-deposited beds. We infer that no more than 3 of the 23 possible tsunamis beds at both sites were deposited following upper plate faulting or submarine landslides independent of megathrust earthquakes. If so, the Semidi segment of the Alaska-Aleutian megathrust near Chirikof Island probably sent high tsunamis southward every 180–270 yr for at least the past 3500 yr.

Alaska↗
USGS2015

An ignimbrite caldera from the bottom up: Exhumed floor and fill of the resurgent Bonanza caldera, Southern Rocky Mountain volcanic field, Colorado

Among large ignimbrites, the Bonanza Tuff and its source caldera in the Southern Rocky Mountain volcanic field display diverse depositional and structural features that provide special insights concerning eruptive processes and caldera development. In contrast to the nested loci for successive ignimbrite eruptions at many large multicyclic calderas elsewhere, Bonanza caldera is an areally isolated structure that formed in response to a single ignimbrite eruption. The adjacent Marshall caldera, the nonresurgent lava-filled source for the 33.9-Ma Thorn Ranch Tuff, is the immediate precursor for Bonanza, but projected structural boundaries of two calderas are largely or entirely separate even though the western topographic rim of Bonanza impinges on the older caldera. Bonanza, source of a compositionally complex regional ignimbrite sheet erupted at 33.12 ± 0.03 Ma, is a much larger caldera system than previously recognized. It is a subequant structure ∼20 km in diameter that subsided at least 3.5 km during explosive eruption of ∼1000 km 3 of magma, then resurgently domed its floor a similar distance vertically. Among its features: (1) varied exposure levels of an intact caldera due to rugged present-day topography—from Paleozoic and Precambrian basement rocks that are intruded by resurgent plutons, upward through precaldera volcanic floor, to a single thickly ponded intracaldera ignimbrite (Bonanza Tuff), interleaved landslide breccia, and overlying postcollapse lavas; (2) large compositional gradients in the Bonanza ignimbrite (silicic andesite to rhyolite ignimbrite; 60%–76% SiO 2 ); (3) multiple alternations of mafic and silicic zones within a single ignimbrite, rather than simple upward gradation to more mafic compositions; (4) compositional contrasts between outflow sectors of the ignimbrite (mainly crystal-poor rhyolite to east, crystal-rich dacite to west); (5) similarly large compositional diversity among postcollapse caldera-fill lavas and resurgent intrusions; (6) brief time span for the entire caldera cycle (33.12 to ca. 33.03 Ma); (7) an exceptionally steep-sided resurgent dome, with dips of 40°–50° on west and 70°–80° on northeast flanks. Some near-original caldera morphology has been erosionally exhumed and remains defined by present-day landforms (western topographic rim, resurgent core, and ring-fault valley), while tilting and deep erosion provide three-dimensional exposures of intracaldera fill, floor, and resurgent structures. The absence of Plinian-fall deposits beneath proximal ignimbrites at Bonanza and other calderas in the region is interpreted as evidence for early initiation of pyroclastic flows, rather than lack of a high eruption column. Although the absence of a Plinian deposit beneath some ignimbrites elsewhere has been interpreted to indicate that abrupt rapid foundering of the magma-body roof initiated the eruption, initial caldera collapse began at Bonanza only after several hundred kilometers of rhyolitic tuff had erupted, as indicated by the minor volume of this composition in the basal intracaldera ignimbrite. Caldera-filling ignimbrite has been largely stripped from the southern and eastern flank of the Bonanza dome, exposing large areas of caldera-floor as a structurally coherent domed plate, bounded by ring faults with locations that are geometrically closely constrained even though largely concealed beneath valley alluvium. The structurally coherent floor at Bonanza contrasts with fault-disrupted floors at some well-exposed multicyclic calderas where successive ignimbrite eruptions caused recurrent subsidence. Floor rocks at Bonanza are intensely brecciated within ∼100 m inboard of ring faults, probably due to compression and crushing of the subsiding floor in proximity to steep inward-dipping faults. Upper levels of the floor are locally penetrated by dike-like crack fills of intracaldera ignimbrite, interpreted as dilatant fracture fills rather than ignimbrite vents. The resurgence geometry at Bonanza has implications for intracaldera-ignimbrite volume; this parameter may have been overestimated at some young calderas elsewhere, with bearing on outflow-intracaldera ratios and times of initial caldera collapse. Such features at Bonanza provide insights for interpreting calderas universally, with respect to processes of caldera collapse and resurgence, inception of subsidence in relation to progression of the ignimbrite eruption, complications with characterizing structural versus topographic margins of calderas, contrasts between intra- versus extracaldera ignimbrite, and limitations in assessing volumes of large caldera-forming eruptions. Bonanza provides a rare site where intact caldera margins and floor are exhumed and exposed, providing valuable perspectives for understanding younger similar calderas in some of the world’s most active and dangerous silicic provinces.

Colorado↗
USGS2019

Syn-collisional exhumation of hot middle crust in the Adirondack Mountains (New York, USA): Implications for extensional orogenesis in the southern Grenville province

Extensional deformation in the lower to middle continental crust is increasingly recognized and shown to have significant impact on crustal architecture, magma emplacement, fluid flow, and ore deposits. Application of the concept of extensional strain to ancient orogenic systems, like the Grenville province of eastern North America, has helped decipher the structural evolution of these regions. The Marcy massif is a ~3000 km2 Mesoproterozoic anorthosite batholith in the Adirondack Mountains (New York, USA) of the southern Grenville province. Bedrock geology mapping at 1:24,000 scale paired with characterization of bedrock exposed by recent landslides provides a glimpse into the structural architecture of the massif and its margin. New data demonstrate granulite- to amphibolite-facies deformational fabrics parallel the margin of the batholith, and that the Marcy massif is draped by a southeast-directed detachment zone. Within the massif, strain is localized into mutually offsetting conjugate shear zones with antithetic kinematic indicators. These relationships indicate that strain was coaxial within the Marcy massif, and that subsimple shear components of strain were partitioned along its margin. In situ U–Th–total Pb monazite analysis shows that deformation around and over the Marcy massif occurred from 1070 to 1060 Ma during granulite-facies metamorphism, and monazite from all samples record evidence for fluid-mediated dissolution reprecipitation from 1050 to 980 Ma. We interpret that rocks cooled isobarically after accretionary orogenesis and emplacement of the anorthosite- mangerite- charnockitegranite plutonic suite at ca. 1160–1140 Ma. Gravitational collapse during the Ottawan phase of the Grenville orogeny initiated along a southeast-directed detachment zone (Marcy massif detachment zone), which accommodated intrusion of the Lyon Mountain Granite Gneiss, and facilitated substantial fluid flow that catalyzed the formation of major ore deposits in the Adirondack Highlands.

New York↗
USGS2021

Latest Quaternary slip rates of the San Bernardino strand of the San Andreas fault, southern California, from Cajon Creek to Badger Canyon

Four new latest Pleistocene slip rates from two sites along the northwestern half of the San Bernardino strand of the San Andreas fault suggest the slip rate decreases southeastward as slip transfers from the Mojave section of the San Andreas fault onto the northern San Jacinto fault zone. At Badger Canyon, offsets coupled with radiocarbon and optically stimulated luminescence (OSL) ages provide three independent slip rates (with 95% confidence intervals): (1) the apex of the oldest dated alluvial fan (~30–28 ka) is right-laterally offset about 300–400 m yielding a slip rate of 13.5 +2.2/-2.5 mm/yr; (2) a terrace riser incised into the northwestern side of this alluvial fan is offset about 280-290 m and was abandoned ~23 ka, yielding a slip rate of 11.9 +0.9/-1.2 mm/yr; and (3) a younger alluvial fan (13-15 ka) has been offset 120–200 m from the same source canyon, yielding a slip rate of 11.8 +4.2/-3.5 mm/yr. These rates are all consistent and result in a preferred, time-averaged rate for the past ~28 ka of 12.8 +5.3/-4.7 mm/yr (95% confidence interval), with an 84% confidence interval of 10-16 mm/yr. At Matthews Ranch, in Pitman Canyon, about 13 km northwest of Badger Canyon, a landslide offset ~ 650 m with a 10Be age of ~47 ka yields a slip rate of 14.5 +9.9/-6.2 mm/yr (95% confidence interval). All of these slip rates for the San Bernardino strand are significantly slower than a previously published rate of 24.5 ± 3.5 mm/yr at the southern end of the Mojave section of the San Andreas fault (Weldon and Sieh, 1985), suggesting that ~12 mm/yr of slip transfers from the Mojave section of the San Andreas fault to the northern San Jacinto fault zone (and other faults) between Lone Pine Canyon and Badger Canyon, with most (if not all) of this slip transfer happening near Cajon Creek. This has been a consistent behavior of the fault for at least the past ~47 ka.

California↗
USGS2017

The Great Acceleration and the disappearing surficial geologic record

The surficial geologic record is the relatively thin veneer of young (<~1 Ma) and mostly unconsolidated sediments that cover portions of Earth’s terrestrial surface (Fig. 1). Once largely ignored as “overburden” by geologists, surficial deposits are now studied to address a wide range of issues related to the sustainability of human societies. Geologists use surficial deposits to determine the frequency and severity of past climatic changes, quantify natural and anthropogenic erosion rates, identify hazards, and calculate recurrence intervals associated with earthquakes, landslides, tsunamis, and volcanic eruptions. Increasingly, however, humans are eradicating the surficial geologic record in many key areas through progressive modification of Earth’s surface.

GSA Today↗
USGS1990

Chapter 22: Changing patterns of extensional tectonics; Overprinting of the basin of the middle and upper Miocene Esmeralda Formation in western Nevada by younger structural basins

The middle and upper Miocene Esmeralda Formation of western Nevada was deposited in a continental basin that crops out over an area of about 2,000 km 2 . The formation consists of thin, westerly derived sedimentary rocks in the western three-quarters of the outcrop area and of thick (3+ km) easterly derived sedimentary rocks in the eastern quarter. Megabreccias along the eastern margin of the basin are interpreted as landslide deposits derived from fault scarps. The basin was probably a half-graben with a major fault or faults on the east side. The position of the basin, its size, and the inferred major syndepositional faults on its east side are all unrelated to present-day topography and the distribution of major faults in the area and indicate a change in paleogeography and structural pattern since the late Miocene. In the eastern part of the basin, this change was accompanied by deformation that includes low-angle-fault detachment of the Esmeralda Formation from underlying amphibolite-grade Late Proterozoic and lower Paleozoic rocks in the Mineral Ridge-Weepah Hills area and uplift, folding, faulting, tilting, and surface exposure of the entire 3+ km thickness of the formation. Middle and upper Miocene basins in the northern part of the Basin and Range province commonly have been attributed to the onset of basin-range tectonism. The tectonic history of the Esmeralda Formation, however, indicates that some of these basins do not occupy the same area, nor are they related to the same syndepositional faults, as present basins. The extensional basin of the Esmeralda Formation can be viewed either as an early structure in an evolving, but kinematically related, extensional terrane in which the distribution of basins and faults changed gradually with time, or as the product of an extensional event kinematically different from that which produced present structures in the same area. The latter hypothesis is favored because of the marked contrasts in the paleogeography and structures associated with the Esmeralda basin compared with modern basins.

Nevada↗
USGS1988

Large-scale bedforms in boulder gravel produced by giant waves in Hawaii

Approximately 105,000 yr ago (based on uranium-series dating), waves in a giant wave train swept up to an elevation of about 375 m on the island of Lanai. The waves deposited the Hulopoe Gravel, which near the present shoreline consists of basalt boulders, coral fragments, and calcareous beachrock slabs, and near the upper limit of the deposit consists of sand and shell fragments. The maximum heights of similar but lower deposits on nearby islands, when adjusted for their estimated subsidence due to volcanic loading during the past 105,000 yr, indicate that the source of the wave was about 50 km southwest of Lanai. We hypothesize that failure and downward movement of the huge Lanai submarine landslide created an ocean disturbance, which produced waves that rushed across the Lanai reef and beach, picked up limestone and lava fragments, and deposited them high on the island as the Hulopoe Gravel. Backwash from the waves stripped soil and rock from the islands and carried much of it to the sea. The Hulopoe Gravel is 8 m thick in a gulch 200 m inland from the Lanai shoreline, where it consists of three beds, successively, 2, 4, and 2 m thick. These beds are considered to have been laid down by successive waves in the wave train. Each bed consists of two units: a lower unit of basalt and limestone boulders, cobbles, and sand, and an upper bimodal unit of large basalt boulders with a pebbly sand matrix. These subunits are inferred to be deposited from the runup and backwash of each wave. At the upper surface of the Hulopoe Gravel, basalt boulders averaging 0.5 m in diameter are arranged in dunelike ridges about 1 m high and 10 m apart. Nearby, where young streams have cut into and exposed the lower beds of the Hulopoe, clasts at the boundaries between the beds are locally imbricated and dip landward. We interpret these features as aspects of torrential flow and crossbedding created during the high-speed backwash of the great waves.

Hawaii↗
USGS1993

Intensity and duration of chemical weathering: An example from soil clays of the southeastern Koolau Mountains, Oahu, Hawaii

Orographic precipitation on the southern flank of the southeastern Koolau Mountains produces a pronounced precipitation gradient. The corresponding gradient in the intensity of the chemical weathering environment provides an opportunity to address the effects of varying chemical weathering intensity on the composition of clay-size weathering products in soils developed on basalt. In addition, little-modified remnants of the constructional surface of the Koolau Volcano, isolated by stream dissection, remain as facets on the southern ends of the parallel ridges of the study area. By comparing clay mineralogy of soils developed on these older geomorphic surfaces with those developed on the younger sharp-crested ridges and steep side slopes, the effects of weathering duration on clay mineralogy can also be addressed. Soil clays in this part of the Koolau Mountains are mineralogically complex; principal phases include smectite, kaolinite, and halloysite, but pure end member phases are uncommon. Rather, most phases contain some amount of mixed layering. Smectite may contain small (<5%) amounts of randomly interstratified halloysite. Similarly, kaolinite commonly contains a small proportion of halloysite interlayers. A complex halloysitic phase shows evidence of interstratification with both smectite and kaolinite. Nonphyllosilicates found in the clay fraction include gibbsite, goethite, rare quartz, and perhaps cristobalite. The gradient in precipitation is reflected in soil clay mineralogy by varying proportions of dominantly smectitic, kaolinitic, and halloysitic phases. In regions of relatively low precipitation (<2,000 mm/yr), soils are dominated by the smectitic and halloysitic phases. With increased precipitation (as much as ∼4,000 mm/yr), kaolinitic and halloysitic phases become the dominant clay minerals, and goethite and gibbsite become increasingly abundant. Older soils developed on geomorphic surfaces representing the original constructional surface of Koolau Volcano are markedly more leached than those from younger landscapes in the same precipitation regime. Although smectite may be present, kaolinite is the dominant phase, and accumulations of Fe and Ti occur in the uppermost soil levels. Enrichment of Zr and Ti in these soils, as compared to concentrations in the original basaltic parent material, indicates that as much as 75% of the parent material has been lost. Thus weathering duration may affect soil clay composition in the same way as weathering intensity. Because smectite and halloysite are expandable clay minerals, their presence in soils may decrease slope stability and influence the nature of slope processes. Soil avalanches occur on steep slopes throughout the study area, whereas slow-moving landslides appear to be restricted to gentler slopes in drier parts of the study area where smectite is abundant. The clay mineralogy of soils thus appears to influence the nature of slope processes in the southeastern Koolau Mountains.

Special Paper of the Geological Society of America↗
USGS2012

Submarine slides north of Puerto Rico and their tsunami potential

New multibeam bathymetry of the entire Puerto Rico trench reveals numerous retrograde slope failures at various scales at the edge of the carbonate platform north of Puerto Rico and the Virgin Islands. The slumped material comprises carbonate blocks, which fail, at least in initial stages, as a coherent rock mass. This, combined with the fact that the edge of the carbonate platform is steeper than most continental slopes, indicates a higher potential for tsunami run-up than along many other U.S. coasts. The style of slope failure appears to be rock falls, slide blocks and debris avalanches. Secondary failure of the failed carbonate products and of the underlying forearc sediments and rocks may lead to debris flows and turbidity flows. Fissures, discovered in the ocean floor near the edge of the platform, indicate that the process is expected to continue in the future. One of the slope failures, the Arecibo amphitheater, previously thought to represent a single giant slide with a volume of 900-1500 cu. km, appears to comprise smaller failures. The expected maximum tsunami run-up on the northern coast of Puerto Rico from one of these slope failures is <20 m, much lower than previously estimated. A recurrence rate of 100s kyr was calculated for the largest observed failures along the edge of the carbonate platform north of Puerto Rico using simple assumptions. Smaller landslides would presumably occur at a higher recurrence rate. Elsewhere around the island, a 20-km wide failure scarp was discovered in the Upper Mona rift and could be associated with the 1918 tsunami and earthquake that hit northwestern Puerto Rico. Large slope failures were also discovered for the first time on the northern side of the Puerto Rico trench. Because Puerto Rico trench slides occur at large water depths (~6000 m), have large horizontal and vertical dimensions, and the directivity from tsunamis emanating from these slides is toward Puerto Rico, they may be of particular concern and may necessitate further study.

Puerto Rico↗
USGS2018

Slope failure and mass transport processes along the Queen Charlotte Fault, southeastern Alaska

The Queen Charlotte Fault defines the Pacific–North America transform plate boundary in western Canada and southeastern Alaska for c. 900 km. The entire length of the fault is submerged along a continental margin dominated by Quaternary glacial processes, yet the geomorphology along the margin has never been systematically examined due to the absence of high-resolution seafloor mapping data. Hence the geological processes that influence the distribution, character and timing of mass transport events and their associated hazards remain poorly understood. Here we develop a classification of the first-order shape of the continental shelf, slope and rise to examine potential relationships between form and process dominance. We found that the margin can be split into six geomorphic groups that vary smoothly from north to south between two basic end-members. The northernmost group (west of Chichagof Island, Alaska) is characterized by concave-upwards slope profiles, gentle slope gradients (<6°) and relatively low along-strike variance, all features characteristic of sediment-dominated siliciclastic margins. Dendritic submarine canyon/channel networks and retrogressive failure complexes along relatively gentle slope gradients are observed throughout the region, suggesting that high rates of Quaternary sediment delivery and accumulation played a fundamental part in mass transport processes. Individual failures range in area from 0.02 to 70 km2 and display scarp heights between 10 and 250 m. Transpression along the Queen Charlotte Fault increases southwards and the slope physiography is thus progressively more influenced by regional-scale tectonic deformation. The southernmost group (west of Haida Gwaii, British Columbia) defines the tectonically dominated end-member: the continental slope is characterized by steep gradients (>20°) along the flanks of broad, margin-parallel ridges and valleys. Mass transport features in the tectonically dominated areas are mostly observed along steep escarpments and the larger slides (up to 10 km2) appear to be failures of consolidated material along the flanks of tectonic features. Overall, these observations highlight the role of first-order margin physiography on the distribution and type of submarine landslides expected to occur in particular morphological settings. The sediment-dominated end-member allows for the accumulation of under-consolidated Quaternary sediments and shows larger, more frequent slides; the rugged physiography of the tectonically dominated end-member leads to sediment bypass and the collapse of uplifted tectonic features. The maximum and average dimensions of slides are an order of magnitude smaller than those of slides observed along other (passive) glaciated margins. We propose that the general patterns observed in slide distribution are caused by the interplay between tectonic activity (long- and short-term) and sediment delivery. The recurrence (<100 years) of M > 7 earthquakes along the Queen Charlotte Fault may generate small, but frequent, failures of under-consolidated Quaternary sediments within the sediment-dominated regions. By contrast, the tectonically dominated regions are characterized by the bypass of Quaternary sediments to the continental rise and the less frequent collapse of steep, uplifted and consolidated sediments.

Alaska↗
USGS2018

Slope failure and mass transport processes along the Queen Charlotte Fault Zone, western British Columbia

Multibeam echosounder (MBES) images, 3.5 kHz seismic-reflection profiles and piston cores obtained along the southern Queen Charlotte Fault Zone are used to map and date mass-wasting events at this transform margin – a seismically active boundary that separates the Pacific Plate from the North American Plate. Whereas the upper continental slope adjacent to and east (upslope) of the fault zone offshore of the Haida Gwaii is heavily gullied, few large-sized submarine landslides in this area are observed in the MBES images. However, smaller submarine seafloor slides exist locally in areas where fluid flow appears to be occurring and large seafloor slides have recently been detected at the base of the steep continental slope just above its contact with the abyssal plain on the Queen Charlotte Terrace. In addition, along the subtle slope re-entrant area offshore of the Dixon Entrance shelf bathymetric data suggest that extensive mass wasting has occurred in the vicinity of an active mud volcano venting gas. We surmise that the relative lack of submarine slides along the upper slope in close proximity to the Queen Charlotte Fault Zone may be the result of seismic strengthening (compaction and cohesion) of a sediment-starved shelf and slope through multiple seismic events.

British Columbia↗
USGS2023

DisasterNet: Causal Bayesian networks with normalizing flows for cascading hazards

Sudden-onset hazards like earthquakes often induce cascading secondary hazards (e.g., landslides, liquefaction, debris flows, etc.) and subsequent impacts (e.g., building and infrastructure damage) that cause catastrophic human and economic losses. Rapid and accurate estimates of these hazards and impacts are critical for timely and effective post-disaster responses. Emerging remote sensing techniques provide pre- and post-event satellite images for rapid hazard estimation. However, hazards and damage often co-occur or colocate with underlying complex cascading geophysical processes, making it challenging to directly differentiate multiple hazards and impacts from satellite imagery using existing single-hazard models. We introduce DisasterNet, a novel family of causal Bayesian networks to model processes that a major hazard triggers cascading hazards and impacts and further jointly induces signal changes in remotely sensed observations. We integrate normalizing flows to effectively model the highly complex causal dependencies in this cascading process. A triplet loss is further designed to leverage prior geophysical knowledge to enhance the identifiability of our highly expressive Bayesian networks. Moreover, a novel stochastic variational inference with normalizing flows is derived to jointly approximate posteriors of multiple unobserved hazards and impacts from noisy remote sensing observations. Integrating with the USGS Prompt Assessment of Global Earthquakes for Response (PAGER) system, our framework is evaluated in recent global earthquake events. Evaluation results show that DisasterNet significantly improves multiple hazard and impact estimation compared to existing USGS products.

Conference Paper↗
USGS2000

Geomorphometry-Diversity in quantitative surface analysis

A widening variety of applications is diversifying geomorphometry ( digital terrain modelling ), the quantitative study of topography. An amalgam of earth science, mathematics, engineering and computer science, the discipline has been revolutionized by the computer manipulation of gridded terrain heights, or digital elevation models (DEMs). Its rapid expansion continues. This article reviews the remarkable diversity of recent morphometric work in 15 selected topics and discusses their significance and prospects. The quantitative analysis of industrial microsurface topography is introduced to the earth science community. The 14 other topics are Internet access to geomorphometry; global DEMs; DEM modelling of channel networks; self-organized criticality; fractal and wavelet analysis; soil resources; landslide hazards; barchan dunes; harvesting wind energy; sea-ice surfaces; sea-floor abyssal hills; Japanese work in morphometry; and the emerging fields of landscape ecology and image understanding. Closing remarks note reasons for the diversity within geomorphometry, speculate on future trends and recommend creating a unified field of surface representation.

Progress in Physical Geography: Earth and Environm↗
USGS2017

The role of initial coherence and path materials in the dynamics of three rock avalanche case histories

Background Rock avalanches are flow-like landslides that can travel at extremely rapid velocities and impact surprisingly large areas. The mechanisms that lead to the unexpected mobility of these flows are unknown and debated. Mechanisms proposed in the literature can be broadly classified into those that rely on intrinsic characteristics of the rock avalanche material, and those that rely on extrinsic factors such as path material. In this work a calibration-based numerical model is used to back-analyze three rock avalanche case histories. The results of these back-analyses are then used to infer factors that govern rock avalanche motion Results Our study has revealed two key insights that must be considered when analyzing rock avalanches. Results from two of the case histories demonstrate the importance of accounting for the initially coherent phase of rock avalanche motion. Additionally, the back-analyzed basal resistance parameters, as well as the best-fit rheology, are different for each case history. This suggests that the governing mechanisms controlling rock avalanche motion are unlikely to be intrinsic. The back-analyzed strength parameters correspond well to those that would be expected by considering the path material that the rock avalanches overran. Conclusion Our results show that accurate simulation of rock avalanche motion must account for the initially coherent phase of movement, and that the mechanisms governing rock avalanche motion are unlikely to be intrinsic to the failed material. Interaction of rock avalanche debris with path materials is the likely mechanism that governs the motion of many rock avalanches.

Geoenvironmental Disasters↗