USGS Science⌕ Search

SEARCH · USGS Science

Results for “One Earth”

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,279 records · Page 71Linked to original sources
USGS2016

Illuminating wildfire erosion and deposition patterns with repeat terrestrial lidar

Erosion following a wildfire is much greater than background erosion in forests because of wildfire-induced changes to soil erodibility and water infiltration. While many previous studies have documented post-wildfire erosion with point and small plot-scale measurements, the spatial distribution of post-fire erosion patterns at the watershed scale remains largely unexplored. In this study lidar surveys were collected periodically in a small, first-order drainage basin over a period of 2 years following a wildfire. The study site was relatively steep with slopes ranging from 17° to > 30°. During the study period, several different types of rain storms occurred on the site including low-intensity frontal storms (2.4 mm h −1 ) and high-intensity convective thunderstorms (79 mm h −1 ). These storms were the dominant drivers of erosion. Erosion resulting from dry ravel and debris flows was notably absent at the site. Successive lidar surveys were subtracted from one another to obtain digital maps of topographic change between surveys. The results show an evolution in geomorphic response, such that the erosional response after rain storms was strongly influenced by the previous erosional events and pre-fire site morphology. Hillslope and channel roughness increased over time, and the watershed armored as coarse cobbles and boulders were exposed. The erosional response was spatially nonuniform; shallow erosion from hillslopes (87% of the study area) contributed 3 times more sediment volume than erosion from convergent areas (13% of the study area). However, the total normalized erosion depth (volume/area) was highest in convergent areas. From a detailed understanding of the spatial locations of erosion, we made inferences regarding the processes driving erosion. It appears that hillslope erosion is controlled by rain splash (for detachment) and overland flow (for transport and quasi-channelized erosion), with the sites of highest erosion corresponding to locations with the lowest roughness. By contrast, in convergent areas we found erosion caused by overland flow. Soil erosion was locally interrupted by immobile objects such as boulders, bedrock, or tree trunks, resulting in a patchy erosion network with increasing roughness over time.

Colorado↗
USGS2014

Three-dimensional imaging, change detection, and stability assessment during the centerline trench levee seepage experiment using terrestrial light detection and ranging technology, Twitchell Island, California, 2012

A full scale field seepage test was conducted on a north-south trending levee segment of a now bypassed old meander belt on Twitchell Island, California, to understand the effects of live and decaying root systems on levee seepage and slope stability. The field test in May 2012 was centered on a north-south trench with two segments: a shorter control segment and a longer seepage test segment. The complete length of the trench area measured 40.4 meters (m) near the levee centerline with mature trees located on the waterside and landside of the levee flanks. The levee was instrumented with piezometers and tensiometers to measure positive and negative porewater pressures across the levee after the trench was flooded with water and held at a constant hydraulic head during the seepage test—the results from this component of the experiment are not discussed in this report. We collected more than one billion three-dimensional light detection and ranging (lidar) data points before, during, and after the centerline seepage test to assess centimeter-scale stability of the two trees and the levee crown. During the seepage test, the waterside tree toppled (rotated 20.7 degrees) into the water. The landside tree rotated away from the levee by 5 centimeters (cm) at a height of 2 m on the tree. The paved surface of the levee crown had three regions that showed subsidence on the waterside of the trench—discussed as the northern, central, and southern features. The northern feature is an elongate region that subsided 2.1 cm over an area with an average width of 1.35 m that extends 15.8 m parallel to the trench from the northern end of the trench to just north of the trench midpoint, and is associated with a crack 1 cm in height that formed during the seepage test on the trench wall. The central subsidence feature is a semicircular region on the waterside of the trench that subsided by as much as 6.2 cm over an area 3.4 m wide and 11.2 m long. The southern feature is an elongate region that has a maximum subsidence of 3.5 cm over an area 0.75 m wide and 8.1 m long and is associated with a number of small fractures in the pavement that are predominately north-south-trending and parallel to the trench. We determined that there was no significant motion of the levee flank during the last week of the seepage test. We also determined biomorphic parameters for the landside tree, such as the 3D positioning on the levee, tree height, levee parallel/perpendicular cross sectional area, and canopy centroid. These biomorphic parameters were requested to support a University of California Berkeley team studying seepage and stability on the levee. A gridded, 2-cm bare-earth digital elevation model of the levee crown and the landside levee flank from the final terrestrial lidar (T-Lidar) survey provided detailed topographic data for future assessment. Because the T-Lidar was not integrated into the project design, other than an initial courtesy dataset to help characterize the levee surface, our ability to contribute to the overall science goals of the seepage test was limited. Therefore, our analysis focused on developing data collection and processing methodology necessary to align ultra high-resolution T-Lidar data (with an average spot spacing 2–3 millimeters on the levee crown) from several instrument setup locations to detect, measure, and characterize dynamic centimeter-scale deformation and surface changes during the seepage test.

