Geology ReportsSearch

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 523 records · Page 29Linked to original sources

Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor

Recent observations of seismic events at various volcanoes suggest that harmonic tremor results from the sustained occurrence of so-called long-period or low-frequency events. Accordingly, we can view the long-period volcanic event as the elementary process of tremor and interpret it as the impulse response of the tremor-generating system. We present a seismic model in which the source of tremor is the acoustic resonance of a fluid-filled volcanic pipe triggered by excess gas pressure. The model consists of three elements, namely, a triggering mechanism, a resonator, and a radiator. For simplicity, we assume a hemispherical trigger, cylindrial resonator, and circular radiator set in a vertical configuration with the trigger capping the top of the pipe and the disk-shaped radiator shutting off its bottom. Considering the simple case of a source buried in a homogeneous half space, we then apply the discrete wave number method to obtain a complete representation of the ground motion response at near and intermediate distances. The results demonstrate that the displacement attributed to the pipe dominates the near-field motion, while that due to the disk is representative of the intermediate and far fields. The trigger itself has a smaller contribution, mainly limited to the field in the proximity of the source. The characteristics displayed by the free surface response evolve from a strong impulsive signature in the immediate vicinity of the epicenter to a well-developed harmonic wave train dominated by Rayleigh waves at larger distances. No clear shear arrival can be detected in the synthetic seismograms. The displacement spectrum reflects the organ-pipe modes of the conduit, and the bandwidth associated with the dominant spectral peak of motion is controlled by the combined losses due to viscous attenuation in the fluid and elastic radiation into the solid. In the case of the cylindrical magma column considered, the radiation loss is proportional to the square of the pipe radius, while the loss related to viscous damping is inversely proportional to the same factor, indicating that the relative importance of the two loss mechanisms is critically dependent on the geometry of the magma reservoir. The relative importance of the pipe and disk elements, likewise, is a function of the conduit cross section. This suggests the possibility of determining the geometry of the source as well as the radiation loss and in situ magma viscosity from a comparison of near- and far-field observations.

Journal of Geophysical Research Solid Earth

Assessment of Mars Pathfinder landing site predictions

Remote sensing data at scales of kilometers and an Earth analog were used to accurately predict the characteristics of the Mars Pathfinder landing site at a scale of meters. The surface surrounding the Mars Pathfinder lander in Ares Vallis appears consistent with orbital interpretations, namely, that it would be a rocky plain composed of materials deposited by catastrophic floods. The surface and observed maximum clast size appears similar to predictions based on an analogous surface of the Ephrata Fan in the Channeled Scabland of Washington state. The elevation of the site measured by relatively small footprint delay-Doppler radar is within 100 m of that determined by two-way ranging and Doppler tracking of the spacecraft. The nearly equal elevations of the Mars Pathfinder and Viking Lander 1 sites allowed a prediction of the atmospheric conditions with altitude (pressure, temperature, and winds) that were well within the entry, descent, and landing design margins. High-resolution (~38 m/pixel) Viking Orbiter 1 images showed a sparsely cratered surface with small knobs with relatively low slopes, consistent with observations of these features from the lander. Measured rock abundance is within 10% of that expected from Viking orbiter thermal observations and models. The fractional area covered by large, potentially hazardous rocks observed is similar to that estimated from model rock distributions based on data from the Viking landing sites, Earth analog sites, and total rock abundance. The bulk and fine-component thermal inertias measured from orbit are similar to those calculated from the observed rock size-frequency distribution. A simple radar echo model based on the reflectivity of the soil (estimated from its bulk density), and the measured fraction of area covered by rocks was used to approximate the quasi-specular and diffuse components of the Earth-based radar echos. Color and albedo orbiter data were used to predict the relatively dust free or unweathered surface around the Pathfinder lander compared to the Viking landing sites. Comparisons with the experiences of selecting the Viking landing sites demonstrate the enormous benefit the Viking data and its analyses and models had on the successful predictions of the Pathfinder site. The Pathfinder experience demonstrates that, in certain locations, geologic processes observed in orbiter data can be used to infer surface characteristics where those processes dominate over other processes affecting the Martian surface layer. Copyright 1999 by the American Geophysical Union.

Journal of Geophysical Research E: Planets

Historic surface slip along the San Andreas Fault near Parkfield, California

