Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Interpretation Research”

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 505 records · Page 28Linked to original sources
USGS2012

An axial view of a metamorphic core complex: Crustal structure of the Whipple and Chemehuevi Mountains, southeastern California

A 135‐km‐long, NW‐SE trending, seismic refraction/wide‐angle reflection profile provides a unique along‐strike view of the crustal structure of a belt of metamorphic core complexes in southeastern California: the Whipple, Chemehuevi, and Sacramento mountains metamorphic core complexes. Interpretation of the seismic data was done by two‐dimensional forward modeling of travel times and amplitudes. The final model consists of (1) a thin (< 1.5 km) veneer of upper plate and fractured lower plate rocks (velocities of 1.5–5.3 km s −1 ) overlying a fairly homogeneous basement with velocities of 6.0 km s −1 ; (2) a localized, high‐velocity (6.4 km s −1 ) body, situated directly beneath the Whipple Mountains; (3) a 6.3–6.4 km s −1 middle crust that is thickest beneath the core complexes; (4) a 6.65±0.15 km s −1 lower crust; (5) crustal thickness of 27 km with a deeper crustal root (3 km) beneath the Whipple Mountains metamorphic core complex; and (6) a P n velocity of 8.0±0.10 km s −1 . The crustal structure that underlies the belt of metamorphic core complexes provides new insights into the processes that control extension in the deep crust. Upper crustal velocities are higher beneath the Whipple Mountains (where velocities increase to 6.4 km s −1 at ∼5 km depth) than beneath the Chemehuevi and Sacramento mountains. In addition, midcrustal discontinuities rise 2–5 km beneath the Whipple complex compared to the other complexes. These observations support greater uplift and a slightly deeper midcrustal origin for the rocks now exposed in the core of the Whipple Mountains compared to rocks in the Chemehuevi and Sacramento mountains. Despite the enhanced uplift and extension in the Whipple Mountains, the crust is thicker here (30 km) than anywhere else along the Colorado River extensional corridor. This may be in part a relic of compressional and magmatic thickening during the Mesozoic. However, we suggest that inflation of the crust during Tertiary extension was the dominant mechanism. Both mantle‐derived magmatism and lateral ductile inflow in the crust are proposed.

Arizona, California↗
USGS2012

Igneous history of the Koyukuk terrane, western Alaska: Constraints on the origin, evolution, and ultimate collision of an accreted island arc terrane

The Koyukuk terrane of western Alaska consists of volcanic, volcaniclastic, and plutonic rocks which range from Late Paleozoic to Early Cretaceous in age. The terrane crops out in a U-shaped belt which is roughly paralleled by outer belts of ultramafic rocks, oceanic plate basalts and cherts, and retrograded blueschist facies rocks of continental protolith. These rocks have been interpreted as components of a volcanic arc terrane that collided with the North American continental margin in Early Cretaceous time. The Koyukuk terrane consists of four time-stratigraphic units: (1) pre-Middle Jurassic basalts, (2) Middle and Late Jurassic granitic rocks, (3) lower Lower Cretaceous volcanic rocks, and (4) upper Lower Cretaceous volcanic rocks. Limited chemical data from the basalts of unit 1 indicate that they were erupted in a nonarc tectonic environment, possibly in an oceanic island or back arc setting. Units 2, 3, and 4 have the characteristics of subduction-related volcanic rocks (i.e., depleted Nb and Ta and enriched alkaline elements, relative to the light rare earth elements). Unit 3 contains tholeiitic, calc-alkaline, and alkaline rocks with chondrite-normalized rare earth element patterns that range from flat (La N /Yb N = 1) to highly light rare earth element enriched (La N /Yb N > 15). The highly alkaline or shoshonitic lavas were erupted toward the end of unit 3 time (Valanginian) during the final stages of arc-continent collision. These alkaline lavas could have been derived by very small degrees of partial melting of a similar source to that of the earlier arc lavas. Unit 4 lavas are also alkaline or shoshonitic, but their incompatible element composition indicates that they were derived from a different source than that of the earlier arc lavas. These late alkaline lavas are chemically similar to crosscutting mid-Cretaceous plutons whose isotopic compositions (Arth et al., this issue ( a )) suggest derivation by partial melting of distinctly older subcontinental lithosphere. We speculate that the parental magmas of unit 4 lavas may also have been derived by partial melting of this subcontinental mantle which was underthrust beneath the Koyukuk arc terrane during the final stage of arc-continent collision.

Journal of Geophysical Research Solid Earth↗
USGS2012

The Yellowstone-Snake River Plain seismic profiling experiment: Crustal structure of the Eastern Snake River Plain