California↗
USGS2024

Seismic attenuation and stress on the San Andreas Fault at Parkfield: Are we critical yet?

The Parkfield transitional segment of the San Andreas Fault (SAF) is characterized by the production of frequent quasi-periodical M6 events that break the very same asperity. The last Parkfield mainshock occurred on 28 September 2004, 38 years after the 1966 earthquake, and after the segment showed a ∼22 years average recurrence time. The main reason for the much longer interevent period between the last two earthquakes is thought to be the reduction of the Coulomb stress from the M6.5 Coalinga earthquake of 2 May 1983, and the M6 Nuñez events of June 11th and 22 July 1983. Plausibly, the transitional segment of the SAF at Parkfield is now in the late part of its seismic cycle and current observations may all be relative to a state of stress close to criticality. However, the behavior of the attenuation parameter in the last few years seems substantially different from the one that characterized the years prior to the 2004 mainshock. A few questions arise: (i) Does a detectable preparation phase for the Parkfield mainshocks exist, and is it the same for all events? (ii) How dynamically/kinematically similar are the quasi-periodic occurrences of the Parkfield mainshocks? (iii) Are some dynamic/kinematic characteristics of the next mainshock predictable from the analysis of current data? (e.g., do we expect the epicenter of the next failure to be co-located to that of 2004?) (iv) Should we expect the duration of the current interseismic period to be close to the 22-year “undisturbed” average value? We respond to the questions listed above by analyzing the non-geometric attenuation of direct S-waves along the transitional segment of the SAF at Parkfield, in the close vicinity of the fault plane, between January 2001 and November 2023. Of particular interest is the preparatory behavior of the attenuation parameter as the 2004 mainshock approached, on both sides of the SAF. We also show that the non-volcanic tremor activity modulates the seismic attenuation in the area, and possibly the seismicity along the Parkfield fault segment, including the occurrence of the mainshocks.

California↗
USGS2015

Aurora painting pays tribute to Civil War's end

This year marks the sesquicentennial anniversary of the end of the American Civil War, a conflict that Abraham Lincoln called a “mighty scourge.” It was one of the most poignant periods in U.S. history, laying bare political, economic, social, and moral divergence between Northern and Southern states. The cause of the divergence that led to war was slavery [e.g., McPherson , 1988, chap. 3]—an institution that, by the 19th century, had been effectively abolished in the North but remained firmly entrenched in the South. War erupted in 1861 after a confederacy of Southern states declared secession from the Union of the United States. When the war finally ended in 1865, the Union had prevailed, and afterward, slavery was abolished throughout the United States. This outcome was obtained at the cost of 750,000 American lives and substantial destruction, especially in the South [e.g., Gugliotta , 2012]. In 1865, the same year the war ended, the American landscape artist Frederic Edwin Church unveiled Aurora Borealis (pictured above), a dramatic and mysterious painting that can be interpreted in terms of 19th century romanticism, scientific philosophy, and Arctic missions of exploration. Aurora Borealis can also be viewed as a restrained tribute to the end of the Civil War—a moving example of how science and current events served as the muses of late romantic artists [e.g., Carr , 1994, p. 277; Avery , 2011; Harvey , 2012].

Eos, Earth and Space Science News↗
USGS2012

Seismic constraints on the nature of lower crustal reflectors beneath the extending Southern Transition Zone of the Colorado Plateau, Arizona

