Geology Reports⌕ Search

Geology topics

Peter J. Haeussler

Publications and source records attributed to Peter J. Haeussler.

At least 109 records · Page 6Linked to original sources

Lidar-revised geologic map of the Wildcat Lake 7.5' quadrangle, Kitsap and Mason Counties, Washington

This map is an interpretation of a 6-ft-resolution (2-m-resolution) lidar (light detection and ranging) digital elevation model combined with the geology depicted on the Geologic Map of the Wildcat Lake 7.5' quadrangle, Kitsap and Mason Counties, Washington (Haeussler and Clark, 2000). Haeussler and Clark described, interpreted, and located the geology on the 1:24,000-scale topographic map of the Wildcat Lake 7.5' quadrangle. This map, derived from 1951 aerial photographs, has 20-ft contours, nominal horizontal resolution of approximately 40 ft (12 m), and nominal mean vertical accuracy of approximately 10 ft (3 m). Similar to many geologic maps, much of the geology in the Haeussler and Clark (2000) map-especially the distribution of surficial deposits-was interpreted from landforms portrayed on the topographic map. In 2001, the Puget Sound lidar Consortium obtained a lidar-derived digital elevation model (DEM) for Kitsap Peninsula including all of the Wildcat Lake 7.5' quadrangle. This new DEM has a horizontal resolution of 6 ft (2 m) and a mean vertical accuracy of about 1 ft (0.3 m). The greater resolution and accuracy of the lidar DEM compared to topography constructed from air photo stereo models have much improved the interpretation of geology in this heavily vegetated landscape, especially the distribution and relative age of some surficial deposits. Many contacts of surficial deposits are adapted unmodified or slightly modified from Haugerud (2009).

Washington↗

Combined effects of tectonic and landslide-generated Tsunami Runup at Seward, Alaska during the Mw 9.2 1964 earthquake

We apply a recently developed and validated numerical model of tsunami propagation and runup to study the inundation of Resurrection Bay and the town of Seward by the 1964 Alaska tsunami. Seward was hit by both tectonic and landslide-generated tsunami waves during the M W "> M W MW 9.2 1964 megathrust earthquake. The earthquake triggered a series of submarine mass failures around the fjord, which resulted in landsliding of part of the coastline into the water, along with the loss of the port facilities. These submarine mass failures generated local waves in the bay within 5 min of the beginning of strong ground motion. Recent studies estimate the total volume of underwater slide material that moved in Resurrection Bay to be about 211 million m 3 (Haeussler et al. in Submarine mass movements and their consequences, pp 269–278, 2007 ). The first tectonic tsunami wave arrived in Resurrection Bay about 30 min after the main shock and was about the same height as the local landslide-generated waves. Our previous numerical study, which focused only on the local landslide-generated waves in Resurrection Bay, demonstrated that they were produced by a number of different slope failures, and estimated relative contributions of different submarine slide complexes into tsunami amplitudes (Suleimani et al. in Pure Appl Geophys 166:131–152, 2009 ). This work extends the previous study by calculating tsunami inundation in Resurrection Bay caused by the combined impact of landslide-generated waves and the tectonic tsunami, and comparing the composite inundation area with observations. To simulate landslide tsunami runup in Seward, we use a viscous slide model of Jiang and LeBlond (J Phys Oceanogr 24(3):559–572, 1994 ) coupled with nonlinear shallow water equations. The input data set includes a high resolution multibeam bathymetry and LIDAR topography grid of Resurrection Bay, and an initial thickness of slide material based on pre- and post-earthquake bathymetry difference maps. For simulation of tectonic tsunami runup, we derive the 1964 coseismic deformations from detailed slip distribution in the rupture area, and use them as an initial condition for propagation of the tectonic tsunami. The numerical model employs nonlinear shallow water equations formulated for depth-averaged water fluxes, and calculates a temporal position of the shoreline using a free-surface moving boundary algorithm. We find that the calculated tsunami runup in Seward caused first by local submarine landslide-generated waves, and later by a tectonic tsunami, is in good agreement with observations of the inundation zone. The analysis of inundation caused by two different tsunami sources improves our understanding of their relative contributions, and supports tsunami risk mitigation in south-central Alaska. The record of the 1964 earthquake, tsunami, and submarine landslides, combined with the high-resolution topography and bathymetry of Resurrection Bay make it an ideal location for studying tectonic tsunamis in coastal regions susceptible to underwater landslides.