The Parkfield Earthquake Prediction Experiment is focusing close attention on the 44-km-long section of the San Andreas fault that last ruptured seismically in 1966 ( M s 6.0). The 20-km-long central segment of the 1966 Parkfield rupture, extending from the mainshock epicenter at Middle Mountain southeastward to Gold Hill, forms a 1- to 2-km salient northeastward away from the dominant N40°W strike. Following the 1966 earthquake afterslip, aseismic slip has been nearly constant. Moderate Parkfield earthquakes have recurred on average every 21 years since 1857, when a great earthquake ( M ≈8) ruptured at least as far north as the southern Parkfield segment. Many measurements of slip have been made near Parkfield since 1966. Nevertheless, much of the history of surface slip remained uncertain, especially the total amount associated with the 1966 event. In 1985 we measured accumulated slip on the four oldest cultural features offset by the fault along the 1966 Parkfield rupture segment. We interpret net slip on each feature as a sum of event slip (sum of coseismic and rapid preseismic and postseismic slip) from Parkfield earthquakes and steady interseismic slip as measured over the last 20 years on nearby alinement arrays, creep meters, and trilateration lines. We assumed for each site that event slip was identical for the 1922, 1934, and 1966 Parkfield events and that long-term average rate of interseismic slip was constant between all events. Two fences on the southern segment, southeast of Gold Hill, indicate event slip of 13 and 15 cm and interseismic slip rate of 0.36 and 0.30 cm/yr since 1959 and 1908, respectively. At these sites, redundant independent data support our assumption that both event and interseismic slip occur uniformly. On the central segment, near Parkfield, both the 1934 and the 1966 ruptures offset a bridge built in 1932. Interseismic slip rate near the bridge has been about 1.1 cm/yr since 1966; thus we deduce an average event slip of 31 cm for the 1934 and 1966 earthquakes. On a parallel fault trace, 1 km to the southwest, slip was about 8 cm in 1966; thus total event slip summed across the entire fault zone near Parkfield was nearly 40 cm. On Middle Mountain, 4 km north of the 1966 mainshock epicenter, an offset fence indicates 17 cm of slip in 1966 and a 2.26-cm/yr interseismic slip rate since circa 1946. Thus the central segment of the 1966 rupture is characterized by much larger event slip (∼40 cm) than both distal segments (∼15 cm). This amount of surface slip per event is about twice what had been previously assumed. Larger 1966 surface slip in the central part of the rupture is geodetically compatible with a coseismic slip of 65±10 cm slip on a narrow, buried asperity between Middle Mountain and Gold Hill that has been inferred from the depth distribution of early aftershocks. Assuming our characteristic surface slip model, one can further deduce a deficit in slip since the great 1857 earthquake. Taking the long-term slip rate as 3.3 cm/yr, the surface slip deficit is 3±0.2 m south of Gold Hill but only 0.3±0.3 m northward from Parkfield.

Journal of Geophysical Research Solid Earth

The Lepanto Cu–Au deposit, Philippines: A fossil hyperacidic volcanic lake complex

Hyperacidic lakes and associated solfatara in active volcanoes are the expression of magmatic gas expansion from source to surface. Here we show for the first time, that the vein system that comprises the ~ 2 Ma high-sulfidation, Lepanto copper–gold deposit in the Mankayan district (Philippines) was associated with a contemporary hyperacidic volcanic lake complex—possibly the first such lake recognized in the geological record. A 15–20‰ difference in sulfur isotopic composition between barite and sulfides and sulfosalts in the vent fumarole encrustations supports the interpretation that SO 2 -rich volcanic gas vented into the base of the lake and marginal to it and ties the mineralization directly to magmatic gas expansion, fracture propagation, and mineralization that occurred through a series of decompression steps within the feeder fracture network. These data confirm that crater lake environments such as Kawah Ijen (Java, Indonesia) provide modern day analogs of the Lepanto and other high sulfidation Cu–Au depositing environments. We also provide extensive analysis of sulfosalt–sulfide reactions during vein formation within the hyperacidic lake complex. Pyrite ± silica deposited first at high temperature followed by enargite that preserves the vapor–solid diffusion of, for example, antimony, tin, and tellurium into the vapor from the crystallizing solid. Subsolidus, intra-crystalline diffusion continued as temperature declined. Pyrite and enargite are replaced by Fe-tennantite in the lodes which initially has low Sb/(Sb + As) atomic ratios around 13.5% close to the ideal tennantite formula, but evolves to higher ratios as crystallization proceeds. Fumarole encrustation clasts and sulfosalts in the lake sediment are more highly evolved with a larger range of trace element substitutions, including antimony. Substitution of especially Zn, Te, Ag, and Sn into tennantite records metal and semi-metal fractionation between the expanding magmatic gas and deposited sulfide sublimates provides a rare insight into the fate of metals and semi-metals in the shallower parts of fracture arrays that feed modern hyperacidic lakes. These data support a growing understanding of the formation of high-sulfidation gold deposits as the consequence of single-phase expansion of gas from magmatic-gas reservoirs beneath the surface of active volcanoes without the intervention of a later aqueous fluid including groundwater. Aggressive sulfide–sulfosalt reactions, including pitting and the almost complete dissolution of earlier minerals, are persistent characteristics of the vein assemblages and precious metals typically occur late in pits or along brittle fractures. These characteristics support a hypothesis of mineral deposition at temperatures of the order of 600 °C in contrast to available fluid inclusion data from enargite that record temperatures following phase transitions in the sulfosalt during the retrograde devolution of the deposit in the presence of groundwater.

