USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Interpretation”

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 1,135 records · Page 63Linked to original sources
USGS2023

Precaldera mafic magmatism at Long Valley, California: Magma-tectonic siting and incubation of the Great Rhyolite System

The iconic volcanic center at Long Valley has released ∼820 km 3 of rhyolite in at least 110 eruptions. From 2.2 Ma until 0.23 Ma, products were exclusively rhyolitic, and ∼ 700 km 3 were high-silica rhyolite severely depleted in Sr, Ba, and Eu. The rhyolitic interval was preceded by an interval from 3.9 to 2.6 Ma with numerous basalt-andesite-dacite eruptions accompanied by no rhyolite at all. We have now mapped the circumcaldera products of this interval, defined 107 eruptive units, characterized them all chemically and petrographically, and dated many by 40 Ar/ 39 Ar. Here we display and describe them by sector around the caldera, interpret the nature of the transcrustal magma system that eventuated in the 35-km-wide Long Valley granite-rhyolite pluton, and analyze regional tectonic factors that did or did not contribute to siting the system. Nine Miocene (12–6 Ma) eruptive units close to Long Valley were followed by a Pliocene flare-up that released >300 mafic eruptions in a SW–NE swath 170 km long, centered across the later site of Long Valley. The basalts and their fractionates are intraplate alkalic products dominated by a continental lithosphere that had long been fluxed by Mesozoic subduction. Tertiary arc volcanism had not impinged on the area of Long Valley. Volumes estimated for the 107 Neogene precaldera eruptive units (only 40 of which exceeded 0.1 km 3 ) total ∼ 27 km 3 ± 50%—only ∼3% of the subsequent volume of rhyolite erupted. Such a volume of high-silica rhyolite with ultra-low Sr and Eu is not a product of partial melting but requires as proximate parent a leucogranitic crystal mush that is itself the upper level of a long-lived plutonic reservoir that extends to the lower crust. The 27 km 3 of Neogene magma that erupted was a small contingent of the mantle-derived basaltic flux needed to energize (and contribute its fractionated melt to) a 30-km-deep compositionally graded crustal column, which culminated in ∼10,000 km 3 of granitoid mush from which 820 km 3 of Quaternary high-silica rhyolitic melt escaped and erupted. Pliocene basaltic eruptions ceased at ∼2.6 Ma, probably because the basaltic flux intensified sufficiently to render the mushy upper crust impenetrable. The 2.6–2.2 Ma quiescent interval represented culmination of thermal activation of the plutonic column and refinement of its leucogranitic mushy upper layer, from which extreme melts escaped for the next 2 Myr. The Pliocene mafic swath crosses the Sierran rangefront fault zone coincident with a left-stepping extensional reentrant that also began developing at ∼3 Ma. Moreover, Long Valley overlies a dextral offset in the initial Sr-isotope 0.706 line, which may reflect the rifted or attenuated edge of Proterozoic crust and mantle lithosphere. Concatenation of these three influences may account for siting of intensified edge-focused magmatism that produced the great Quaternary pluton.

California↗
USGS2012

Earthquake depths and the relation to strain accumulation and stress near strike-slip faults in southern California

In this paper I report and attempt to interpret several observations about the depth distribution of earthquakes in the region of southern California containing the southern San Andreas and San Jacinto fault zones. These observations are as follows: (1) Earthquakes in the major fault zones are predominantly deep (maximum, 22 km; characteristic, 11–18 km). (2) Earthquakes in the crustal blocks bounding the fault zones are predominantly shallow (maximum, 17 km; characteristic, 1–6 km). (3) In the San Jacinto fault zone, maximum earthquake depths correlate with surface heat flow. These relations together with focal mechanisms, geodetic strain measurements, and fault zone models are consistent with the following ideas: (1) Interseismic plate motion is accommodated by aseismic slip along an extension of the major fault zone below a brittle zone that is locked between large earthquakes. (2) The aseismic slip in a narrow fault zone in the brittle-plastic transition region concentrates strain at the base of the brittle fault zone. (3) Deep earthquakes occur in the lower part of the brittle fault zone due to stick-slip failure of highly stressed patches. (4) Background earthquakes and aftershocks that occur several kilometers deeper than large earthquake hypocenters suggest that a zone of mixed slip behavior may exist between the stable sliding (deep) and stick-slip (shallow) regions of the fault zone. This zone of mixed slip may fail predominantly by stable sliding, analogous to the central San Andreas and Calaveras faults of California. Seismic precursors to large earthquakes may occur in this zone, and large earthquakes may nucleate near the upper edge of this zone. Furthermore, the difference in seismicity between the San Jacinto and southern San Andreas faults suggests that the nature of this mixed zone may evolve as total displacement in the fault zone increases. (5) Shear stress may be less in the crustal blocks than in the deep brittle fault zones and generally at a level sufficient to cause brittle failure only shallow in the crustal blocks. (6) In the stress field produced by plate motion and slip in the deep fault zone, the upper brittle fault zone is not oriented favorably for shear failure. Lack of shallow earthquakes in the fault zones and the predominance of shallow earthquakes on favorably oriented fractures in the adjacent crustal blocks suggest that either stress in the upper brittle fault zone is relatively low or the upper fault zone is effectively strong due to its orientation.

California↗
USGS2010

Seismicity and fluid geochemistry at Lassen Volcanic National Park, California: Evidence for two circulation cells in the hydrothermal system