Alaska↗

Structural evolution of an arc-basin: The Gravina Belt in central southeastern Alaska

The upper Middle Jurassic-Lower Cretaceous Gravina belt lies along the eastern margin of the Alexander terrane in southeastern Alaska. This group of turbidites and mafic to intermediate volcanic rocks was deformed during mid to Late Cretaceous time during the closing of a basin of unknown size between the Alexander terrane on the west and the Stikine terrane to the east. Therefore structures of Gravina belt rocks largely reflect the final accretion and subsequent transport of the Alexander terrane. Six geologic transects across the central Gravina belt (southern Mitkof Island to northern Douglas Island) define a structural history that includes (1) syndepositional soft-sediment folding and faulting, possibly in conjunction with block tilting and extension; (2) tight to isoclinal folding or thrust faulting, with formation of a slaty cleavage (S 1 ) striking ∼330°; (3) local coaxial refolding with formation of crenulation cleavage (S 2 ); (4) development of domainal crenulation folds and cleavage (S 3 ) oriented at a large angle to the margin of the belt; (5) intrusion of tonalitic plutons around 90 Ma; and (6) right-lateral strike-slip displacement on faults oriented ∼330°. Finite strain measurements on sedimentary rocks suggest the belt was at least twice its present width, normal to the foliation, before deformation. Subhorizontal margin-parallel fold axes, margin-parallel slaty cleavage, and margin-perpendicular stretching lineations suggest orthogonal contraction of the Gravina basin, assuming that oblique plate convergence was not decoupled along strike-slip faults. After contractional deformation, strike-slip faults indicate dextral displacement (probably of the order of several tens of kilometers) of the Alexander terrane with respect to the terranes to the east. Domainal crenulation folds and cleavages at a high angle to the margin of the belt suggest that the change in convergence directions occurred while the rocks (presently at the surface) could still plastically deform. This kinematic interpretation of structures is consistent with changes in plate motions [Engebretson et al., 1985]. Before 100 Ma, the convergence directions between the Kula and North America plates were at a high angle to the continental margin, whereas after 100 Ma, convergence directions were at a small angle to the continental margin. In addition, after 100 Ma, the Kula-North America, and not the Farallon-North America, convergence direction is most compatible with a N-S principal paleostress orientation derived from inversion of strike-slip fault data. This relationship suggests that it may have been the Kula plate that drove northward transport of the Alexander terrane along the margin of North America.

Alaska↗

Historic and paleo-submarine landslide deposits imaged beneath Port Valdez, Alaska: Implications for tsunami generation in a glacial fiord

During the 1964 M9.2 great Alaskan earthquake, submarine-slope failures resulted in the generation of highly destructive tsunamis at Port Valdez, Alaska. A high-resolution, mini-sparker reflection profiler was used to image debris lobes, which we attribute to slope failures that occurred both during and prior to the 1964 megathrust event. In these reflection profiles, debris lobe deposits are indicated by acoustically opaque units that are separated by undisturbed parallel-layered reflectors. Near-surface debris lobes attributed to the 1964 earthquake include: (1) a debris lobe over 30 m thick that emanates from the fiord-head delta in eastern Port Valdez; and (2) debris flow lobes incorporating large, intact blocks up to 40 m high in western Port Valdez, off the Shoup Glacier moraine. In addition to the near-surface debris lobes, we imaged at least five additional debris lobe deposits buried beneath the 1964 deposit. The debris lobe directly beneath the 1964 deposit has a similar thickness and spatial distribution as the 1964 deposit. However, the older, deeper, debris lobes are thinner, less extensive, and separated by thinner sequences of parallel-layered reflectors. Glacier retreat and concomitant build-up of the fiord-head delta combined with longer time intervals between megathrust events may have resulted in more extensive delta failures and thus thicker debris lobes through time.

