Search USGSSearch

SEARCH · Search USGS

Results for “Stacks Journal”

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 55 records · Page 3Linked to original sources

Chemical variations in Yellowknife Bay formation sedimentary rocks analyzed by ChemCam on board the Curiosity rover on Mars

The Yellowknife Bay formation represents a ~5 m thick stratigraphic section of lithified fluvial and lacustrine sediments analyzed by the Curiosity rover in Gale crater, Mars. Previous works have mainly focused on the mudstones that were drilled by the rover at two locations. The present study focuses on the sedimentary rocks stratigraphically above the mudstones by studying their chemical variations in parallel with rock textures. Results show that differences in composition correlate with textures and both manifest subtle but significant variations through the stratigraphic column. Though the chemistry of the sediments does not vary much in the lower part of the stratigraphy, the variations in alkali elements indicate variations in the source material and/or physical sorting, as shown by the identification of alkali feldspars. The sandstones contain similar relative proportions of hydrogen to the mudstones below, suggesting the presence of hydrous minerals that may have contributed to their cementation. Slight variations in magnesium correlate with changes in textures suggesting that diagenesis through cementation and dissolution modified the initial rock composition and texture simultaneously. The upper part of the stratigraphy (~1 m thick) displays rocks with different compositions suggesting a strong change in the depositional system. The presence of float rocks with similar compositions found along the rover traverse suggests that some of these outcrops extend further away in the nearby hummocky plains.

Journal of Geophysical Research E: Planets

Topobathymetric elevation model development using a new methodology: Coastal National Elevation Database

During the coming decades, coastlines will respond to widely predicted sea-level rise, storm surge, and coastalinundation flooding from disastrous events. Because physical processes in coastal environments are controlled by the geomorphology of over-the-land topography and underwater bathymetry, many applications of geospatial data in coastal environments require detailed knowledge of the near-shore topography and bathymetry. In this paper, an updated methodology used by the U.S. Geological Survey Coastal National Elevation Database (CoNED) Applications Project is presented for developing coastal topobathymetric elevation models (TBDEMs) from multiple topographic data sources with adjacent intertidal topobathymetric and offshore bathymetric sources to generate seamlessly integrated TBDEMs. This repeatable, updatable, and logically consistent methodology assimilates topographic data (land elevation) and bathymetry (water depth) into a seamless coastal elevation model. Within the overarching framework, vertical datum transformations are standardized in a workflow that interweaves spatially consistent interpolation (gridding) techniques with a land/water boundary mask delineation approach. Output gridded raster TBDEMs are stacked into a file storage system of mosaic datasets within an Esri ArcGIS geodatabase for efficient updating while maintaining current and updated spatially referenced metadata. Topobathymetric data provide a required seamless elevation product for several science application studies, such as shoreline delineation, coastal inundation mapping, sediment-transport, sea-level rise, storm surge models, and tsunami impact assessment. These detailed coastal elevation data are critical to depict regions prone to climate change impacts and are essential to planners and managers responsible for mitigating the associated risks and costs to both human communities and ecosystems. The CoNED methodology approach has been used to construct integrated TBDEM models in Mobile Bay, the northern Gulf of Mexico, San Francisco Bay, the Hurricane Sandy region, and southern California.

Journal of Coastal Research

Multiple stacked plays to drive continued Permian development

The exploitation of unconventional (continuous) resource plays with horizontal drilling and fracturing techniques continues to spur renewed development in mature producing basins in the US. The Permian basin of west Texas has emerged as a particularly strong prospect due to the potential for multiple stacked resource plays in formations such as the Wolfcamp shale and Spraberry and Bone Spring formations. Production began in the Permian basin in the early 1920s, and more than 116,000 wells have been drilled in Midland basin alone. 1 2 In 2017, total Permian basin oil production climbed to nearly 1.7 million b/d. 3 The US Geological Survey (USGS) conducts geologic-based assessments of undiscovered, technically recoverable domestic hydrocarbon resources. The USGS recently assessed the Wolfcamp shale in the Midland basin portion of the Permian basin for continuous (unconventional) oil and gas resources, and completed a reassessment of the Spraberry formation for both continuous and conventional oil and gas. 4 5 Combined, the assessments of the Spraberry formation and underlying Wolfcamp shale in Midland basin comprise the largest domestic continuous oil assessment in the contiguous US (Fig. 1).