Seismic refraction profiles recorded along the eastern Snake River Plain (ESRP) in southeastern Idaho during the 1978 Yellowstone-Snake River Plain cooperative seismic profiling experiment are interpreted to infer the crustal velocity and attenuation (Q-1) structure of the ESRP. Travel-time and synthetic seismogram modeling of a 250 km reversed refraction profile as well as a 100 km detailed profile indicate that the crust of the ESRP is highly anomalous. Approximately 3 to 6 km of volcanic rocks (with some interbedded sediments) overlie an upper-crustal layer (compressional velocity ≅6.1 km/s) which thins southwestward along the ESRP from a thickness of 10 km near Island Park Caldera to 2 to 3 km beneath the central and southwestern portions of the ESRP. An intermediate-velocity (≅6.5 km/s) layer extends from ≅10 to ≅20 km depth. A thick (≅22 km) lower crust of compressional velocity 6.8 km/s, a total crustal thickness of ≅42 km, and a P n velocity of ≅7.9 km/s is observed in the ESRP, similar to the western Snake River Plain and the Rocky Mountains Provinces. High attenuation is evident on the amplitude corrected seismic data due to low-Q values in the volcanic rocks (Q p = 20 to 200) and throughout the crust (Q p = 160 to 300). Based on these characteristics of the crustal structure and volcanic-age progression data, it is suggested that the ESRP has resulted from an intensive period of intrusion of mantle-derived basaltic magma into the upper crust generating explosive silicic volcanism and associated regional uplift and caldera collapse. This activity began about 15 m.y. ago in southwestern Idaho and has migrated northeast to its present position at Yellowstone. Subsequent cooling of the intruded upper crust results in the 6.5 km/s velocity intermediate layer. Crustal subsidence and periodic basaltic volcanism as represented by the ESRP complete the sequence of crustal evolution.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Correlation of the Deccan and Rajahmundry Trap lavas: Are these the longest and largest lava flows on Earth?

We propose that the Rajahmundry Trap lavas, found near the east coast of peninsular India , are remnants of the longest lava flows yet recognized on Earth (∼ 1000 km long). These outlying Deccan-like lavas are shown to belong to the main Deccan Traps. Several previous studies have already suggested this correlation, but have not demonstrated it categorically. The exposed Rajahmundry lavas are interpreted to be the distal parts of two very-large-volume pāhoehoe flow fields, one each from the Ambenali and Mahabaleshwar Formations of the Wai Sub-group in the Deccan Basalt Group. Eruptive conditions required to emplace such long flows are met by plausible values for cooling and eruption rates, and this is shown by applying a model for the formation of inflated pāhoehoe sheet flow lobes. The model predicts flow lobe thicknesses similar to those observed in the Rajahmundry lavas. For the last 400 km of flow, the lava flows were confined to the pre-existing Krishna valley drainage system that existed in the basement beyond the edge of the gradually expanding Deccan lava field, allowing the flows to extend across the subcontinent to the eastern margin where they were emplaced into a littoral and/or shallow marine environment. These lavas and other individual flow fields in the Wai Sub-group may exceed eruptive volumes of 5000 km 3 , which would place them amongst the largest magnitude effusive eruptive units yet known. We suggest that the length of flood basalt lava flows on Earth is restricted mainly by the size of land masses and topography. In the case of the Rajahmundry lavas, the flows reached estuaries and the sea, where their advance was perhaps effectively terminated by cooling and/or disruption. However, it is only during large igneous province basaltic volcanism that such huge volumes of lava are erupted in single events, and when the magma supply rate is sufficiently high and maintained to allow the formation of very long lava flows. The Rajahmundry lava fields were emplaced around 65 Ma during the later times of Deccan volcanism, probably just after the K/T environmental crisis. However, many lava-forming eruptions of similar magnitude and style straddled the K/T boundary.

Rajahmundry Trap↗
USGS2022

Olivine and glass chemistry record cycles of plumbing system recovery after summit collapse events at Kīlauea Volcano, Hawai‘i

The eruptive activity of Kīlauea Volcano (Hawai‘i) in the past 2500 years has alternated between centuries-long periods dominated either by explosive or effusive eruptions. The onset of explosive periods appears to be marked by caldera collapse events at the volcano's summit accompanied by draining of Kīlauea's magmatic plumbing system. Here we leverage >1800 olivine forsterite (Fo) contents, >900 glass MgO contents, and estimated magma supply rates from the past six centuries to describe the relationships between summit collapse and the composition of erupted material. On a first order basis, the major element chemistry of the centuries-long eruptive periods largely originates from fundamental differences between fractional crystallization of shallowly stored magmas during high-supply effusive-dominated periods versus little evolution of mafic recharge magmas during low-supply explosive-dominated periods. The modern effusive period (1820s-present) is dominated by relatively evolved olivine forsterite contents (Fo81–82) for Kīlauea, which is interpreted to reflect a buffered crustal reservoir system in which shallow storage and fractional crystallization control the composition of magmas. In contrast, olivine crystals from the explosive Keanakāko‘i Tephra (1500 - early 1800s C.E.) are dominated by higher olivine forsterite contents (Fo89) which are interpreted to reflect more primitive compositions, are correlated with glass MgO compositions extending to high values (e.g.,11.0 wt%), and show signs of magma mixing (zoned olivine, bimodal Fo populations). These signatures reflect a disrupted reservoir system in which high-MgO recharge melts mix with melts left over from draining of the shallow (<5 km) magma plumbing. Superimposed on these explosive-effusive periods are three decades- to centuries long periods of progressively evolving olivine and glass compositions. Eruptions that occur after caldera collapse in ~1500C.E. and smaller scale crater collapse events in 1790 (inferred) and 1924 have heterogeneous olivine populations dominated by ≥Fo88 and typically high MgO glasses. These compositions reflect inefficient mixing of stored and primitive recharge magmas after the disruption of the shallow plumbing system. After these collapses, olivine Fo and glass MgO subsequently evolve to <Fo82 and <7.0 wt% compositions, reflecting the recovery of the crustal plumbing system to an end-member system state characterized by efficient mixing of recharge and stored magmas that serve to buffer the shallow magma reservoirs. These evolved signatures suggest that a mature and buffered reservoir system may be a key condition for significant disruptions of volcanic plumbing systems. Plumbing system recovery is slower following large-scale caldera collapse (hundreds of years) compared to recovery following smaller crater collapse (tens of years), which may be modulated by differences in magma supply rates. Following the 2018 crater collapse olivine populations have high-Fo but glasses are low MgO, suggesting that this collapse might have disrupted shallow magma pathways but not strongly impacted the reservoir(s). Ultimately, olivine and glass major element chemistry record the impacts of caldera and smaller but significant summit crater collapse events at Kīlauea and could be used to provide a framework for better characterizing long-term volcano evolution in Hawai‘i and shield volcanoes elsewhere.