Seismic analysis and geochemical interpretations provide evidence that two separate hydrothermal cells circulate within the greater Lassen hydrothermal system. One cell originates south to SW of Lassen Peak and within the Brokeoff Volcano depression where it forms a reservoir of hot fluid (235–270°C) that boils to feed steam to the high-temperature fumarolic areas, and has a plume of degassed reservoir liquid that flows southward to emerge at Growler and Morgan Hot Springs. The second cell originates SSE to SE of Lassen Peak and flows southeastward along inferred faults of the Walker Lane belt (WLB) where it forms a reservoir of hot fluid (220–240°C) that boils beneath Devils Kitchen and Boiling Springs Lake, and has an outflow plume of degassed liquid that boils again beneath Terminal Geyser. Three distinct seismogenic zones (identified as the West , Middle , and East seismic clusters) occur at shallow depths (<6 km) in Lassen Volcanic National Park, SW to SSE of Lassen Peak and adjacent to areas of high-temperature (≤161°C) fumarolic activity (Sulphur Works, Pilot Pinnacle, Little Hot Springs Valley, and Bumpass Hell) and an area of cold, weak gas emissions (Cold Boiling Lake). The three zones are located within the inferred Rockland caldera in response to interactions between deeply circulating meteoric water and hot brittle rock that overlies residual magma associated with the Lassen Volcanic Center. Earthquake focal mechanisms and stress inversions indicate primarily N–S oriented normal faulting and E–W extension, with some oblique faulting and right lateral shear in the East cluster. The different focal mechanisms as well as spatial and temporal earthquake patterns for the East cluster indicate a greater influence by regional tectonics and inferred faults within the WLB. A fourth, deeper (5–10 km) seismogenic zone (the Devils Kitchen seismic cluster) occurs SE of the East cluster and trends NNW from Sifford Mountain toward the Devils Kitchen thermal area where fumarolic temperatures are ≤123°C. Lassen fumaroles discharge geothermal gases that indicate mixing between a N 2 -rich, arc-type component and gases derived from air-saturated meteoric recharge water. Most gases have relatively weak isotopic indicators of upper mantle or volcanic components, except for gas from Sulphur Works where δ 13 C–CO 2 , δ 34 S–H 2 S, and δ 15 N–N 2 values indicate a contribution from the mantle and a subducted sediment source in an arc volcanic setting.

California↗
USGS2000

Sea-floor environments within Long Island Sound: A regional overview

Modern sea-floor sedimentary environments within the glaciated, topographically complex Long Island Sound estuary have been interpreted and mapped from an extensive collection of sidescan sonographs, bottom samples, and video-camera observations together with supplemental bathymetric, marine-geologic, and bottom-current data. Four categories of environments are present that reflect the dominant long-term processes of erosion or nondeposition; coarsegrained bedload transport; sediment sorting and reworking; and fine-grained deposition. (1) Environments of erosion or nondeposition contain exposures of glacial drift, coarse lag deposits, and possibly bedrock and include sediments which range from boulder fields to gravelly coarse-to-medium sands. (2) Environments of coarse-grained bedload transport are mantled by sand ribbons and sand waves and contain mostly coarse-to-fine sands with only small amounts of mud. (3) Environments of sediment sorting and reworking comprise both uniform and heterogeneous sediment types and contain variable amounts of fine sand and mud. (4) Environments of fine-grained deposition are blanketed by muds and sandy muds. The patchy distribution of sedimentary environments within Long Island Sound reflects both regional and local changes in bottom processes. Regional changes are primarily the result of a strong, east-to-west decreasing gradient of bottom tidal-current speeds, coupled with the net (westward) estuarine bottom drift. The regional current regime has produced a westward succession of environments along the basin floor beginning with erosion or nondeposition at the narrow eastern entrance to the Sound, changing to an extensive area of coarse-grained bedload transport, passing into a contiguous band of sediment sorting, and ending with broad areas of fine-grained deposition in the central and western Sound. However, local changes in processes are superimposed on the regional conditions within the central and western parts of the basin and along the nearshore margins. Within the central and western basin, localized sedimentary environments are produced where the bottom flow is enhanced by, and interacts with, the bottom topography, whereas along the nearshore margins, they variously reflect wave-produced currents, the irregular bathymetry, the indented shoreline, and the proximal supply of sediments. Results from this study (1) confirm the high trapping efficiency of fine-grained sediments in the Sound, (2) suggest that fine-grained sediments accumulate at an average (regional) rate of 0.08 g/cm 2 /y , and (3) indicate that the postglacial delta in the eastern Sound was a significant source of fine-grained sediments now buried beneath depositional areas.

Connecticut, New York↗
USGS2008

Methanogenic pathways of coal-bed gas in the Powder River Basin, United States: The geologic factor

Coal-bed gas of the Tertiary Fort Union and Wasatch Formations in the Powder River Basin in Wyoming and Montana, U.S. was interpreted as microbial in origin by previous studies based on limited data on the gas and water composition and isotopes associated with the coal beds. To fully evaluate the microbial origin of the gas and mechanisms of methane generation, additional data for 165 gas and water samples from 7 different coal-bed methane-bearing coal-bed reservoirs were collected basinwide and correlated to the coal geology and stratigraphy. The C1/(C2 + C3) ratio and vitrinite reflectance of coal and organic shale permitted differentiation between microbial gas and transitional thermogenic gas in the central part of the basin. Analyses of methane ??13C and ??D, carbon dioxide ??13C, and water ??D values indicate gas was generated primarily from microbial CO2 reduction, but with significant gas generated by microbial methyl-type fermentation (aceticlastic) in some areas of the basin. Microbial CO2 reduction occurs basinwide, but is generally dominant in Paleocene Fort Union Formation coals in the central part of the basin, whereas microbial methyl-type fermentation is common along the northwest and east margins. Isotopically light methane ??13C is distributed along the basin margins where ??D is also depleted, indicating that both CO2-reduction and methyl-type fermentation pathways played major roles in gas generation, but gas from the latter pathway overprinted gas from the former pathway. More specifically, along the northwest basin margin gas generation by methyl-type fermentation may have been stimulated by late-stage infiltration of groundwater recharge from clinker areas, which flowed through highly fractured and faulted coal aquifers. Also, groundwater recharge controlled a change in gas composition in the shallow Eocene Wasatch Formation with the increase of nitrogen and decrease of methane composition of the coal-bed gas. Other geologic factors, such as burial, thermal and maturation history, lateral and vertical continuity, and coalification of the coal beds, also played a significant role in controlling methanogenic pathways and provided new perspectives on gas evolution and emplacement. The early-stage gas produced by CO2 reduction has mixed with transitional thermogenic gas in the deeper, central parts of the Powder River Basin to form 'old' gas, whereas along the basin margins the overprint of gas from methyl-type fermentation represents 'new' gas. Thus, a clear understanding of these geologic factors is necessary to relate the microbiological, biogeochemical, and hydrological processes involved in the generation of coal-bed gas.

