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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

68 records · Page 4Linked to original sources

USGS2016

Mihi Breccia: A stack of lacustrine sediments and subaqueous pyroclastic flows within the Taupo Volcanic Zone, New Zealand

The Taupo Volcanic Zone (TVZ), New Zealand, encompasses a wide variety of arc-related strata, although most of its small-volume (non-caldera-forming) eruptions are poorly-exposed and extensively hydrothermally altered. The Mihi Breccia is a stratigraphic sequence consisting of interbedded rhyolitic pyroclastic flows and lacustrine sediments with eruption ages of 281 ± 18 to at least 239 ± 6 ka (uncertainties at 2σ). In contrast to other small-volume rhyolitic eruptions within the TVZ, Mihi Breccia is relatively well-exposed within the Paeroa fault block, and contains minimal hydrothermal alteration. Pyroclastic flow characteristics and textures including: 1) breadcrusted juvenile clasts, 2) lack of welding, 3) abundant ash-rich matrix, 4) lack of fiamme and eutaxitic textures, 5) lack of thermal oxidation colors, 6) lack of cooling joints, 7) exclusive lacustrine sediment lithic clasts, and 8) interbedding with lacustrine sediments, all indicating that Mihi Breccia strata originated in a paleo-lake system. This ephemeral paleo-lake system is inferred to have lasted for > 50 kyr (based on Mihi Breccia age constraints), and referred to as Huka Lake. Mihi Breccia pyroclastic flow juvenile clast geochemistry and petrography correspond with similar-aged (264 ± 8, 263 ± 10, and 247 ± 4 ka) intra-caldera rhyolite domes filling the Reporoa caldera (source of the 281 ± 81 Kaingaroa Formation ignimbrite). These exposed intra-caldera rhyolite domes (as well as geophysically inferred subsurface domes) are proposed to be source vents for the Mihi Breccia pyroclastic flows. Soft-sediment deformation associated with Mihi Breccia strata indicate either seismic shock, rapid sediment loading during pyroclastic flow emplacement, or both. Thus, the Mihi Breccia reflects a prolonged series of subaqueous rhyolite dome building and associated pyroclastic flows, accompanied by seismic activity, emplaced into a large paleo-lake system within the TVZ.

Taupo Volcanic Zone↗

USGS2018

Soil property and class maps of the conterminous United States at 100-meter spatial resolution

With growing concern for the depletion of soil resources, conventional soil maps need to be updated and provided at finer and finer resolutions to be able to support spatially explicit human–landscape models. Three US soil point datasets—the National Cooperative Soil Survey Characterization Database, the National Soil Information System, and the Rapid Carbon Assessment dataset—were combined with a stack of over 200 environmental datasets and gSSURGO polygon maps to generate complete coverage gridded predictions at 100-m spatial resolution of six soil properties (percentage of organic C, total N, bulk density, pH, and percentage of sand and clay) and two US soil taxonomic classes (291 great groups [GGs] and 78 modified particle size classes [mPSCs]) for the conterminous United States. Models were built using parallelized random forest and gradient boosting algorithms as implemented in the ranger and xgboost packages for R. Soil property predictions were generated at seven standard soil depths (0, 5, 15, 30, 60, 100, and 200 cm). Prediction probability maps for US soil taxonomic classifications were also generated. Cross validation results indicated an out-of-bag classification accuracy of 60% for GGs and 66% for mPSCs; for soil properties, RMSE for leave-location-out cross-validation was 0.74 ( R 2 = 0.68), 17.8 wt% ( R 2 = 0.57), 12 wt% ( R 2 = 0.46), 3.63 wt% ( R 2 = 0.41), 0.2 g cm −3 ( R 2 = 0.42), and 0.27 wt% ( R 2 = 0.39) for pH, percent sand and clay, weight percentage of organic C, bulk density, and weight percentage of total N, respectively. Nine independent validation datasets were used to assess prediction accuracies for soil class models, and results ranged between 24 and 58% and between 24 and 93% for GG and mPSC prediction accuracies, respectively. Although mapping accuracies were variable and likely lower than gSSURGO in some areas, this modeling approach can enable easier integration of soil information with spatially explicit models compared with multicomponent map units.

Soil Science Society of America Journal↗

USGS2024

Pliocene–Pleistocene warm-water incursions and water mass changes on the Ross Sea continental shelf (Antarctica) based on foraminifera from IODP Expedition 374

