USGS Science⌕ Search

SEARCH · USGS Science

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,333 records · Page 74Linked to original sources

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

The structure and volume of large geysers in Yellowstone National Park, USA and the mineralogy and chemistry of their silica sinter deposits

Siliceous sinter is formed by biogenic and abiogenic opal deposition around hot springs and geysers. Using Structure-from-Motion photogrammetry we generated three-dimensional models of Giant and Castle Geysers from the Upper Geyser Basin of Yellowstone National Park. We use these models to calculate an approximate mass of sinter for each (~2 and ~ 5 kton, respectively) and estimate a range of plausible long-term deposition rates for Castle Geyser (470 to 940 kg·yr −1 ). We estimate ~2% of the silica discharged from Castle Geyser is deposited as sinter in the cone and proximal terraces. We collected 15 sinter samples following the stratigraphy of each geyser from an older terrace to a younger cone and examined them using a variety of analytical methods. We find that young opaline sinter with a water content of <12 wt% (from loss on ignition) contains higher concentrations of major and trace elements, notably As, Sb, Rb, Ga and Cs, relative to older dehydrated sinter. Rare earth element (REE) concentrations in sinter are 2–3 orders of magnitude higher than in the thermal water from which they are deposited. Sinter deposits are enriched in light REE, Gd and Yb when normalized to concentrations in thermal water and enriched in Eu, Tm, and Yb when normalized to the underlying rhyolite. Sinter samples with the highest REE concentrations are also enriched in organic material, implying either microbial uptake of REE, or that organic molecules are efficient ligands that form metal complexes.

Wyoming↗

Development, structure, and behavior of a perched lava channel at Kīlauea Volcano, Hawaiʻi, during 2007

Channelized lava flows are commonly produced during the early stages of basaltic eruptions. These channels usually maintain their morphology until the eruption ends or discharge is diverted. In some instances, narrower channels can roof over, developing into lava tubes. We report here on a channelized flow erupted at Kīlauea volcano in 2007 that evolved into a “perched lava channel” composed of a string of interconnected, elongate lava pools, forming a lava channel/lava pool hybrid. The lava channel, which had a time-averaged discharge rate of ~3–9 m 3 /s, initially fed a series of flow branches that exhibited cooling-limited and volume-limited controls on flow length, sometimes with each process controlling a different morphological aspect of a single flow branch. The perched lava channel grew vertically primarily by overplating of the channel levees from frequent overflows, forming a compound flow field. This vertical growth only occurred when the distal end of the channel was blocked. When levee failure at the distal end of the channel caused the lava level in the channel to drop below the levee rim, no vertical growth occurred. Seeps of spiny lava and slabby pāhoehoe were common, erupting from uplift scarps on the channel levees, apparently fed by sills from denser, relatively crystal-rich material filling the bottom of the channel. We infer that lava in the channel was stratified in vesicularity and velocity, with foamy, vesicular, faster-moving lava at the top of the lava stream and denser, relatively outgassed, slower-moving lava filling the bottom of the channel. The channel levees were unstable, failing on several occasions, perhaps triggered by the levee seeps. The appearance of seeps, therefore, is one way of assessing the collapse potential of similar perched lava structures.

Hawaii↗

Mapping the resistivity structure of Walker Ridge 313 in the Gulf of Mexico using the marine CSEM method