Hawaii↗
USGS1999

Chronologic model and transgressive-regressive signatures in the late neocene siliciclastic foundation (long key formation) of the Florida keys

Recent drilling of continuous cores in southernmost Florida has documented a thick unit of upper Neogene siliciclastics subjacent to surficial shallow-water Quaternary carbonates exposed on islands of the Florida Keys. The siliciclastics comprise the Long Key Formation and were identified in two cores collected from the middle and upper Florida Keys. A chronologic model based on new planktic foraminiferal biochronology and strontium-isotope chronology suggests the timing of siliciclastic deposition and provides a basis for regional correlation. The chronologic model, supplemented by vertical trends in quartz grain size, pattern of planktic menardiiform coiling direction, and paleoenvironmental interpretations of benthic foraminiferal assemblages, shows that the Long Key Formation contains three intervals (I-III) of varying thickness, grain-size composition, and paleowater depth. Interval I is uppermost Miocene. The quartz grains in Interval I fine upward from basal very coarse sand to fine and very fine sand. Benthic foraminifera indicate an upward shift from an outer-shelf to inner-shelf depositional environment. Interval II, deposited during the late early to early late Pliocene, contains reworked upper Miocene siliciclastics and faunas. In the upper Keys, quartz grains in Interval II range from very coarse sand that fines upward to very fine sand and then coarsens to very coarse and medium sand. In situ benthic faunas indicate an upward shift from outer-shelf to inner-shelf deposition. In the middle Keys, Interval II is different, with the quartz grains ranging primarily from medium to very fine sand. In situ benthic taxa indicate deposition on an inner shelf. In both the middle and upper Keys, the upper Pliocene siliciclastics of Interval III contain quartz grains ranging from very coarse to very fine sands that were deposited on an inner shelf. A sequence boundary between Interval I and Interval II is suggested by; an abrupt shift in the strontium-isotope chemostratigraphy; coarsening in quartz grain size above the boundary; an abrupt landward shift in depositional facies in the upper Keys core; and a distinct variation in the predominant coiling direction of the menardiiform planktic foraminifera, from fluctuating dextral-sinistral to dextral in the upper Keys core. Successive siliciclastic infilling, likely associated with eustatic sea-level change and current redeposition, formed a foundation for subsequent carbonate deposition. Deep-sea biostratigraphic techniques, integrated with ages derived from strontium-isotope chemostratigraphy, can be successfully applied to coastal-margin sequences, even though a depauperate suite of faunal markers is common.

Journal of Sedimentary Research↗
USGS1996

Crustal structure of the Colorado Plateau, Arizona: Application of new long-offset seismic data analysis techniques

The Colorado Plateau is a large crustal block in the southwestern United States that has been raised intact nearly 2 km above sea level since Cretaceous marine sediments were deposited on its surface. Controversy exists concerning the thickness of the plateau crust and the source of its buoyancy. Interpretations of seismic data collected on the plateau vary as to whether the crust is closer to 40 or 50 km thick. A thick crust could support the observed topography of the Colorado Plateau isostatically, while a thinner crust would indicate the presence of an underlying low-density mantle. This paper reports results on long-offset seismic data collected during the 1989 segment of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that extended from the Transition Zone into the Colorado Plateau in northwest Arizona. We apply two new methods to analyze long-offset data that employ finite difference travel time calculations: (1) a first-arrival time inverter to find upper crustal velocity structure and (2) a forward-modeling technique that allows the direct use of the inverted upper crustal solution in modeling secondary reflected arrivals. We find that the crustal thickness increases from 30 km beneath the metamorphic core complexes in the southern Basin and Range province to about 42 km beneath the northern Transition Zone and southern Colorado Plateau margin. We observe some crustal thinning (to ∼37 km thick) and slightly higher lower crustal velocities farther inboard; beneath the Kaibab uplift on the north rim of the Grand Canyon the crust thickens to a maximum of 48 km. We observe a nonuniform crustal thickness beneath the Colorado Plateau that varies by ∼15% and corresponds approximately to variations in topography with the thickest crust underlying the highest elevations. Crustal compositions (as inferred from seismic velocities) appear to be the same beneath the Colorado Plateau as those in the Basin and Range province to the southwest, implying that the plateau crust represents an unextended version of the Basin and Range. Some of the variability in crustal structure appears to correspond to preserved lithospheric discontinuities that date back to the Proterozoic Era.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Precise U‐Pb ages of Duluth Complex and related mafic intrusions, northeastern Minnesota: Geochronological insights to physical, petrogenetic, paleomagnetic, and tectonomagmatic processes associated with the 1.1 Ga Midcontinent Rift System