Alaska↗

Numerical study of tsunami generated by multiple submarine slope failures in Resurrection Bay, Alaska, during the M W 9.2 1964 earthquake

We use a viscous slide model of Jiang and LeBlond (1994) coupled with nonlinear shallow water equations to study tsunami waves in Resurrection Bay, in south-central Alaska. The town of Seward, located at the head of Resurrection Bay, was hit hard by both tectonic and local landslide-generated tsunami waves during the M W 9.2 1964 earthquake with an epicenter located about 150 km northeast of Seward. Recent studies have estimated the total volume of underwater slide material that moved in Resurrection Bay during the earthquake to be about 211 million m 3 . Resurrection Bay is a glacial fjord with large tidal ranges and sediments accumulating on steep underwater slopes at a high rate. Also, it is located in a seismically active region above the Aleutian megathrust. All these factors make the town vulnerable to locally generated waves produced by underwater slope failures. Therefore it is crucial to assess the tsunami hazard related to local landslide-generated tsunamis in Resurrection Bay in order to conduct comprehensive tsunami inundation mapping at Seward. We use numerical modeling to recreate the landslides and tsunami waves of the 1964 earthquake to test the hypothesis that the local tsunami in Resurrection Bay has been produced by a number of different slope failures. We find that numerical results are in good agreement with the observational data, and the model could be employed to evaluate landslide tsunami hazard in Alaska fjords for the purposes of tsunami hazard mitigation. ?? Birkh??user Verlag, Basel 2009.

Conference Paper↗

The role of ridge subduction in determining the geochemistry and Nd–Sr–Pb isotopic evolution of the Kodiak batholith in southern Alaska

The Paleocene Kodiak batholith, part of the Sanak&ndash;Baranof belt of Tertiary near-trench intrusive rocks, forms an elongate body (~ 150 km long) that transects Kodiak Island from SW to NE. The batholith consists of three zones (Southern, Central, and Northern) of kyanite-, muscovite-, and garnet-bearing biotite tonalite and granodiorite and less abundant granite that intruded an accretionary prism (Kodiak Formation, and Ghost Rocks Formation). Small and likely coeval bodies (Northern, Western, and Eastern satellite groups) of quartz gabbro, diorite, tonalite, granodiorite, and leucogranite flank the batholith. The batholith is calc-alkalic, has an aluminum saturation index of > 1.1, FeO t /(FeO t + MgO) ~ 0.65 (at SiO 2 = 65 wt.%), and increases in SiO 2 (~ 61 wt.%&ndash;73 wt.%) and decreases in TiO 2 (~ 0.9 wt.%&ndash;0.3 wt.%) from SW to NE. As a group, the granitic rocks have light REE-enriched chondrite-normalized patterns with small or no negative Eu anomalies, primitive mantle-normalized negative anomalies for Nb and Ti, and positive anomalies for Pb. Small to large negative anomalies for Th are also distinctive. The quartz gabbros and diorites are generally characterized by generally flat to light REE chondrite-normalized patterns (no Eu anomalies), and mantle-normalized negative anomalies for Nb, Ti, and P. Pb isotopic compositions ( 206 Pb/ 204 Pb = 18.850&ndash;18.960; 207 Pb/ 204 Pb = 15.575&ndash;15.694; 208 Pb/ 204 Pb = 38.350&ndash;39.039) are intermediate between depleted mantle and average continental crust. The Southern zone and a portion of the Central zone are characterized by negative &epsilon; Nd values of &minus; 3.7 to &minus; 0.3 and T DM ages ranging from ~ 838 Ma to 1011 Ma. Other granitic rocks from the Central and Northern zones have higher &epsilon; Nd values of &minus; 0.4 to + 4.7 and younger T DM ages of ~ 450 to 797 Ma. Granitic and mafic plutons from the Eastern satellites show a wide range of &epsilon; Nd values of &minus; 2.7 to + 6.4, and T DM ages from 204 Ma to 2124 Ma. 87 Sr/ 86 Sr values of the Southern and Central zones overlap and tend to be slightly more radiogenic ( 87 Sr/ 86 Sr > 0.70426) than the Northern zone ( 87 Sr/ 86 Sr < 0.70472). 206 Pb/ 204 Pb values increase slightly from the Southern and Central zones toward the Northern zone. There is no clear correlation of the major or trace elements with &epsilon; Nd , Pb or Sr isotopic values. Kodiak Formation and the Ghost Rocks Formation overlap the isotopic compositions (e.g., 206 Pb/ 204 Pb = 18.978 to 19.165, 87 Sr/ 86 Sr of 0.705715 to 0.707118, and &epsilon; Nd of &minus; 6.7 to &minus; 1.5 at 59 Ma) and T DM values (959 to 1489 Ma) of the batholith. Production of large volumes of granitic rocks in the Sanak&ndash;Baranof belt, and particularly on Kodiak Island, reflects a sequence of processes that includes underplating of mantle-derived mafic (possibly from the mantle wedge) and intermediate rocks under the accretionary flysch, interlayering of mantle-derived and flyschoid rocks, and partial melting of the mixed lithologic assemblages. Limited degrees of fractional crystallization or assimilation and fractional crystallization influenced compositions of the granitic rocks. The contribution of mantle-derived rocks that resided in the accretionary prism for only a short period of time prior to partial melting likely exceeds 40% (up to 80%). The balance (60 to 20%) is from a recently recycled crustal component represented by the Kodiak Formation. This type of progressive intracrustal melting from mixed sources controlled the geochemical character of the batholith and is most consistent with the hypothesis that the granitic rocks are associated with a slab-window produced by collision of a spreading oceanic center and a subduction zone and migration beneath the accretionary prism.