We determine the reflection polarity and exploit variations in P and S wave reflectivity and P wave amplitude versus offset (AVO) to constrain the origin of lower crustal reflectivity observed on new three-component seismic data recorded across the structural transition of the Colorado Plateau. The near vertical incidence reflection data were collected by Stanford University in 1989 as part of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that traversed the Arizona Transition Zone of the Colorado Plateau. The results of independent waveform modeling methods are consistent with much of the lower crustal reflectivity resulting from thin, high-impedance layers. The reflection polarity of the cleanest lower crustal events is positive, which implies that these reflections result from high-velocity contrasts, and the waveform character indicates that the reflectors are probably layers less than or approximately equal to 200 m thick. The lower crustal events are generally less reflective to incident S waves than to P waves, which agrees with the predicted behavior of high-velocity mafic layering. Analysis of the P wave AVO character of lower crustal reflections demonstrates that the events maintain a constant amplitude with offset, which is most consistent with a mafic-layering model. One exception is a high-amplitude (10 dB above background) event near the base of lower crustal reflectivity which abruptly decreases in amplitude at increasing offsets. The event has a pronounced S wave response, which along with its negative AVO trend is a possible indication of the presence of fluids in the lower crust. The Arizona Transition Zone is an active but weakly extended province, which causes us to discard models of lower crustal layering resulting from shearing because of the high degree of strain required to create such layers. Instead, we favor horizontal basaltic intrusions as the primary origin of high-impedance reflectors based on (1) The fact that most xenoliths in eruptive basalts of the Transition Zone are of mafic igneous composition, (2) indications that a pulse of magmatic activity crossed the Transition Zone in the late Tertiary period, and (3) the high regional heat flow observed in the Transition Zone. The apparent presence of fluids near the base of the reflective zone may indicate a partially molten intrusion. We present a mechanism by which magma can be trapped and be induced to intrude horizontally at rheologic contrasts in extending crust.

Arizona↗
USGS1998

Suppression of large earthquakes by stress shadows: A comparison of Coulomb and rate-and-state failure

Stress shadows generated by California's two most recent great earthquakes (1857 Fort Tejon and 1906 San Francisco) substantially modified 19th and 20th century earthquake history in the Los Angeles basin and in the San Francisco Bay area. Simple Coulomb failure calculations, which assume that earthquakes can be modeled as static dislocations in an elastic half-space, have done quite well at approximating how long the stress shadows, or relaxing effects, should last and at predicting where subsequent large earthquakes will not occur. There has, however, been at least one apparent exception to the predictions of such simple models. The 1911 M >6.0 earthquake near Morgan Hill, California, occurred at a relaxed site on the Calaveras fault. We examine how the more complex rate-and-state friction formalism based on laboratory experiments might have allowed the 1911 earthquake. Rate-and-state time-to-failure calculations are consistent with the occurrence of the 1911 event just 5 years after 1906 if the Calaveras fault was already close to failure before the effects of 1906. We also examine the likelihood that the entire 78 years of relative quiet (only four M ≥6 earthquakes) in the bay area after 1906 is consistent with rate-and-state assumptions, given that the previous 7 decades produced 18 M ≥6 earthquakes. Combinations of rate-and-state variables can be found that are consistent with this pattern of large bay area earthquakes, assuming that the rate of earthquakes in the 7 decades before 1906 would have continued had 1906 not occurred. These results demonstrate that rate-and-state offers a consistent explanation for the 78-year quiescence and the 1911 anomaly, although they do not rule out several alternate explanations.

California↗
USGS2007

Paleomagnetic and mineral magnetic constraints on Zn-Pb ore genesis in the Pend Oreille Mine, Metaline district, Washington, USA

Zinc-lead mineralization in the Metaline mining district of northeastern Washington, USA, is hosted by the Cambrian Metaline Formation and is classified into Yellowhead-type (YO) and Josephine-type (JO) ore based on texture and mineralogy. Paleomagnetic results are reported for four Cambrian Metaline Formation sites, one Ordovician Ledbetter slate site, 12 YO and 13 JO (including two breccia sites) mineralization sites in the Pend Oreille Mine, and eight sites from the nearby Cretaceous Kaniksu granite batholith. Thermal and alternating field step demagnetization, saturation isothermal remanence analysis, and synthetic specimen tests show that the remanence in the host carbonates and Zn-Pb mineralization is carried mostly by pseudosingle (PSD) to single domain (SD) pyrrhotite and mostly by PSD to SD magnetite in the Kaniksu granite. Based on thermomagnetic measurements, sphalerite and galena concentrates and tailings from the mine's mill contain hexagonal and monoclimc pyrrhotite. The postfolding characteristic remanent magnetization (ChRM), known thermal data, and paleoarc method of dating suggest that the Zn-Pb mineralization carries a primary chemical remanent magnetization (CRM), and Metaline Formation carbonates a secondary CRM that were acquired during the Middle Jurassic (166 ??6 Ma) during the waning stages of the Nevadan orogeny. A paleomagnetic breccia test favours a solution-collapse origin for the Josephine breccia. Finally, the Kaniksu paleopole is concordant with the North American Cretaceous reference paleopole, suggesting the Kootenay terrane has not been rotated since emplacement of the batholith at ???94 Ma. ?? 2007 NRC Canada.