Texas

Hydrological, geomorphological, and chemical effects of Hurricane Andrew on coastal marshes of Louisiana

Hydrological conditions are reported for the hurricane-induced storm surge at marsh sites in Terrebonne and Barataria Basins. There was a 6 hour or greater delay in the arrival of the storm surge at interior marsh sites compared to coastal locations. A water surge of about 1.5 m with a salinity of 10–15‰ was recorded at an intermediate salinity marsh site in Terrebonne Basin where significant geomorphological changes occurred with lateral movement and stacking of marsh vegetation and sediment, erosion, and the conversion of wetlands to open water. A residual salt wedge attributed to the storm surge was detected in sediments at this intermediate marsh site about 55 days after the hurricane. In the pore water, at 30 cm in depth, there were increases in alkalinity, ammonia, phosphate, and manganese, but not iron or sulfide, although there was a significant depletion of sulfur with respect to chloride. The increases in concentration of these species in the pore water suggests that there was an increase in organic matter decomposition. Although increases in sulfide phytotoxicity were not obvious, the residual salinity of pore water may influence wetland loss through longer-term impacts on macrophyte production. Storm-induced sediment deposition, which will help ameliorate wetland loss, was observed at a salt marsh site in Terrebonne Basin.

Louisiana

3-D P- and S-wave velocity structure and low-frequency earthquake locations in the Parkfield, California region

To refine the 3-D seismic velocity model in the greater Parkfield, California region, a new data set including regular earthquakes, shots, quarry blasts and low-frequency earthquakes (LFEs) was assembled. Hundreds of traces of each LFE family at two temporary arrays were stacked with time–frequency domain phase weighted stacking method to improve signal-to-noise ratio. We extend our model resolution to lower crustal depth with LFE data. Our result images not only previously identified features but also low velocity zones (LVZs) in the area around the LFEs and the lower crust beneath the southern Rinconada Fault. The former LVZ is consistent with high fluid pressure that can account for several aspects of LFE behaviour. The latter LVZ is consistent with a high conductivity zone in magnetotelluric studies. A new Vs model was developed with S picks that were obtained with a new autopicker. At shallow depth, the low Vs areas underlie the strongest shaking areas in the 2004 Parkfield earthquake. We relocate LFE families and analyse the location uncertainties with the NonLinLoc and tomoDD codes. The two methods yield similar results.

California

Semiautomated estimates of directivity and related source properties of small to moderate southern California earthquakes using second seismic moments

We develop a semiautomated method for estimating with second seismic moments the directivity, rupture area, duration, and centroid velocity of earthquakes. The method is applied to 41 southern California earthquakes with magnitude in the range 3.5–5.2 and provides stable results for 28 events. Apparent source time functions (ASTFs) of P and S phases are derived using deconvolution with three stacked empirical Green's functions (seGf). The use of seGf suppresses nongeneric source effects, improves the focal mechanism correspondence to the analyzed earthquakes, and typically allows inclusion of 5 to 15 more ASTFs compared with analysis using a single eGf. Most analyzed earthquakes in the Trifurcation area of the San Jacinto Fault have directivities toward the northwest, while events around Cajon Pass and San Gabriel Mountain tend to propagate toward the southeast. These results are generally consistent with predictions for dynamic rupture on bimaterial interfaces associated with the imaged velocity contrasts in the area. The second moment inversions also provide constraints on the upper and lower bounds of rupture areas in our data set. Stress drops and uncertainties are estimated for elliptical ruptures using the derived characteristic rupture length and width. The semiautomated second moment method with seGfs can be used for routine application to moderate earthquakes in locations with good station coverage.