Gulf of Alaska↗

Surface rupture map of the 2002 M7.9 Denali fault earthquake, Alaska: Digital data

The November 3, 2002, Mw7.9 Denali Fault earthquake produced about 340 km of surface rupture along the Susitna Glacier Thrust Fault and the right-lateral, strike-slip Denali and Totschunda Faults. Digital photogrammetric methods were primarily used to create a 1:500-scale, three-dimensional surface rupture map, and 1:6,000-scale aerial photographs were used for three-dimensional digitization in ESRI's ArcMap GIS software, using Leica's StereoAnalyst plug in. Points were digitized 4.3 m apart, on average, for the entire surface rupture. Earthquake-induced landslides, sackungen, and unruptured Holocene fault scarps on the eastern Denali Fault were also digitized where they lay within the limits of air photo coverage. This digital three-dimensional fault-trace map is superior to traditional maps in terms of relative and absolute accuracy, completeness, and detail and is used as a basis for three-dimensional visualization. Field work complements the air photo observations in locations of dense vegetation, on bedrock, or in areas where the surface trace is weakly developed. Seventeen km of the fault trace, which broke through glacier ice, were not digitized in detail due to time constraints, and air photos missed another 10 km of fault rupture through the upper Black Rapids Glacier, so that was not mapped in detail either.

Alaska↗

Studies by the U.S. Geological Survey in Alaska, 2007

The collection of papers that follow continues the series of U.S. Geological Survey (USGS) investigative reports in Alaska under the broad umbrella of the geologic sciences. This series represents new and sometimes-preliminary findings that are of interest to Earth scientists in academia, government, and industry; to land and resource managers; and to the general public. The reports presented in Studies by the U.S. Geological Survey in Alaska cover a broad spectrum of topics from various parts of the State, serving to emphasize the diversity of USGS efforts to meet the Nation's needs for Earth-science information in Alaska. This professional paper is one of a series of 'online only' versions of Studies by the U.S. Geological Survey in Alaska, reflecting the current trend toward disseminating research results on the World Wide Web with rapid posting of completed reports.

