SEARCH · Geology Reports

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,261 records · Page 70Linked to original sources

USGS2012

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

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

Journal of Geophysical Research Solid Earth↗

USGS2012

Seismic refraction measurements within the Peninsular terrane, south central Alaska

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

Alaska↗

USGS2005

Temporal patterns of Northern Goshawk nest area occupancy and habitat: A retrospective analysis

We studied occupancy and habitat associations of Northern Goshawks (Accipiter gentilis) at nest areas in south-central Oregon in 1992-94. We surveyed 51 pre-1992 nest areas (i.e., historical breeding areas first discovered during 1973-91) for goshawks and used aerial-photograph interpretation to document forest cover conditions and changes over time between areas that were occupied by goshawks and those where we did not detect goshawks (no-response sites). We also surveyed for new nests during 1992-94. Of 38 occupied nests first found in 1992-94 (i.e., post-1992 nest areas), 86% (33/38) were in mid-aged (mean stand DBH 23-53 cm, <15 trees/ha >53 cm DBH) or late (???15 trees/ha >53 cm DBH; mean stand DBH >53 cm) closed (>50% canopy closure) structural-stage forest. Occupancy of historical (pre-1992) nest areas by goshawks was 29% (15/51). Of 46 pre-1992 nest areas that we examined for habitat change, 15 were occupied by goshawks in 1994 and had more mid-aged closed and late closed forest in 12-, 24-, 52-, 120-, and 170-ha circular areas centered on nest locations than did 31 no-response areas. There was no difference in the amount of late closed and mid-aged closed forest in pre-1992 nest areas compared with occupied pre-1992 nest areas. A logistic regression model for all occupied nest areas confirmed that late closed and mid-aged closed forest variables were important indicators of forest conditions that supported breeding pairs. Goshawks were more likely to persist in the historical nest areas that had about 50% of mature and older closed-canopy forest within the 52ha scale. We recommend retaining existing late closed, late open, and mid closed structure within 52ha scale of the nest site. Moreover, late closed and mid closed structure combined should not fall below 50% within the 52-ha scale and should exceed 40% within the 170-ha scale surrounding the nest site. To optimize conditions for breeding goshawks, we recommend retaining large trees (>53 cm DBH) to help preserve stand integrity, maintain closed canopies, and provide connectivity to alternative nest sites within nest areas. ?? 2005 The Raptor Research Foundation, Inc.

Journal of Raptor Research↗

USGS2006

Isotope geochemistry of early Kilauea magmas from the submarine Hilina bench: The nature of the Hilina mantle component

Submarine lavas recovered from the Hilina bench region, offshore Kilauea, Hawaii Island provide information on ancient Kilauea volcano and the geochemical components of the Hawaiian hotspot. Alkalic lavas, including nephelinite, basanite, hawaiite, and alkali basalt, dominate the earliest stage of Kilauea magmatism. Transitional basalt pillow lavas are an intermediate phase, preceding development of the voluminous tholeiitic subaerial shield and submarine Puna Ridge. Most alkalic through transitional lavas are quite uniform in Sr–Nd–Pb isotopes, supporting the interpretation that variable extent partial melting of a relatively homogeneous source was responsible for much of the geochemical diversity of early Kilauea magmas ( Sisson et al., 2002 ). These samples are among the highest 206 Pb/ 204 Pb known from Hawaii and may represent melts from a distinct geochemical and isotopic end-member involved in the generation of most Hawaiian tholeiites. This end-member is similar to the postulated literature Kea component, but we propose that it should be renamed Hilina, to avoid confusion with the geographically defined Kea-trend volcanoes. Isotopic compositions of some shield-stage Kilauea tholeiites overlap the Hilina end-member but most deviate far into the interior of the isotopic field defined by magmas from other Hawaiian volcanoes, reflecting the introduction of melt contributions from both “Koolau” (high 87 Sr/ 86 Sr, low 206 Pb/ 204 Pb) and depleted (low 87 Sr/ 86 Sr, intermediate 206 Pb/ 204 Pb) source materials. This shift in isotopic character from nearly uniform, end-member, and alkalic, to diverse and tholeiitic corresponds with the major increase in Kilauea's magmatic productivity. Two popular geodynamic models can account for these relations: (1) The upwelling mantle source could be concentrically zoned in both chemical/isotopic composition, and in speed/extent of upwelling, with Hilina (and Loihi) components situated in the weakly ascending margins and the Koolau component in the interior. The depleted component could be refractory and spread throughout or scavenged from the overlying lithosphere. (2) The Hilina (and Loihi) components could be a more fertile material (lower melting temperature) spread irregularly throughout the Hawaiian source in a matrix of more refractory depleted and Koolau compositions. Modest upwelling along the leading hotspot margin melts the fertile domains predominantly, while the refractory matrix also partially melts in the more vigorously upwelling hotspot interior, diluting the Hilina and Loihi components and yielding voluminous isotopically diverse tholeiitic magmas.