International Ocean Discovery Program (IODP) Expedition 374 sailed to the Ross Sea in 2018 to reconstruct paleoenvironments, track the history of key water masses, and assess model simulations that show warm-water incursions from the Southern Ocean led to the loss of marine-based Antarctic ice sheets during past interglacials. IODP Site U1523 (water depth 828 m) is located at the continental shelf break, northeast of Pennell Bank on the southeastern flank of Iselin Bank, where it lies beneath the Antarctic Slope Current (ASC). This site is sensitive to warm-water incursions from the Ross Sea Gyre and modified Circumpolar Deep Water (mCDW) today and during times of past warming climate. Multiple incursions of subpolar or temperate planktic foraminifera taxa occurred at Site U1523 after 3.8 Ma and prior to ∼ 1.82 Ma. Many of these warm-water taxa incursions likely represent interglacials of the latest Early Pliocene and Early Pleistocene, including Marine Isotope Stage (MIS) Gi7 to Gi3 ( ∼ 3.72–3.65 Ma), and Early Pleistocene MIS 91 or 90 ( ∼ 2.34–2.32 Ma) and MIS 77–67 ( ∼ 2.03–1.83 Ma) and suggest warmer-than-present conditions and less ice cover in the Ross Sea. However, a moderately resolved age model based on four key events prohibits us from precisely correlating with Marine Isotope Stages established by the LR04 Stack; therefore, these correlations are best estimates. Diatom-rich intervals during the latest Pliocene at Site U1523 include evidence of anomalously warm conditions based on the presence of subtropical and temperate planktic foraminiferal species in what likely correlates with interglacial MIS G17 ( ∼ 2.95 Ma), and a second interval that likely correlates with MIS KM3 ( ∼ 3.16 Ma) of the mid-Piacenzian Warm Period. Collectively, these multiple incursions of warmer-water planktic foraminifera provide evidence for polar amplification during super-interglacials of the Pliocene and Early Pleistocene. Higher abundances of planktic and benthic foraminifera during the Mid- to Late Pleistocene associated with interglacials of the MIS 37–31 interval ( ∼ 1.23–1.07 Ma), MIS 25 ( ∼ 0.95 Ma), MIS 15 ( ∼ 0.60 Ma), and MIS 6–5e transition ( ∼ 0.133–0.126 Ma) also indicate a reduced ice shelf and relatively warm conditions, including multiple warmer interglacials during the Mid-Pleistocene Transition (MPT). A decrease in sedimentation rate after ∼ 1.78 Ma is followed by a major change in benthic foraminiferal biofacies marked by a decrease in Globocassidulina subglobosa and a decrease in mud ( < 63 µ m) after ∼ 1.5 Ma. Subsequent dominance of Trifarina earlandi biofacies beginning during MIS 15 ( ∼ 600 ka) indicate progressive strengthening of the Antarctic Slope Current along the shelf edge of the Ross Sea during the mid to Late Pleistocene. A sharp increase in foraminiferal fragmentation after the MPT ( ∼ 900 ka) and variable abundances of T. earlandi indicate higher productivity, a stronger but variable ASC during interglacials, and/or corrosive waters, suggesting changes in water masses entering (mCDW) and exiting (High Salinity Shelf Water or Dense Shelf Water) the Ross Sea since the MPT.

Journal of Micropalaeontology↗

USGS2001

The geology of a part of Acadia and the nature of the Acadian orogeny across Central and Eastern Maine