Professional Paper↗

Quaternary uplift history of Wingham Island, south-central Alaska

Marine terraces cut into Pleistocene deposits on Wingham Island in the Gulf of Alaska provide new constraints on the position of sea level, ice thickness and total glacioisostatic rebound at the end of the Last Glacial Maximum. A radiocarbon age of 13.9 +- 0.15 ka on the most prominent terrace is coincident with the end of meltwater pulse 1A, possibly suggestive of a link between changes in relative sea level and terrace formation. Isostatic modeling suggests a local ice thickness of 600 to 700 m with high (~10 cm/yr) initial rates of postglacial rebound. In addition to the unique ties to meltwater pulse 1A, the timing of emergence for Wingham Island following the Last Glacial Maximum has implications for the early migration of humans into North America.

Alaska↗

Detrital zircon geochronology of Cretaceous and Paleogene strata across the south-central Alaskan convergent margin

Ages of detrital zircons are reported from ten samples of Lower Cretaceous to Paleogene metasandstones and sandstones from the Chugach Mountains, Talkeetna Mountains, and western Alaska Range of south-central Alaska. Zircon ages are also reported from three igneous clasts from two conglomerates. The results bear on the regional geology, stratigraphy, tectonics, and mineral resource potential of the southern Alaska convergent margin. Chugach Mountains - The first detrital zircon data are reported here from the two main components of the Chugach accretionary complex - the inboard McHugh Complex and the outboard Valdez Group. Detrital zircons from sandstone and two conglomerate clasts of diorite were dated from the McHugh Complex near Anchorage. This now stands as the youngest known part of the McHugh Complex, with an inferred Turonian (Late Cretaceous) depositional age no older than 91-93 Ma. The zircon population has probability density peaks at 93 and 104 Ma and a smattering of Early Cretaceous and Jurassic grains, with nothing older than 191 Ma. The two diorite clasts yielded Jurassic U-Pb zircon ages of 179 and 181 Ma. Together, these findings suggest a Mesozoic arc as primary zircon source, the closest and most likely candidate being the Wrangellia composite terrane. The detrital zircon sample from the Valdez Group contains zircons as young as 69 and 77 Ma, consistent with the previously assigned Maastrichtian to Campanian (Late Cretaceous) depositional age. The zircon population has peaks at 78, 91, 148, and 163 Ma, minor peaks at 129, 177, 330, and 352 Ma, and no concordant zircons older than Devonian. A granite clast from a Valdez Group conglomerate yielded a Triassic U-Pb zircon age of 221 Ma. Like the McHugh Complex, the Valdez Group appears to have been derived almost entirely from Mesozoic arc sources, but a few Precambrian zircons are also present. Talkeetna Mountains - Detrital zircons ages were obtained from southernmost metasedimentary rocks of the Talkeetna Mountains (schist of Hatcher Pass) and, immediately to the south, the northernmost sedimentary sequence of the Matanuska forearc basin (Arkose Ridge Formation). Detrital zircons from the Paleogene Arkose Ridge Formation are as young as 61 and 70 Ma; the population is dominated by a single Late Cretaceous peak at 76 Ma; the oldest zircon is 181 Ma. Sedimentological evidence clearly shows that the conglomeratic Arkose Ridge Formation was derived from the Talkeetna Mountains; our detrital zircon data support this inference. Zircons dated at ca. 90 Ma in the Arkose Ridge sample suggest that buried or unmapped plutons of this age may exist in the Talkeetnas. This is a particularly interesting age as it corresponds to the age of the supergiant Pebble gold-molybdenum-copper porphyry prospect near Iliamna and suggests a new area of prospectivity for Pebble-type deposits. The schist of Hatcher Pass, which was previously assigned a Jurassic depositional age, yielded surprisingly young Late Cretaceous detrital zircons, the youngest at 75 Ma. The probability density curve has four Cretaceous peaks from 76 to 102 Ma, a pair of Late Jurassic peaks at 155 and 166 Ma, three Early Jurassic to Late Triassic peaks at 186, 197, and 213 Ma, minor Carboniferous peaks at 303 and 346 Ma, and a minor Paleoproterozoic peak at 1828 Ma. The schist of Hatcher Pass was largely derived from Mesozoic arc sources, most likely the Wrangellia composite terrane, with some contribution from one or more older, inboard sources, probably including the Yukon-Tanana terrane. We postulate that the schist of Hatcher Pass represents metamorphosed rocks of the Valdez Group that were subducted and then exhumed along the Chugach terrane's 'backstop' during Paleogene transtension. Western Alaska Range - Six detrital zircon samples were collected from a little studied belt of turbidites in Tyonek quadrangle on strike with the Kahiltna assemblage of the central Alaska Range.

