Search NASA⌕ Search

Engineering topics

Jeffrey T. Freymueller

Publications and source records attributed to Jeffrey T. Freymueller.

Quantifying mass flows at Mt. Cleveland, Alaska between 2001 and 2020 using satellite photogrammetry

Measuring eruption volume provides constraints on the magma supply rate and plumbing systems and therefore is a critical component for monitoring volcanoes. We use ArcticDEM—a large collection of time-dependent digital elevation models (DEMs) derived from satellite stereo-photogrammetry — to construct a first-of-its-kind measurement of the volumes of recent mass flows at Mount Cleveland, Alaska. We quantify the volume of the products of the 2001 eruption (the largest in the past half-century) as (54.8 ± 0.5) × 10 -6 m -3 covering a total area of 5.09 km -2 . The total volume of material loss at the summit crater is (0.67 ± 0.02) × 10 -6 m -3 , which is likely caused by later explosions and collapses of the shallow magma chamber. The total eruptive volume between 2017 and 2020 is (0.086 ± 0.002) × 10 -6 m -3 . Elevation changes associated with lahars are variable. On the upper northern slopes of the volcano, the lahar channels were almost exclusively erosive, suggesting that lahars originating at the summit eroded and entrained loose materials high on the cone. In general, lahar deposits were thickest near their toes, with some reaching ~20 m thickness.

Chunli Dai↗

Cascading Rupture of a Megathrust

Understanding variability in the size and location of large earthquakes along subduction margins is crucial for evaluating seismic and tsunami hazards. We present a coseismic slip model for the 2021 M8.2 Chignik earthquake and investigate the relationship of this earthquake to previous major events in the Alaska Peninsula region and to interseismic coupling. Stress changes from the 2020 M7.8 Simeonof event triggered the Chignik event, and together, the earthquakes partially filled an unruptured section along a 3000-km subduction margin that has experienced a series of ruptures along almost its entire length over the past century. Variations in coupling and structural characteristics may make the region more prone to nucleating M7 to M8 events rather than larger M > 8.5 earthquakes. Stress changes and rupture areas suggest that the two recent earthquakes may be part of an 80-year-long rupture cascade and may have advanced seismic hazard in the region.

Julie L. Elliott↗

Stress Promotion of the 1958 Mw∼7.8 Fairweather Fault Earthquake and Others in Southeast Alaska by Glacial Isostatic Adjustment and Inter-earthquake Stress Transfer

We assess how recent glacial isostatic adjustment (GIA) and inter-earthquake stress transfer modulated the state of stress on major faults in Southeast Alaska, and how these effects may have influenced recent large and moderate earthquakes. The Glacier Bay Icefield has lost >3,000 km3 of ice since ca. 1770, with ice thinning up to 1.5 km, and the rest of southeast Alaska has been deglaciating since ca. 1900. The resulting GIA response in the solid earth produces present-day surface uplift rates of up to ∼4 cm/yr, among the fastest on Earth. This rapid deformation occurs directly alongside and atop the Fairweather Fault, which accommodates 4–5 cm/yr of right-lateral motion through slip in frequent large earthquakes, recently including a Mw∼7.8 earthquake in 1958. We find that the 1958 earthquake nucleated very close to the location of maximum GIA-induced Coulomb stress increase on the Fairweather Fault between 1770 and 1958 (estimated at ∼0.2–0.6 MPa), on the only section of fault where GIA-induced stress changes had involved both right-lateral shearing and unclamping. The September 10, 1899 Mw∼8.1 Yakutat Bay earthquake also promoted failure on at least the northwest part of the eventual 1958 rupture. In turn, the 1899 and 1958 earthquakes and GIA promoted failure along the St. Elias compressional margin, the site of a Mw∼7.4 event in 1979. While rapid tectonic loading is very likely the main driver of earthquakes in southeast Alaska, we estimate that 23 of 30 instrumentally constrained Mw ≥ 5.0 earthquakes in this region were also promoted by post-1770 GIA.

Chris Rollins↗

Characterization of the 2008 Phreatomagmatic Eruption of Okmok From ArcticDEM and InSAR: Deposition, Erosion, and Deformation

Okmok volcano, one of the most active volcanoes in the Aleutian Islands, erupted between 12 July and 19 August 2008 and produced a new tephra cone (Ahmanilix) up to 253 m high. A coeruptive deposit thickness map and the posteruptive elevation change rate map are generated for this phreatomagmatic eruption using the new, high-resolution (2 m) digital elevation model (DEM) time series provided by ArcticDEM. The total dense rock equivalent (DRE) volume over the proximal deposit field classified from the elevation change map is estimated as 0.147 ± 0.001 km 3 DRE (bulk volume of 0.453 ± 0.003 km 3 ). Combined with the distal deposit volume derived from field measurements, the total DRE is 0.25 km 3 (bulk volume of 0.76 km 3 ). The estimated erosion rate is up to −15 ± 3 m/year along the flanks of Ahmanilix, and the redeposition rate is about 1 ± 0.3 m/year at the base and 6 ± 1 m/year within the center of Ahmanilix. The surface elevation change also reveals a large ground surface collapse that disrupted the old Cone D lava bench, with a total collapse area of 0.15 km2 and a volume of (7.1 ± 0.4) × 10−3 km 3 . We also produce a coeruptive deformation map from satellite interferometric synthetic aperture radar (InSAR) using topography from both ArcticDEM and the Shuttle Radar Topography Mission (SRTM) and compare the results. Finally, based on historical eruption volumes, the magma supply rate at Okmok is estimated as (1.77 ± 0.1) × 10−3 km 3 /year.

Chunli Dai↗