Benguet

Debris-flow origin for the Simud/Tiu deposit on Mars

A late Hesperian smooth plains deposit on Mars interpreted as a debris flow extends more than 2000 km from Hydraotes Chaos, through Simud and Tiu Valles, and into Chryse Planitia. The Simud/Tiu deposit widens out to >1000 km and embays streamlined landforms and knobs made up of sedimentary and perhaps volcanic deposits that were carved by earlier channeling activity. Morphologic features of the Simud/Tiu deposit observed in Viking and Pathfinder images are generally consistent with a debris-flow origin, but some of the deposit's salient features are not readily explained by catastrophic flooding or ice flow. Internal depressions appear to be bounded by linear scarps along flow margins where differential shearing may have occurred and in areas where flow spreading may have produced zones of extensional breakup and thinning within the flow. Possible flow lobes within the deposit may have formed by successive flow surges within the flow unit. The Pathfinder landing site is on the Simud/Tiu deposit, and the observations there are consistent with debris flow. The low, longitudinal ridges at the site may have formed by clast interactions as the flow ground to a halt. Imbricated, planar rocks on the ridges, such as in the Rock Garden, also may have been emplaced by debris or ice flow. However, stream energy calculations at Ares Vallis and channel geology indicate that flooding probably was incapable of emplacing the meter-size boulders observed at the Pathfinder site. Dewatering of pressurized zones in the debris flow or underlying material may be responsible for mud eruptions that formed a couple of patches of low pancakelike shields up to 5 km in diameter and for probable water flows that formed two small rille channels a few kilometers long. Local irregular grooves may be cracks that resulted from later desiccation and contraction of the flow material. The debris-flow unit apparently coalesced from outflows of water-fluidized debris originating from beneath chaotic and hummocky terrains within and along the margins of Simud and Tiu Valles. The deposit is onlapped from the north by another flow deposit originating from Acidalia Planitia. If the Simud/Tiu debris flow had entered a standing body of water, a turbidity current may have arisen from the debris flow and then backflowed over the debris flow to account for the Acidalia deposit.

Journal of Geophysical Research E: Planets

Seafloor environments within the Boston Harbor- Massachusetts Bay sedimentary system: A regional synthesis

Modern seafloor sedimentary environments within the glaciated, topographically complex Boston Harbor and Massachusetts Bay area have been interpreted and mapped from an extensive collection of sidescan sonar records and supplemental marine geologic data. Three categories of environments are present that reflect the dominant long-term processes of erosion or nondeposition, deposition, and sediment reworking. (1) Environments of erosion or nondeposition comprise exposures of bedrock, glacial drift, coarse lag deposits, and possibly coastal plain rocks that contain sediments (where present) ranging from boulder fields to gravelly sands and occur in areas of relatively strong currents. (2) Environments of deposition contain fine-grained sediments ranging from muddy sands to muds that have accumulated in areas of predominantly weak bottom currents. (3) Environments of sediment reworking contain patches with textures ranging from sandy gravels to muds that have been produced by a combination of erosion and deposition in areas with variable bottom currents. The distribution of sedimentary environments across the Boston Harbor-Massachusetts Bay area is extremely patchy. Locally, this patchiness is due either to modifications of bottom-current strength (caused by the irregular topography and differences in water depth) or to small-scale changes in the supply of fine-grained sediments. Regional patchiness, however, reflects differences in geologic and oceanographic conditions among the estuarine, inner shelf, and basinal parts of the sedimentary system. The estuarine part of the system (Boston Harbor) is a depositional trap for fine-grained sediments because it is protected from large waves, has generally weak and variable tidal currents, and receives a large supply of fine grained detritus from natural and anthropogenic sources. The inner shelf, on the other hand, is largely an area of erosion or nondeposition due to sediment removal and redistribution during past sea-level changes, to sediment resuspension and winnowing by modern waves and currents, and to an inadequate supply of fine-grained sediments. The basinal part of the system (Stellwagen Basin) is mainly a tranquil depositional environment in which fine-grained sediments from several potential sources settle through the water column and accumulate under weak bottom currents. This study indicates areas within the Boston Harbor-Massachusetts Bay sedimentary system where fine-grained sediments and associated contaminants are likely to be either moved or deposited. It also provides a guide to the locations and variability of benthic habitats.

