Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Structure”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 1,207 records · Page 67Linked to original sources
USGS2018

Ghost forests of Marco Island: Mangrove mortality driven by belowground soil structural shifts during tidal hydrologic alteration

Land use changes often create in situ stress and eventual mortality in mangroves as unsuitable hydroperiods are created through tidal flow alterations. Here, we document mangrove forest and soil structural changes within transects established in tidally restricted areas on Marco Island (Collier County, Florida, USA), which has broad swaths of dead-standing or unhealthy mangroves (“ghost forests”). Transects (N = 4) were arranged to include full canopy, transitional, and open canopy (dead) forests, and compared to nearby reference forests. Aboveground and belowground carbon (C) stocks (Total C) ranged from 288 to 304 Mg C ha −1 on full canopy, transitional, and reference sites, which did not differ from each other. However, Total C was lower for dead sites (233 Mg C ha −1 ) dictated entirely by differences in aboveground C (live and dead trees, downed wood); no differences were found among forest condition in belowground C stocks. This belowground C has been persistent in the soil for 85 years since initial tidal restriction. Nevertheless, hydrologic rehabilitation has the potential to increase total C stocks on dead sites by 70–110 Mg C ha −1 . Collapse of the soil surface by 6–8 cm just under the active root zone in chronically stressed mangroves was evident within the bulk density profiles from transitional versus dead sites, suggesting that surface elevation loss as root turnover ceases may work correlatively with chronic stressors (anoxia, P limitation) to affect rapid mortality of forests over short periods of time many years after stress initiation. Hydrologic rehabilitation of stressed or denuded mangroves must also include an understanding of how these soil processes might be re-established.

Florida↗
USGS2024

Prototyping structured decision making for water resource management in the San Francisco Bay-Delta

A structured decision making (SDM) approach can help evaluate tradeoffs between conservation and human-benefit objectives by fostering communication and knowledge transfer among stakeholders, decision makers, and the public. However, the process is iterative and completing the full process may take years. It can be difficult to initiate an SDM effort when problems seem insurmountable. Occasionally, SDM may not even be the best or correct approach for addressing the conservation problem at hand. We describe the implementation of an SDM process to help inform difficult decisions related to competing objectives. We convened a diverse stakeholder group from the largest estuary in the western United States; the San Francisco Bay and Sacramento-San Joaquin Delta (Bay-Delta). The stakeholder group consisted of representatives from local, state, and federal agencies, non-profit organizations, and recreational fishers. The stakeholder group agreed on a problem statement and identified four priority objectives related to Chinook salmon, delta smelt, water availability and reliability, and agricultural water use. Furthermore, they proposed 14 candidate management actions to achieve their objectives. The group then used existing quantitative models and data to evaluate trade-offs in proposed management actions to identify areas of agreement of proposed candidate actions. The clear communication of the problem statement and objectives among the stakeholder group, along with evaluation of tradeoffs and uncertainty via decision-support models suggest that a full SDM approach may work in the Bay-Delta. We further communicate lessons learned during our implementation of SDM to help guide future SDM efforts in the region and elsewhere.

California↗
USGS2006

Crustal structure across the Altyn Tagh Range at the northern margin of the Tibetan Plateau and tectonic implications

We present new seismic refraction/wide-angle-reflection data across the Altyn Tagh Range and its adjacent basins. We find that the crustal velocity structure, and by inference, the composition of the crust changes abruptly beneath the Cherchen fault, i.e., ∼100 km north of the northern margin of the Tibetan plateau. North of the Cherchen fault, beneath the Tarim basin, a platform-type crust is evident. In contrast, south the Cherchen fault the crust is characterized by a missing high-velocity lower-crustal layer. Our seismic model indicates that the high topography (∼3 km) of the Altyn Tagh Range is supported by a wedge-shaped region with a seismic velocity of 7.6–7.8 km/s that we interpret as a zone of crust–mantle mix. We infer that the Altyn Tagh Range formed by crustal-scale strike-slip motion along the North Altyn Tagh fault and northeast–southwest contraction over the range. The contraction is accommodated by (1) crustal thickening via upper-crustal thrusting and lower-crustal flow (i.e., creep), and (2) slip-parallel (SW-directed) underthrusting of only the lower crust and mantle of the eastern Tarim basin beneath the Altyn Tagh Range.