California

Exhumation history of the NW Indian Himalaya revealed by fission track and 40Ar/39Ar ages

New fission track and Ar/Ar geochronological data provide time constraints on the exhumation history of the Himalayan nappes in the Mandi (Beas valley) &mdash; Tso Morari transect of the NW Indian Himalaya. Results from this and previous studies suggest that the SW-directed North Himalayan nappes were emplaced by detachment from the underthrusted upper Indian crust by 55 Ma and metamorphosed by ca. 48&ndash;40 Ma. The nappe stack was subsequently exhumed to shallow upper crustal depths (<10 km) by 40&ndash;30 Ma in the Tso Morari dome (northern section of the transect) and by 30&ndash;20 Ma close to frontal thrusts in the Baralacha La region. From the Oligocene to the present, exhumation continued slowly. Metamorphism started in the High Himalayan nappe prior to the Late Oligocene. High temperatures and anatexis of the subducting upper Indian crust engendered the buoyancy-driven ductile detachment and extrusion of the High Himalayan nappe in the zone of continental collision. Late extrusion of the High Himalayan nappe started about 26 Ma ago, accompanied by ductile extensional shearing in the Zanskar shear zone in its roof between 22 and 19 Ma concomitant with thrusting along the basal Main Central Thrust to the south. The northern part of the nappe was then rapidly exhumed to shallow depth (<10 km) between 20 and 6 Ma, while its southern front reached this depth at 10&ndash;5 Ma.

Journal of Asian Earth Sciences

Geometry of the décollement below eastern Bangladesh and implications for seismic hazard

Eastern Bangladesh sits on the seismically active Chittagong-Myanmar fold and thrust belt (CMFB), a north-trending accretionary wedge on the eastern side of the India-Eurasia collision. Earthquakes on the basal décollement and associated thrusts within the CMFB present a hazard to this densely populated region. In this study, we interpret 28 seismic reflection profiles from both published and unpublished sources to constrain the depth of the basal décollement. To convert profiles from the time domain to the depth domain, we integrate sonic log and seismic stacking velocity data to generate time-velocity relationships for different parts of the CMFB. Our analysis reveals that the décollement is ∼9 km deep in northeast and southeast Bangladesh, but shallows to ∼5 km in east-central Bangladesh. The décollement has an area of 7.25 × 10 4 km 2 (∼150 × 450 km), making it capable of an Mw 8.5 earthquake. However, the warped geometry of this fault might act as a rupture barrier were a large earthquake to occur on the décollement. Our combined velocity and fault model lay the groundwork for future studies to address seismic segmentation, ground shaking, and rupture modeling in the CMFB. Finally, we use our compiled data set to analyze the evolution of fold kinematics in the CMFB. We observe that folding style and failure mode varies, from mainly ductile deformation in the foreland to mainly brittle in the hinterland. The dual-failure modes within the CMFB support the hypothesis that a region with ductile deformation may still be capable of seismic behavior.

Journal of Geophysical Research: Solid Earth

Probing the volcanic-plutonic connection and the genesis of crystal-rich rhyolite in a deeply dissected supervolcano in the Nevada Great Basin: Source of the late Eocene Caetano Tuff