International Journal of Coal Geology↗
USGS2017

Magmatic degassing, lava dome extrusion, and explosions from Mount Cleveland volcano, Alaska, 2011–2015: Insight into the continuous nature of volcanic activity over multi-year timescales

Mount Cleveland volcano (1730 m) is one of the most active volcanoes in the Aleutian arc, Alaska, but heightened activity is rarely accompanied by geophysical signals, which makes interpretation of the activity difficult. In this study, we combine volcanic gas emissions measured for the first time in August 2015 with longer-term measurements of thermal output and lava extrusion rates between 2011 and 2015 calculated from MODIS satellite data with the aim to develop a better understanding of the nature of volcanic activity at Mount Cleveland. Degassing measurements were made in the month following two explosive events (21 July and 7 August 2015) and during a period of new dome growth in the summit crater. SO 2 emission rates ranged from 400 to 860 t d − 1 and CO 2 /SO 2 ratios were < 3, consistent with the presence of shallow magma in the conduit and the observed growth of a new lava dome. Thermal anomalies derived from MODIS data from 2011 to 2015 had an average repose time of only 4 days, pointing to the continuous nature of volcanic activity at this volcano. Rapid increases in the cumulative thermal output were often coincident with visual confirmation of dome growth or accumulations of tephra in the crater. The average rate of lava extrusion calculated for 9 periods of rapid increase in thermal output was 0.28 m 3 s − 1 , and the total volume extruded from 2011 to 2015 was 1.9–5.8 Mm 3 . The thermal output from the lava extrusion events only accounts for roughly half of the thermal budget, suggesting a continued presence of shallow magma in the upper conduit, likely driven by convection. Axisymmetric dome morphology and occasional drain back of lava into the conduit suggests low-viscosity magmas drive volcanism at Mount Cleveland. It follows also that only small overpressures can be maintained given the small domes and fluid magmas, which is consistent with the low explosivity of most of Mount Cleveland's eruptions. Changes between phases of dome growth and explosive activity are somewhat unpredictable and likely result from plugs that are related to the dome obtaining a critical dimension, or from small variations in the magma ascent rate that lead to crystallization-induced blockages in the upper conduit, thereby reducing the ability of magma to degas. We suggest the small magma volumes, slow ascent rates, and low magma viscosity lead to the overall lack of anomalous geophysical signals prior to eruptions, and that more continuous volcanic degassing measurements might lead to more successful eruption forecasting at this continuously-active open-vent volcano.

Alaska↗
USGS2010

Depositional setting, petrology and chemistry of Permian coals from the Paraná Basin: 2. South Santa Catarina Coalfield, Brazil

In Brazil economically important coal deposits occur in the southern part of the Paraná Basin, where coal seams occur in the Permian Rio Bonito Formation, with major coal development in the states of Rio Grande de Sul and Santa Catarina. The current paper presents results on sequence stratigraphic interpretation of the coal-bearing strata , and petrological and geochemical coal seam characterization from the South Santa Catarina Coalfield, Paraná Basin. In terms of sequence stratigraphic interpretation the precursor mires of the Santa Catarina coal seams formed in an estuarine-barrier shoreface depositional environment , with major peat accumulation in a high stand systems tract (Pre-Bonito and Bonito seams), a lowstand systems tract (Ponta Alta seam, seam A, seam B) and a transgressive systems tract (Irapuá, Barro Branco and Treviso seams). Seam thicknesses range from 1.70 to 2.39 m, but high proportions of impure coal (coaly shale and shaley coal), carbonaceous shale and partings reduce the net coal thickness significantly. Coal lithoypes are variable, with banded coal predominant in the Barro Branco seam, and banded dull and dull coal predominantly in Bonito and Irapuá seams, respectively. Results from petrographic analyses indicate a vitrinite reflectance range from 0.76 to 1.63 %Rrandom (HVB A to LVB coal). Maceral group distribution varies significantly, with the Barro Branco seam having the highest vitrinite content (mean 67.5 vol%), whereas the Irapuá seam has the highest inertinite content (33.8 vol%). Liptinite mean values range from 7.8 vol% (Barro Branco seam) to 22.5 vol% (Irapuá seam). Results from proximate analyses indicate for the three seams high ash yields (50.2 – 64.2 wt.%). Considering the International Classification of in-Seam Coals, all samples are in fact classified as carbonaceous rocks (> 50 wt.% ash). Sulfur contents range from 3.4 to 7.7 wt.%, of which the major part occurs as pyritic sulfur. Results of X-ray diffraction indicate the predominance of quartz and kaolinite (also pyrite). Gypsum , gibbsite , jarosite and calcite were also identified in some samples. Feldspar was noted but is rare. The major element distribution in the three seams (coal basis) is dominated by SiO 2 (31.3 wt.%, mean value), Al 2 O 3 (14.5 wt.%, mean value) and Fe 2 O 3 (6.9 wt.%, mean value). Considering the concentrations of trace elements that are of potential environmental hazards the Barro Branco, Bonito and Irapuá seams (coal base) are significantly enriched in Co (15.7 ppm), Cr (54.5 ppm), Li (59.3 ppm), Mn (150.4 ppm), Pb (58.0 ppm) and V (99.6 ppm), when compared to average trace elements contents reported for U. S. coals. Hierarchical cluster analysis identified, based on similarity levels, three groups of major elements and seven groups of trace elements. Applying discriminant analyses using trace and major element distribution, it could be demonstrated that the three seams from Santa Catarina show distinct populations in the discriminant analyses plots, and also differ from the coals of Rio Grande do Sul analyzed in a previous study.

Santa Catarina↗
USGS2016

RMT focal plane sensitivity to seismic network geometry and faulting style