Earth and Planetary Science Letters↗
USGS2006

Crustal structure and tectonics of the northern part of the Southern Granulite Terrane, India

Deep seismic reflection studies investigating the exposed Archean lower continental crust of the Southern Granulite Terrane, India, yield important constraints on the nature and evolution of the deep crust, including the formation and exhumation of granulites. Seismic reflection images along the Kuppam–Bhavani profile reveal a band of reflections that dip southward from 10.5 to 15.0 s two-way-time (TWT), across a distance of 50 km. The bottom of these reflections beneath the Dharwar craton is interpreted as the Moho. Further south, another reflection band dipping northward is observed. These bands of reflectivity constitute a divergent reflection fabric that converges at the Moho boundary observed at the Mettur shear zone. Reflection fabrics that intersect at a steep angle are interpreted as a collisional signature due to the convergence of crustal blocks, which we infer resulted in crustal thickening and the formation of granulites. Anomalous gravity and magnetic signatures are also observed across the Mettur shear zone. The gravity model derived from the Bouguer gravity data corroborates seismic results. The tectonic regime and seismic reflection profiles are combined in a 3-D representation that illustrates our evidence for paleo-subduction at a collision zone. The structural dissimilarities and geophysical anomalies suggest that the Mettur shear zone is a suture between the Dharwar craton in the north and another crustal block in the south. This study contributes significantly to our understanding of the operation of Archean plate tectonics, here inferred to involve collision and subduction. Furthermore, it provides an important link between the Gondwanaland and global granulite evolution occurring throughout the late Archean.

Earth and Planetary Science Letters↗
USGS2019

A comparison of age- and size-structured assessment models applied to a stock of cisco in Thunder Bay, Ontario

Stock assessments are critical to modern fisheries management, supporting the calculation of key reference variables used to make informed management decisions. However, there is still considerable uncertainty as to which class of assessment models is appropriate to use under different circumstances. A common class of models used when age data are available are statistical catch-at-age assessment (SCAA) models, which track annual cohorts through time. When age data are unavailable, as is often the case in invertebrate fisheries where the lack of a bony structure such as otoliths makes aging difficult, statistical catch-at-size assessment (SCSA) models are more often employed, tracking fish or invertebrates through time by size-classes rather than ages. Do SCAA models actually perform better than SCSA models when age data are available, or is this just an assumption we make in fisheries research and management? We examined this question by evaluating the effectiveness of both SCAA and SCSA models in characterizing cisco, Coregonus artedi , population dynamics in Thunder Bay, Ontario. Both models were fit using an integrated framework with multiple sources of data including hydroacoustic estimates of spawning stock, fishery-dependent and -independent age/length compositions, and harvest data. Our results suggest that for cisco in Thunder Bay, data-limitations related to lack of size-composition data over the size range for which cisco growth is rapid resulted in difficulty estimating relative year-class strength within a SCSA. This led to parameter confounding and ultimately the inability to estimate natural mortality within a SCSA. This hampered the utility of a SCSA model in comparison with a SCAA model when age-composition data were available.

Ontario↗
USGS2019

Surface water connectivity controls fish food web structure and complexity across local- and meta-food webs in Arctic Coastal Plain lakes

The need for theories that address food web assembly and complexity over multiple spatial scales are critical to understanding their stability and persistence. In a meta-food web – an integrated network of local food webs – spatial heterogeneity in physical processes may have profound effects on food web function and energy flow. In the Arctic, surface water connectivity plays a vital role in determining fish assemblage composition, and potentially, food web structure. We examined lentic food web complexity associated with heterogeneity in surface water connectivity among Arctic lakes at the at the local scale, by contrasting lakes over a stream-lake connectivity gradient, and at the regional scale, by contrasting two locations with different surface water conditions (i.e., wet and dry) on the Arctic Coastal Plain of Alaska. Among lakes and across locations, increased hydrologic connectivity between streams and lakes increased the number of fish species and increased the complexity of the food web. The interaction of the region’s hydrologic connectivity, local stream-lake connections, and the trophic niches of relevant fish species produced integrated, complex meta-food webs. Fully understanding mechanisms that support meta-food web stability are crucial when assessing future changes to Arctic stream-lake networks and the function and persistence of aquatic food webs.

Food Webs↗
USGS2020