Canadian Journal of Earth Sciences↗
USGS2025

An automated compositing method for producing annual clear images from Landsat Collection 2 for annual NLCD production

Quality image input is fundamental to the quality of derived land cover products. Substantial time and effort are usually required to prepare images. Here, we present a novel and streamlined compositing algorithm that ingests Landsat Collection 2 Analysis Ready Data (ARD) and outputs cloud-free and gap-free composite imagery, which can be directly used for classification. This method leverages and improves the previous National Land Cover Database (NLCD) Virtual Median Value Point (VMVP) compositing method, the first part of the image preparation for NLCD 2019 operational production. The NLCD 2019 image preparation approach includes a second part, a residual cloud and cloud shadow detection and gap-filling method, to produce final cloud-free and gap-free composite imagery. The second part requires one clear reference image for each target year. Additional reference images are needed for producing reasonable observations for perennial ice/snow areas because Pixel QA (Quality Assessment) from ARD has difficulties differentiating ice/snow areas from clouds. Unlike the NLCD 2019 image preparation approach, our new compositing method, which is referred to as Automated VMVP (AVMVP), uses Landsat ARD as the only input and does not require reference images and extra steps. In this method, we developed new spectral filter criteria coupled with counts of clear observations using Pixel QA to identify potential cloud and cloud shadow observations on initially selected observations from the NLCD VMVP compositing algorithm. We also automate “gap-filling” using clear observations retrieved from a maximum of ±2 years around the target year when needed. Finally, a percentile-filtered compositing method was developed for the perennial ice/snow areas. All these steps are streamlined, pixel-based, and directly run on Landsat Collection 2 ARD. We have run successful tests on the conterminous United States (CONUS). Composite images derived from our innovative method were used to produce the CONUS Annual NLCD Collection 1 product suite that covers the period from 1985 to 2023.

conterminous United States↗
USGS2018

An updated method for estimating landslide‐event magnitude

Summary statistics derived from the frequency–area distribution (FAD) of inventories of triggered landslides allows for direct comparison of landslides triggered by one event (e.g. earthquake, rainstorm) with another. Such comparisons are vital to understand links between the landslide‐event and the environmental characteristics of the area affected. This could lead to methods for rapid estimation of landslide‐event magnitude, which in turn could lead to estimates of the total triggered landslide area. Previous studies proposed that the FAD of landslides follows an inverse power‐law, which provides the basis to model the size distribution of landslides and to estimate landslide‐event magnitude ( mLS ), which quantifies the severity of the event. In this study, we use a much larger collection of earthquake‐induced landslide (EQIL) inventories ( n =45) than previous studies to show that size distributions are much more variable than previously assumed. We present an updated model and propose a method for estimating mLS and its uncertainty that better fits the observations and is more reproducible, robust, and consistent than existing methods. We validate our model by computing mLS for all of the inventories in our dataset and comparing that with the total landslide areas of the inventories. We show that our method is able to estimate the total landslide area of the events in this larger inventory dataset more successfully than the existing methods.

Earth Surface Processes and Landforms↗
USGS2018

The influence of neap-spring tidal variation and wave energy on sediment flux in salt marsh tidal creeks

Sediment flux in marsh tidal creeks is commonly used to gage sediment supply to marshes. We conducted a field investigation of temporal variability in sediment flux in tidal creeks in the accreting tidal marsh at China Camp State Park adjacent to northern San Francisco Bay. Suspended-sediment concentration (SSC), velocity, and depth were measured near the mouths of two tidal creeks during three six-to-ten-week deployments: two in winter and one in summer. Currents, wave properties and SSC were measured in the adjacent shallows. All deployments spanned the largest spring tides of the season. Results show that tidally-averaged suspended-sediment flux (SSF) in the tidal creeks decreased with increasing tidal energy, and SSF was negative (bayward) for tidal cycles with maximum water surface elevation above the marsh plain. Export during the largest spring tides dominated the cumulative SSF measured during the deployments. During ebb tides following the highest tides, velocities exceeded 1 m/s in the narrow tidal creeks, resulting in negative tidally-averaged water flux, and mobilizing sediment from the creek banks or bed. Storm surge also produced negative SSF. Tidally-averaged SSF was positive in wavey conditions with moderate tides. Spring-tide sediment export was about 50% less at a station 130 m further up the tidal creek than at the creek mouth. The negative tidally-averaged water flux near the creek mouth during spring tides indicates that in the lower marsh, some of the water flooding directly across the bay--marsh interface drains through the tidal creeks, and suggests that this interface may be a pathway for sediment supply to the lower marsh as well.