The zone of Acadian collision between the Medial New England and Composite Avalon terranes is well preserved in Maine. A transect from northwest (Rome) to southeast (Camden) crosses the eastern part of Medial New England comprising the Central Maine basin, Liberty-Orrington thrust sheet, and Fredericton trough, and the western part of Composite Avalon, including the Graham Lake, Clarry Hill, and Clam Cove thrust sheets. U-Pb geochronology of events before, during, and after the Acadian orogeny helps elucidate the nature and distribution of tectonostratigraphic belts in this zone and the timing of some Acadian events in the Northern Appalachians. The Central Maine basin consists of sedimentary and volcanic rocks of Middle Ordovician (∼470 to ∼460 Ma) age overlain with probable conformity by latest Ordovician(?) through earliest Devonian marine rift and flysch sedimentary rocks; these are intruded by weakly to undeformed plutonic rocks of Early and Middle Devonian age (∼399–378 Ma). The Fredericton trough consists of Early Silurian gray pelite and sandstone to earliest Late Silurian calcareous turbidite, deformed and variably metamorphosed prior to the emplacement of Late Silurian (∼422 Ma) and Early to Late Devonian (∼418 to ∼368 Ma) plutons. The Liberty-Orrington thrust sheet consists of Cambrian(?)-Ordovician (>∼474 to ∼469 Ma and younger) clastic sedimentary and volcanic rocks intruded by highly deformed Late Silurian (∼424 to ∼422 Ma) and Devonian (∼418 to ∼389 Ma) plutons, possibly metamorphosed in Late Silurian time (prior to ∼417 Ma), and metamorphosed to amphibolite facies in Early to Middle Devonian time (∼400 to ∼381 Ma). The Graham Lake thrust sheet contains possible Precambrian rocks, Cambrian sedimentary rocks with a volcanic unit dated at ∼503 Ma, and Ordovician rocks with possible Caradocian Old World fossils, metamorphosed and deformed in Silurian time and intruded by mildly to undeformed Late Silurian (∼421 Ma) and Late Devonian (∼371 to ∼368 Ma) plutons. The Clarry Hill thrust sheet consists of poorly studied, highly metamorphosed Cambrian (?) rocks. The Clam Cove thrust sheet contains highly deformed Precambrian limestone, shale, sandstone, and conglomerate, metamorphosed to epidote amphibolite facies and intruded by a mildly deformed pluton dated at ∼421 Ma. Metamorphism, deformation, and voluminous intrusive igneous activity of Silurian age are common to both the most southeastern parts of Medial New England and the thrust sheets of Composite Avalon. In contrast to Medial New England, the thrust sheets of Composite Avalon show only modest effects of Devonian deformation and metamorphism. Regional stratigraphic relations, paleontologic findings, and U-Pb geochronology suggest that the Graham Lake, Clarry Hill, and Clam Cove thrust sheets are far-traveled allochthons that were widely separated from Medial New England in the Silurian. One hundred nine (109) new U-Pb analyses of zircon, monazite, and sphene from 25 samples of metamorphosed, stratified, and intrusive igneous rocks are used to decipher the history of events along the transect. In our view, many of the structures within Medial New England were formed during the closing of the back-arc ocean during the waning stages of the Taconian orogeny. These features include the precursors to the upright folds found in Silurian sedimentary rocks, the west-facing Liberty-Orrington thrust sheet, and the Silurian plutons of Medial New England. Acadian tectonic features, including flysch sedimentation, igneous activity, deformation with nappe emplacement, and metamorphism record the progressive loading of Medial New England by a stack of thrust nappes emplaced in latest Silurian to Middle Devonian time. The Acadian orogeny is a prolonged event, lasting from earliest Late Silurian to the Late Devonian, whose evolution involved: (1) convergence between Medial New England and Composite Avalon along an east-dipping subduction zone from earliest Late Silurian to Early Devonian time; (2) collision and concurrent delamination of lithospheric mantle beneath Medial New England in Early Devonian time resulting in deformation, high-grade metamorphism, and intrusive igneous activity in the most eastern part of Medial New England and the western parts of Composite Avalon; (3) Early to Middle Devonian northwest-migrating penetrative deformations of the Acadian Main stage, including northwest-directed thrusting and recumbent folding followed by tightening of folds possibly produced in the waning stages of the Taconian orogeny and forming folded wave trains of isoclinal folds (∼419 to ∼404 Ma), (4) asymmetrical folds produced by east-west shortening (∼399 to ∼380 Ma); and (5) final westward emplacement of Composite Avalon thrust sheets onto Medial New England in Late Devonian time (∼380 to ∼371 Ma). The present boundary between rocks of Medial New England and Composite Avalon at the surface is the Sennebec Pond fault, a high-angle fault that cuts the Graham Lake and Clarry Hill thrusts, and is intruded by the Mt. Waldo pluton (∼371), one of several Late Devonian plutons (∼371 to ∼367 Ma) that mark the end of the Acadian orogeny in coastal Maine. The thrust sheets of Composite Avalon carried rocks of peri-Gondwanan affinity and with possible Old World fauna many kilometers to the west over Medial New England, thus potentially hiding the original suture at depth under the Gulf of Maine.

Maine↗

USGS2022

High-rate very-long-period seismicity at Yasur volcano, Vanuatu: source mechanism and decoupling from surficial explosions and infrasound