Quantifying 40 years of rockfall activity in Yosemite Valley with historical Structure-from-Motion photogrammetry and terrestrial laser scanning

Rockfalls and rockslides are often dominant geomorphic processes in steep bedrock landscapes, but documenting their occurrence can be challenging, requiring frequent monitoring and well resolved spatial data. Repeat application of remote sensing methods such as Terrestrial Laser Scanning (TLS) and Structure-from-Motion (SfM) photogrammetry can detect even very small rockfalls, but typically these acquisitions span only years and may not record rockfall activity representative of longer-term rates of cliff erosion. Inventory databases can extend rockfall records, but are commonly incomplete and prone to observation bias. We employed TLS and SfM on two adjacent cliffs (El Capitan and Middle Brother) in Yosemite Valley, integrating semi-annual data collections from 2010 to 2017 with “historical” (archival) SfM models derived from oblique photographs taken in 1976. Comparing the 1976 SfM models against more recent data allows for more accurate and precise rockfall detection and volume measurement over a 40-year period. Change detection indicates that 235 rockfalls occurred from the two cliffs, more than twice as many events as are recorded in Yosemite's inventory database. Although individual rockfall volumes reported in the inventory database vary from those measured by SfM-TLS, reported cumulative volumes are similar to measured volumes, likely because the large-volume events that account for most of the cumulative volume tend to be widely observed and well-documented. Volume-frequency relationships indicate that the cliffs erode predominantly by less frequent, larger-volume rockfalls, at rates of 0.9 to 1.7 mm/yr. Our study demonstrates how integrated SfM and TLS measurements, especially utilizing SfM models derived from historical imagery, allow detection and quantification of rockfalls spanning several decades, complementing and improving inventory databases, informing rockfall hazard assessment, and providing longer-term rates of cliff erosion.

California↗
USGS2016

Helminth community structure in two species of arctic-breeding waterfowl

Climate change is occurring rapidly at high latitudes, and subsequent changes in parasite communities may have implications for hosts including wildlife and humans. Waterfowl, in particular, harbor numerous parasites and may facilitate parasite movement across broad geographic areas due to migratory movements. However, little is known about helminth community structure of waterfowl at northern latitudes. We investigated the helminth communities of two avian herbivores that breed at high latitudes, Pacific black brant ( Branta bernicla nigricans ), and greater white-fronted geese ( Anser albifrons ), to examine effects of species, geographic area, age, and sex on helminth species richness, aggregation, prevalence, and intensity. We collected 83 and 58 black brant and white-fronted geese, respectively, from Arctic and Subarctic Alaska July-August 2014. We identified 10 known helminth species ( Amidostomum anseris , Amidostomum spatulatum , Drepanidotaenia lanceolata , Epomidiostomum crami , Heterakis dispar , Notocotylus attenuatus , Tetrameres striata , Trichostrongylus tenuis , Tschertkovilepis setigera , and Wardoides nyrocae ) and 1 previously undescribed trematode. All geese sampled were infected with at least one helminth species. All helminth species identified were present in both age classes and species, providing evidence of transmission at high latitudes and suggesting broad host susceptibility. Also, all but one helminth species were present at both sites, suggesting conditions are suitable for transmission across a large latitudinal/environmental gradient. Our study provides important baseline information on avian parasites that can be used to evaluate the effects of a changing climate on host-parasite distributions.

Alaska↗
USGS2022

Genome-wide genetic diversity may help identify fine-scale genetic structure among lake whitefish spawning groups in Lake Erie

In Lake Erie, lake whitefish Coregonus clupeaformis supported lucrative fisheries before populations were decimated by overfishing and water quality degradation. In recent years, there has been a renewed interest in lake whitefish and management of the fishery they support. Lake whitefish spawn on several reefs throughout Lake Erie, but the relative recruitment dynamics and contributions of spawning groups to the fishery are not well understood. Modern high-throughput sequencing approaches offer new opportunities to census population diversity and to identify subtle differences among closely related populations. We used high-throughput sequencing data to evaluate the genetic structure and diversity of lake whitefish collected opportunistically across broad spatial scales in Lake Erie. Using RAD-capture (Rapture), we sequenced and genotyped individuals (N = 88) from the west, central, and east basin of Lake Erie at 120,268 single nucleotide polymorphisms (SNPs). Lake whitefish from Niagara and Crib Reefs (west basin) diverged from the three collections. Interestingly, these were the only lake whitefish collected during the act of spawning (late November), and all other fish were collected pre-spawn (August-early November). These results suggest that some lake whitefish spawning reefs may be reproductively isolated, though definition of these groups into stocks will require more intentional sampling during the act of spawning.