Modern tectonic studies often use regional moment tensors (RMTs) to interpret the seismotectonic framework of an earthquake or earthquake sequence; however, despite extensive use, little existing work addresses RMT parameter uncertainty. Here, we quantify how network geometry and faulting style affect RMT sensitivity. We examine how data-model fits change with fault plane geometry (strike and dip) for varying station configurations. We calculate the relative data fit for incrementally varying geometries about a best-fitting solution, applying our workflow to real and synthetic seismograms for both real and hypothetical station distributions and earthquakes. Initially, we conduct purely observational tests, computing RMTs from synthetic seismograms for hypothetical earthquakes and a series of well-behaved network geometries. We then incorporate real data and station distributions from the International Maule Aftershock Deployment (IMAD), which recorded aftershocks of the 2010 M W 8.8 Maule earthquake, and a set of regional stations capturing the ongoing earthquake sequence in Oklahoma and southern Kansas. We consider RMTs computed under three scenarios: (1) real seismic records selected for high data quality; (2) synthetic seismic records with noise computed for the observed source-station pairings and (3) synthetic seismic records with noise computed for all possible station-source pairings. To assess RMT sensitivity for each test, we observe the ‘fit falloff’, which portrays how relative fit changes when strike or dip varies incrementally; we then derive the ranges of acceptable strikes and dips by identifying the span of solutions with relative fits larger than 90 per cent of the best fit. For the azimuthally incomplete IMAD network, Scenario 3 best constrains fault geometry, with average ranges of 45° and 31° for strike and dip, respectively. In Oklahoma, Scenario 3 best constrains fault dip with an average range of 46°; however, strike is best constrained by Scenario 1, with a range of 26°. We draw two main conclusions from this study. (1) Station distribution impacts our ability to constrain RMTs using waveform time-series; however, in some tectonic settings, faulting style also plays a significant role and (2) increasing station density and data quantity (both the number of stations and the number of individual channels) does not necessarily improve RMT constraint. These results may be useful when organizing future seismic deployments (e.g. by concentrating stations in alignment with anticipated nodal planes), and in computing RMTs, either by guiding a more rigorous data selection process for input data or informing variable weighting among the selected data (e.g. by eliminating the transverse component when strike-slip mechanisms are expected).

Geophysical Journal International↗
USGS1989

Multiple hydrothermal and metamorphic events in the Kidd Creek volcanogenic massive sulphide deposit, Timmins, Ontario: evidence from tourmalines and chlorites

Tourmaline and chlorite are the principal ferromagnesian silicate minerals in the Kidd Creek massive sulphide deposit. Tourmaline is most common in sphalerite-rich peripheral margins of the chalcopyrite stringer zone. Within the north orebody, samples typically contain <1% tourmaline, but small areas (hand-specimen scale) may have 10–20%. Chlorite is more widely distributed and in places constitutes 30–50% of rock volumes. Associated assemblages may include quartz, sulphides (principally chalcopyrite, sphalerite, and (or) pyrite), carbonate, albite, sericite, and rare fluorite, allanite, or zoisite(?).The tourmalines and chlorites record a series of multiple hydrothermal and metamorphic events. Paragenetic studies suggest that tourmaline was deposited during several discrete stages of mineralization, as evidenced by brecciation and cross-cutting relationships. Most of the tourmalines have two concentric growth zones defined by different colours (green, brown, blue, yellow). Some tourmalines also display pale discordant rims that cross-cut and embay the inner growth zones and polycrystalline, multiple-extinction domains. Late sulphide veinlets (chalcopyrite, pyrrhotite) transect the inner growth zones and pale discordant rims of many crystals. The concentric growth zones are interpreted as primary features developed by the main ore-forming hydrothermal system, whereas the discordant rims, polycrystalline domains, and cross-cutting sulphide veinlets reflect post-ore metamorphic processes.Detailed electron microprobe analyses of tourmalines show a wide compositional range, from Fe-rich dravite nearly to end-member schorl, with Fe/(Fe + Mg) ratios varying from 0.33 to 0.92; only minor amounts of Ca are present, yielding uniformly high Na/(Na + Ca) ratios of 0.84–0.99. Two sets of chemical zoning trends are identified in the tourmalines, involving systematic changes in Fe/(Fe + Mg), Na/(Na + Ca), Al, and Ti that are believed to reflect internal coupled substitutions (e.g.,  + Ti = Na + Al) and local mineral equilibria (e.g., tourmaline–chlorite). Analyses of the pale discordant reaction rims show consistent depletion of Fe, Ca, and Ti, presumably by fluid–solid reactions during post-ore metamorphism.Chlorites also show an extensive range in composition, from ripidolite nearly to end-member daphnite, with Fe/(Fe + Mg) ratios of 0.43–0.98 and Si cation values of 5.00–5.39. Chlorites from the fringes of the footwall stringer zone have narrow compositional ranges, whereas chlorites near footwall rhyolite sills in the core of the stringer zone display major variations in Fe/(Fe + Mg) ratios, including one sample with a range of 0.68–0.95. The former group of chlorites has Fe/(Fe + Mg) ratios that correlate well with those of coexisting tourmalines (exclusive of late reaction rims). Data for the latter group, in contrast, fall off equilibrium K D curves, indicating that the tourmalines and chlorites within these samples are not in chemical equilibrium. The chlorites are believed to have been altered (overprinted) by Fe-rich hydrothermal fluids apparently generated during intrusion of the rhyolite sills. The tourmalines, however, are unaffected and retain primary chemical signatures.Variations in mineral proportions and mineral chemistry within the deposit mainly depend on fluctuations in temperature, pH, water/rock ratios, and amounts of entrained seawater. The major proposed control is mixing between high-temperature, Fe-rich end-member hydrothermal fluids and cold, Mg-rich entrained seawater. Fe/(Fe + Mg) variations in footwall tourmalines (and equilibrium chlorites) are believed to largely reflect the progressive infiltration of Mg-rich seawater into the margins and top of the hydrothermal system. The more Fe-rich compositions of Kidd Creek tourmalines relative to those from sediment-hosted massive sulphide deposits (e.g., Sullivan, British Columbia) may be related to the preferential generation of end-member hydrothermal fluids in proximal volcanic environments like that at Kidd Creek.

Ontario↗
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↗
USGS2007

Fault locking, block rotation and crustal deformation in the Pacific Northwest