California↗
USGS2012

Re–Os geochronology of the lacustrine Green River Formation: Insights into direct depositional dating of lacustrine successions, Re–Os systematics and paleocontinental weathering

Lacustrine sedimentary successions provide exceptionally high-resolution records of continental geological processes, responding to tectonic, climatic and magmatic influences. These successions are therefore essential for correlating geological and climatic phenomena across continents and furthermore the globe. Producing accurate geochronological frameworks within lacustrine strata is challenging because the stratigraphy is often bereft of biostratigraphy and directly dateable tuff horizons. The rhenium–osmium (Re–Os) geochronometer is a well-established tool for determining precise and accurate depositional ages of marine organic-rich rocks. Lake systems with stratified water columns are predisposed to the preservation of organic-rich rocks and thus should permit direct Re–Os geochronology of lacustrine strata. We present Re–Os systematics from one of the world's best documented lacustrine systems, the Eocene Green River Formation, providing accurate Re–Os depositional dates that are supported by Ar–Ar and U–Pb ages of intercalated tuff horizons. Precision of the Green River Formation Re–Os dates is controlled by the variation in initial 187 Os/ 188 Os and the range of 187 Re/ 188 Os ratios, as also documented in marine systems. Controls on uptake and fractionation of Re and Os are considered to relate mainly to depositional setting and the type of organic matter deposited, with the need to further understand the chelating precursors of Re and Os in organic matter highlighted. In addition to geochronology, the Re–Os data records the 187Os/188Os composition of lake water (1.41–1.54) at the time of deposition, giving an insight into continental runoff derived from weathering of the geological hinterland of the Green River Formation. Such insights enable us to evaluate fluctuations in continental climatic, tectonic and magmatic processes and provide the ability for chemostratigraphic correlation combined with direct depositional dates. Furthermore, initial 187Os/188Os values can be used as a diagnostic tool to distinguish between lacustrine and marine depositional settings when compared to known oceanic 187 Os/ 188 Os values.

Colorado;Utah;Wyoming↗
USGS1989

Style of extensional tectonism during rifting, Red Sea and Gulf of Aden

Models describing the development of the Red Sea and the Gulf of Aden, prior to the present periods of sea-floor spreading, include those that use block faulting on steep normal faults, uniform diffuse shear in continental crust, simple shear on large detachment faults that cut the entire lithosphere, combinations involving detachment faults/ductile deformation/plutonic inflation, and ones that minimize the role of mechanical extension in favor of an earlier stage of sea-floor spreading. Geologic and geophysical studies from the Arabian continental margin in the southern Red Sea and LANDSAT analysis of the northern Somalia margin in the Gulf of Aden suggest that the early continental rifts were long narrow features that formed by extension on closely spaced normal faults above moderate- to shallow-dipping detachments with break-away zones defining one rift flank and root zones under the opposing rift flank. The rift flanks presently form the opposing continental margins across each ocean basin. The detachment on the Arabian margin dips gently to the west, with a breakaway zone now eroded above the deeply dissected terrain of the Arabian escarpment. The Arabian detachment projects westward to middle crustal levels beneath the sediment of the southern Red Sea coastal plain. Strata in the upper plate dip as steeply as 60° to the west, and the beds are repeated by numerous planar and listric normal faults that dip to the east. Most of the faults truncate downward at the detachment. Thus, the upper plate is highly extended and the rocks in its eastern part have been translated about 20 km westward and 21/2- to 5-km downward relative to the rest of Arabia. A prominent detachment surface, with a north dip, is evident in northernmost Somalia where it breaks away north of the Somalian escarpment in an otherwise undeformed section of cratonic strata of Jurassic to Eocene age. The upper plate of the Somalian detachment consists of a highly faulted collage of the cratonic strata. This fault projects to middle crustal levels in the opposing Arabian margin to the northeast.

Journal of African Earth Sciences↗
USGS2009