Yasur volcano, Vanuatu is a continuously active open-vent basaltic-andesite stratocone with persistent and long-lived eruptive activity. We present results from a seismo-acoustic field experiment at Yasur, providing locally dense broad-band seismic and infrasonic network coverage from 2016 July 27 to August 3. We corroborate our seismo-acoustic observations with coincident video data from cameras deployed at the crater and on an unoccupied aircraft system (UAS). The waveforms contain a profusion of signals reflecting Yasur’s rapidly occurring and persistent explosive activity. The typical infrasonic signature of Yasur explosions is a classic short-duration and often asymmetric explosion waveform characterized by a sharp compressive onset and wideband frequency content. The dominant seismic signals are numerous repetitive very-long-period (VLP) signals with periods of ∼2–10 s. The VLP seismic events are ‘high-rate’, reoccurring near-continuously throughout the data set with short interevent times (∼20–60 s). We observe variability in the synchronization of seismic VLP and acoustic sources. Explosion events clearly delineated by infrasonic waveforms are underlain by seismic VLPs. However, strong seismic VLPs also occur with only a weak infrasonic expression. Multiplet analysis of the seismic VLPs reveals a systematic progression in the seismo-acoustic source decoupling. The same dominant seismic VLP multiplet occurs with and without surficial explosions and infrasound, and these transitions occur over a timescale of a few days during our field campaign. We subsequently employ template matching, stacking, and full-waveform inversion to image the source mechanism of the dominant VLP multiplet. Inversion of the dominant VLP multiplet stack points to a composite source consisting of either a dual-crack (plus forces) or pipe-crack (plus forces) mechanism. The derived mechanisms correspond to a point-source directly beneath the summit vents with centroid depths in the range ∼900–1000 m below topography. All mechanisms suggest a northeast trending crack dipping relatively shallowly to the northwest and indicate a VLP source centroid and mechanism controlled by a stable structural geologic feature beneath Yasur. We interpret the results in the framework of gas slug ascent through the conduit responsible for Yasur explosions. The VLP mechanism and timing with infrasound (when present) are explained by a shallow-buffered top-down model in which slug ascent is relatively aseismic until reaching the base of a shallow section. Slug disruption in this shallow zone triggers a pressure disturbance that propagates downward and couples at the conduit base (VLP centroid). If the shallow section is open, an explosion propagates to the surface, producing infrasound. In the case of (the same multiplet) VLPs occurring without surficial explosions and weak or no infrasound, the decoupling of the dominant VLPs at ∼900–1000 m depth from surficial explosions and infrasound strongly indicates buffering of the terminal slug ascent. This buffering could be achieved by a variety of conditions at or directly beneath the vents, such as a high-viscosity layer of crystal-rich magma, a debris cap from backfill, a foam layer, or a combination of these. The dominant VLP at Yasur captured by our experiment has a source depth and mechanism separated from surface processes and is stable over time.

Geophysical Journal International↗

USGS2016

Tomographic Rayleigh-wave group velocities in the Central Valley, California centered on the Sacramento/San Joaquin Delta

If shaking from a local or regional earthquake in the San Francisco Bay region were to rupture levees in the Sacramento/San Joaquin Delta then brackish water from San Francisco Bay would contaminate the water in the Delta: the source of fresh water for about half of California. As a prelude to a full shear-wave velocity model that can be used in computer simulations and further seismic hazard analysis, we report on the use of ambient noise tomography to build a fundamental-mode, Rayleigh-wave group velocity model for the region around the Sacramento/San Joaquin Delta in the western Central Valley, California. Recordings from the vertical component of about 31 stations were processed to compute the spatial distribution of Rayleigh wave group velocities. Complex coherency between pairs of stations were stacked over 8 months to more than a year. Dispersion curves were determined from 4 to about 18 seconds. We calculated average group velocities for each period and inverted for deviations from the average for a matrix of cells that covered the study area. Smoothing using the first difference is applied. Cells of the model were about 5.6 km in either dimension. Checkerboard tests of resolution, which is dependent on station density, suggest that the resolving ability of the array is reasonably good within the middle of the array with resolution between 0.2 and 0.4 degrees. Overall, low velocities in the middle of each image reflect the deeper sedimentary syncline in the Central Valley. In detail, the model shows several centers of low velocity that may be associated with gross geologic features such as faulting along the western margin of the Central Valley, oil and gas reservoirs, and large cross cutting features like the Stockton arch. At shorter periods around 5.5s, the model’s western boundary between low and high velocities closely follows regional fault geometry and the edge of a residual isostatic gravity low. In the eastern part of the valley, the boundaries of the low velocity zone and gravity anomaly are better aligned at longer periods (around 10.5s) suggesting that the eastern edge of the gravity low is associated with deeper structure. There is a strong correspondence between a low in gravity near the Kirby Hills fault and low velocities from the ambient noise tomography. At longer periods, higher velocities creep in from the east and narrow the overall dimension defined by the lower velocities. Overall, there is a strong correspondence between the shape and location of low velocities in the Rayleigh wave velocity images, and geological and geophysical features.

California↗

USGS2004

The global short-period wavefield modelled with a Monte Carlo seismic phonon method