We interpret Global Positioning System (GPS) measurements in the northwestern United States and adjacent parts of western Canada to describe relative motions of crustal blocks, locking on faults and permanent deformation associated with convergence between the Juan de Fuca and North American plates. To estimate angular velocities of the oceanic Juan de Fuca and Explorer plates and several continental crustal blocks, we invert the GPS velocities together with seafloor spreading rates, earthquake slip vector azimuths and fault slip azimuths and rates. We also determine the degree to which faults are either creeping aseismically or, alternatively, locked on the block-bounding faults. The Cascadia subduction thrust is locked mainly offshore, except in central Oregon, where locking extends inland. Most of Oregon and southwest Washington rotate clockwise relative to North America at rates of 0.4–1.0 ° Myr −1 . No shear or extension along the Cascades volcanic arc has occurred at the mm/yr level during the past decade, suggesting that the shear deformation extending northward from the Walker Lane and eastern California shear zone south of Oregon is largely accommodated by block rotation in Oregon. The general agreement of vertical axis rotation rates derived from GPS velocities with those estimated from palaeomagnetic declination anomalies suggests that the rotations have been relatively steady for 10–15 Ma. Additional permanent dextral shear is indicated within the Oregon Coast Range near the coast. Block rotations in the Pacific Northwest do not result in net westward flux of crustal material—the crust is simply spinning and not escaping. On Vancouver Island, where the convergence obliquity is less than in Oregon and Washington, the contractional strain at the coast is more aligned with Juan de Fuca—North America motion. GPS velocities are fit significantly better when Vancouver Island and the southern Coast Mountains move relative to North America in a block-like fashion. The relative motions of the Oregon, western Washington and Vancouver Island crustal blocks indicate that the rate of permanent shortening, the type that causes upper plate earthquakes, across the Puget Sound region is 4.4 ± 0.3 mm yr −1 . This shortening is likely distributed over several faults but GPS data alone cannot determine the partitioning of slip on them. The transition from predominantly shear deformation within the continent south of the Mendocino Triple Junction to predominantly block rotations north of it is similar to changes in tectonic style at other transitions from shear to subduction. This similarity suggests that crustal block rotations are enhanced in the vicinity of subduction zones possibly due to lower resisting stress.

Geophysical Journal International↗
USGS2001

Anomalously high b-values in the South Flank of Kilauea volcano, Hawaii: Evidence for the distribution of magma below Kilauea's East rift zone

The pattern of b-value of the frequency-magnitude relation, or mean magnitude, varies little in the Kaoiki-Hilea area of Hawaii, and the b-values are normal, with b = 0.8 in the top 10 km and somewhat lower values below that depth. We interpret the Kaoiki-Hilea area as relatively stable, normal Hawaiian crust. In contrast, the b-values beneath Kilauea's South Flank are anomalously high (b = 1.3-1.7) at depths between 4 and 8 km, with the highest values near the East Rift zone, but extending 5-8 km away from the rift. Also, the anomalously high b-values vary along strike, parallel to the rift zone. The highest b-values are observed near Hiiaka and Pauahi craters at the bend in the rift, the next highest are near Makaopuhi and also near Puu Kaliu. The mildest anomalies occur adjacent to the central section of the rift. The locations of the three major and two minor b-value anomalies correspond to places where shallow magma reservoirs have been proposed based on analyses of seismicity, geodetic data and differentiated lava chemistry. The existence of the magma reservoirs is also supported by magnetic anomalies, which may be areas of dike concentration, and self-potential anomalies, which are areas of thermal upwelling above a hot source. The simplest explanation of these anomalously high b-values is that they are due to the presence of active magma bodies beneath the East Rift zone at depths down to 8 km. In other volcanoes, anomalously high b-values correlate with volumes adjacent to active magma chambers. This supports a model of a magma body beneath the East Rift zone, which may widen and thin along strike, and which may reach 8 km depth and extend from Kilauea's summit to a distance of at least 40 km down rift. The anomalously high b-values at the center of the South Flank, several kilometers away from the rift, may be explained by unusually high pore pressure throughout the South Flank, or by anomalously strong heterogeneity due to extensive cracking, or by both phenomena. The major b-value anomalies are located SSE of their parent reservoirs, in the direction of motion of the flank, suggesting that magma reservoirs leave an imprint in the mobile flank. We hypothesize that the extensive cracking may have been acquired when the anomalous parts of the South Flank, now several kilometers distant from the rift zone, were generated at the rift zone near persistent reservoirs. Since their generation, these volumes may have moved seaward, away from the rift, but earthquakes occurring in them still use the preexisting complex crack distribution. Along the decollement plane at 10 km depth, the b-values are exceptionally low (b = 0.5), suggesting faulting in a more homogeneous medium. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗
USGS2008

Trace fossils and sedimentary facies from a Late Cambrian‐Early Ordovician tide‐dominated shelf (Santa Rosita Formation, northwest Argentina): Implications for ichnofacies models of shallow marine successions