Hawai'i↗

USGS2005

The Circumpolar Arctic vegetation map

Question: What are the major vegetation units in the Arctic, what is their composition, and how are they distributed among major bioclimate subzones and countries? Location: The Arctic tundra region, north of the tree line. Methods: A photo-interpretive approach was used to delineate the vegetation onto an Advanced Very High Resolution Radiometer (AVHRR) base image. Mapping experts within nine Arctic regions prepared draft maps using geographic information technology (ArcInfo) of their portion of the Arctic, and these were later synthesized to make the final map. Area analysis of the map was done according to bioclimate subzones, and country. The integrated mapping procedures resulted in other maps of vegetation, topography, soils, landscapes, lake cover, substrate pH, and above-ground biomass. Results: The final map was published at 1:7 500 000 scale map. Within the Arctic (total area = 7.11 x 106 km 2), about 5.05 ?? 106 km2 is vegetated. The remainder is ice covered. The map legend generally portrays the zonal vegetation within each map polygon. About 26% of the vegetated area is erect shrublands, 18% peaty graminoid tundras, 13% mountain complexes, 12% barrens, 11% mineral graminoid tundras, 11% prostrate-shrub tundras, and 7% wetlands. Canada has by far the most terrain in the High Arctic mostly associated with abundant barren types and prostrate dwarf-shrub tundra, whereas Russia has the largest area in the Low Arctic, predominantly low-shrub tundra. Conclusions: The CAVM is the first vegetation map of an entire global biome at a comparable resolution. The consistent treatment of the vegetation across the circumpolar Arctic, abundant ancillary material, and digital database should promote the application to numerous land-use, and climate-change applications and will make updating the map relatively easy. ?? IAVS; Opulus Press.

Journal of Vegetation Science↗

USGS2011

Lower crustal relaxation beneath the Tibetan Plateau and Qaidam Basin following the 2001 Kokoxili earthquake

In 2001 November a magnitude 7.8 earthquake ruptured a 400 km long portion of the Kunlun fault, northeastern Tibet. In this study, we analyse over five years of post-seismic geodetic data and interpret the observed surface deformation in terms of stress relaxation in the thick Tibetan lower crust. We model GPS time-series (first year) and InSAR line of sight measurements (years two to five) and infer that the most likely mechanism of post-seismic stress relaxation is time-dependent distributed creep of viscoelastic material in the lower crust. Since a single relaxation time is not sufficient to model the observed deformation, viscous flow is modelled by a lower crustal Burgers rheology, which has two material relaxation times. The optimum model has a transient viscosity 9 × 10 17 Pa s, steady-state viscosity 1 × 10 19 Pa s and a ratio of long term to Maxwell shear modulus of 2:3. This model gives a good fit to GPS stations south of the Kunlun Fault, while displacements at stations north of the fault are over-predicted. We attribute this asymmetry in the GPS residual to lateral heterogeneity in rheological structure across the southern margin of the Qaidam Basin, with thinner crust/higher viscosities beneath the basin than beneath the Tibetan Plateau. Deep afterslip localized in a shear zone beneath the fault rupture gives a reasonable match to the observed InSAR data, but the slip model does not fit the earlier GPS data well. We conclude that while some localized afterslip likely occurred during the early post-seismic phase, the bulk of the observed deformation signal is due to viscous flow in the lower crust. To investigate regional variability in rheological structure, we also analyse post-seismic displacements following the 1997 Manyi earthquake that occurred 250 km west of the Kokoxili rupture. We find that viscoelastic properties are the same as for the Kokoxili area except for the transient viscosity, which is 5 × 10 17 Pa s. The viscosities estimated for the Manyi and Kokoxili areas are consistent with constraints obtained from other earthquakes in the northwest and south central parts of the Tibetan Plateau.