Diverse lavas from closely spaced volcanoes drawing from a common parent: Emmons Lake Volcanic Center, Eastern Aleutian Arc

Emmons Lake Volcanic Center (ELVC) on the lower Alaskan Peninsula is one of the largest and most diverse volcanic centers in the Aleutian Arc. Since the Middle Pleistocene, eruption of ~ 350 km 3 of basalt through rhyolite has produced a 30 km, arc front chain of nested calderas and overlapping stratovolcanoes. ELVC has experienced as many as five major caldera-forming eruptions, the most recent, at ~ 27 ka, produced ~ 50 km 3 of rhyolitic ignimbrite and ash fall. These violent silicic events were interspersed with less energetic, but prodigious, outpourings of basalt through dacite. Holocene eruptions are mostly basaltic andesite to andesite and historically recorded activity includes over 40 eruptions within the last 200 yr, all from Pavlof volcano, the most active site in the Aleutian Arc. Geochemical and geophysical observations suggest that although all ELVC eruptions derive from a common clinopyroxene + spinel + plagioclase fractionating high-aluminum basalt parent in the lower crust, magma follows one of two closely spaced, but distinct paths to the surface. Under the eastern end of the chain, magma moves rapidly and cleanly through a relatively young (~ 28 ka), hydraulically connected dike plexus. Steady supply, short magma residence times, and limited interaction with crustal rocks preserve the geochemistry of deep crustal processes. Below the western part of the chain, magma moves haltingly through a long-lived (~ 500 ka) and complex intrusive column in which many generations of basaltic to andesitic melts have mingled and fractionated. Buoyant, silicic melts periodically separate from the lower parts of the column to feed voluminous eruptions of dacite and rhyolite. Mafic lavas record a complicated passage through cumulate zones and hydrous silicic residues as manifested by disequilibrium phenocryst textures, incompatible element enrichments, and decoupling of REEs and HFSEs ratios. Such features are absent in mafic lavas from the younger part of the chain, highlighting the importance of plumbing architecture and longevity in creating petrologic diversity. Supplemental Data include 156 major element (XRF) and 128 trace element (ICP-MS) whole-rock analyses, 23 new 40 Ar/ 39 Ar ages, a generalized geologic map with associated unit descriptions and field photographs, and photomicrographs of key petrographic features.

Alaska↗
USGS2003

Geochronology and eruptive history of the Katmai volcanic cluster, Alaska Peninsula

In the Katmai district of the Alaska Peninsula, K–Ar and 40 Ar/ 39 Ar ages have been determined for a dozen andesite–dacite stratocones on the arc front and for 11 rear-arc volcanoes, 10 of which are monogenetic. Tied to mapping and stratigraphic studies, our dating emphasized proximal basal lavas that rest on basement rocks, in order to estimate ages of inception of each polygenetic cone. Oldest among arc-front cones is Alagogshak Volcano (690–43 ka), succeeded in the Holocene by the active Mount Martin cone. Mount Mageik consists of four overlapping subedifices, basal lavas of which give ages of 93, 71, and 59 ka, and Holocene. The three small prehistoric cones of Trident Volcano yield ages of 143, 101–58, and 44 ka. Falling Mountain and Mount Cerberus, dacite domes near the 1912 Novarupta vent, are related compositionally to the Trident group and give ages of 70 ka and 114 ka. Mount Katmai, which underwent caldera collapse in 1912, consists of two subedifices that overlapped in space and time, and is the only arc-front center here to include basalt and rhyolite; one cone began by 90 ka, the other by 47 ka. Snowy Mountain also consists of two contiguous cones, which started around 200 and 171 ka, respectively, the younger remaining active into the Holocene. Devils Desk, the only mafic cone on the arc front, was short-lived at about 245 ka. In the rear-arc, (1) Mount Griggs produced mafic-to-silicic andesite in several episodes between 292 ka and the Holocene; (2) the Savonoski River cluster includes a Pliocene dacite dome and five small mafic cones (390–88 ka); (3) Gertrude Creek cone (49.8% SiO 2 ) yields an age of 500 ka; and (4) the Saddlehorn Creek cluster includes five Pliocene basalt-to-andesite remnants. Eruptive volumes were reconstructed, permitting estimates of average eruption rates for edifice lifetimes. Since the mid Pleistocene, total volume erupted along the arc front here is 210±47 km 3 and in the rear-arc 39±6 km 3 , of which Mount Griggs alone accounts for 35±5 km 3 . Most productive has been Mount Katmai at 70±18 km 3 , yielding a rate of ∼1 km 3 /kyr, followed by Mount Mageik (0.33 km 3 /kyr) and Mount Griggs (0.3 km 3 /kyr since 50 ka).