At high frequencies (∼1 Hz), much of the seismic energy arriving at teleseismic distances is not found in the main phases (e.g. P , PP , S , etc.) but is contained in the extended coda that follows these arrivals. This coda results from scattering off small-scale velocity and density perturbations within the crust and mantle and contains valuable information regarding the depth dependence and strength of this heterogeneity as well as the relative importance of intrinsic versus scattering attenuation. Most analyses of seismic coda to date have concentrated on S -wave coda generated from lithospheric scattering for events recorded at local and regional distances. Here, we examine the globally averaged vertical-component, 1-Hz wavefield (>10° range) for earthquakes recorded in the IRIS FARM archive from 1990 to 1999. We apply an envelope-function stacking technique to image the average time–distance behavior of the wavefield for both shallow (≤50 km) and deep (≥500 km) earthquakes. Unlike regional records, our images are dominated by P and P coda owing to the large effect of attenuation on PP and S at high frequencies. Modelling our results is complicated by the need to include a variety of ray paths, the likely contributions of multiple scattering and the possible importance of P -to- S and S -to- P scattering. We adopt a stochastic, particle-based approach in which millions of seismic phonons are randomly sprayed from the source and tracked through the Earth. Each phonon represents an energy packet that travels along the appropriate ray path until it is affected by a discontinuity or a scatterer. Discontinuities are modelled by treating the energy normalized reflection and transmission coefficients as probabilities. Scattering probabilities and scattering angles are computed in a similar fashion, assuming random velocity and density perturbations characterized by an exponential autocorrelation function. Intrinsic attenuation is included by reducing the energy contained in each particle as an appropriate function of traveltime. We find that most scattering occurs in the lithosphere and upper mantle, as previous results have indicated, but that some lower-mantle scattering is likely also required. A model with 3 to 4 per cent rms velocity heterogeneity at 4-km scale length in the upper mantle and 0.5 per cent rms velocity heterogeneity at 8-km scale length in the lower mantle (with intrinsic attenuation of Q α = 450 above 200 km depth and Q α = 2500 below 200 km) provides a reasonable fit to both the shallow- and deep-earthquake observations, although many trade-offs exist between the scale length, depth extent and strength of the heterogeneity.

Geophysical Journal International↗

USGS2013

Very long period conduit oscillations induced by rockfalls at Kilauea Volcano, Hawaii

Eruptive activity at the summit of Kilauea Volcano, Hawaii, beginning in 2010 and continuing to the present time is characterized by transient outgassing bursts accompanied by very long period (VLP) seismic signals triggered by rockfalls from the vent walls impacting a lava lake in a pit within the Halemaumau pit crater. We use raw data recorded with an 11-station broadband network to model the source mechanism of signals accompanying two large rockfalls on 29 August 2012 and two smaller average rockfalls obtained by stacking over all events with similar waveforms to improve the signal-to-noise ratio. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. We apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of Green's functions representing the seismometer response to translation and tilt ground motions. This method enables a robust description of the source mechanism over the period range 1–1000 s. The VLP signals associated with the rockfalls originate in a source region ∼1 km below the eastern perimeter of the Halemaumau pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks including an east striking crack (dike) dipping 80° to the north, intersecting a north striking crack (another dike) dipping 65° to the east. Each rockfall is marked by a similar step-like inflation trailed by decaying oscillations of the volumetric source, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes induced by the rock mass impacting the top of the lava column. Assuming a simple lumped parameter representation of the shallow magmatic system, the observed pressure and volume variations can be modeled with the following attributes: rockfall volume (200–4500 m 3 ), length of magma column (120–210 m), diameter of pipe connecting the Halemaumau pit crater to the subjacent dike system (6 m), average thickness of the two underlying dikes (3–6 m), and effective magma viscosity (30–210 Pa s). Most rockfalls occur during episodes of sustained deflation of the Kilauea summit. The mass loss rate in the shallow magmatic system is estimated to be 1400–15,000 kg s −1 based on measurements of the temporal variation of VLP period in the two large rockfalls that occurred on 29 August 2012.

Hawaii↗

USGS1994

Earthquake classification, location, and error analysis in a volcanic environment: implications for the magmatic system of the 1989-1990 eruptions at redoubt volcano, Alaska

Determination of the precise locations of seismic events associated with the 1989-1990 eruptions of Redoubt Volcano posed a number of problems, including poorly known crustal velocities, a sparse station distribution, and an abundance of events with emergent phase onsets. In addition, the high relief of the volcano could not be incorporated into the hypoellipse earthquake location algorithm. This algorithm was modified to allow hypocenters to be located above the elevation of the seismic stations. The velocity model was calibrated on the basis of a posteruptive seismic survey, in which four chemical explosions were recorded by eight stations of the permanent network supplemented with 20 temporary seismographs deployed on and around the volcanic edifice. The model consists of a stack of homogeneous horizontal layers; setting the top of the model at the summit allows events to be located anywhere within the volcanic edifice. Detailed analysis of hypocentral errors shows that the long-period (LP) events constituting the vigorous 23-hour swarm that preceded the initial eruption on December 14 could have originated from a point 1.4 km below the crater floor. A similar analysis of LP events in the swarm preceding the major eruption on January 2 shows they also could have originated from a point, the location of which is shifted 0.8 km northwest and 0.7 km deeper than the source of the initial swarm. We suggest this shift in LP activity reflects a northward jump in the pathway for magmatic gases caused by the sealing of the initial pathway by magma extrusion during the last half of December. Volcano-tectonic (VT) earthquakes did not occur until after the initial 23-hour-long swarm. They began slowly just below the LP source and their rate of occurrence increased after the eruption of 01:52 AST on December 15, when they shifted to depths of 6 to 10 km. After January 2 the VT activity migrated gradually northward; this migration suggests northward propagating withdrawal of magma from a plexus of dikes and/or sills located in the 6 to 10 km depth range. Precise relocations of selected events prior to January 2 clearly resolve a narrow, steeply dipping, pencil-shaped concentration of activity in the depth range of 1-7 km, which illuminates the conduit along which magma was transported to the surface. A third event type, named hybrid, which blends the characteristics of both VT and LP events, originates just below the LP source, and may reflect brittle failure along a zone intersecting a fluid-filled crack. The distribution of hybrid events is elongated 0.2-0.4 km in an east-west direction. This distribution may offer constraints on the orientation and size of the fluid-filled crack inferred to be the source of the LP events. ?? 1994.