Alaska↗

Combined high-resolution LIDAR topography and multibeam bathymetry for northern Resurrection Bay, Seward, Alaska

A new Digital Elevation Model was created using the best available high-resolution topography and multibeam bathymetry surrounding the area of Seward, Alaska. Datasets of (1) LIDAR topography collected for the Kenai Watershed Forum, (2) Seward harbor soundings from the U.S. Army Corp of Engineers, and (3) multibeam bathymetry from the National Oceanic and Atmospheric Administration contributed to the final combined product. These datasets were placed into a common coordinate system, horizontal datum, vertical datum, and data format prior to being combined. The projected coordinate system of Universal Transverse Mercator Zone 6 North American Datum of 1927 was used for the horizontal coordinates. Z-values in meters were referenced to the tidal datum of Mean High Water. Gaps between the datasets were interpolated to create the final seamless 5-meter grid covering the area of interest around Seward, Alaska.

Alaska↗

Submarine slope failures near Seward, Alaska, during the M9.2 1964 earthquake

Following the 1964 M9.2 megathrust earthquake in southern Alaska, Seward was the only town hit by tsunamis generated from both submarine landslides and tectonic sources. Within 45 seconds of the start of the earthquake, a 1.2-km-long section of waterfront began sliding seaward, and soon after, ~6-8-m high waves inundated the town. Studies soon after the earthquake concluded that submarine landslides along the Seward waterfront generated the tsunamis that occurred immediately after the earthquake. We analyze pre- and post-earthquake bathymetry data to assess the location and extent of submarine mass failures and sediment transport. New NOAA multibeam bathymetry shows the morphology of the entire fjord at 15 m resolution. We also assembled all older soundings from smooth sheets for comparison to the multibeam dataset. We gridded the sounding data, applied corrections for coseismic subsidence, post-seismic rebound, unrecovered co-seismic subsidence, sea-level rise (vertical datum shift), and measurement errors. The difference grids show changes resulting from the 1964 earthquake. We estimate the total volume of slide material to be about 211 million m3. Most of this material was transported to a deep, flat area, which we refer to as &ldquo;the bathtub&rdquo;, about 6 to 13 km south of Seward. Sub-bottom profiling of the bathtub shows an acoustically transparent unit, which we interpret as a sediment flow deposit resulting from the submarine landslides. The scale of the submarine landslides and the distance over which sediment was transported is much larger than previously appreciated.

Alaska↗

Reassessment of seismically induced, tsunamigenic submarine slope failures in Port Valdez, Alaska, USA

The M9.2 Alaska earthquake of 1964 caused major damage to the port facilities and town of Valdez, most of it the result of submarine landslide and the consequent tsunamis. Recent bathymetric multibeam surveys, high-resolution subbottom profiles, and dated sediment cores in Port Valdez supply new information about the morphology and character of the landslide deposits. A comparison of pre- and post-earthquake bathymetry provides an estimate of the net volume of landslide debris deposited in the basin and the volume of sediment removed from the source region. Landslide features include (1) large blocks (up to 40-m high) near the location of the greatest tsunamiwave runup (~50 m), (2) two debris lobes associated with the blocks, (3) a series of gullies, channels and talus, near the fjord-head delta and badly damaged old town of Valdez, and (4) the front of a debris lobe that flowed half-way down the fjord from the east end.