Late Cenozoic faulting and large-magnitude extension in the Great Basin of the western USA has created locally deep windows into the upper crust, permitting direct study of volcanic and plutonic rocks within individual calderas. The Caetano caldera in north–central Nevada, formed during the mid-Tertiary ignimbrite flare-up, offers one of the best exposed and most complete records of caldera magmatism. Integrating whole-rock geochemistry, mineral chemistry, isotope geochemistry and geochronology with field studies and geologic mapping, we define the petrologic evolution of the magmatic system that sourced the >1100 km 3 Caetano Tuff. The intra-caldera Caetano Tuff is up to ∼5 km thick, composed of crystal-rich (30–45 vol. %), high-silica rhyolite, overlain by a smaller volume of comparably crystal-rich, low-silica rhyolite. It defies classification as either a monotonous intermediate or crystal-poor zoned rhyolite, as commonly ascribed to ignimbrite eruptions. Crystallization modeling based on the observed mineralogy and major and trace element geochemistry demonstrates that the compositional zonation can be explained by liquid–cumulate evolution in the Caetano Tuff magma chamber, with the more evolved lower Caetano Tuff consisting of extracted liquids that continued to crystallize and mix in the upper part of the chamber following segregation from a cumulate-rich, and more heterogeneous, source mush. The latter is represented in the caldera stratigraphy by the less evolved upper Caetano Tuff. Whole-rock major, trace and rare earth element geochemistry, modal mineralogy and mineral chemistry, O, Sr, Nd and Pb isotope geochemistry, sanidine Ar–Ar geochronology, and zircon U–Pb geochronology and trace element geochemistry provide robust evidence that the voluminous caldera intrusions (Carico Lake pluton and Redrock Canyon porphyry) are genetically equivalent to the least evolved Caetano Tuff and formed from magma that remained in the lower chamber after ignimbrite eruption and caldera collapse. Thus, the Caetano Tuff contradicts models for the mutually exclusive origins of voluminous volcanic and plutonic magmas in the upper crust. Crystal-scale O isotope data indicate that the Caetano Tuff is one of the most 18 O-enriched rhyolites in the Great Basin (δ 18 O magma = 10·2 ± 0·2‰), supporting anatexis of local metasedimentary basement crust. Metapelite xenoliths in the Carico Lake pluton and ubiquitous xenocrystic zircons in the Caetano Tuff provide constraints for the anatexis process; these data point to shallow (<15 km) dehydration melting of a protolith similar to the Proterozoic McCoy Creek Group siliciclastic sediments in eastern Nevada, projected beneath Caetano in fault-stacked shelf sediments that were thickened during Mesozoic crustal shortening. Mean zircon U–Pb ages for different stratigraphic levels of the intra-caldera Caetano Tuff are 34·2–34·5 Ma, 0·2–0·5 Myr older than the caldera sanidine 40 Ar/ 39 Ar age of 34·00 ± 0·03 Ma, documenting protracted duration of assembly and homogenization of isotopically diverse upper crustal melts, followed by crystallization and zonation to generate the Caetano Tuff magma chamber. Sanidine rims in the least evolved Caetano Tuff and in the Carico Lake pluton and Redrock Canyon porphyry have sharply zoned Ba domains that point to crystal growth during magmatic recharge events. The recharge magma is inferred to have been compositionally similar to the Caetano Tuff magma, with increased Ba resulting from remelting of Ba-rich sanidine cumulates. Mush reactivation to generate the Caetano Tuff eruption was sufficiently rapid to preserve compositional gradients in the intracaldera ignimbrite, calling into question models that predict homogeneity as a prerequisite for remobilizing crystal-rich ignimbrite magmas.

Journal of Petrology

An integrated mudstone facies classification scheme and revised interpretation of the sedimentologic processes driving carbon burial in the Cenomanian–Turonian Greenhorn Formation, Colorado, U.S.A.