Tibetan Plateau and Qaidam Basin↗

USGS2007

Cleats and their relation to geologic lineaments and coalbed methane potential in Pennsylvanian coals in Indiana

Cleats and fractures in Pennsylvanian coals in southwestern Indiana were described, statistically analyzed, and subsequently interpreted in terms of their origin, relation to geologic lineaments, and significance for coal permeability and coalbed gas generation and storage. These cleats can be interpreted as the result of superimposed endogenic and exogenic processes. Endogenic processes are associated with coalification (i.e., matrix dehydration and shrinkage), while exogenic processes are mainly associated with larger-scale phenomena, such as tectonic stress. At least two distinct generations of cleats were identified on the basis of field reconnaissance and microscopic study: a first generation of cleats that developed early on during coalification and a second generation that cuts through the previous one at an angle that mimics the orientation of the present-day stress field. The observed parallelism between early-formed cleats and mapped lineaments suggests a well-established tectonic control during early cleat formation. Authigenic minerals filling early cleats represent the vestiges of once open hydrologic regimes. The second generation of cleats is characterized by less prominent features (i.e., smaller apertures) with a much less pronounced occurrence of authigenic mineralization. Our findings suggest a multistage development of cleats that resulted from tectonic stress regimes that changed orientation during coalification and basin evolution. The coals studied are characterized by a macrocleat distribution similar to that of well-developed coalbed methane basins (e.g., Black Warrior Basin, Alabama). Scatter plots and regression analyses of meso- and microcleats reveal a power-law distribution between spacing and cleat aperture. The same distribution was observed for fractures at microscopic scale. Our observations suggest that microcleats enhance permeability by providing additional paths for migration of gas out of the coal matrix, in addition to providing access for methanogenic bacteria. The abundance, distribution, and orientation of cleats control coal fabric and are crucial features in all stages of coalbed gas operations (i.e., exploration and production). Understanding coal fabric is important for coal gas exploration as it may be related to groundwater migration and the occurrence of methanogenic bacteria, prerequisite to biogenic gas accumulations. Likewise, the distribution of cleats in coal also determines pathways for migration and accumulation of thermogenic gas generated during coalification. ?? 2007 Elsevier B.V. All rights reserved.

International Journal of Coal Geology↗

USGS2013

Reliability of fish size estimates obtained from multibeam imaging sonar

Multibeam imaging sonars have considerable potential for use in fisheries surveys because the video-like images are easy to interpret, and they contain information about fish size, shape, and swimming behavior, as well as characteristics of occupied habitats. We examined images obtained using a dual-frequency identification sonar (DIDSON) multibeam sonar for Atlantic sturgeon Acipenser oxyrinchus oxyrinchus , striped bass Morone saxatilis , white perch M. americana , and channel catfish Ictalurus punctatus of known size (20–141 cm) to determine the reliability of length estimates. For ranges up to 11 m, percent measurement error (sonar estimate – total length)/total length × 100 varied by species but was not related to the fish's range or aspect angle (orientation relative to the sonar beam). Least-square mean percent error was significantly different from 0.0 for Atlantic sturgeon ( x̄ = −8.34, SE = 2.39) and white perch ( x̄ = 14.48, SE = 3.99) but not striped bass ( x̄ = 3.71, SE = 2.58) or channel catfish ( x̄ = 3.97, SE = 5.16). Underestimating lengths of Atlantic sturgeon may be due to difficulty in detecting the snout or the longer dorsal lobe of the heterocercal tail. White perch was the smallest species tested, and it had the largest percent measurement errors (both positive and negative) and the lowest percentage of images classified as good or acceptable. Automated length estimates for the four species using Echoview software varied with position in the view-field. Estimates tended to be low at more extreme azimuthal angles (fish's angle off-axis within the view-field), but mean and maximum estimates were highly correlated with total length. Software estimates also were biased by fish images partially outside the view-field and when acoustic crosstalk occurred (when a fish perpendicular to the sonar and at relatively close range is detected in the side lobes of adjacent beams). These sources of bias are apparent when files are processed manually and can be filtered out when producing automated software estimates. Multibeam sonar estimates of fish size should be useful for research and management if these potential sources of bias and imprecision are addressed.

Journal of Fish and Wildlife Management↗

USGS2025