Massachusetts Bay

Mars Global Surveyor Thermal Emission Spectrometer experiment: Investigation description and surface science results

The Thermal Emission Spectrometer (TES) investigation on Mars Global Surveyor (MGS) is aimed at determining (1) the composition of surface minerals, rocks, and ices; (2) the temperature and dynamics of the atmosphere; (3) the properties of the atmospheric aerosols and clouds; (4) the nature of the polar regions; and (5) the thermophysical properties of the surface materials. These objectives are met using an infrared (5.8- to 50-μm) interferometric spectrometer, along with broadband thermal (5.1- to 150-μm) and visible/near-IR (0.3- to 2.9-μm) radiometers. The MGS TES instrument weighs 14.47 kg, consumes 10.6 W when operating, and is 23.6×35.5×40.0 cm in size. The TES data are calibrated to a 1-σ precision of 2.5 −6 ×10 −8 W cm −2 sr −1 /cm −1 , 1.6×10 −6 W cm −2 sr −1 , and ∼0.5 K in the spectrometer, visible/near-IR bolometer, and IR bolometer, respectively. These instrument subsections are calibrated to an absolute accuracy of ∼4×10 −8 W cm −2 sr −1 /cm −1 (0.5 K at 280 K), 1–2%, and ∼1–2 K, respectively. Global mapping of surface mineralogy at a spatial resolution of 3 km has shown the following: (1) The mineralogic composition of dark regions varies from basaltic, primarily plagioclase feldspar and clinopyroxene, in the ancient, southern highlands to andesitic, dominated by plagioclase feldspar and volcanic glass, in the younger northern plains. (2) Aqueous mineralization has produced gray, crystalline hematite in limited regions under ambient or hydrothermal conditions; these deposits are interpreted to be in-place sedimentary rock formations and indicate that liquid water was stable near the surface for a long period of time. (3) There is no evidence for large-scale (tens of kilometers) occurrences of moderate-grained (>50-μm) carbonates exposed at the surface at a detection limit of ∼10%. (4) Unweathered volcanic minerals dominate the spectral properties of dark regions, and weathering products, such as clays, have not been observed anywhere above a detection limit of ∼10%; this lack of evidence for chemical weathering indicates a geologic history dominated by a cold, dry climate in which mechanical, rather than chemical, weathering was the significant form of erosion and sediment production. (5) There is no conclusive evidence for sulfate minerals at a detection limit of ∼15%. The polar region has been studied with the following major conclusions: (1) Condensed CO 2 has three distinct end-members, from fine-grained crystals to slab ice. (2) The growth and retreat of the polar caps observed by MGS is virtually the same as observed by Viking 12 Martian years ago. (3) Unique regions have been identified that appear to differ primarily in the grain size of CO 2 ; one south polar region appears to remain as black slab CO 2 ice throughout its sublimation. (4) Regional atmospheric dust is common in localized and regional dust storms around the margin and interior of the southern cap. Analysis of the thermophysical properties of the surface shows that (1) the spatial pattern of albedo has changed since Viking observations, (2) a unique cluster of surface materials with intermediate inertia and albedo occurs that is distinct from the previously identified low-inertia/bright and high-inertia/dark surfaces, and (3) localized patches of high-inertia material have been found in topographic lows and may have been formed by a unique set of aeolian, fluvial, or erosional processes or may be exposed bedrock.

Journal of Geophysical Research E: Planets

Three-dimensional models of deformation near strike-slip faults