Journal of Volcanology and Geothermal Research↗

USGS2026

Subsurface structure across the Tacoma Basin, Washington State, using trans-dimensional Bayesian inversion of fundamental mode spatial autocorrelation data

Spatial autocorrelation (SPAC), the azimuthal average of the normalized cross-correlation between equidistant station pairs deployed in a 2-D array, is widely used to image the subsurface structure. However, the rigorous estimate of subsurface structure and its uncertainties as a function of depth using SPAC data is challenging due to the nonlinear relation between the SPAC data and Earth structure as well as the trade-off between depth and velocity. Additionally, data noise is strongly correlated due to data processing (e.g. filtering, stacking from multiple time segments and azimuthal averaging). Most studies do not account for the correlated noise and fix the ratio of compressional-wave velocity ( V P ) to shear-wave velocity ( V s ) (i.e. V P / V s ratio) and the number of layers, both of which are typically unknown. To address these challenges, we develop a hierarchical trans-dimensional Bayesian inversion of fundamental mode of SPAC data that properly accounts for the correlated data noise, samples the V P / V s ratio and relaxes the number of layers (i.e. model parametrization) to be unknown in the inversion. We further examine the limitation of using only fundamental modes in the inversion. Our synthetic experiments show that the inversion recovers an incorrect model unless we sample the correlated noise and V P / V s ratio in the inversion. The inversion is then applied to SPAC data acquired at 19 sites across the Tacoma basin in Washington State to characterize the V s and the time-averaged V s over 30-m depth ( V s 30 ). Our results show that the V s 30 varies from ∼200 to 800 m s −1 . The V s 30 within the basin is higher in the middle and lower on the east and west sides. We find that these V s 30 values vary with geologic unit. The uncertainties for V s 30 are within 20 m s −1 in average except for the most eastern site TB28. Additionally, the uncertainties are greater for deeper depths beneath most of the sites as the sensitivity decreases as a function of depth. The Vs structure as a function of depth is also complex beneath some sites, possibly because the SPAC curves are affected by higher order Rayleigh modes that are not considered in the inversion. To better constrain the deeper V s structure, V s 30 and/or other average measures of V s over depth, additional constraints from complementary data, such as ellipticity or geologic data are needed. Moreover, our synthetic experiments show that higher order modes can have significant effect in the inversion results, particularly when there is a low-velocity layer.

Washington↗

USGS2021

Deepwater debrites and linked megaturbidites in confined basins: An example from the Onnuri Basin, East Sea of Korea

We analyzed data from seven piston cores, multi-channel seismic-reflection (MCS) and chirp profiles, and multibeam echosounder (MBES) data to study the distribution, emplacement time, sedimentary facies, and depositional processes of sediment-gravity-flow deposits in the Onnuri Basin, a confined basin in the East Sea. These data reveal that debris flows have traveled ca. 30 km downslope, forming a seismic facies consisting of stacked, wedge-shaped, transparent units separated by high-amplitude continuous reflectors. Analysis of piston cores shows three distinct sedimentary units, throughout the basin. The lowest unit, I, is a debrite containing numerous mud clasts of varying size and color distributed in a mud-rich matrix; it is absent over elevated basinal highs or ridges, such as the Onnuri Ridge, suggesting that local topography controls its distribution. The debrite forms a recognizable acoustically transparent layer on subbottom chirp profiles (av. 7 m thick), covers approximately 500 km 2 , and has an estimated volume of ∼ 3.5 km 3 . The overlying unit, II, contains normally graded beds composed of massive sand, laminated and cross-laminated sand and silt, and a thick cap of structureless mud. This unit is interpreted to be a megaturbidite deposited from turbidity currents that originated from the flow transformation of debris flows on the upper continental slope. The megaturbidite covers the entire basin (at least 650 km 2 ), and has an average thickness of 2.8 m (maximum thickness of 4.35 m), and comprises a volume of 1.8 km 3 . Variations in grain size and sedimentary structures suggest that the megaturbidite was deposited by progressively waning flows that reflected off basin flanks and ridges. The thick (up to 3.65 m) structureless mud cap further indicates deposition in a confined basin. The sharp basal contact, together with the lack of hemipelagic sediments between debrite and overlying megaturbidite, suggest that both were deposited during the same flow event, likely to have originated from a single catastrophic slope failure. Collapsing slide material evolved into a debris flow, from which a turbidite formed by dilution of the debris flow. Radiocarbon dates suggest that the slope failure occurred about 13–11 ka, a time when sea level was ca. 50 m lower than at the present day. Hemipelagic sediments in the topmost unit, III-2, above the megaturbidite indicate that the basin has been stable since ca. 11 ka. We provide robust evidence that submarine slope failures evolve downslope into slides, debris flows, and finally, thick megaturbidites. This contribution highlights the importance of seafloor morphology on the distribution and stratigraphy of submarine flows in confined basins.