Greater sage-grouse seasonal habitat associations: A review and considerations for interpretation and management applications

Habitat features needed by wildlife can change in composition throughout the year, particularly in temperate ecosystems, leading to distinct seasonal spatial-use patterns. Studies of species-habitat associations therefore often focus on understanding relationships within discrete seasonal periods with common goals of prediction (e.g., habitat mapping) and inference (e.g., interpreting model coefficients). Across the range of the greater sage-grouse ( Centrocercus urophasianus ) of western North America, the increasing use of high-frequency tracking devices has led to a surge in habitat association studies covering multiple temporal periods and spatial extents. We reviewed the literature for seasonal habitat association studies corresponding to the second and third orders of selection (Johnson 1980). Our objectives were to summarize the methodological approaches used to estimate habitat associations to aid understanding in cross-study comparisons and identify common habitat features reported as selected or avoided within different seasonal periods. We reviewed 71 second- and third-order studies published from 2007–2023 that evaluated covariates collected in a geographic information system (GIS) and modeled probability of selection or intensity of use. The most common study design evaluated a single level of selection (third-order) and was multi-scale (i.e., covariates were measured at varying spatial grains). The most common model approach estimated habitat associations using resource selection functions (RSFs) fit with logistic regression. Studies mostly focused on the breeding periods and winter, but all seasons throughout the annual cycle were covered. There was clear support for selection of sagebrush and avoidance of trees and rugged terrain across seasons, and strong selection of mesic conditions in summer. However, habitat associations for most covariates were mixed, with proportionally equivalent selection and avoidance reported, even within the same seasons. Different factors hampered cross-study comparisons, including variation in study design, but additional contributors likely included important context-dependent habitat associations, such as functional responses to changing habitat availability. We suggest collaborative studies leveraging multiple datasets can help improve seasonal habitat inference by removing the effects of variable study designs.

Journal of Wildlife Management↗

USGS2005

Thermal observations of gas pistoning at Kilauea Volcano

Data acquired by three continuously recording thermal infrared thermometers situated on the north rim of Pu'u'O' o Crater at Kilauea Volcano during 2002 revealed episodes of periodic thermal pulses originating from a degassing vent on the crater floor. These thermal pulses are interpreted as gas release (jetting events) associated with gas pistoning, a mechanism observed previously at both Mauna Ulu and Pu'u'O' o. During a 35-day-long period spanning June and July 2002, gas pistoning was frequently the dominant mode of gas release, with as many as several hundred pulses occurring in uninterrupted series. On other days, degassing alternated between periods of quasi-continuous gas jetting and intervals of gas pistoning that contained a few to a few dozen pulses. Characteristic time intervals between pistoning events ranged from 2 up to 7 min. We identify three types of pistoning. Type 1 involves emission of lava, followed by gas jetting and drain back; type 2 is the same but the elevated position of the vent does not allow postjet drain back; and type 3 involves gas jetting only with no precursory lava flow. To explain gas pistoning, we apply a model whereby a stagnant cap of degassed magma develops in the conduit below the vent. Gas bubbles rise through the magma column and collect under the cap. The collective buoyancy of these bubbles pushes the cap upward. When the cap reaches the surface, it erupts from the vent as a lava flow. Unloading of the conduit magma in this way results in an abrupt pressure drop (i.e., the overburden felt by the bubbles is reduced), causing explosive gas expansion in the form of gas jetting from the vent. This terminates the event and lava drains back into the conduit to start the cycle anew. In the case where there is no surface lava emission or drain back, the cap instead pushes into and spreads out within a subsurface cavity. Again, this unloads the conduit magma and terminates in explosive gas release. Once gas is expelled, lava in the cavity is free to drain back. We hypothesize that pistoning is a stable mode of degassing for low-viscosity basaltic magmas with appropriate conduit geometries and volatile supply rates. Copyright 2005 by the American Geophysical Union.

Hawaii↗

USGS1990

Maceral and palynomorph facies from two tertiary peat-forming environments in the Powder River Basin, U.S.A.