Michigan, New York, Ohio, Ontario, Pennsylvania↗
USGS2017

Structured decision making as a conservation tool for recovery planning of two endangered salamanders

At least one-third of all amphibian species face the threat of extinction, and current amphibian extinction rates are four orders of magnitude greater than background rates. Preventing extirpation often requires both ex situ (i.e., conservation breeding programs) and in situ strategies (i.e., protecting natural habitats). Flatwoods salamanders ( Ambystoma bishopi and A. cingulatum ) are protected under the U.S. Endangered Species Act. The two species have decreased from 476 historical locations to 63 recently extant locations (86.8% loss). We suggest that recovery efforts are needed to increase populations and prevent extinction, but uncertainty regarding optimal actions in both ex situ and in situ realms hinders recovery planning. We used structured decision making (SDM) to address key uncertainties regarding both captive breeding and habitat restoration, and we developed short-, medium-, and long-term goals to achieve recovery objectives. By promoting a transparent, logical approach, SDM has proven vital to recovery plan development for flatwoods salamanders. The SDM approach has clear advantages over other previous approaches to recovery efforts, and we suggest that it should be considered for other complex decisions regarding endangered species.

Journal for Nature Conservation↗
USGS2003

Reinterpretation of the stratigraphy and structure of the Rancho Las Norias area, central Sonora, Mexico

New geologic mapping and fossil data in the vicinity of Rancho Las Norias, 30 km east of Hermosillo, Sonora, Mexico, show that rocks previously mapped as Precambrian instead are Paleozoic. Previous geologic maps of the Rancho Las Norias area show northeast-directed, southwest-dipping reverse or thrust faults deforming both Precambrian and Paleozoic rocks. The revised stratigraphy requires reinterpretation of some of these faults as high-angle normal or oblique-slip faults and the elimination of other faults. We agree with earlier geologic map interpretations that compressional structures have affected the Paleozoic rocks in the area, but our mapping suggests that the direction of compression is from southeast to northwest. Published by Elsevier Ltd.

Journal of South American Earth Sciences↗
USGS2013

Seismic structure of the crust and uppermost mantle of South America and surrounding oceanic basins

We present a new set of contour maps of the seismic structure of South America and the surrounding ocean basins. These maps include new data, helping to constrain crustal thickness, whole-crustal average P-wave and S-wave velocity, and the seismic velocity of the uppermost mantle ( P n and S n ). We find that: (1) The weighted average thickness of the crust under South America is 38.17 km (standard deviation, s.d. ±8.7 km), which is ∼1 km thinner than the global average of 39.2 km (s.d. ±8.5 km) for continental crust. (2) Histograms of whole-crustal P-wave velocities for the South American crust are bi-modal, with the lower peak occurring for crust that appears to be missing a high-velocity (6.9–7.3 km/s) lower crustal layer. (3) The average P-wave velocity of the crystalline crust ( P cc ) is 6.47 km/s (s.d. ±0.25 km/s). This is essentially identical to the global average of 6.45 km/s. (4) The average P n velocity beneath South America is 8.00 km/s (s.d. ±0.23 km/s), slightly lower than the global average of 8.07 km/s. (5) A region across northern Chile and northeast Argentina has anomalously low P- and S-wave velocities in the crust. Geographically, this corresponds to the shallowly-subducted portion of the Nazca plate (the Pampean flat slab first described by Isacks et al., 1968), which is also a region of crustal extension. (6) The thick crust of the Brazilian craton appears to extend into Venezuela and Colombia. (7) The crust in the Amazon basin and along the western edge of the Brazilian craton may be thinned by extension. (8) The average crustal P-wave velocity under the eastern Pacific seafloor is higher than under the western Atlantic seafloor, most likely due to the thicker sediment layer on the older Atlantic seafloor.

Bulletin of the Seismological Society of America↗
USGS2005

Petrologic constraints on the thermal structure of the Cascades arc