The Santa Rosita Formation is one the most widely distributed lower Paleozoic units of northwest Argentina. At the Quebrada del Salto Alto section, east of Purmamarca, Jujuy Province, it is represented by four sedimentary facies: thick‐bedded planar cross‐stratified quartzose sandstones (A), thin‐bedded planar cross‐stratified quartzose sandstones and mudstones (B), wave‐rippled sandstones and bioturbated mudstones (C), and black and greenish gray shales (D). Paleocurrent data, sandstone architecture, and sedimentary structures from faciès A and B indicate bipolar/bimodal paleoflows, suggesting the action of tidal currents. The succession is interpreted as that of a tide‐dominated shelf, with only secondary influence of wave processes. Trace fossils are restricted to facies B and C. The Cruziana ichnocoenosis is preserved on the soles of thin‐bedded planar cross‐stratified quartzose sandstones (faciès B). This ichnocoenosis consists of Conostichus isp., Cruziana omanica, C. semiplicata, C. cf. tortworthi, Cruziana isp. Helminthopsis abeli, Monomorphichnus bilinearis, M. multilineatus, Palaeophycus tubularis, Rusophycus carbonarias, R. latus , and R. isp. The occurrence of Cruziana semiplicata , C. omanica, C. cf. tortworthi , and Rusophycus latus supports a Late Cambrian‐Tremadoc age. Slabbing of Cruziana shows complex interactions between biologic and sedimentologic processes, and suggests a predominance of exhumed traces, washed out and recast by tractive sand deposition. Sandstone soles are densely packed with biogenic structures and exhibit distinctive clusters of Rusophycus isp. that most likely represent trilobite nesting burrows. The Cruziana ichnocoenosis records the resident fauna of a protected, lower intertidal to subtidal interbar setting. The Skolithos ichnocoenosis is represented by high to low density vertical burrows of Skolithos linearis , which extend downwards to the quartzose sandstone soles of faciès B and cross the Cruziana ichnocoenosis. The Skolithos ichnocoenosis represents colonization by suspension‐feeding organisms following a major change in environmental conditions, related to the migration of lower intertidal to subtidal sandwaves. The Planolites ichnocoenosis consists exclusively of Planolites montanus within mudstones overlying wave‐rippled sandstones (facies C). The Planolites ichnocoenosis records opportunistic colonization by inf aunal deposit feeders that mined the organic‐rich fine‐grained sediment during the waning phase of storms that scoured organic detritus from the sea bottom. The section records, from base to top, a Cruziana‐Skolithos ichnofacies zone, a Skolithos ichnofacies zone and an unbioturbated zone typified by the thick‐bedded cross‐stratified quartzose sandstone (fades A). This trend reflects progressively higher energy conditions linked to the establishment of a large sand wave complex. The presence of a mixed Cruziana‐Skolithos ichnofacies in the lower interval reflects changes in substrate and energy levels, rather than water depth. Accordingly, contrasting ichnocoenoses from interbars (Cruziana) and sandwaves (Skolithos) must be considered an example of ichnofacies controlled by local parameters instead of general bathymétrie trends. Conversely, the vertical replacement of the Cruziana ichnofacies by the Skolithos ichnofacies towards the middle interval of the section reflects the environmental changes associated with the transition between the intertidal and subtidal zones. As overall tidal energy increases from supratidal to subtidal settings, the Skolithos ichnofacies tends to occur seaward of the Cruziana ichnofacies in tide‐dominated shallow marine environments. Therefore, onshore‐offshore ichnofacies replacement in tide‐dominated shallow seas is opposite to that in wave‐dominated marine settings.

Ichnos: An International Journal for Plant and Ani↗
USGS2012

Basin and range crustal and upper mantle structure, northwest to central Nevada

We present an interpretation of the crustal and uppermost mantle structure of the Basin and Range of northwestern Nevada based on seismic refraction/wide‐angle reflection, near‐vertical reflection, and gravity data. In comparison to most previous estimates, we find that the crust is somewhat thicker (32–36 km versus 22–30 km), and the uppermost mantle velocity is somewhat higher (8.0 km/s versus 7.3–7.9 km/s). Along our transects, the crust is thinnest (32 km) in the Carson Sink‐Buena Vista Valley region and increases by 2–4 km to the west and east, respectively. There is considerable complexity throughout the crust where velocities range from of 2.5 km/s at the surface to 7.4 km/s in the lowermost crust. Variations in velocity and structure of the upper crustal layers reveal apparent basement velocity depressions (areas of lower velocities extending up to 10 km in depth) that underlie some surface ranges as well as the basins. The middle crust rises from about 20 km beneath central Nevada to within 12 km of the surface beneath the area of thinnest crust and is characterized by a modest (∼0.1 km/s) change in velocity and low‐velocity gradients. These midcrustal layers mark the onset of high crustal reflectivity and the apparent limiting depth to which Basin and Range faults can be traced in near‐vertical reflection profiles, suggesting that these midcrustal layers represent the transition between the brittle and ductile zones of the crust. The lower crust is more structurally complex, with layers thickening and thinning in a systematic manner with the upper crustal layers; generally, where there are velocity depressions in the upper crust, the lower crust is thickest and shallowest. The geometry of these lower crustal layers (derived from refraction modeling) coincides with changes in the crustal reflectivity, determined from the Consortium of Continental Reflection Profiling reflection data. The lower crustal layer is unusually high in velocity (7.4 km/s) and is likely the layer identified as mantle in some previous studies. We do not identify the 7.4 km/s layer as mantle because (1) there is an underlying layer with a velocity (8.0 km/s) that is more consistent with the worldwide average velocity for the upper mantle, and (2) the 7.4 km/s layer does not correspond to the “reflection” Moho. Gravity modeling and comparison to existing seismic models show a general consensus in many aspects with respect to crustal structure. This new model forms the basis for speculation on some of the processes associated with rifting of the Basin and Range Province. One such process, lithospheric magmatism, is inferred from the strong attenuation of transmitted seismic waves, which occurs at the same interface at which high‐amplitude, bright spot reflections originate. Unlike previous models, the overall structure and velocity of the crust and uppermost mantle of our new model are similar to other regions worldwide which have undergone high degrees of extension.

Nevada↗
USGS2016

Erosional and depositional history of the Atlantic passive margin as recorded in detrital zircon fission-track ages and lithic detritus in Atlantic Coastal plain sediments