Standardizing facies descriptions has proven key to integrating interpretations of depositional processes and environments from sedimentologic observations with geochemistry data for mudstone lithologies. Because of their fine-grained nature, high degree of compaction, and heterogeneous composition, standardizing methods for mudstone descriptions has proven difficult, but it is critical to formulating meaningful interpretations of the processes that govern the accumulation of organic-rich lithologies and their role in both petroleum systems and the global carbon cycle. In this study, we have developed a modified facies classification scheme for mudstone lithologies that incorporates sedimentologic and compositional observation at the hand-sample and thin-section scales with geochemical measurements, including bulk organic and inorganic geochemistry, to characterize these rocks and their variability more completely for improved interpretations of depositional environments during a low-order sea-level transgression. The facies described in this study are of the Cenomanian–Turonian Greenhorn Formation in the USGS #1 Portland Core drilled in Fremont County, Colorado. Strata of the Greenhorn Formation span Oceanic Anoxic Event 2 (OAE-2) and the preceding interval. Lithologies range from organic-rich argillaceous mudstones with varied sedimentary structures to organic-lean, highly bioturbated limestones. Six facies were identified, each differentiated by varied sedimentary structures and geochemical composition. These facies occur in a predictable stratigraphic stacking pattern that represents a low-order sea-level transgression with interpreted depositional environments ranging from terrigenous-dominated pro-delta and muddy continental shelf at the base of the interval to pelagic offshore marine at the top of the Greenhorn Formation. Though the facies are consistent with previous interpretations of depositional environments at this locale in the Cretaceous Western Interior Seaway during the Greenhorn cyclothem, the sedimentary processes governing the accumulation of organic-rich strata that have defined this interval are significantly revised. Variability in the proximity and intensity of bottom currents driven by storms and geostrophic flows were key to the accumulation of each facies, with significant sediment transport occurring even through deposition in the most oxygen-depleted bottom waters. The methodology and interpretations provided here are now being employed to basin-scale predictions of organic enrichment utilizing calibrated petrophysical methods. The approach and results from this study improve understanding of how organic and inorganic carbon was sequestered during perturbations to the global carbon cycle associated with events such as OAE-2.

Colorado

Crystal structure and chemistry of lithium-bearing trioctahedral micas-3T

Chemical analyses and crystal structure refinements were performed on lithian siderophyllite-3T crystals from granitic pegmatites of the anorogenic Pikes Peak batholith (Colorado) to characterize the crystal chemistry and relations with trioctahedral lithium-bearing micas showing different stacking sequences. Chemical data show that the studied samples fall on the siderophyllite-polylithionite join, closer to the siderophyllite end-member. Single-crystal X-ray refinements were carried out on three samples (two of which were taken from core and rim of the same crystal) in space-group P3 1 12 (the agreement factor, Robs, varies between 0.034 and 0.036). Mean bond distances and mean electron counts of M1, M2 and M3 octahedral sites indicate an ordered cation distribution with M1 and M3 positions substantially larger than M2. In the sample with the largest iron content, the M2 mean electron count increases as well as the mean distance, whereas remains smaller than or . The tetrahedral cation-oxygen atom mean distances range from 1.614 to 1.638 A and from 1.663 to 1.678 A for T1 and T2 sites, respectively, being consistent with Al 3+ enrichment in the T2 sites. The tetrahedral rotation angle, α, is generally small (3.1 ≤ α ≤ 4.6) and decreases with siderophyllite content. As Fe increases, the T1 tetrahedron becomes flatter (112.4 ≤ t1 ≤ 110.5??), whereas T2 tetrahedron distortion appears unchanged (110.7 ≤ T2 ≤ 110.9).

European Journal of Mineralogy

Using open-science workflow tools to produce SCEC CyberShake physics-based probabilistic seismic hazard models