The differences between the depositional settings of the Smith and Anderson subbituminous coal beds (Paleocene, central Powder River Basin, U.S.A.) are interpreted on the basis of their petrographic composition and palynologic assemblages. The Smith coal bed is relatively thin (< 5 m ) "> (< 5m) and has a high degree of thickness variation (0–5 m) over short distances (< 8 km ) "> (< 8km) , the result of deposition in an anastomosed fluvial environment where numerous lakes and small channels limited the extent of peat deposits. Although plants related to the living genus Glyptostrobus dominated the swamp, peat-forming plant communities contained mixes of ancestral species of Platanus and Ulmus as well as vegetation of lower stature such as Spagnum . Some intervals within the coal bed contain anomalously high concentrations of cell walls (humotelinite) and cell fillings (corpohuminite), some of which came from locally abundant, decay-resistant Pinaceae (or pinaceous) vegetation. Raised areas of the peat swamp characterized by Sphagnum were also favorable for the accumulation of carbonized plant components (inertinites). Because the peat of the Anderson coal bed formed on top of thick (> 50 m ) "> (> 50m) sandstone bodies of an abandoned meander-belt complex, the coal bed is generally thicker (> 7 m ) "> (> 7m) and more widespread (> 15 km ) "> (> 15km) in extent than the Smith coal bed. The sands provided a relatively stable, poorly compactable platform that was favorable to the growth of large, arborescent vegetation, such as the dominant ancestral Glyptostrobus , as well as ancestral Nyssa, Carya , and Betulaceae in a well-drained but moist swamp environment. The stability of the peat-forming environment resulted in a raised peat deposit of relatively uniform paleoflora and peat composition. In the thicker areas of the Anderson coal bed, the upward increase in carbonized plant components indicates a progressively drier or better-drained swamp environment. Intervals within the coal bed that overlie or are lateral to crevasse-splay deposits contain a high concentration of pollen attributable to Pterocarya and an absence of carbonized plant remains, an indication that ancestral Pterocarya preferred a water-saturated environment close to the edge of the swamp where detrial influx occurred.

International Journal of Coal Geology↗

USGS2012

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

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

California↗

USGS2007

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

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

Journal of Geophysical Research B: Solid Earth↗

USGS2022

Guidelines for volcano-observatory operations during crises: Recommendations from the 2019 Volcano Observatory Best Practices meeting

In November 2019, the fourth meeting on Volcano Observatory Best Practices workshop was held in Mexico City as a series of talks, discussions, and panels. Volcanologists from around the world offered suggestions for ways to optimize volcano-observatory crisis operations. By crisis, we mean unrest that may or may not lead to eruption, the eruption itself, or its aftermath, all of which require analysis and communications by the observatory. During a crisis, the priority of the observatory should be to acquire, process, analyze, and interpret data in a timely manner. A primary goal is to communicate effectively with the authorities in charge of civil protection. Crisis operations should rely upon exhaustive planning in the years prior to any actual unrest or eruptions. Ideally, nearly everything that observatories do during a crisis should be envisioned, prepared, and practiced prior to the actual event. Pre-existing agreements and exercises with academic and government collaborators will minimize confusion about roles and responsibilities. In the situation where planning is unfinished, observatories should prioritize close ties and communications with the land and civil-defense authorities near the most threatening volcanoes. To a large extent, volcanic crises become social crises, and any volcano observatory should have a communication strategy, a lead communicator, regular status updates, and a network of colleagues outside the observatory who can provide similar messaging to a public that desires consistent and authoritative information. Checklists permit tired observatory staff to fulfill their duties without forgetting key communications, data streams, or protocols that need regular fulfilment (Bretton et al. 2018; Newhall et al. 2020). Observatory leaders need to manage staff workload to prevent exhaustion and ensure that expertise is available as needed. Event trees and regular group discussions encourage multi-disciplinary thinking, consideration of disparate viewpoints, and documentation of all group decisions and consensus. Though regulations, roles and responsibilities differ around the world, scientists can justify their actions in the wake of an eruption if they document their work, are thoughtful and conscientious in their deliberations, and carry out protocols and procedures developed prior to volcanic unrest. This paper also contains six case studies of volcanic eruptions or observatory actions that illustrate some of the topics discussed herein. Specifically, we discuss Ambae (Vanuatu) in 2017–2018, Kīlauea (USA) in 2018, Etna (Italy) in 2018, Bárðarbunga (Iceland) in 2014, Cotopaxi (Ecuador) in 2015, and global data sharing to prepare for eruptions at Nyiragongo (Democratic Republic of Congo).

Hawaii↗

USGS2016

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

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

Journal of Volcanology and Geothermal Research↗