Comparison of fission-track (FT) ages of detrital zircons recovered from Atlantic Coastal Plain sediments to FT ages of zircons from bedrock in source terranes in the Appalachians provides a key to understanding the provenance of the sediments and, in turn, the erosional and depositional history of the Atlantic passive margin. In Appalachian source terranes, the oldest zircon fission-track (ZFT) ages from bedrock in the western Appalachians (defined for this paper as the Appalachian Plateau, Valley and Ridge, and far western Blue Ridge) are notably older than the oldest ages from bedrock in the eastern Appalachians (Piedmont and main part of the Blue Ridge). The age difference is seen both in ZFT sample ages and in individual zircon grain ages and reflects differences in the thermotectonic history of the rocks. In the east, ZFT data indicate that the rocks cooled from temperatures high enough to partially or totally reset ZFT ages during the Paleozoic and (or) Mesozoic. The majority of the rocks are interpreted to have cooled through the ZFT closure temperature (∼235 °C) at various times during the late Paleozoic Alleghanian orogeny. In contrast, most of the rocks sampled in the western Appalachians have never been heated to temperatures high enough to totally reset their ZFT ages. Reflecting their contrasting thermotectonic histories, nearly 80 percent of the sampled western rocks yield one or more zircon grains with very old FT ages, in excess of 800 Ma; zircon grains yielding FT ages this old have not been found in rocks in the Piedmont and main part of the Blue Ridge. The ZFT data suggest that the asymmetry of zircon ages of exposed bedrock in the eastern and western Appalachians was in evidence by no later than the Early Cretaceous and probably by the Late Triassic. Detrital zircon suites from sands collected in the Atlantic Coastal Plain provide a record of detritus eroded from source terranes in the Appalachians during the Mesozoic and Cenozoic. In Virginia and Maryland, sands of Early Cretaceous through late early Oligocene age do not yield any old zircons comparable in age to the old zircons found in bedrock in the western Appalachians. Very old zircons yielding FT ages >800 Ma are only encountered in Coastal Plain sands of middle early Miocene and younger age. Miocene and younger fluvial-deltaic deposits associated with the major mid-Atlantic Coastal Plain rivers that now head in the western Appalachians (the Hudson, Delaware, Susquehanna, Potomac, James, and Roanoke) contain abundant clasts of fossiliferous chert and quartzite and other distinctive rock types derived from Paleozoic rocks of the western Appalachians. These distinctive clasts have not been reported in older Coastal Plain sediments. The ZFT and lithic detritus data indicate that the drainage divide for one or more east-flowing mid-Atlantic rivers migrated west into the western Appalachians, and the river(s) began transporting western Appalachian detritus to the Atlantic Coastal Plain, sometime between the late early Oligocene and middle early Miocene. By no later than late middle Miocene most if not all of the major rivers that now head west of the Blue Ridge were transporting western Appalachian detritus to the Coastal Plain. Prior to the drainage divide migrating into the western Appalachians, the ZFT data are consistent with the dominant source of Atlantic Coastal Plain sediments being detritus from the Piedmont and main part of the Blue Ridge, with possible input from distant volcanic sources. The ZFT data suggest that the rapid increase in the rate of siliciclastic sediment accumulation in middle Atlantic margin offshore basins that peaked in the middle Miocene and produced almost 30 percent of the total volume of post-rift siliciclastic sediments in the offshore basins began in the early Miocene when Atlantic river(s) gained access to the relatively easily eroded Paleozoic sedimentary rocks of the western Appalachians.

American Journal of Science↗
USGS1982

Application of automated image analysis to coal petrography

The coal petrologist seeks to determine the petrographic characteristics of organic and inorganic coal constituents and their lateral and vertical variations within a single coal bed or different coal beds of a particular coal field. Definitive descriptions of coal characteristics and coal facies provide the basis for interpretation of depositional environments, diagenetic changes, and burial history and determination of the degree of coalification or metamorphism. Numerous coal core or columnar samples must be studied in detail in order to adequately describe and define coal microlithotypes, lithotypes, and lithologic facies and their variations. The large amount of petrographic information required can be obtained rapidly and quantitatively by use of an automated image-analysis system (AIAS). An AIAS can be used to generate quantitative megascopic and microscopic modal analyses for the lithologic units of an entire columnar section of a coal bed. In our scheme for megascopic analysis, distinctive bands 2 mm or more thick are first demarcated by visual inspection. These bands consist of either nearly pure microlithotypes or lithotypes such as vitrite/vitrain or fusite/fusain, or assemblages of microlithotypes. Megascopic analysis with the aid of the AIAS is next performed to determine volume percentages of vitrite, inertite, minerals, and microlithotype mixtures in bands 0.5 to 2 mm thick. The microlithotype mixtures are analyzed microscopically by use of the AIAS to determine their modal composition in terms of maceral and optically observable mineral components. Megascopic and microscopic data are combined to describe the coal unit quantitatively in terms of (V) for vitrite, (E) for liptite, (I) for inertite or fusite, (M) for mineral components other than iron sulfide, (S) for iron sulfide, and (VEIM) for the composition of the mixed phases (X i ) i = 1,2, etc. in terms of the maceral groups vitrinite V, exinite E, inertinite I, and optically observable mineral content M. The volume percentage of each component present is indicated by a subscript. For example, a lithologic unit was determined megascopically to have the composition (V) 13 (I) 1 (S) 1 (X 1 ) 83 (X 2 ) 2 . After microscopic analysis of the mixed phases, this composition was expressed as (V) 13 (I) 1 (S) 1 (V 63 E 19 I 14 M 4 ) 83 (V 67 E 11 I 13 M 9 ) 2 . Finally, these data were combined in a description of the bulk composition as V 67 E 16 I 13 M 3 S 1 . An AIAS can also analyze textural characteristics and can be used for quick and reliable determination of rank (reflectance). Our AIAS is completely software based and incorporates a television (TV) camera that has optimum response characteristics in the range of reflectance less than 5%, making it particularly suitable for coal studies. Analysis of the digitized signal from the TV camera is controlled by a microprocessor having a resolution of 64 gray levels between full illumination and dark current. The processed image is reconverted for display on a TV monitor screen, on which selection of phases or features to be analyzed is readily controlled and edited by the operator through use of a lightpen. We expect that automated image analysis, because it can rapidly provide a large amount of pertinent information, will play a major role in the advancement of coal petrography.

International Journal of Coal Geology↗
USGS2008

The evolution of volcano-hosted geothermal systems based on deep wells from Karaha-Telaga Bodas, Indonesia