Onnuri Basin↗

USGS2004

Application of ground-penetrating radar, digital optical borehole images, and cores for characterization of porosity hydraulic conductivity and paleokarst in the Biscayne aquifer, southeastern Florida, USA

This paper presents examples of ground-penetrating radar (GPR) data from two study sites in southeastern Florida where karstic Pleistocene platform carbonates that comprise the unconfined Biscayne aquifer were imaged. Important features shown on resultant GPR profiles include: (1) upward and lateral qualitative interpretative distribution of porosity and hydraulic conductivity; (2) paleotopographic relief on karstic subaerial exposure surfaces; and (3) vertical stacking of chronostratigraphic high-frequency cycles (HFCs). These characteristics were verified by comparison to rock properties observed and measured in core samples, and identified in digital optical borehole images. Results demonstrate that an empirical relation exists between measured whole-core porosity and hydraulic conductivity, observed porosity on digital optical borehole images, formation conductivity, and GPR reflection amplitudes-as porosity and hydraulic conductivity determined from core and borehole images increases, formation conductivity increases, and GPR reflection amplitude decreases. This relation allows for qualitative interpretation of the vertical and lateral distribution of porosity and hydraulic conductivity within HFCs. Two subtidal HFCs in the uppermost Biscayne aquifer have significantly unique populations of whole-core porosity values and vertical hydraulic conductivity values. Porosity measurements from one cycle has a median value about two to three times greater than the values from the other HFC, and median values of vertical hydraulic-conductivity about three orders of magnitude higher than the other HFC. The HFC with the higher porosity and hydraulic conductivity values is shown as a discrete package of relatively low-amplitude reflections, whereas the HFC characterized by lower porosity and hydraulic-conductivity measurements is expressed by higher amplitude reflections. Porosity and hydraulic-conductivity values measured from whole-core samples, and vuggy porosity identified on digital borehole images from shallowing-upward, peritidal HFCs show that the highest porosity occurs at the base of the cycles, moderate porosity at the middle of the cycles, and lowest porosity occurs at the top of cycles. Hydraulic conductivity is also highest at the base of the peritidal cycles and lowest in the middle to upper parts of cycles. This change in porosity and hydraulic conductivity from bottom to top is visible as an upward variation in reflection amplitude on GPR profiles-lowest amplitudes at the base and highest at the cycle tops. This study demonstrates that GPR can be used to show the qualitative distribution of porosity and hydraulic conductivity within a cycle-stratigraphic framework composed of carbonate HFCs. The distribution of porosity and hydraulic conductivity within HFCs is related to depositional textures. The upward and lateral patterns of the rock facies within the HFCs can be translated to geophysical-log properties and radar facies configurations that could aid in interpretation and prediction of ground-water flow through a carbonate aquifer. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Applied Geophysics↗

USGS2012

Pillow basalts of the Angayucham terrane: Oceanic plateau and island crust accreted to the Brooks Range