The Statewide (formerly Southern) California Earthquake Center (SCEC) conducts multidisciplinary earthquake system science research that aims to develop predictive models of earthquake processes, and to produce accurate seismic hazard information that can improve societal preparedness and resiliency to earthquake hazards. As part of this program, SCEC has developed the CyberShake platform, which calculates physics-based probabilistic seismic hazard analysis (PSHA) models for regions with high-quality seismic velocity and fault models. The CyberShake platform implements a sophisticated computational workflow that includes over 15 individual codes written by 6 developers. These codes are heterogeneous, ranging from short-running high-throughput serial CPU codes to large, long-running, parallel GPU codes. Additionally, CyberShake simulation campaigns are computationally extensive, typically producing tens of terabytes of meaningful scientific data and metadata over several months of around-the-clock execution on leadership-class supercomputers. To meet the needs of the CyberShake platform, we have developed an extreme-scale workflow stack, including the Pegasus Workflow Management System, HTCondor, Globus, and custom tools. We present this workflow software stack and identify how the CyberShake platform and supporting tools enable us to meet a variety of challenges that come with large-scale simulations, such as automated remote job submission, data management, and verification and validation. This platform enabled us to perform our most recent simulation campaign, CyberShake Study 22.12, from December 2022 to April 2023. During this time, our workflow tools executed approximately 32,000 jobs, and used up to 73% of the Summit system at Oak Ridge Leadership Computing Facility. Our workflow tools managed about 2.5 PB of total temporary and output data, and automatically staged 19 million output files totaling 74 TB back to archival storage on the University of Southern California's Center for Advanced Research Computing systems, including file-based relational data and large binary files to efficiently store millions of simulated seismograms. CyberShake extreme-scale workflows have generated simulation-based probabilistic seismic hazard models that are being used by seismological, engineering, and governmental communities.

Frontiers Earth Science Journal

Geology and geochronology of the Spirit Mountain batholith, southern Nevada: Implications for timescales and physical processes of batholith construction

The Spirit Mountain batholith (SMB) is a ∼ 250 km2 composite silicic intrusion located within the Colorado River Extensional Corridor in southernmost Nevada. Westward tilting of 40–50° has exposed a cross-section from the roof through deep levels of the batholith. Piecemeal construction is indicated by zircon geochronology, field relations, and elemental geochemistry. Zircon U/Pb data (SHRIMP) demonstrates a ∼ 2 million year (17.4–15.3 Ma) history for the SMB. Individual samples contain zircons with ages that span the lifetime of the batholith, suggesting recycling of extant zircon into new magma pulses. Field relations reveal several distinct intrusive episodes and suggest a common injection geometry of stacked horizontal sheets. The largest unit of the SMB is a gradational section (from roof downward) of high-silica leucogranite through coarse granite into foliated quartz monzonite. Solidification of this unit spans most of the history of the batholith. The 25 km × 2 km leucogranite was emplaced incrementally as subhorizontal sheets over most or all of the history of this section, suggesting repeated fractional crystallization and melt segregation events. The quartz monzonite and coarse granite are interpreted to be cumulate residuum of this fractionation. Age data from throughout this gradational unit show multiple zircon populations within individual samples. Subsequent distinct intrusions that cut this large unit, which include minor populations of zircons that record the ages of earlier events in construction of the batholith, preserve a sheeted, sill-on-sill geometry. We envision the SMB to have been a patchwork of melt-rich, melt-poor, and entirely solid zones throughout its active life. Preservation of intrusion geometries and contacts depended on the consistency of the host rock. Zircons recycled into new pulses of magma document remobilization of previously emplaced crystal mush, suggesting the mechanisms by which evidence for initial construction of the batholith became blurred.

Nevada

Long period seismicity and very long period infrasound driven by shallow magmatic degassing at Mount Pagan, Mariana Islands

Long period (LP) seismicity and very long period infrasound (iVLP) were recorded during continuous degassing from Mount Pagan, Mariana Islands, in July 2013 to January 2014. The frequency content of the LP and iVLP events and delay times between the two arrivals were remarkably stable and indicate nearly co-located sources. Using phase-weighted stacking over similar events to dampen noise, we find that the LP source centroid is located 60&thinsp;m below and 180&thinsp;m west of the summit vent. The moment tensor reveals a volumetric source modeled as resonance of a subhorizontal sill intersecting a dike. We model the seismoacoustic wavefields with a coupled earth-air 3-D finite difference code. The ratios of pressure to velocity measured at the infrasound arrays are an order of magnitude larger than the synthetic ratios, so the iVLP is not the result of LP energy transmitting into the atmosphere at its epicenter. Based on crater shape and dimensions determined by structure from motion, we model the iVLP as acoustic resonance of an exponential horn. The source of the continuous plume from gas analysis is shallow magmatic degassing, which repeatedly pressurized the dike-sill portion of the conduit over the 7&thinsp;months of observation. Periodic gas release caused the geologically controlled sill to partially collapse and resonate, while venting of gas at the surface triggered resonance in the crater. LP degassing only accounts for ~12% of total degassing, indicating that most degassing is relatively aseismic and that multiple active pathways exist beneath the vent.