Alaska↗

Holocene slip rate and revised characteristic earthquake parameters for the western segment of the Castle Mountain fault, Alaska

The western segment of the Castle Mountain fault poses a significant seismic hazard to the most populated region of south-central Alaska. We identify a previously unrecognized margin of a postglacial outwash channel that is offset right laterally 36 ± 4 m across the western segment of the Castle Mountain fault. This offset occurred after glaciers withdrew from the lowland 11,300–15,380 cal yr b.p. and after outwash channel margins were cut and stabilized 11,210–13,470 cal yr b.p. Using these ages and the measured separation, we obtain a maximum slip rate of 3.0 ± 0.6 mm yr −1 and a minimum slip rate of 2.8 ± 0.7 mm yr −1 . These are the first lateral slip rates for the Castle Mountain fault established by a field measurement. Based on timing of the most recent earthquake, 670 ± 60 yr b.p. , the Castle Mountain fault could have accumulated an average single-event slip of about 1.9 m (extremes range from 1.3 to 2.6 m). The fault consists of two segments; a surface-rupturing earthquake likely will be limited to the 62-km-long western segment. Area-magnitude regression calculations suggest that such an earthquake on the western Castle Mountain fault would have a moment magnitude of 6.9 to 7.3.

Alaska↗

3D Visualization of Earthquake Focal Mechanisms Using ArcScene

We created a new tool, 3D Focal Mechanisms (3DFM), for viewing earthquake focal mechanism symbols three dimensionally. This tool operates within the Environmental Systems Research Institute (ESRI&reg;) GIS software ArcScene&reg; 9.x. The program requires as input a GIS point dataset of earthquake locations containing strike, dip, and rake values for a nodal plane of each earthquake. Other information, such as depth and magnitude of the earthquake, may also be included in the dataset. By default for each focal point, 3DFM will create a black and white sphere or &ldquo;beach ball&rdquo; that is oriented based on the strike, dip, and rake values. If depth values for each earthquake are included, the focal symbol will also be placed at its appropriate location beneath the Earth's surface. In addition to the default settings, there are several other options in 3DFM that can be adjusted. The appearance of the symbols can be changed by (1) creating rings around the fault planes that are colored based on magnitude, (2) showing only the fault planes instead of a sphere, (3) drawing a flat disc that identifies the primary nodal plane, (4) or by displaying the null, pressure, and tension axes. The size of the symbols can be changed by adjusting their diameter, scaling them based on the magnitude of the earthquake, or scaling them by the estimated size of the rupture patch based on earthquake magnitude. It is also possible to filter the data using any combination of the strike, dip, rake, magnitude, depth, null axis plunge, pressure axis plunge, tension axis plunge, or fault type values of the points. For a large dataset, these filters can be used to create different subsets of symbols. Symbols created by 3DFM are stored in graphics layers that appear in the ArcScene&reg; table of contents. Multiple graphics layers can be created and saved to preserve the output from different symbol options.

Data Series↗

Emplacement of the Kodiak batholith and slab-window migration

The Kodiak batholith is one of the largest, most elongate intrusive bodies in the forearc Sanak-Baranof plutonic belt located in southern Alaska. This belt is interpreted to have formed during the subduction of an oceanic spreading center and the associated migration of a slab window. Individual plutons of the Kodiak batholith track the location and evolution of the underlying slab window. Six U/Pb zircon ages from the axis of the batholith exhibit a northeastward-decreasing age progression of 59.2 &plusmn; 0.2 Ma at the southwest end to 58.4 &plusmn; 0.2 Ma at the northeast tip. The trench-parallel rate of age progression is within error of the average slab-window migration rate for the entire Sanak-Baranof belt (~19 cm/yr). Structural relationships, U/Pb ages, and a model of new gravity data indicate that magma from the Kodiak batholith ascended 5-10 km as a northeastward-younging series of 1-8-km-diameter viscoelastic diapirs. Individual plutons ascended by multiple emplacement mechanisms including downward flow, collapse of wall rock, stoping, and diking. Stokes flow xenolith calculations suggest ascent rates of 5-100 m/yr and an effective magmatic viscosity of 107-108 Pa s. Pre-existing structural or lithologic heterogeneities did not dominantly control the location of the main batholith. Instead, its location was determined by migration of the slab window at depth.