USGS2008

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

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

Journal of Volcanology and Geothermal Research↗

USGS2012

Interpretation of changes in water level accompanying fault creep and implications for earthquake prediction

Quantitative calculations for the effect of a fault creep event on observations of changes in water level in wells provide an approach to the tectonic interpretation of these phenomena. For the pore pressure field associated with an idealized creep event having an exponential displacement versus time curve, an analytic expression has been obtained in terms of exponential-integral functions. The pore pressure versus time curves for observation points near the fault are pulselike; a sharp pressure increase (or decrease, depending on the direction of propagation) is followed by more gradual decay to the normal level after the creep event. The time function of the water level change may be obtained by applying the filter—derived by A. G. Johnson and others to determine the influence of atmospheric pressure on water level—to the analytic pore pressure versus time curves. The resulting water level curves show a fairly rapid increase (or decrease) and then a very gradual return to normal. The results of this analytic model do not reproduce the steplike changes in water level observed by Johnson and others. If the procedure used to obtain the water level from the pore pressure is correct, these results suggest that steplike changes in water level are not produced by smoothly propagating creep events but by creep events that propagate discontinously, by changes in the bulk properties of the region around the well, or by some other mechanism. In addition, simplistic calculations show that significant pressure field variations and water level changes near the surface may be expected to accompany a propagating creep event on a buried fault. Water level changes of as much as several meters may be expected at the surface for a creep event having a dislocation amplitude of l m on a semi-infinite dislocation surface extending downward from a depth of 10 km and propagating horizontally at a rate of 10 km/day. The maximum near-surface effect should be observed at a horizontal distance from the fault about equal to the depth of the top of the dislocation surface. These results are consistent with the observations made in China of large water level changes preceding large earthquakes, if some sort of aseismic creep event at depth precedes these earthquakes.

Journal of Geophysical Research Solid Earth↗

USGS1996

Depositional controls on coal distribution and quality in the Eocene Brunner Coal Measures, Buller Coalfield, South Island, New Zealand

The Buller Coalfield on the West Coast of the South Island, New Zealand, contains the Eocene Brunner Coal Measures. The coal measures unconformably overlie Paleozoic-Cretaceous basement rocks and are conformably overlain by, and laterally interfinger with, the Eocene marine Kaiata Formation. This study examines the lithofacies frameworks of the coal measures in order to interpret their depositional environments. The lower part of the coal measures is dominated by conglomeratic lithofacies that rest on a basal erosional surface and thicken in paleovalleys incised into an undulating peneplain surface. These lithofacies are overlain by sandstone, mudstone and organic-rich lithofacies of the upper part of the coal measures. The main coal seam of the organic-rich lithofacies is thick (10-20 m), extensive, locally split, and locally absent. This seam and associated coal seams in the Buller Coalfield are of low- to high-volatile bituminous rank (vitrinite reflectance between 0.65% and 1.75%). The main seam contains a variable percentage of ash and sulphur. These values are related to the thickening and areal distribution of the seam, which in turn, were controlled by the nature of clastic deposition and peat-forming mire systems, marine transgression and local tidal incursion. The conglomeratic lithofacies represent deposits of trunk and tributary braided streams that rapidly aggraded incised paleovalleys during sea-level stillstands. The main seam represents a deposit of raised mires that initially developed as topogenous mires on abandoned margins of inactive braidbelts. Peat accumulated in mires as a response to a rise in the water table, probably initially due to gradual sea-level rise and climate, and the resulting raised topography served as protection from floods. The upper part of the coal measures consists of sandstone lithofacies of flu vial origin and bioturbated sandstone, mudstone and organic-rich lithofacies, which represent deposits of paralic (deltaic, barrier shoreface, tidal and mire) and marine environments. The fluvial sandstone lithofacies accumulated in channels during a sea-level stillstand. The channels were infilled by coeval braided and meandering streams prior to transgression. Continued transgression, ranging from tidal channel-estuarine incursions to widespread but uneven paleoshoreline encroachment, accompanied by moderate basin subsidence, is marked by a stacked, back-stepping geometry of bioturbated sandstone and marine mudstone lithofacies. Final retrogradation (sea-level highstand) is marked by backfilling of estuaries and by rapid landward deposition of the marine Kaiata Formation in the late Eocene.

International Journal of Coal Geology↗