Primitive late Cenozoic basaltic lavas from the Cascades volcanic arc near latitude 46°N comprise two distinct compositional groups. Group I includes samples with low Ba/Nb (<20) and other compositional similarities to oceanic island and MORB lavas from within-plate settings. In contrast, Group II exhibits enrichment of Ba and large-ion lithophile elements (LILE) and depletion of Nb and high-field strength elements (HFSE) as seen commonly in calcalkalic lavas from other volcanic arcs. Lavas of both groups are widely distributed across the transect, and Group I lavas are found as much as 30–40 km trenchward of stratovolcanoes that define the High Cascades ‘volcanic front (VF)’. The most primitive lavas are sparsely porphyritic, have elevated Ni, Cr, and Mg#, high calculated magmatic temperatures (1200–1300 °C), and lack evidence of shallow (crustal level) storage and crystallization. Compositions of parental liquids were calculated for each primitive sample on the premise of Fe–Mg equilibrium with mantle peridotite. Assuming that such magmas ascended rapidly from accumulation zones in the mantle, we estimate P and T of segregation. We infer that (a) Group I magmas ascended from systematically greater depths (∼50–70 km) than Group II (∼30–50 km), implying the possible existence of compositional stratification in the mantle wedge; (b) Group I basalts show the least evidence for slab-derived contributions in their sources despite their apparently greater segregation depths (approaching the locus of the Cascadia slab beneath the frontal arc region); (c) Group II lavas with the strongest slab compositional signature have temperatures far exceeding the wet peridotite solidus at high pressure; and (d) the inferred thermal structure of the mantle wedge is very warm, implying a significant component of mantle upwelling and convection. Group I lavas are interpreted as decompression melts from this mantle, and their compositions suggest that their source was little modified by slab-derived contributions. We speculate that melting to produce Group II magmas occurs in the shallow mantle, possibly in response to heating by hot ascending Group I magmas. If true, it seems unlikely that the slab-like signal in Group II lavas can be attributed to modern slab inputs; rather, we postulate that this signature may reflect melting of lithospheric mantle domains containing a ‘stored’ slab-derived component inherited from earlier stages of Cascadia subduction. This scenario differs from the standard paradigm for subduction zones (SZs), and stresses the importance of convecting asthenospheric mantle in driving arc magmatism, particularly in warm subduction zones where slab fluid contributions likely are minimal. In contrast, because tectonic conditions in more typical volcanic arcs favor subduction of cooler, less dehydrated oceanic lithosphere, slab-derived fluids may promote extensive flux-melting in the wedge. Such melts may dominate the magmatic output and mask wedge contributions. The Cascade arc thus provides rarely afforded insights into arc magma genesis.

British Columbia, California, Oregon, Washington↗
USGS2007

Three-dimensional P-wave velocity structure derived from local earthquakes at the Katmai group of volcanoes, Alaska

The three-dimensional P-wave velocity structure beneath the Katmai group of volcanoes is determined by inversion of more than 10,000 rays from over 1000 earthquakes recorded on a local 18 station short-period network between September 1996 and May 2001. The inversion is well constrained from sea level to about 6 km below sea level and encompasses all of the Katmai volcanoes; Martin, Mageik, Trident, Griggs, Novarupta, Snowy, and Katmai caldera. The inversion reduced the average RMS travel-time error from 0.22 s for locations from the standard one-dimensional model to 0.13 s for the best three-dimensional model. The final model, from the 6th inversion step, reveals a prominent low velocity zone (3.6–5.0 km/s) centered at Katmai Pass and extending from Mageik to Trident volcanoes. The anomaly has values about 20–25% slower than velocities outboard of the region (5.0–6.5 km/s). Moderately low velocities (4.5–6.0 km/s) are observed along the volcanic axis between Martin and Katmai Caldera. Griggs volcano, located about 10 km behind (northwest of) the volcanic axis, has unremarkable velocities (5.0–5.7 km/s) compared to non-volcanic regions. The highest velocities are observed between Snowy and Griggs volcanoes (5.5–6.5 km/s). Relocated hypocenters for the best 3-D model are shifted significantly relative to the standard model with clusters of seismicity at Martin volcano shifting systematically deeper by about 1 km to depths of 0 to 4 km below sea level. Hypocenters for the Katmai Caldera are more tightly clustered, relocating beneath the 1912 scarp walls. The relocated hypocenters allow us to compare spatial frequency-size distributions ( b -values) using one-dimensional and three-dimensional models. We find that the distribution of b is significantly changed for Martin volcano, which was characterized by variable values (0.8 < b < 2.0) with standard locations and more uniform values (0.8 < b < 1.2) after relocation. Other seismic clusters at Mageik (1.2 < b < 2.2), Trident (0.5 < b < 1.5) and Katmai Caldera (0.8 < b < 1.8) had stable b -values indicating the robustness of the observations. The strong high b -value region at Mageik volcano is mainly associated with an earthquake swarm in October, 1996 that possibly indicates a shallow intrusion or influx of gas. The new velocity and spatial b -value results, in conjunction with prior gravity (Bouguer anomalies up to − 40 mgal) and interferometry (several cm uplift) data, provide strong evidence in favor of partially molten rock at shallow depths beneath the Mageik–Katmai–Novarupta region. Moderately low velocities beneath Martin and Katmai suggest that old, mostly solidified intrusions exist beneath these volcanoes. Higher relative velocities beneath the Griggs and Snowy vents suggest that no magma is resident in the shallow crust beneath these volcanoes.

