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.

85 records · Page 5Linked to original sources

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↗

USGS2012

Seismic reflection images of the crust of the northern part of the Chugach terrane, Alaska: Results of a survey for the Trans-Alaska Crustal Transect (TACT)

Deep crustal seismic reflection data show strong reflections from the middle and lower crust of the convergent continental margin near the eastern end of the Aleutian trench. These data were collected across the Border Ranges fault system, a major suture zone that separates the Peninsular and Chugach tectonostratigraphic terranes. The Chugach terrane consists of metamorphosed, strongly deformed flysch, subduction melange as well as blueschist- and greenschist-facies rocks. Even so, highly reflective rocks lie at great depth within or beneath this terrane. The shallow part of the seismic section (0–5 s) reveals only a few reflections, which undulate and interweave, indicating that some of them originate outside of the plane of the seismic section. Other undulating events could reveal early Cenozoic granitic plutons or antiformal stacks of rock imbricated along thrust faults. In contrast to this shallow reflection pattern, three reflection bands contain most of the events on the middle part of the seismic section (5–12 s or 12–34 km). We interpret divergent reflections within the uppermost band (5–5.5 s) as gently north dipping thrust faults, and subparallel reflections truncated at the boundaries of the upper and middle bands may image a metamorphic foliation or thrust faults. If the truncated reflections represent thrust faults, then the reflection geometry suggests that subparallel roof and floor thrust faults bound imbricated rocks, forming duplex structures. The upper and middle bands are separated by an area on the seismic section that reveals few reflections. This area correlates approximately in depth with a high velocity (7.5 km/s) layer evident in refraction models. The poorly reflective, high velocity unit might be volcanic rocks like those exposed in the structurally lowest part of the Chugach terrane. If so, then the bottom of the poorly reflective rocks indicates the base of the Chugach terrane. Rocks that cause the upper and middle reflection bands and the intervening, poorly reflective volcanic rocks may form a midcrustal shear zone that is about 10 km thick. The deepest band of subparallel reflections (30–35 km) correlates closely in depth with the top of the Wadati-Benioff zone associated with the underthrusting Pacific plate. The reflections probably stem from the decollement that separates lithospheric plates within the subduction zone. This reflection band increases abruptly in thickness, possibly indicating local underplating of subducted sediment. No events from the oceanic Moho have been recognized. Neither the Border Ranges fault system, the suture between the Chugach and Peninsular terranes, nor deep crustal layers of the Peninsular terrane are evident in these seismic reflection data.

Journal of Geophysical Research Solid Earth↗

USGS2002

Late quaternary evolution of the Orinoco Delta, Venezuela

The modern Orinoco Delta is the latest of a series of stacked deltas that have infilled the Eastern Venezuelan Basin (EVB) since the Oligocene. During the late Pleistocene sea-level lowstand (20,000 to 16,000 yrs BP), bedrock control points at the position of the present delta apex prevented the river channel from incising as deeply as many other major river systems. Shallow seismic data indicate that the late Pleistocene Orinoco incised into the present continental shelf, where it formed a braided-river complex that transported sediment to a series of shelf-edge deltas. As sea level rose from 16,000 to 9,500 yrs BP, the Orinoco shoreline shifted rapidly landward, causing shallow-marine waves and currents to form a widespread transgressive sand unit. Decelerating sea-level rise and a warmer, wetter climate during the early Holocene (9,500 to 6,000 yrs BP) induced delta development within the relatively quiet-water environment of the EVB embayment. Sea level approached its present stand in the middle Holocene (6,000 to 3,000 yrs BP), and the Orinoco coast prograded, broadening the delta plain and infilling the EVB embayment. Significant quantities of Amazon sediment began to be transported to the Orinoco coast by littoral currents. Continued progradation in the late Holocene caused the constriction at Boca de Serpientes to alter nearshore and shelf hydrodynamics and subdivide the submarine delta into two distinct areas: the Atlantic shelf and the Gulf of Paria. The increased influence of littoral currents along the coast promoted mudcape development. Because most of the water and sediment were transported across the delta plain through the Rio Grande distributary in the southern delta, much of the central and northwestern delta plain became sediment starved, promoting widespread accumulation of peat deposits. Human impacts on the delta are mostly associated with the Volca??n Dam on Can??o Manamo. However, human activities have had relatively little effect on the delta processes and environments.