We use three-dimensional elastic models to help guide the kinematic interpretation of crustal deformation associated with strike-slip faults. Deformation of the brittle upper crust in the vicinity of strike-slip fault systems is modeled with the assumption that upper crustal deformation is driven by the relative plate motion in the upper mantle. The driving motion is represented by displacement that is specified on the bottom of a 15-km-thick elastic upper crust everywhere except in a zone of finite width in the vicinity of the faults, which we term the “shear zone.” Stress-free basal boundary conditions are specified within the shear zone. The basal driving displacement is either pure strike slip or strike slip with a small oblique component, and the geometry of the fault system includes a single fault, several parallel faults, and overlapping en echelon faults. We examine the variations in deformation due to changes in the width of the shear zone and due to changes in the shear strength of the faults. In models with weak faults the width of the shear zone has a considerable effect on the surficial extent and amplitude of the vertical and horizontal deformation and on the amount of rotation around horizontal and vertical axes. Strong fault models have more localized deformation at the tip of the faults, and the deformation is partly distributed outside the fault zone. The dimensions of large basins along strike-slip faults, such as the Rukwa and Dead Sea basins, and the absence of uplift around pull-apart basins fit models with weak faults better than models with strong faults. Our models also suggest that the length-to-width ratio of pull-apart basins depends on the width of the shear zone and the shear strength of the faults and is not constant as previously suggested. We show that pure strike-slip motion can produce tectonic features, such as elongate half grabens along a single fault, rotated blocks at the ends of parallel faults, or extension perpendicular to overlapping en echelon faults, which can be misinterpreted to indicate a regional component of extension. Zones of subsidence or uplift can become wider than expected for transform plate boundaries when a minor component of oblique motion is added to a system of parallel strike-slip faults.

Journal of Geophysical Research B: Solid Earth

The competition between thermal contraction and differentiation in the stress history of the Moon

The scarcity of both extension and compression features on the Moon strongly constrains the history of the lunar radius—to variations of less than ±1 km over the past 3.8 Gyr. This limit has traditionally been interpreted as requiring a delicate balance between thermal contraction of the near‐surface and expansion of a substantial cold interior region. Recent theories of lunar origin (e.g., giant impact), in contrast, favor a “hot” initial state. We propose that a reconciliation may be possible by taking account of the volume change Δ V / V | d due to differentiation. We calculate STP densities based on simplified normative mineralogies for a suite of estimates of the bulk lunar composition, of primary lunar basalt, and of the residuum left when the maximum amount of the latter is extracted from the former. Typically Δ V / V | d ≃ 2 to 5%—an expansion equivalent to heating by ∼10 3 K. Provided the timing of differentiation is correct, one might offset the cooling of a magma ocean as much as 630 km deep by differentiation of the remainder of the Moon (which need not start much below the solidus temperature). A large but not impossible amount of gabbroic melt production is implied: ∼100 times the volume of mare basalts known to have been extruded. We do not address the detailed genetic relationship of this melt to the basalts observed on the lunar surface but point out that it need not have reached the surface directly or even have entered the crust in order for the expansion to have occurred. To assess the timing of melt formation, we investigate a simple conductive lunar thermal model which takes account of both Δ V / V | d and thermal contraction. Our initial state is characterized by a central temperature T c and a depth Z 0 above which the material (derived from the magma ocean) is already at the solidus and is not suceptible to volume changes upon further differentiation. We find a range of models satisfying the limits on radius increase and decrease. The hottest has T c = 1210 K, Z 0 = 400 km; without Δ V / V | d , we would need a larger or colder (or both) core, e.g., T c ≲ 700 K for Z 0 = 200–400 km, in agreement with previous investigators. Our modeling thus lends credence to the idea that the Moon could have been initially ≳50% molten (with the remainder relatively close to the solidus) and yet experienced little volume change over the last 3.8 Gyr.

Journal of Geophysical Research B: Solid Earth

Late Quaternary stratigraphy and depositional history of the Long Island Sound basin