Mount Pagan

Tectonic setting of the Yukon-Koyukuk basin and its borderlands, western Alaska

The Yukon-Koyukuk basin of western Alaska is composed of an arcuate belt of Jurassic and Early Cretaceous subduction-related volcanic and plutonic rocks (Koyukuk terrane) flanked by deep subbasins filled with mid-Cretaceous terrigenous sedimentary rocks. The basin is bordered on three sides by metamorphosed Proterozoic and Paleozoic continental rocks (Seward, Arctic Alaska, and Ruby terranes) and is separated from the metamorphic borderlands by a narrow, highly tectonized belt of oceanic crust and mantle rocks (composite Angayucham-Tozitna terrane). The oceanic and mantle rocks, which dip inward beneath the basin and are thrust outward onto the borderlands, are divided into three separate thrust panels: (1) a structurally lowest panel (Slate Creek) composed of phyllite and metagraywacke of probable Devonian age, (2) a middle panel (Narvak) composed of imbricated basalt, chert, and gabbro of Devonian to Early Jurassic age, and (3) a structurally highest panel (Kanuti) composed of gabbro and peridotite of probable Middle and Late Jurassic age. The three thrust panels appear to represent a reversely stacked sequence that progresses from continental slope deposits in the lower panel to cumulus and mantle peridotites in the upper. Metamorphic mineral K-Ar ages from garnet amphibolite on the sole of the upper panel suggest that the upper panel was emplaced on the middle panel in the Middle to Late Jurassic. Subsequent accretion of all three panels to the continental rocks of the borderlands occurred in the latest Jurassic and Early Cretaceous, synchronous with arc volcanism within the basin. Arc volcanism waned and the accretionary phase ended in the middle of Early Cretaceous time. Uplift and erosion of the metamorphic borderlands and the obducted oceanic rocks began in late Early Cretaceous and was accompanied by the rapid filling of two flanking subbasins with turbiditic sediments. In the latest Early Cretaceous and early Late Cretaceous, shallow marine and nonmarine conglomerates were deposited around the margins of the basin, and a prograding delta was built out from the southeast margin of the basin across the turbiditic subbasins and the remnant volcanic arc. In the Late Cretaceous, western Alaska was subjected to strong east-west compression which severely deformed both the Yukon-Koyukuk basin and the borderlands. Several widespread magmatic episodes in the mid- and Late Cretaceous and in early Tertiary transgress the basin boundaries and stitch together the accreted arc and oceanic terranes and the continental borderlands.

Journal of Geophysical Research Solid Earth

Eruption dynamics at Mount St. Helens imaged from broadband seismic waveforms: Interaction of the shallow magmatic and hydrothermal systems