Journal of Coastal Research↗

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↗

USGS1996

Trace fossils from Jurassic lacustrine turbidites of the Anyao Formation (Central China) and their environmental and evolutionary significance

The Lower Jurassic Anyao Formation crops out near Jiyuan city, western Henan Province, central China. It is part of the infill of the nonmarine early Mesozoic Jiyuan‐Yima Basin. In the Jiyuan section, this unit is about 100 m thick and consists of laterally persistent, thin and thick‐bedded turbidite sandstones and mudstones displaying complete and base‐or top‐absent Bouma sequences, and thick‐bedded massive sandstones. The Anyao Formation records sedimentation within a lacustrine turbidite system developed in a pull‐apart basin. Processes involved include high and low density turbidity currents, sometimes affected by liquefaction or fluidization. Facies analysis suggests that this succession is formed by stacked aggradational turbidite lobes. The absence of thick mudstone packages indicates that background sedimentation was subordinate to high frequency turbidite deposition. The Anyao Formation hosts a moderately diverse ichnofauna preserved as hypichnial casts on the soles of thin‐bedded turbidite sandstones. The ichnofauna consists of Cochlichnus anguineus, Hel‐minthoidichnites tenuis, Helminthopsis abeli, H. hieroglyphica, Mono‐morphichnus lineatus, Paracanthorhaphe togwunia, Tuberculichnus vagans, Vagorichnus anyao, tiny grazing trails, and irregularly branching burrows. Vagorichnus anyao occurs not only as a discrete trace, but also as a compound ichnotaxon intergrading with Gordia marina and Tuberculichnus vagans. Both predepositional and post‐depositional traces are present on the soles of turbidites. This ichnofauna comprises both feeding and grazing traces produced by a deposit‐feeding lacustrine benthic biota. Crawling traces are rare. Although certain ichnofossils (e.g. V. anyao, P. togwunia) show overall similarities with deep‐sea agrichnia, they differ in reflecting remarkably less specialized feeding strategies, displaying overcrossing between specimens (and to a lesser extent, self‐crossing), and in the case of V’ anyao recording post‐turbidite burrowing activity. The development of less specialized strategies than those displayed by deep‐marine ichnofaunas may be related to less stable conditions, typical of lake settings. Oxyenation, energy, sedimentation rate (both event and background), food supply, soft‐sediment deformation and erosion rate have mainly influenced trace‐fossil distribution. Turbidity currents would have ensured oxygen (as well as food) supply to deep lake settings, thus allowing the establishment of a moderately diverse biota. Biogenic structures were mostly confined to the outer, low energy areas. High sedimentation rates and strong erosion precluded preservation of ichnofossils in inner lobe settings. The Anyao ichnofauna is of significance in furthering knowledge of the colonization of deep lakes throughout the stratigraphie record and in identification of additional nonmarine ichnofacies. The ana‐lyzed ichnofauna resembles late Paleozoic lacustrine assemblages described from different localities around the world and is regarded as a Mesozoic example of the Mermia ichnofacies. However, when compared with Paleozoic assemblages, the Anyao ichnofauna shows a clear dominance of burrows over surface trails, deeper burrowing penetration, larger size, and presence of relatively more complex structures. The high burrow/surface trail ratio may be indicative of lower preservation potential in the latter, thus reflecting a tap‐honomic overprint. In contrast to the Paleozoic examples, the establishment of a relatively well‐developed lacustrine infauna may have precluded preservation of surface trails. Burrower activity probably obliterated biogenic structures formed close to the sediment‐water interface.

Ichnos: An International Journal for Plant and Ani↗

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↗