The stratigraphy of Late Quaternary geologic units beneath Long Island Sound (LIS) is interpreted from 3,500 km of high-resolution, seismic-reflection profiles supplemented by vibracore data. Knowledge gained from onshore regional geologic studies and previous offshore investigations is also incorporated in these interpretations. Glacial deposits overlie and nearly fill an Inner Lowland which is floored by crystalline rocks in the north and bounded to the south by an irregular, north-facing cuesta scarp in Cretaceous Coastal-Plain strata. The thick glacial section consists largely of sediment that was deposited in glacial Lake Connecticut during the retreat of the last (late Wisconsin) ice sheet. Ice-marginal lacustrine fan deposits and submerged extensions of the southeastern Connecticut recessional moraines mark positions of the ice margin during systematic northward retreat. When the ice margin stood at positions just north of the present shoreline of Connecticut (17.6-16.5 ka), glacial Lake Connecticut was slightly larger than Long Island Sound is today. Ice-marginal deltas and near-ice marginal fluviodeltaic deposits occur along coastal Connecticut. Delta-building was diachronous; progressively younger deltas record northwesterly retreat of the ice margin and lowering lake levels. Distally from the deltas, varied lake-clays (up to 150 m thick) drape the underlying topography. The lake spillway cut across the Harbor Hill-Roanoke Point-Fishers Island-Charlestown end moraine at its lowest point just west of Fishers Island (The Race). Lake levels fell as the spillway notch was erosionally deepened, and gradually the lakebed became subaerially exposed. A largely preserved channel system cuts the lake-deposit surface; fluvial channel-fill sediments are overlain by estuarine sediments in the channel system. Perhaps as early as 15 ka, but no later than 13 ka, the sea entered the Basin through the channel at The Race and quickly transgressed westward via the lakebed channel system. As the sea rose further and spread throughout the Basin, a wave-cut surface was formed. Extensive marine deltaic deposits, up to 40 m thick, overlie the wave-cut marine unconformity in the east-central part of the Basin. These sediments were derived from erosion of the drained bed of Lake Hitchcock, in the Connecticut valley to the north, and were deposited in a -40-m relative sea at about 12.4 ka. Where quiet waters prevail, marine mud generally less than 15 m thick blankets the older deposits of the Basin. Elsewhere, especially in eastern LIS, tidal currents are actively reworking and transporting glacial and postglacial deposits.

Connecticut, New York

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

Scientific results of the Hydrate-01 Stratigraphic Test Well Program, Western Prudhoe Bay Unit, Alaska North Slope

The United States Department of Energy, the MH21-S Research Consortium of Japan, and the United States Geological Survey are collaborating to enable gas hydrate scientific drilling and extended-duration reservoir response testing on the Alaska North Slope. To feasibly execute such a test, a location is required that is accessible from existing roads and gravel pads and that can be occupied without disrupting ongoing industry operations. A review of potential locations meeting these criteria determined the likely occurrence of gas hydrate in two fine-grained marginal-marine sands of Tertiary age in the vicinity of the inactive “Kuparuk State 7-11-12” exploration pad in the western Prudhoe Bay Unit (PBU). Existing well and seismic data for that site were insufficient to preclude the potential for free gas occurrence within the deeper (and most prospective) target sand. Therefore, with support from the PBU Working Interest Owners, Alaska Department of Natural Resources, and Petrotechnical Resources Alaska, the Hydrate-01 Stratigraphic Test Well (STW) was drilled in December 2018 to confirm the suitability of the site for future gas hydrate scientific testing. The Hydrate-01 well was successfully drilled to −3290 ft (1003 m) subsea vertical depth at a bottom hole location of approximately 900 ft (∼275 m) east of the surface location. The drilling program featured acquisition of a full suite of logging while drilling data, the collection of side-wall pressure cores, and the installation of distributed temperature and distributed acoustic sensor fiber-optic cables. The log data acquired confirmed the occurrence of gas hydrate at high saturation in two target sands. Integrated evaluation of log and sidewall core data provide petrophysical and geomechanical property information that allow for potential reservoir response to depressurization to be simulated. The deeper “B1 sand” is deemed to be most favorable for reservoir response testing as a result of confirmed gas hydrate occurrence in sediments of high intrinsic permeability, location within 100 ft (30 m) of the base of gas hydrate stability, and minimal risk for direct communication with permeable water-bearing (hydrate-free) zones. The shallower “D1 sand” provides a secondary target that is differentiated by colder in situ temperatures and the interpreted direct hydraulic communication to a lower section of non-hydrate-bearing, water-saturated sand. The Hydrate-01 log data also confirm the occurrence of at least one sub-seismic fault in close proximity to the B1 sand reservoir. To better image the distribution of the gas-hydrate-bearing reservoir sections and associated faults, a three-dimensional (3D) vertical seismic profile was conducted in early 2019 using the distributed acoustic sensors installed as part of the Hydrate-01 STW completion. Detailed two-dimensional (2D) and 3D geologic models have been constructed to enable numerical simulations to inform the planning for potential future scientific tests of reservoir response to depressurization at the site.

Alaska

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

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

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

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

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

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