Precise resolution of the timing of igneous activity is crucial to understanding the dynamic processes associated with continental rifting. Mafic intrusive rocks constitute a major portion of the exposed 1100 Ma (Keweenawan) Midcontinent Rift system in northeastern Minnesota; however, prior to this study, geochronological data were insufficient to allow rigorous interpretation of intrusive histories and their relationships to extrusive suites. Eight anorthositic and gabbroic intrusives were chosen to represent both the temporal and spatial ranges of plutonic activity that formed the Duluth Complex and related mafic intrusions. U‐Pb isotopic analyses from zircons and baddeleyites result in U‐Pb concordant ages with little or no ambiguity introduced by inherited components, Pb loss or common Pb. The earliest Keweenawan plutonism exposed in Minnesota occurs along the northeastern flank of the Duluth Complex as a series of layered gabbros (Nathan's layered series) emplaced at 1106.9 ± 0.6 Ma. This sequence of gabbro sheets shares temporal, spatial, and compositional similarities with the nearby Logan sills in Ontario. Four Duluth Complex anorthositic and troctolitic series samples from widely separated areas have unresolvable ages between 1099.3 ± 0.3 and 1098.6 ± 0.5 Ma, indicating a very short duration for peak intrusive activity (0.5–1 m.y.). The unresolvable ages between anorthositic and troctolitic plutons suggest that these two magma series are more closely related than previously modeled and argue strongly for the need to reexamine their fundamental petrogenetic relationships. These dates also imply that the major reverse‐to‐normal magnetic polarity switch, used throughout the rift system as an important correlation tool, occurred prior to 1099 Ma. This date is several million years earlier than previously suspected and emphasizes the need for further paleomagnetic and geochronological data from the overlying volcanics. Much of the hypabyssal intrusive suite within the volcanic pile overlying Duluth Complex plutons may be significantly younger than the main pulse of plutonic activity. Two hypabyssal bodies, the Sonju Lake intrusion and gabbro at Silver Bay, were emplaced at 1096.1 ± 0.8 Ma and 1095.8 ± 1.2 Ma, respectively. Dates reported here and in previous studies support the concept of episodic tectonomagmatic rift development where magmatism was apparently concentrated in episodes of short duration (<1–3 m.y.) interspersed with longer hiatuses (2–8 m.y.).

Michigan, Minnesota, Ontario, Wisconsin↗
USGS2012

Active faulting and deformation of the Coalinga anticline as interpreted from three-dimensional velocity structure and seismicity

This work gives a clear picture of the geometry of aftershock seismicity in a large thrust earthquake. Interpretation of hypocenters and fault plane solutions, from the 1983 Coalinga, Coast Range California, earthquake sequence, in combination with the three-dimensional velocity structure shows that the active faulting beneath the fold primarily consists of a set of southwest dipping thrusts uplifting blocks of higher-velocity material. Above the main listric blind thrust there is a conjugate fault, steeply northeast dipping, that provides the western limit of the aftershocks within the Coalinga Anticline and that corresponds in location and spatial extent with the adjacent Pleasant Valley syncline. The character of the seismicity varies with the degree of previous deformation on each section of the anticline. Where the previous uplift was largest, the shallow seismicity shows secondary faulting on either side of the fold with orientations that correspond to the preexisting geologic structure. Diffuse seismicity characterizes the area with the least previous deformation. The mainshock rupture terminated where the fold trend was no longer uniform but had competing north and west trending features. The upward extent of the mainshock rupture ended at the approximate boundary between Franciscan and Great Valley Sequence rocks. Above that depth the main thrust appears to splay into a steeper segment and a near-horizontal segment. Thus the extent of rupture area is limited by the area of uniform structural orientation and by the variation in the type of material. With the three-dimensional velocity model each individual hypocenter moved slightly (0–2 km) in accord with the details of the surrounding velocity structure, so that secondary features in the seismicity pattern are more detailed than with a local one-dimensional model and station corrections. The overall character of the fault plane solutions was not altered by the three-dimensional model, but the more accurate ray paths did result in distinct changes. In particular, the mainshock has a fault plane dipping 30° southwest instead of the 23° obtained with the one-dimensional model.

Journal of Geophysical Research Solid Earth↗
USGS2012

Seismic refraction measurements within the Peninsular terrane, south central Alaska

We present an interpretation of crustal seismic refraction data from the Peninsular terrane, one of the many exotic terranes that have been accreted to the continental margin of southern Alaska in the past 200 m.y. A seismic refraction line was collected along the Glenn Highway in the Copper River Basin of south central Alaska in 1984 and 1985, as part of the U.S. Geological Survey Trans‐Alaska Crustal Transect (TACT) program. P wave velocities of 2.7–3.5 km/s and thicknesses of 1–2 km characterize post‐Lower Jurassic sedimentary rocks that underlie most of the seismic refraction line. An average crustal velocity structure includes the following five velocity divisions. Beneath the sedimentary rocks lie 1–2 km of 4.0–4.6 km/s materials, correlating with andesitic volcaniclastic sedimentary rocks and lava flows of the Lower Jurassic Talkeetna Formation. Below these rocks, seismic velocity increases rapidly, from 5.0 to 6.1 km/s, in 2–3 km. At 7–8 km depth, velocity jumps to 6.3 km/s and increments to 6.6 km/s by 10–12 km depth. Velocities increase from 6.8 to 7.0 km/s between 12 to 20 km depth. At about 22 km depth, a jump in velocity from 7.0 to 7.4 km/s is inferred but is poorly resolved. Depth to the Moho discontinuity could not be determined from our data. The absence of clear PmP reflections may indicate that Moho is deeper than 40 km. Data from two offset shotpoints northeast of the line and within the Wrangellia terrane constrain the deep structure transition between Peninsular and Wrangellia terranes. The 6.3–6.6 km/s material thickens to the northeast, toward the suture between Peninsular and Wrangellia terranes, but southwest of its mapped trace at the West Fork fault. Peninsular terrane crustal structure appears dissimilar to that of continental interiors. It is similar to velocity structures determined for accreted island arc fragments in California, such as the basement of the Great Valley and the Klamath Mountains.

Alaska↗
USGS2022

Mars Science Laboratory CheMin data from the Glen Torridon region and the significance of lake-groundwater interactions in interpreting mineralogy and sedimentary history

The Glen Torridon (GT) region is positioned in terrains with strong clay mineral signatures, as inferred from orbital spectroscopy. The GT campaign confirmed orbital distinctions with in situ measurements by the Mars Science Laboratory rover, Curiosity , and the CheMin X-ray diffraction instrument with of some of the highest clay mineral abundances to date. Additionally, GT is unique because of distinct phase identifications for the first time by CheMin, including: (i) Fe-carbonates, and (ii) a novel peak in the XRD patterns of some GT samples, with an interplanar spacing of at 9.2 Å. Fe-carbonates have been previously suggested from other instruments onboard, but this is the first definitive reporting by CheMin of multiple samples with Fe-carbonate. This new phase has never been observed in Gale crater with the CheMin instrument and may be a new mineral for Mars, but discrete identification still remains enigmatic because no single phase on Earth is able to account the mineralogical, geochemical, and sedimentological constraints in the GT region. Here, we modeled XRD profiles and propose an interstratified clay mineral, specifically greenalite-minnesotaite (G-M), as a reasonable candidate. The coexistence of Fe-carbonate and Fe-rich clay minerals in the GT samples, supports a conceptual model of a lacustrine groundwater mixing environment in the ancient Gale crater lake. Groundwater interaction with percolating lake waters in the sediments is common in terrestrial lacustrine settings, and the diffusion of two distinct water bodies within the subsurface can create a geochemical gradient and unique mineral front in the sediments. Ultimately, the proximity to this mixing zone controlled the secondary minerals preserved in the Jura, Knockfarril Hill, and Glasgow members of the Mt. Sharp group exposed in GT.

Journal of Geophysical Research E: Planets↗
USGS1975

Stratigraphy, conodont dating, and paleotectonic interpretation of the type Milligen Formation (Devonian), Wood River area, Idaho

The Milligen Formation at and near its type locality in the Wood River area is considerably older than and unrelated to rocks of Early Mississippian age called Milligen Formation in the Lost River Range and other ranges of east-central Idaho. Conodont faunas were found in limestones of a thin upper member of the sparsely fossiliferous marine Milligen Formation in its principal reference section at Milligen Gulch, at Fisher Canyon, and near Bellevue, Idaho. The faunas include indigenous conodonts here assigned to the early Late Devonian (early Frasnian) Lower and Middle Polygnathus asymmetricus Zones, and reworked conodonts derived from several Middle and Early Devonian conodont zones. An underlying much thicker argillite member of the Milligen contains fewer limestones, but a thin encrinite interbed near the middle of the member yielded early Middle Devonian (Eifelian) conodonts. This lower member probably represents most of Middle and Early Devonian time. Although its base is nowhere exposed in the Wood River area, the Milligen is inferred to have been deposited on the Silurian Trail Creek Formation, which crops out just to the east in the Pioneer Mountains. The age of the Milligen is therefore wholly Devonian and the highest fossiliferous beds are no younger than early Late Devonian. Reworked Middle and Early Devonian conodonts in limestone turbidites of the upper member of the Milligen Formation are identical to conodonts found in shelf (miogeosynclinal) carbonate rocks farther east. A postulated eastern source for the turbidites is supported by new data on the distribution, thickness, and tectonic facies of Devonian rocks that suggest the presence of a Late Devonian ridge on the continental shelf east of the Milligen depositional area. The Milligen Formation was intensely folded and was emergent during most of the Mississippian time when it formed part of the Antler Highlands, which shed flysch sediments eastward into the Copper basin. The Wood River Formation of Pennsylvanian and Permian age was then deposited over a subdued topography on the Milligen Formation. The Hailey Conglomerate Member at the base of the Wood River filled many irregularities in the surface. This depositional contact later was largely destroyed and the contact between the Milligen and Wood River is now a regional thrust fault at most places.

Idaho↗
USGS2012

Seismic guided waves trapped in the fault zone of the Landers, California, earthquake of 1992

A mobile seismic array of seven stations was deployed at 11 sites along the fault trace of the M 7.4 Landers earthquake of June 28, 1992, with a maximum offset of 1 km from the trace. We found a distinct wave train with a relatively long period following the S waves that shows up only when both the stations and the events are close to the fault trace. This wave train is interpreted as a seismic guided wave trapped in a low‐velocity fault zone. To study the distribution of amplitude of the guided waves with distance from the fault trace and also their attenuation with travel distance along the fault zone, we eliminated source and recording site effects by the coda normalization method. The normalized amplitudes of guided waves show a spectral peak at 3–4 Hz, which decays sharply with distance from the fault trace. Spectral amplitudes at high frequencies (8–15 Hz) show an opposite trend, increasing with distance from the fault trace. The normalized amplitudes of guided waves at 3–4 Hz also show a systematic decrease with hypocentral distance along the fault zone, from which we infer an apparent Q of 50. In order to confirm the existence of the guided waves, a dense array of 31 stations was deployed at one of the 11 sites. The resultant records revealed unequivocal evidence for the existence of guided waves associated with the fault zone. By modeling the waveforms as S waves trapped in a low‐velocity waveguide sandwiched between two homogeneous half‐spaces with velocity V s = 3.0 km/s, we infer a waveguide width of about 180 m, a shear velocity of 2.0–2.2 km/s, and a Q of ∼50. Hypocenters of aftershocks with clear guided waves show a systematic distribution both laterally and with depth delineating the extent of the low‐velocity fault zone in three dimensions. We find that the zone extends to a depth of at least 10 km. This zone apparently continues to the south across the Pinto Mountain fault because guided waves are observed at stations north of the Pinto Mountain fault for earthquakes with epicenters south of it. On the other hand, the zone appears to be discontinuous at the fault bend located about 20 km north of the mainshock epicenter; guided waves were observed for stations and epicenters which are located on the same sides of the fault bend but not for those on the opposite sides.

California↗
USGS2007

Earthquake-by-earthquake fold growth above the Puente Hills blind thrust fault, Los Angeles, California: Implications for fold kinematics and seismic hazard

Boreholes and high-resolution seismic reflection data collected across the forelimb growth triangle above the central segment of the Puente Hills thrust fault (PHT) beneath Los Angeles, California, provide a detailed record of incremental fold growth during large earthquakes on this major blind thrust fault. These data document fold growth within a discrete kink band that narrows upward from ???460 m at the base of the Quaternary section (200-250 m depth) to <150 m at 2.5 m depth, with most growth during the most recent folding event occurring within a zone only ???60 m wide. These observations, coupled with evidence from petroleum industry seismic reflection data, demonstrate that most (>82% at 250 m depth) folding and uplift occur within discrete kink bands, thereby enabling us to develop a paleoseismic history of the underlying blind thrust fault. The borehole data reveal that the youngest part of the growth triangle in the uppermost 20 m comprises three stratigraphically discrete growth intervals marked by southward thickening sedimentary strata that are separated by intervals in which sediments do not change thickness across the site. We interpret the intervals of growth as occurring after the formation of now-buried paleofold scarps during three large PHT earthquakes in the past 8 kyr. The intervening intervals of no growth record periods of structural quiescence and deposition at the regional, near-horizontal stream gradient at the study site. Minimum uplift in each of the scarp-forming events, which occurred at 0.2-2.2 ka (event Y), 3.0-6.3 ka (event X), and 6.6-8.1 ka (event W), ranged from ???1.1 to ???1.6 m, indicating minimum thrust displacements of ???2.5 to 4.5 m. Such large displacements are consistent with the occurrence of large-magnitude earthquakes (Mw > 7). Cumulative, minimum uplift in the past three events was 3.3 to 4.7 m, suggesting cumulative thrust displacement of ???7 to 10.5 m. These values yield a minimum Holocene slip rate for the PHT of ???0.9 to 1.6 mm/yr. The borehole and seismic reflection data demonstrate that dip within the kink band is acquired incrementally, such that older strata that have been deformed by more earthquakes dip more steeply than younger strata. Specifically, strata dip 0.4?? at 4 m depth, 0.7?? at 20 m depth, 8?? at 90 m, 16?? at 110 m, and 17?? at 200 m. Moreover, structural restorations of the borehole data show that the locus of active folding (the anticlinal active axial surface) does not extend to the surface in exactly the same location from earthquake to earthquake. Rather, that the axial surfaces migrate from earthquake to earthquake, reflecting a component of fold growth by kink band migration. The incremental acquisition of bed dip in the growth triangle may reflect some combination of fold growth by limb rotation in addition to kink band migration, possibly through a component of trishear or shear fault bend folding. Alternatively, the component of limb rotation may result from curved hinge fault bend folding, and/or the mechanical response of loosely consolidated granular sediments in the shallow subsurface to folding at depth. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2013

Merapi 2010 eruption—Chronology and extrusion rates monitored with satellite radar and used in eruption forecasting

Despite dense cloud cover, satellite-borne commercial Synthetic Aperture Radar (SAR) enabled frequent monitoring of Merapi volcano's 2010 eruption. Near-real-time interpretation of images derived from the amplitude of the SAR signals and timely delivery of these interpretations to those responsible for warnings, allowed satellite remote sensing for the first time to play an equal role with in situ seismic, geodetic and gas monitoring in guiding life-saving decisions during a major volcanic crisis. Our remotely sensed data provide an observational chronology for the main phase of the 2010 eruption, which lasted 12 days (26 October–7 November, 2010). Unlike the prolonged low-rate and relatively low explosivity dome-forming and collapse eruptions of recent decades at Merapi, the eruption began with an explosive eruption that produced a new summit crater on 26 October and was accompanied by an ash column and pyroclastic flows that extended 8 km down the flanks. This initial explosive event was followed by smaller explosive eruptions on 29 October–1 November, then by a period of rapid dome growth on 1–4 November, which produced a summit lava dome with a volume of ~ 5 × 10 6 m 3 . A paroxysmal VEI 4 magmatic eruption (with ash column to 17 km altitude) destroyed this dome, greatly enlarged the new summit crater and produced extensive pyroclastic flows (to ~ 16 km radial distance in the Gendol drainage) and surges during the night of 4–5 November. The paroxysmal eruption was followed by a period of jetting of gas and tephra and by a second short period (12 h) of rapid dome growth on 6 November. The eruption ended with low-level ash and steam emissions that buried the 6 November dome with tephra and continued at low levels until seismicity decreased to background levels by about 23 November. Our near-real-time commercial SAR documented the explosive events on 26 October and 4–5 November and high rates of dome growth (> 25 m 3 s − 1 ). An event tree analysis for the previous 2006 Merapi eruption indicated that for lava dome extrusion rates > 1.2 m 3 s − 1 , the probability of a large (1872-scale) eruption was ~ 10%. Consequently, the order-of-magnitude greater rates in 2010, along with the explosive start of the eruption on 26 October, the large volume of lava accumulating at the summit by 4 November, and the rapid and large increases in seismic energy release, deformation and gas emissions were the basis for warnings of an unusually large eruption by the Indonesian Geological Agency's Center for Volcanology and Geologic Hazard Mitigation (CVGHM) and their Volcano Research and Technology Development Center (BPPTK) in Yogyakarta — warnings that saved thousands of lives.

Merapi Volcano↗
USGS2016

The chemistry and isotopic composition of waters in the low-enthalpy geothermal system of Cimino-Vico Volcanic District, Italy

Geothermal energy exploration is based in part on interpretation of the chemistry, temperature, and discharge rate of thermal springs. Here we present the major element chemistry and the δD, δ 18 O, 87 Sr/ 86 Sr and δ 11 B isotopic ratio of groundwater from the low-enthalpy geothermal system near the city of Viterbo in the Cimino-Vico volcanic district of west-Central Italy. The geothermal system hosts many thermal springs and gas vents, but the resource is still unexploited. Water chemistry is controlled by mixing between low salinity,HCO 3 -rich fresh waters (<24.2°C) flowing in shallow volcanic rocks and SO 4 -rich thermal waters (25.3°C to 62.2°C) ascending from deep, high permeability Mesozoic limestones. The (equivalent) SO 4 /Cl (0.01–0.02), Na/Cl (2.82–5.83) and B/Cl ratios (0.02–0.38) of thermal waters differs from the ratios in other geothermal systems from Central Italy, probably implying a lack of hydraulic continuity across the region. The δ 18 O (−6.6‰ to −5.9‰) and δD (−40.60‰ to −36.30‰) isotopic composition of spring water suggest that the recharge area for the geothermal system is the summit region of Mount Cimino. The strontium isotope ratios ( 87 Sr/ 86 Sr) of thermal waters (0.70797–0.70805) are consistent with dissolution of the Mesozoic evaporite-carbonate units that constitute the reservoir, and the ratios of cold fresh waters mainly reflect shallow circulation through the volcanic cover and some minor admixture (<10%) of thermal water as well. The boron isotopic composition (δ 11 B) of fresh waters (−5.00 and 6.12‰) is similar to that of the volcanic cover, but the δ 11 B of thermal waters (−8.37‰ to −4.12‰) is a mismatch for the Mesozoic reservoir rocks and instead reflects dissolution of secondary boron minerals during fluid ascent through flysch units that overlie the reservoir. A slow and tortuous ascent enhances extraction of boron but also promotes conductive cooling, partially masking the heat present in the reservoir. Overall data from this study is consistent with previous studies that concluded that the geothermal system has a large energy potential.

Journal of Volcanology and Geothermal Research↗
USGS2008

Using amphibole phenocrysts to track vapor transfer during magma crystallization and transport: An example from Mount St. Helens, Washington

In order to evaluate and further constrain models for volatile movement and vapor enrichment of magma stored at shallow levels, amphibole phenocrysts from 2004–2005 Mount St. Helens dacite were analyzed for major and selected trace elements (Li, Cu, Zn, Mn, and REE) and Li isotopes. Several recent studies have examined fluid-mobile trace element abundances in phencryst phases and melt inclusions as a means of tracking volatile movement within subvolcanic magmatic systems, and high Li contents in plagioclase phenocrysts from 1980 and 2004 Mount St. Helens dacites have been interpreted as evidence that shallow magma was fluxed by a Li-bearing vapor phase prior to eruption. In amphibole phenocrysts, Zn and Mn behave compatibly, correlating to FeO ⁎ and Al 2 O 3 , and show no systematic change with time. In contrast, Li and Cu abundances in amphibole vary by up to 3 orders of magnitude (7.6–1140 μg/g and 1.7 to 94 μg/g, respectively), and do not generally correlate with either major or trace elements. However, they do correlate moderately well ( R 2 = 0.54, >> 95% confidence) with each other and show systematic temporal variations that are opposite to those observed for plagioclase, precluding a simple 1-step diffusion model for Li enrichment. We propose a Diffusion-Crystallization Multi-Stage (DCMS) model to explain the temporal variations and co-variations of Li and Cu. In early erupted dacite (October–December 2004) profiles of Li isotopes in conjunction with measured 7 Li intensities and core-to-rim increases in Li concentration are characteristic of Li diffusion into the amphiboles, consistent with prior models of plagioclase enrichment. In amphiboles from 2005 dacite, average Li and Cu concentrations are high (∼ 260–660 μg/g and ∼ 29–45 μg/g, respectively) and in contrast to amphiboles from earlier-erupted dacite, correlate weakly with Al 2 O 3 wt.%. Amphibole Al 2 O 3 concentrations are an indicator of pressure, with high-Al amphiboles crystallizing at higher pressures, and we suggest that Li and Cu are partitioned into a fluid phase during ascent and crystallization of the magma so that amphiboles crystallizing at lower pressure have correspondingly lower Li and Cu concentrations. However, low Li and Cu in amphiboles from the dacite at the start of the eruption also require crystallization from a low Li–Cu bearing melt or residence times long enough for amphiboles to re-equilibrate with a Li–Cu depleted melt. Estimated residence times suggest that amphiboles in early dacite could have been present since the end of the 1980–1986 eruptive episode at Mount St. Helens.

Journal of Volcanology and Geothermal Research↗
USGS2012

The Steens Mountain (Oregon) geomagnetic polarity transition, 2. Field intensity variations and discussion of reversal models

We carried out an extensive paleointensity study of the 15.5±0.3 m.y. Miocene reversed‐to‐normal polarity transition recorded in lava flows from Steens Mountain (south central Oregon). One hundred eighty‐five samples from the collection whose paleodirectional study is reported by Mankinen et al. (this issue) were chosen for paleointensity investigations because of their low viscosity index, high Curie point and reversibility, or near reversibility, of the strong field magnetization curve versus temperature. Application of the Thellier stepwise double heating method was very successful, yielding 157 usable paleointensity estimates corresponding to 73 distinct lava flows. After grouping successive lava flows that did not differ significantly in direction and intensity, we obtained 51 distinguishable, complete field vectors of which 10 are reversed, 28 are transitional, and 13 are normal. The record is complex, quite unlike that predicted by simple flooding or standing nondipole field models. It begins with an estimated several thousand years of reversed polarity with an average intensity of 31.5±8.5 μT, about one third lower than the expected Miocene intensity. This difference is interpreted as a long‐term reduction of the dipole moment prior to the reversal. When site directions and intensities are considered, truly transitional directions and intensities appear almost at the same time at the beginning of the transition, and they disappear simultaneously at the end of the reversal. Large deviations in declination occur during this approximately 4500±1000 year transition period that are compatible with roughly similar average magnitudes of zonal and nonzonal field components at the site. The transitional intensity is generally low, with an average of 10.9±4.9 μT for directions more than 45° away from the dipole field and a minimum of about 5 μT. The root‐mean‐square of the three field components X , Y , and Z are of the same order of magnitude for the transitional field and the historical nondipole field at the site latitude. However, a field intensity increase to pretransitional values occurs when the field temporarily reaches normal directions, which suggests that dipolar structure could have been briefly regenerated during the transition in an aborted attempt to reestablish a stationary field. Changes in the field vector are progressive but jerky, with at least two, and possibly three, large swings at astonishingly high rates. Each of those transitional geomagnetic impulses occurs when the field intensity is low (less than 10 μT) and is followed by an interval of directional stasis during which the magnitude of the field increases greatly. For the best documented geomagnetic impulse the rapid directional change corresponds to a vectorial intensity change of 6700±2700 nT yr −1 , which is about 15–50 times larger than the maximum rate of change of the nondipole field observed during the last centuries. The occurrence of geomagnetic impulses seems to support reversal models assuming an increase in the level of turbulence within the liquid core during transitions. The record closes with an estimated several thousand years of normal polarity with an average intensity of 46.7±20.1 μT, agreeing with the expected Miocene value. However, the occurrence of rather large and apparently rapid intensity fluctuations accompanied by little change in direction suggests that the newly reestablished dipole was still somewhat unstable.

Journal of Geophysical Research B: Solid Earth↗