The Angayucham Mountains (north margin of the Yukon-Koyukuk province) are made up of an imbricate stack of four to eight east-west trending, steeply dipping, fault slabs composed of Paleozoic (Devonian to Mississippean), Middle to Late Triassic, and Early Jurassic oceanic upper crustal rocks (pillow basalt, subordinate diabase, basaltic tuff, and radiolarian chert). Field relations and geochemical characteristics of the basaltic rocks suggest that the fault slabs were derived from an oceanic plateau or island setting and were emplaced onto the Brooks Range continental margin. The basalts are variably metamorphosed to prehnite-pumpellyite and low-greenschist facies. Major element analyses suggest that many are hypersthene-normative olivine tholeiites. Classification based on immobile trace elements confirms the tholeiitic character of most of the basalts but suggests that some had primary compositions transitional to alkali basalt. Although field and petrographic features of the basalts are similar, trace element characteristics allow definition of geographically distinct suites. A central outcrop belt along the crest of the mountains is made up of basalt with relatively flat rare earth element (REE) patterns. This belt is flanked to the north and south by LREE (light rare earth element)-enriched basalts. Radiolarian and conodont ages from interpillow and interlayered chert and limestone indicate that the central belt of basalts is Triassic in age, the southern belt is Jurassic in age, and the northern belt contains a mixture of Paleozoic and Mesozoic ages. Data for most of the basalts cluster in the “within-plate basalt” fields of trace element discriminant diagrams; none have trace-element characteristics of island arc basalt. The Triassic and Jurassic basalts are geochemically most akin to modern oceanic plateau and island basalts. Field evidence also favors an oceanic plateau or island setting. The great composite thickness of pillow basalt probably resulted from obduction faulting, but the lack of fault slabs of gabbro or peridotite suggests that obduction faults did not penetrate below oceanic layer 2, a likely occurrence if layer 2 were anomalously thick, as in the vicinity of an oceanic island. The presence of basaltic tuff interbeds indicates proximity to an explosive basaltic eruptive center. The juxtaposition of submarine basalts of differing chemical affinity and age, adjacent to higher-grade Paleozoic metamorphic rocks of the Brooks Range to the north, may be explained by obduction of internally complex (thickened) oceanic crust formed in an ocean plateau setting. Emplacement and rotation of thrust plates to steep attitudes occurred during accretion of the Brooks Range passive margin, probably beginning in the Late to Middle Jurassic.

Journal of Geophysical Research Solid Earth↗

USGS2012

Source parameters of microearthquakes on an interplate asperity off Kamaishi, NE Japan over two earthquake cycles

We have estimated the source parameters of interplate earthquakes in an earthquake cluster off Kamaishi, NE Japan over two cycles of M~ 4.9 repeating earthquakes. The M~ 4.9 earthquake sequence is composed of nine events that occurred since 1957 which have a strong periodicity (5.5 ± 0.7 yr) and constant size ( M 4.9 ± 0.2), probably due to stable sliding around the source area (asperity). Using P - and S -wave traveltime differentials estimated from waveform cross-spectra, three M~ 4.9 main shocks and 50 accompanying microearthquakes ( M 1.5–3.6) from 1995 to 2008 were precisely relocated. The source sizes, stress drops and slip amounts for earthquakes of M 2.4 or larger were also estimated from corner frequencies and seismic moments using simultaneous inversion of stacked spectral ratios. Relocation using the double-difference method shows that the slip area of the 2008 M~ 4.9 main shock is co-located with those of the 1995 and 2001 M~ 4.9 main shocks. Four groups of microearthquake clusters are located in and around the mainshock slip areas. Of these, two clusters are located at the deeper and shallower edge of the slip areas and most of these microearthquakes occurred repeatedly in the interseismic period. Two other clusters located near the centre of the mainshock source areas are not as active as the clusters near the edge. The occurrence of these earthquakes is limited to the latter half of the earthquake cycles of the M~ 4.9 main shock. Similar spatial and temporal features of microearthquake occurrence were seen for two other cycles before the 1995 M 5.0 and 1990 M 5.0 main shocks based on group identification by waveform similarities. Stress drops of microearthquakes are 3–11 MPa and are relatively constant within each group during the two earthquake cycles. The 2001 and 2008 M~ 4.9 earthquakes have larger stress drops of 41 and 27 MPa, respectively. These results show that the stress drop is probably determined by the fault properties and does not change much for earthquakes rupturing in the same area. The occurrence of microearthquakes in the interseismic period suggests the intrusion of aseismic slip, causing a loading of these patches. We also found that some earthquakes near the centre of the mainshock source area occurred just after the earthquakes at the deeper edge of the mainshock source area. These seismic activities probably indicate episodic aseismic slip migrating from the deeper regions in the mainshock asperity to its centre during interseismic periods. Comparison of the source parameters for the 2001 and 2008 main shocks shows that the seismic moments (1.04 x 10 16 Nm and 1.12 x 10 16 Nm for the 2008 and 2001 earthquakes, respectively) and source sizes (radius = 570 m and 540 m for the 2008 and 2001 earthquakes, respectively) are comparable. Based on careful phase identification and hypocentre relocation by constraining the hypocentres of other small earthquakes to their precisely located centroids, we found that the hypocentres of the 2001 and 2008 M~ 4.9 events are located in the southeastern part of the mainshock source area. This location does not correspond to either episodic slip area or hypocentres of small earthquakes that occurred during the earthquake cycle.

Kamaishi↗