Alaska↗
USGS1988

The surface integral approach to Radarclinometry

Because radarclinometry is fundamentally describable in terms of a nonlinear, first-order, partial differential equation, one expects that it can, in principle, be carried out by direct deterministic integration beginning at a given threshold profile along the azimuthal coordinate. Such a boundary condition could be provided by the altimetry profile obtained on a preceding or succeeding orbital revolution of the radar-bearing spacecraft. Notwithstanding the mismatched resolutions of the radar altimeter and the radar imaging system as planned for the Megallan mission to Venus, there are fundamental considerations, not involving system noise, that influence the possibility of success of this approach. From the topographic map of the Lake Champlain West quadrangle in the Adirondack Mountains of the U.S., a radar image is synthesized. Radarclinometry, in surface integral form, recaptures the topographic map when the applicable radar reflectance function is weakly variable over the range of application, but it diverges beyond a certain point for nominally variable reflectance functions. The effect can be understood by using results from the "shape-from-shading" literature. (This literature is produced by a group within the artificial intelligence community who have been independently attacking, for all practical purposes, photoclinometry, except that they have not given primacy to images of terrain.) The ubiquity of the instability suggests that the value of the surface integral approach is much in doubt. ?? 1988 Kluwer Academic Publishers.

Earth, Moon and Planets↗
USGS2014

Source and progression of a submarine landslide and tsunami: The 1964 Great Alaska earthquake at Valdez

Like many subduction zone earthquakes, the deadliest aspects of the 1964 M = 9.2 Alaska earthquake were the tsunamis it caused. The worst of these were generated by local submarine landslides induced by the earthquake. These caused high runups, engulfing several coastal towns in Prince William Sound. In this paper, we study one of these cases in detail, the Port Valdez submarine landslide and tsunami. We combine eyewitness reports, preserved film, and careful posttsunami surveys with new geophysical data to inform numerical models for landslide tsunami generation. We review the series of events as recorded at Valdez old town and then determine the corresponding subsurface events that led to the tsunami. We build digital elevation models of part of the pretsunami and posttsunami fjord‐head delta. Comparing them reveals a ~1500 m long region that receded 150 m to the east, which we interpret as the primary delta landslide source. Multibeam imagery and high‐resolution seismic reflection data identify a ~400 m wide chute with hummocky deposits at its terminus, which may define the primary slide path. Using these elements we run hydrodynamic models of the landslide‐driven tsunamis that match observations of current direction, maximum inundation, and wave height at Valdez old town. We speculate that failure conditions at the delta front may have been influenced by manmade changes in drainage patterns as well as the fast retreat of Valdez and other glaciers during the past century.

Alaska↗
USGS2020

Seismic and geodetic progression of the 2018 summit caldera collapse of Kīlauea Volcano

The 2018 eruption of Kīlauea volcano, Hawaiʻi, resulted in a major collapse of the summit caldera along with an effusive eruption in the lower East Rift Zone. The caldera collapse comprised 62 highly similar collapse cycles of strong ground deformation and earthquake swarms that ended with a magnitude 5 collapse event and one partial cycle that did not end with a collapse event. We analyzed geodetic and seismic data to better understand how the caldera collapse progressed over 3 months of activity, focusing on the cyclical activity. We identified 3 main phases of collapse: initial ring-fault activation and small explosions (Phase 1), an eastward shift in activity and freeing of the central piston (Phase 2), and a recoupling of the piston to the reservoir followed by relatively steady behavior until the eruption’s end (Phase 3). Additionally, we observed geodetic evidence of tangential motion from the localization of the main ring fault (Phase 2) and the formation of a major peripheral ring fault on the eastern side of the collapse caldera during Phase 3. Both geodetic and seismic parameters suggest that the collapse may have had an eastward-component of motion after the ring fault system had formed. The cyclical seismic and geodetic parameters show no obvious signs that the collapse was coming to an end, with the only notable change being a significant increase in the ratio of cyclical displacement to co-collapse displacement observed during the last complete cycle on GNSS stations outside the caldera region.

Hawaii↗
USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Journal of Geophysical Research B: Solid Earth↗