SEARCH · Geology Reports

Results for “Journal 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,189 records · Page 66Linked to original sources

USGS1996

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

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

Journal of Geophysical Research B: Solid Earth↗

USGS2012

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

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

Michigan, Minnesota, Ontario, Wisconsin↗

USGS2012

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

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

Journal of Geophysical Research Solid Earth↗

USGS2012

Seismic refraction measurements within the Peninsular terrane, south central Alaska

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

Alaska↗

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↗

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↗

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↗

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↗

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↗

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↗

USGS2021

Late Holocene environmental change in Celestun Lagoon, Yucatan, Mexico

Epikarst estuary response to hydroclimate change remains poorly understood, despite the well-studied link between climate and karst groundwater aquifers. The influence of sea-level rise and coastal geomorphic change on these estuaries obscures climate signals, thus requiring careful development of paleoenvironmental histories to interpret the paleoclimate archives. We used foraminifera assemblages, carbon stable isotope ratios (δ 13 C) and carbon:nitrogen (C:N) mass ratios of organic matter in sediment cores to infer environmental changes over the past 5300 years in Celestun Lagoon, Yucatan, Mexico. Specimens (> 125 µm) from modern core top sediments revealed three assemblages: (1) a brackish mangrove assemblage of agglutinated Miliammina and Ammotium taxa and hyaline Haynesina (2) an inner-shelf marine assemblage of Bolivina , Hanzawaia , and Rosalina, and (3) a brackish assemblage dominated by Ammonia and Elphidium . Assemblages changed along the lagoon channel in response to changes in salinity and vegetation, i.e. seagrass and mangrove. In addition to these three foraminifera assemblages, lagoon sediments deposited since 5300 cal yr BP are comprised of two more assemblages, defined by Archaias and Laevipeneroplis, which indicate marine Thalassia seagrasses, and Trichohyalus, which indicates restricted inland mangrove ponds. Our data suggest that Celestun Lagoon displayed four phases of development: (1) an inland mangrove pond (5300 BP) (2) a shallow unprotected coastline with marine seagrass and barrier island initiation (4900 BP) (3) a protected brackish lagoon (3000 BP), and (4) a protected lagoon surrounded by mangroves (1700 BP). Stratigraphic (temporal) changes in core assemblages resemble spatial differences in communities across the modern lagoon, from the southern marine sector to the northern brackish region. Similar temporal patterns have been reported from other Yucatan Peninsula lagoons and from cenotes (Nichupte, Aktun Ha), suggesting a regional coastal response to sea level rise and climate change, including geomorphic controls (longshore drift) on lagoon salinity, as observed today. Holocene barrier island development progressively protected the northwest Yucatan Peninsula coastline, reducing mixing between seawater and rain-fed submarine groundwater discharge. Superimposed on this geomorphic signal, assemblage changes that are observed reflect the most severe regional wet and dry climate episodes, which coincide with paleoclimate records from lowland lake archives (Chichancanab, Salpeten). Our results emphasize the need to consider coastal geomorphic evolution when using epikarst estuary and lagoon sediment archives for paleoclimate reconstruction and provide evidence of hydroclimate changes on the Yucatan Peninsula.

Yucatan, Celestun Lagoon↗

USGS2016

Use of structured decision-making to explicitly incorporate environmental process understanding in management of coastal restoration projects: Case study on barrier islands of the northern Gulf of Mexico

Coastal ecosystem management typically relies on subjective interpretation of scientific understanding, with limited methods for explicitly incorporating process knowledge into decisions that must meet multiple, potentially competing stakeholder objectives. Conversely, the scientific community lacks methods for identifying which advancements in system understanding would have the highest value to decision-makers. A case in point is barrier island restoration, where decision-makers lack tools to objectively use system understanding to determine how to optimally use limited contingency funds when project construction in this dynamic environment does not proceed as expected. In this study, collaborative structured decision-making (SDM) was evaluated as an approach to incorporate process understanding into mid-construction decisions and to identify priority gaps in knowledge from a management perspective. The focus was a barrier island restoration project at Ship Island, Mississippi, where sand will be used to close an extensive breach that currently divides the island. SDM was used to estimate damage that may occur during construction, and guide repair decisions within the confines of limited availability of sand and funding to minimize adverse impacts to project objectives. Sand was identified as more limiting than funds, and unrepaired major breaching would negatively impact objectives. Repairing minor damage immediately was determined to be generally more cost effective (depending on the longshore extent) than risking more damage to a weakened project. Key gaps in process-understanding relative to project management were identified as the relationship of island width to breach formation; the amounts of sand lost during breaching, lowering, or narrowing of the berm; the potential for minor breaches to self-heal versus developing into a major breach; and the relationship between upstream nourishment and resiliency of the berm to storms. This application is a prototype for using structured decision-making in support of engineering projects in dynamic environments where mid-construction decisions may arise; highlights uncertainty about barrier island physical processes that limit the ability to make robust decisions; and demonstrates the potential for direct incorporation of process-based models in a formal adaptive management decision framework.

Mississippi↗