Journal of Volcanology and Geothermal Research↗
USGS2014

Seismicity and seismic structure at Okmok Volcano, Alaska

Okmok volcano is an active volcanic caldera located on the northeastern portion of Umnak Island in the Aleutian arc, with recent eruptions in 1997 and 2008. The Okmok area had ~900 locatable earthquakes between 2003 and June 2008, and an additional ~600 earthquakes from the beginning of the 2008 eruption to mid 2009, providing an adequate dataset for seismic tomography. To image the seismic velocity structure of Okmok, we apply waveform cross-correlation using bispectrum verification and double-difference tomography to a subset of these earthquakes. We also perform P-wave attenuation tomography using a spectral decay technique. We examine the spatio-temporal characteristics of seismicity in the opening sequence of the 2008 eruption to investigate the path of magma migration during the establishment of a new eruptive vent. We also incorporate the new earthquake relocations and three-dimensional (3D) velocity model with first-motion polarities to compute focal mechanisms for selected events in the 2008 pre-eruptive and eruptive periods. Through these techniques we obtain precise relocations, a well-constrained 3D P-wave velocity model, and a marginally resolved S-wave velocity model. We image a main low Vp and Vs anomaly directly under the caldera consisting of a shallow zone at 0–2 km depth connected to a larger deeper zone that extends to about 6 km depth. We find that areas of low Qp are concentrated in the central to southwestern portion of the caldera and correspond fairly well with areas of low Vp. We interpret the deeper part of the low velocity anomaly (4–6 km depth) beneath the caldera as a magma body. This is consistent with results from ambient noise tomography and suggests that previous estimates of depth to Okmok's magma chamber based only on geodetic data may be too shallow. The distribution of events preceding the 2008 eruption suggest that a combination of overpressure in the zone surrounding the magma chamber and the introduction of new material from below were jointly responsible for the explosive eruption. Magma escaping from the top of the main magma chamber likely reacted with both a smaller shallow pod of magma and groundwater on its way up below the Cone D area. The earthquakes in the 2008 pre-eruptive and eruptive periods are found to have a mixture of strike-slip, oblique normal, and oblique thrust mechanisms, with a dominant P-axis orientation that is nearly perpendicular to the regional tectonic stress. This may indicate that the stresses related to magmatic activity locally dominated regional tectonic forces during this time period.

Alaska↗
USGS2017

Early postcaldera rhyolite and structural resurgence at Long Valley Caldera, California