The current eruption at Mount St. Helens is characterized by dome building and shallow, repetitive, long-period (LP) earthquakes. Waveform cross-correlation reveals remarkable similarity for a majority of the earthquakes over periods of several weeks. Stacked spectra of these events display multiple peaks between 0.5 and 2 Hz that are common to most stations. Lower-amplitude very-long-period (VLP) events commonly accompany the LP events. We model the source mechanisms of LP and VLP events in the 0.5-4 s and 8-40 s bands, respectively, using data recorded in July 2005 with a 19-station temporary broadband network. The source mechanism of the LP events includes: 1) a volumetric component modeled as resonance of a gently NNW-dipping, steam-filled crack located directly beneath the actively extruding part of the new dome and within 100 m of the crater floor and 2) a vertical single force attributed to movement of the overlying dome. The VLP source, which also includes volumetric and single-force components, is 250 m deeper and NNW of the LP source, at the SW edge of the 1980s lava dome. The volumetric component points to the compression and expansion of a shallow, magma-filled sill, which is subparallel to the hydrothermal crack imaged at the LP source, coupled with a smaller component of expansion and compression of a dike. The single-force components are due to mass advection in the magma conduit. The location, geometry and timing of the sources suggest the VLP and LP events are caused by perturbations of a common crack system.

Journal of Geophysical Research B: Solid Earth

Using a modified time-reverse imaging technique to locate low-frequency earthquakes on the San Andreas Fault near Cholame, California

We present a new method to locate low-frequency earthquakes (LFEs) within tectonic tremor episodes based on time-reverse imaging techniques. The modified time-reverse imaging technique presented here is the first method that locates individual LFEs within tremor episodes within 5 km uncertainty without relying on high-amplitude P-wave arrivals and that produces similar hypocentral locations to methods that locate events by stacking hundreds of LFEs without having to assume event co-location. In contrast to classic time-reverse imaging algorithms, we implement a modification to the method that searches for phase coherence over a short time period rather than identifying the maximum amplitude of a superpositioned wavefield. The method is independent of amplitude and can help constrain event origin time. The method uses individual LFE origin times, but does not rely on a priori information on LFE templates and families.We apply the method to locate 34 individual LFEs within tremor episodes that occur between 2010 and 2011 on the San Andreas Fault, near Cholame, California. Individual LFE location accuracies range from 2.6 to 5 km horizontally and 4.8 km vertically. Other methods that have been able to locate individual LFEs with accuracy of less than 5 km have mainly used large-amplitude events where a P-phase arrival can be identified. The method described here has the potential to locate a larger number of individual low-amplitude events with only the S-phase arrival. Location accuracy is controlled by the velocity model resolution and the wavelength of the dominant energy of the signal. Location results are also dependent on the number of stations used and are negligibly correlated with other factors such as the maximum gap in azimuthal coverage, source&ndash;station distance and signal-to-noise ratio.

California

The crustal thickness of Australia

We investigate the crustal structure of the Australian continent using the temporary broadband stations of the Skippy and Kimba projects and permanent broadband stations. We isolate near-receiver information, in the form of crustal P-to-S conversions, using the receiver function technique. Stacked receiver functions are inverted for S velocity structure using a Genetic Algorithm approach to Receiver Function Inversion (GARFI). From the resulting velocity models we are able to determine the Moho depth and to classify the width of the crust-mantle transition for 65 broadband stations. Using these results and 51 independent estimates of crustal thickness from refraction and reflection profiles, we present a new, improved, map of Moho depth for the Australian continent. The thinnest crust (25 km) occurs in the Archean Yilgarn Craton in Western Australia; the thickest crust (61 km) occurs in Proterozoic central Australia. The average crustal thickness is 38.8 km (standard deviation 6.2 km). Interpolation error estimates are made using kriging and fall into the range 2.5–7.0 km. We find generally good agreement between the depth to the seismologically defined Moho and xenolith-derived estimates of crustal thickness beneath northeastern Australia. However, beneath the Lachlan Fold Belt the estimates are not in agreement, and it is possible that the two techniques are mapping differing parts of a broad Moho transition zone. The Archean cratons of Western Australia appear to have remained largely stable since cratonization, reflected in only slight variation of Moho depth. The largely Proterozoic center of Australia shows relatively thicker crust overall as well as major Moho offsets. We see evidence of the margin of the contact between the Precambrian craton and the Tasman Orogen, referred to as the Tasman Line.

Journal of Geophysical Research B: Solid Earth