Temperature and pressure surveys, fluid samples, and petrologic analyses of rock samples from deep drill holes at the Karaha - Telaga Bodas geothermal field on the volcanic ridge extending northward from Galunggung Volcano, West Java, have provided a unique opportunity to characterize the evolution of an active volcano-hosted geothermal system. Wells up to 3 km in depth have encountered temperatures as high as 353??C and a weakly altered granodiorite that intruded to within 2 to 3 km of the surface. The intrusion is shallowest beneath the southern end of the field where an acid lake overlies a nearly vertical low resistivity structure (<10 ohm-m) defined by magnetotelluric measurements. This structure is interpreted to represent a vapor-dominated chimney that provides a pathway to the surface for magmatic gases. Four distinct hydrothermal mineral assemblages document the evolution of the geothermal system and the transition from liquid- to vapor-dominated conditions. The earliest assemblage represents the initial liquid-dominated system generated during emplacement of the granodiorite between 5910 ?? 76 and 4200 ?? 150 y BP. Tourmaline, biotite, actinolite, epidote and clay minerals were deposited contemporaneously at progressively greater distances from the intrusive contact (assemblage 1). At 4200 ?? 150 y BP, flank collapse and the formation of the volcano's crater, Kawah Galunggung, resulted in catastrophic decompression and boiling of the hydrothermal fluids. This event initiated development of the modern vapor-dominated regime. Chalcedony and then quartz were deposited as the early low salinity liquids boiled (assemblage 2). Both vapor- and liquid-rich fluid inclusions were trapped in the quartz crystals. Liquid-rich fluid inclusions from the southern part of the field record salinities ranging from 0 to 26 weight percent NaCl- CaCl2 equivalent and locally contain fluorite daughter crystals. We suggest, based on temperature-salinity relationships and evidence of boiling, that these fluids were progressively concentrated as steam was lost from the system. However, mixing with fluids derived from the underlying intrusion or generated during the formation of acid SO4 water on the vapor-dominated chimney margins could have contributed to the observed salinities. As pressures declined, CO2- and SO4-rich steam-heated water drained downward, depositing anhydrite and calcite (assemblage 3) in the fractures, limiting further recharge. Fluid inclusions with salinities up to 31 weight percent NaCl equivalent were trapped in these minerals as the descending water vaporized. The final assemblage is represented by precipitates of NaCl, KCl and FeClx deposited on rock surfaces in portions of the vapor-dominated zone that boiled dry. Vapor-dominated conditions extend over a distance of at least 10 km and to depths below sea level. Deep wells drilled into the underlying liquid-dominated reservoir in the northern and central part of the volcanic ridge produce low salinity fluids representing recent recharge of meteoric and steam-heated water. The evolution of volcanic-hosted vapor-dominated geothermal systems can be described by a five stage model. Stage 1 involves the formation of an over-pressured liquid-dominated geothermal system soon after magmatic intrusion. In Stages 2 and 3, pressures progressively decrease, and a curtain of steam-heated water surrounding a magmatic vapor-dominated chimney at 350??C and 14 ?? 2 MPa develops. The relatively low pressure near the base of the chimney causes liquid inflow adjacent to the intrusion and the development of a secondary marginal vapor-dominated zone. In Stage 4, the magmatic vapor discharge from the intrusion becomes small, vapor pressure declines, and the secondary vapor-dominated zone expands above the intrusion. In Stage 5, the vapor-dominated zone floods because heat from the intrusion is insufficient to boil all liquid inflow. A more common, liquid-dominated volcanic-hosted system the

American Journal of Science↗
USGS2012

The mechanics and three-dimensional internal structure of active magmatic systems: Kilauea volcano, Hawaii

Interpretation of abundant seismic data suggests that Kilauea's primary conduit within the upper mantle is concentrically zoned to about 34-km depth. This zoned structure is inferred to contain a central core region of relatively higher permeability, surrounded by numerous dikes that are in intermittent hydraulic communication with each other and with the central core. During periods of relatively high magma transport, the entire cross section of the conduit is utilized. During periods of relatively low to moderate transport, however, only the central core is active. As the conduit penetrates the oceanic crust and enters the volcanic shield, it simultaneously supplies the deeper sections of the rift zones (6-to 10-km depth) and the roots of the summit reservoir with picritic magma. The rift zones at depth are inferred to be almost wholly molten and to possess a high degree of fluid continuity from Heiheiahulu in the East Rift Zone, 45 km westward through the roots of the summit magma reservoir, and well into the Southwest Rift Zone. Higher in the shield, the subcaldera magma reservoir and the shallow rift zones occupy the 2-to 4-km depth interval. Summit-differentiated olivine tholeiite (ρ ≈ 2.62 g cm −3 ) is periodically injected laterally along a horizon of neutral buoyancy within the rift zones, where the density of the magma is just balanced by the in situ density of the shield (Ryan, 1987 a , b ). Deep rift zone intrusions push seaward the deep tectonic blocks of the volcano's south flank. Shallow rift intrusions build a sheeted dike complex, inferred to be in isostatic equilibrium with the higher-density deep rift cores below. General finite element analyses are presented for the deformation and stress fields surrounding such dikes in the horizontal and vertical planes. The dike tip in two and three dimensions is surrounded by a tubular core of tensile (σ 1 , σ 2 ) and shear stress (τ max ). The displacement field is characterized by counterrotating cells on either side of the dike tip which, in vertical orientation, produce the characteristic subsidence above the dike complex, with uplift on either side, forming a ridge-trough-ridge structure. A finite element model of Kilauea's shield computes the displacement fields and principal stress (σ 1 ) distributions resulting from intrusive activity on each or both of the rift zones. Within the summit region, tensile stress lobes produced by the three-dimensional upward extension of the intrusions superpose constructively to produce calderawide regimes of tensile stress, conducive to caldera development. Parametric studies of (1) intrusion in the East Rift Zone only, (2) intrusion in the Southwest Rift Zone only, and (3) intrusion in both rift zones demonstrate their unique kinematic contributions. For case 1, the caldera undergoes a counterclockwise rotation (torque up state) conducive to the development of rightstepping en echelon eruptive fissures, as exemplified by the August 14, 1971, eruption. For case 2, the caldera undergoes a clockwise rotation (torque down state) conducive to the development of left-stepping eruptive fissures, as occurred during the December 31, 1974, eruption. For case 3, the caldera substructure is driven due southward, producing the southward migration of the upper portions of the summit magma reservoir.

Journal of Geophysical Research Solid Earth↗