A marine controlled source electromagnetic (CSEM) campaign was carried out in the Gulf of Mexico to further develop marine electromagnetic techniques in order to aid the detection and mapping of gas hydrate deposits. Marine CSEM methods are used to obtain an electrical resistivity structure of the subsurface which can indicate the type of substance filling the pore space, such as gas hydrates which are more resistive. Results from the Walker Ridge 313 study (WR 313) are presented in this paper and compared with the Gulf of Mexico Gas Hydrate Joint Industry Project II (JIP2) logging while drilling (LWD) results and available seismic data. The hydrate, known to exist within sheeted sand deposits, is mapped as a resistive region in the two dimensional (2D) CSEM inversion models. This is consistent with the JIP2 LWD resistivity results. CSEM inversions that use seismic horizons provide more realistic results compared to the unconstrained inversions by providing sharp boundaries and architectural control on the location of the resistive and conductive regions in the CSEM model. The seismic horizons include: 1) the base of the gas hydrate stability zone (BGHSZ), 2) the top of salt, and 3) the top and bottom of a fine grained marine mud interval with near vertical hydrate filled fractures, to constrain the CSEM inversion model. The top of salt provides improved location for brines, water saturated salt, and resistive salt. Inversions of the CSEM data map the occurrence of a ‘halo’ of conductive brines above salt. The use of the BGHSZ as a constraint on the inversion helps distinguish between free gas and gas hydrate as well as gas hydrate and water saturated sediments.

Louisiana↗

Magmatic-vapor expansion and the formation of high-sulfidation gold deposits: Structural controls on hydrothermal alteration and ore mineralization

High-sulfidation copper–gold lode deposits such as Chinkuashih, Taiwan, Lepanto, Philippines, and Goldfield, Nevada, formed within 1500 m of the paleosurface in volcanic terranes. All underwent an early stage of extensive advanced argillic silica–alunite alteration followed by an abrupt change to spatially much more restricted stages of fracture-controlled sulfide–sulfosalt mineral assemblages and gold–silver mineralization. The alteration as well as ore mineralization stages of these deposits were controlled by the dynamics and history of syn-hydrothermal faulting. At the Sulfate Stage, aggressive advanced argillic alteration and silicification were consequent on the in situ formation of acidic condensate from magmatic vapor as it expanded through secondary fracture networks alongside active faults. The reduction of permeability at this stage due to alteration decreased fluid flow to the surface, and progressively developed a barrier between magmatic-vapor expansion constrained by the active faults and peripheral hydrothermal activity dominated by hot-water flow. In conjunction with the increased rock strength resulting from alteration, subsequent fault-slip inversion in response to an increase in compressional stress generated new, highly permeable fractures localized by the embrittled, altered rock. The new fractures focused magmatic-vapor expansion with much lower heat loss so that condensation occurred. Sulfide Stage sulfosalt, sulfide, and gold–silver deposition then resulted from destabilization of vapor phase metal species due to vapor decompression through the new fracture array. The switch from sulfate to sulfide assemblages is, therefore, a logical consequence of changes in structural permeability due to the coupling of alteration and fracture dynamics rather than to changes in the chemistry of the fluid phase at its magmatic source.

Ore Geology Reviews↗

Structural stability of coprecipitated natural organic matter and ferric iron under reducing conditions

The objective was to assess the interaction of Fe coprecipitated with dissolved organic matter (DOM) and its effect on Fe (hydr)oxide crystallinity and DOM retention under abiotic reducing conditions. A Fe-based coagulant was reacted with DOM from an agricultural drain and the resulting precipitate (floc) was exposed to S(-II) and Fe(II). Solution concentrations of Fe(II/III) and DOM were monitored, floc crystallinity was determined using X-ray diffraction, and the composition and distribution of functional groups were assessed using scanning transmission X-ray microscopy (STXM) and near edge X-ray absorption fine structure (NEXAFS) spectroscopy. Results indicate coprecipitation of Fe(III) with DOM forms a non-crystalline floc that withstands crystallization regardless of change in pH, Fe:DOM ratio and type of reductant added. There was no evidence that exposure to reducing conditions led to release of DOM from the floc, indicating that coprecipitation with complex natural DOM in aquatic environments may stabilize Fe (hydr)oxides against crystallization upon reaction with reduced species and lead to long term sequestration of the DOM. STXM analysis identified spatially distinct regions with remarkable functional group purity, contrary to the model of DOM as a relatively uniform complex polymer lacking identifiable organic compounds. Polysaccharide-like OM was strongly and directly correlated with the presence of Fe but showed different Fe binding strength depending on the presence of carboxylic acid functional groups, whereas amide and aromatic functional groups were inversely correlated with Fe content.

Organic Geochemistry↗