After the 767-ka caldera-forming eruption of 650 km 3 of rhyolite magma as the Bishop Tuff, 90–100 km 3 of similar rhyolite erupted in the west-central part of Long Valley caldera in as many as 40 batches spread over the 110,000-year interval from ~ 750 ka to ~ 640 ka. Centrally, this Early Rhyolite (ER) is as thick as 622 m, but it spread radially to cover much of the caldera floor, where half its area is now concealed by post-ER sediments and lavas. At least 75% of the ER is aphyric rhyolite tuff. Drillholes encountered ~ 22 (altered) ER lava flows intercalated in the pyroclastic pile, and another 11 units of (largely fresh) ER lava are exposed on the caldera's resurgent dome and at Lookout Mountain. Exposed units have been distinguished, mapped, studied petrographically and chemically, and radioisotopically dated; each is described in detail. Their phenocryst contents range from 0 to 2.5 wt%. All the phyric units have plagioclase, orthopyroxene, and ilmenite; most have biotite and rare tiny magnetite, and a few contain rare zircon. The compositional range of fresh obsidians is narrow—74.3–75.0% SiO 2 , 1.21–1.37% FeO*, and 5.12–5.26% K 2 O, but wider variations in Ti, Ba, Sr, and Zr permit distinction of individual units and eruptive groups. The limited chemical and petrographic variability shown by so many ER batches released episodically for ~ 110,000 years suggests a thermally buffered and well-stirred reservoir. The ER central area, where ER eruptions had taken place, was uplifted ~ 400 m to form a structural dome ~ 10 km in diameter. Most of the inflation is attributable to 10 sills of ER that intrude the Bishop Tuff beneath the uplift, but other processes potentially contributing to resurgence are also considered. As shown by erratics of Mesozoic rocks ice-rafted from the Sierra Nevada and dropped on ER lavas, much of the ER had erupted early enough and at low enough elevation to be inundated by the intracaldera lake and was only later lifted by the resurgence that also raised clusters of the erratics hundreds of meters higher than any shoreline. Most of the uplift was over by ~ 570 ka, but dome-crossing faults that exhibit normal throw of 10–30 m cut lavas as young as 175–125 ka. For most elements, chemical ranges of the ER lie within those of the zoned Bishop Tuff, which had erupted earlier from the same place. Only Ba, Zr, Hf, and Eu/Eu* extend to ranges outside those of the Bishop Tuff, nominally to less evolved compositions. Initial 87 Sr/ 86 Sr values of ER are likewise within the range of the Bishop Tuff, but ER ratios of 143 Nd/ 144 Nd and 206 Pb/ 204 Pb extend beyond those of the Bishop Tuff to values slightly more influenced by upper-crustal contributions. FeTi-oxide geothermometry yields 752°–844 °C for ER, compared to 700°–820 °C for the Bishop Tuff. ER fO 2 values are 0.5–1.0 log units more reduced than those of the T–fO 2 array of the Bishop Tuff. The postcaldera reduction may reflect reaction with graphite from the black lithics of Paleozoic graphitic metapelite so abundant in the Bishop Tuff. Much of the pumice emplaced during the later half of the Bishop Tuff eruption has 10–25 wt% phenocrysts, dominantly quartz and sanidine, but the 100 km 3 of ER has only 0–2.5 wt% and completely lacks quartz and sanidine. Postcaldera processes, including mixing, volatile ascent, and crystal resorption, as well as potential contaminants and magmatic inputs, are all considered.

California↗
USGS2020

High-resolution imaging of hydrothermal heat flux using optical and thermal Structure-from-Motion photogrammetry

Quantifying hydrothermal heat flux at meter-scale resolution over N0.25 km2 is required to bridge in-situ heat flux and satellite-based measurements. We advance a methodology that blends ground-based daytime optical and nighttime thermal infrared (TIR) imagery using Structure-from-Motion photogrammetry to map radiant hydrothermal heat flux over these scales at sites with low signal-to-noise ratios that would otherwise be difficult to characterize using, for example, unmanned aerial systems. The improved method uses a computerized telescopic mount to relocate and align daytime optical acquisitions with nighttime TIR imagery, thereby enabling TIR acquisition from multiple camera orientations positioned throughout a study area. This facilitates mapping of thermal features at sites of varying size and complexity and helps to ameliorate topographic occlusion effects and geometric distortions that can bias radiant hydrothermal heat flux estimates derived from the resulting orthorectified thermal maps. We assessed detection thresholds of this method at three sites across central California, which range in size, topography, and heat flux conditions. We found that blending of optical and thermal acquisitions successfully detected anomalous heat flux, even in cases where temperatures were only slightly greater than background. This approach might be applied to a variety of volcanic and hydrothermal systems to quantify the spatial distribution of heat flux, and how this may relate to factors such as the distribution of ground fractures and lava flow rheology.

Journal of Volcanology and Geothermal Research↗