Geological Society of America Bulletin↗

Imaging the transition from Aleutian subduction to Yakutat collision in central Alaska, with local earthquakes and active source data

In southern and central Alaska the subduction and active volcanism of the Aleutian subduction zone give way to a broad plate boundary zone with mountain building and strike-slip faulting, where the Yakutat terrane joins the subducting Pacific plate. The interplay of these tectonic elements can be best understood by considering the entire region in three dimensions. We image three-dimensional seismic velocity using abundant local earthquakes, supplemented by active source data. Crustal low-velocity correlates with basins. The Denali fault zone is a dominant feature with a change in crustal thickness across the fault. A relatively high-velocity subducted slab and a low-velocity mantle wedge are observed, and high Vp/Vs beneath the active volcanic systems, which indicates focusing of partial melt. North of Cook Inlet, the subducted Yakutat slab is characterized by a thick low-velocity, high-Vp/Vs, crust. High-velocity material above the Yakutat slab may represent a residual older slab, which inhibits vertical flow of Yakutat subduction fluids. Alternate lateral flow allows Yakutat subduction fluids to contribute to Cook Inlet volcanism and the Wrangell volcanic field. The apparent northeast edge of the subducted Yakutat slab is southwest of the Wrangell volcanics, which have adakitic composition consistent with melting of this Yakutat slab edge. In the mantle, the Yakutat slab is subducting with the Pacific plate, while at shallower depths the Yakutat slab overthrusts the shallow Pacific plate along the Transition fault. This region of crustal doubling within the shallow slab is associated with extremely strong plate coupling and the primary asperity of the Mw 9.2 great 1964 earthquake. Copyright 2006 by the American Geophysical Union.

Alaska↗

Deformation driven by subduction and microplate collision: Geodynamics of Cook Inlet basin, Alaska

Late Neogene and younger deformation in Cook Inlet basin is caused by dextral transpression in the plate margin of south-central Alaska. Collision and subduction of the Yakutat microplate at the northeastern end of the Aleutian subduction zone is driving the accretionary complex of the Chugach and Kenai Mountains toward the Alaska Range on the opposite side of the basin. This deformation creates belts of fault-cored anticlines that are prolific traps of hydrocarbons and are also potential sources for damaging earthquakes. The faults dip steeply, extend into the Mesozoic basement beneath the Tertiary basin fill, and form conjugate flower structures at some localities. Comparing the geometry of the natural faults and folds with analog models created in a sandbox deformation apparatus suggests that some of the faults accommodate significant dextral as well as reverse-slip motion. We develop a tectonic model in which dextral shearing and horizontal shortening of the basin is driven by microplate collision with an additional component of thrust-type strain caused by plate subduction. This model predicts temporally fluctuating stress fields that are coupled to the recurrence intervals of large-magnitude subduction zone earthquakes. The maximum principal compressive stress is oriented east-southeast to east-northeast with nearly vertical least compressive stress when the basin's lithosphere is mostly decoupled from the underlying subduction megathrust. This stress tensor is compatible with principal stresses inferred from focal mechanisms of earthquakes that occur within the crust beneath Cook Inlet basin. Locking of the megathrust between great magnitude earthquakes may cause the maximum principal compressive stress to rotate toward the northwest. Moderate dipping faults that strike north to northeast may be optimally oriented for rupture in the ambient stress field, but steeply dipping faults within the cores of some anticlines are unfavorably oriented with respect to both modeled and observed stress fields, suggesting that elevated fluid pressure may be required to trigger fault rupture. ?? 2006 Geological Society of America.

Geological Society of America Bulletin↗