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.

SourceResearch connectionFrom yearThrough year

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

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

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↗

USGS2021

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↗

USGS2022

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↗

USGS2024

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↗

USGS2011

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↗

USGS2012

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↗

USGS2015

From patterns to causal understanding: Structural equation modeling (SEM) in soil ecology

In this perspectives paper we highlight a heretofore underused statistical method in soil ecological research, structural equation modeling (SEM). SEM is commonly used in the general ecological literature to develop causal understanding from observational data, but has been more slowly adopted by soil ecologists. We provide some basic information on the many advantages and possibilities associated with using SEM and provide some examples of how SEM can be used by soil ecologists to shift focus from describing patterns to developing causal understanding and inspiring new types of experimental tests. SEM is a promising tool to aid the growth of soil ecology as a discipline, particularly by supporting research that is increasingly hypothesis-driven and interdisciplinary, thus shining light into the black box of interactions belowground.

Pedobiologia↗

USGS2024

Utilization of stochastic ground motion simulations for scenario-based performance assessment of geo-structures

Probabilistic seismic performance assessments of engineered structures can be highly sensitive to the seismic input excitation and its variability. In the present study, the scenario-based performance assessment recommended by Federal Emergency Management Agency (FEMA) P-58 guidelines is adopted to estimate seismic fragility of concrete dams for various seismic hazard scenarios. Due to the scarcity of recorded ground motions and thereby their poor representation of uncertainties, stochastic ground motion simulation methods are utilized to obtain the required input excitations. Moreover, to understand the uncertainty in ground motion simulation models, two broadband stochastic simulation models are used to generate input excitations representing six seismic hazard scenarios defined by earthquake magnitude, source-to-site distance, and soil conditions. Optimal intensity measure parameters for each scenario are identified using a systematic procedure that considers criteria such as efficiency, practicality, proficiency, sufficiency, and hazard compatibility. Fragility curves and surfaces are derived using the cloud analysis technique, taking into account various damage measures and limit state functions. The study finds that the derived fragility curves are particularly sensitive to the selection of earthquake scenarios, the choice of records, and the methods used to calculate fragility curves, with less sensitivity observed to different engineering demand parameters. Given this sensitivity, particularly to ground motion selection, the study highlights the necessity of incorporating both model-to-model variability (epistemic uncertainty) and record-to-record variability (aleatory uncertainty), alongside the established material and modeling uncertainties, in the probabilistic seismic assessment.

Reliability Engineering and System Safety (RESS)↗

USGS2023

Historical Structure from Motion (HSfM): Automated processing of historical aerial photographs for long-term topographic change analysis

Precisely measuring the Earth’s changing surface on decadal to centennial time scales is critical for many science and engineering applications, yet long-term records of quantitative landscape change are often temporally and geographically sparse. Archives of scanned historical aerial photographs provide an opportunity to augment these records with accurate elevation measurements that capture the historical state of the Earth surface. Structure from Motion (SfM) photogrammetry workflows produce high-quality digital elevation models (DEMs) and orthoimage mosaics from these historical images, but time-intensive tasks like manual image preprocessing (e.g., fiducial marker identification) and ground control point (GCP) selection impede processing at scale. We developed an automated method to process historical images and generate self-consistent time series of high-resolution (0.5–2 m) DEMs and orthomosaics, without manual GCP selection. The method relies on SfM to correct camera interior and exterior orientation and a robust multi-stage co-registration approach using modern reference terrain datasets for geolocation refinement. We demonstrate the method using scanned images from the North American Glacier Aerial Photography (NAGAP) archive collected between 1967 and 1997. We present results for two sites with variable photo acquisition geometry and overlap — Mount Baker and South Cascade Glacier in Washington State, USA. The automated method corrects initial camera position errors of several kilometers and produces accurately georeferenced, high-resolution DEMs and orthoimages, regardless of camera configuration, acquisition geometry, terrain characteristics, and reference DEM properties. The average RMS reprojection error following bundle adjustment optimization was 0.67 px (0.15 m) for the 261 images contributing to 10 final DEM mosaics between 1970 and 1992 at Mount Baker, and 0.65 px (0.13 m) for the 243 images contributing to 18 individual DEMs between 1967 and 1997 at South Cascade Glacier. The relative accuracy of elevation values in the historical time series stacks was 0.68 m at Mount Baker and 0.37 m at South Cascade Glacier. Our products have reduced systematic error and improved accuracy compared to DEM products generated using SfM with manual GCP selection. Final elevation change measurement precision was ∼ "> ∼ 0.7–1.0 m over a 30-year period, enabling the study of processes with rates as low as ∼ "> ∼ 1-3 cm/yr. Our results demonstrate the potential of this scalable method to rapidly process archives of historical imagery and deliver new quantitative insights on long-term geodetic change and Earth surface processes.

Washington↗

USGS2020

Spatial variability of phytoplankton in a shallow tidal freshwater system reveals complex controls on abundance and community structure

Estuaries worldwide are undergoing changes to patterns of aquatic productivity because of human activities that alter flow, impact sediment delivery and thus the light field, and contribute nutrients and contaminants like pesticides and metals. These changes can influence phytoplankton communities, which in turn can alter estuarine food webs. We used multiple approaches-including high-resolution water quality mapping, synoptic sampling, productivity and nitrogen uptake rates, Lagrangian parcel tracking, enclosure experiments and bottle incubations-over a short time period to take a “spatial snapshot” of conditions in the northern region of the San Francisco Estuary (California, USA) to examine how environmental drivers like light availability, nutrients, water residence time, and contaminants affect phytoplankton abundance and community attributes like size distribution, taxonomic structure, and nutrient uptake rates. Zones characterized by longer residence time (15–60 days) had higher chlorophyll-a concentrations (9 ± 4 µg L −1 ) and were comprised primarily of small phytoplankton cells (<5 µm, 74 ± 8%), lower ammonium concentrations (1 ± 0.8 µM), higher nitrate uptake rates, and higher rates of potential carbon productivity. Conversely, zones characterized by shorter residence time (1–14 days) had higher ammonium concentration (13 ± 5 µM) and lower chlorophyll-a concentration (5 ± 1 µg L −1 ) with diatoms making up a larger percent contribution. Longer residence time, however, did not result in the accumulation of large (>5 µm) cells considered important to pelagic food webs. Rather, longer residence time zones had a phytoplankton community comprised primarily of small cells, particularly picocyanobacteria that made up 38 ± 17% of the chlorophyll-a – nearly double the concentration seen in shorter residence time zones (22 ± 7% picocyanobacterial of chlorophyll-a). Our results suggest that water residence time in estuaries may have an effect as large or larger than that experimentally demonstrated for light, contaminants, or nutrients.

California↗

USGS2006

Crustal structure of the northeastern margin of the Tibetan plateau from the Songpan-Ganzi terrane to the Ordos basin

The 1000-km-long Darlag–Lanzhou–Jingbian seismic refraction profile is located in the NE margin of the Tibetan plateau. This profile crosses the northern Songpan-Ganzi terrane, the Qinling-Qilian fold system, the Haiyuan arcuate tectonic region, and the stable Ordos basin. The P-wave and S-wave velocity structure and Poisson's ratios reveal many significant characteristics in the profile. The crustal thickness increases from northeast to southwest. The average crustal thickness observed increases from 42 km in the Ordos basin to 63 km in the Songpan-Ganzi terrane. The crust becomes obviously thicker south of the Haiyuan fault and beneath the West-Qinlin Shan. The crustal velocities have significant variations along the profile. The average P-wave velocities for the crystalline crust vary between 6.3 and 6.4 km/s. Beneath the Songpan-Ganzi terrane, West-Qinling Shan, and Haiyuan arcuate tectonic region P-wave velocities of 6.3 km/s are 0.15 km/s lower than the worldwide average of 6.45 km/s. North of the Kunlun fault, with exclusion of the Haiyuan arcuate tectonic region, the average P-wave velocity is 6.4 km/s and only 0.5 km/s lower than the worldwide average. A combination of the P-wave velocity and Poisson's ratio suggests that the crust is dominantly felsic in composition with an intermediate composition at the base. A mafic lower crust is absent in the NE margin of the Tibetan plateau from the Songpan-Ganzi terrane to the Ordos basin. There are low velocity zones in the West-Qinling Shan and the Haiyuan arcuate tectonic region. The low velocity zones have low S-wave velocities and high Poisson's ratios, so it is possible these zones are due to partial melting. The crust is divided into two layers, the upper and the lower crust, with crustal thickening mainly in the lower crust as the NE Tibetan plateau is approached. The results in the study show that the thickness of the lower crust increases from 22 to 38 km as the crustal thickness increases from 42 km in the Ordos basin to 63 km in the Songpan-Ganzi terrane south of the Kunlun fault. Both the Conrad discontinuity and Moho in the West-Qinling Shan and in the Haiyuan arcuate tectonic region are laminated interfaces, implying intense tectonic activity. The arcuate faults and large earthquakes in the Haiyuan arcuate tectonic region are the result of interaction between the Tibetan plateau and the Sino–Korean and Gobi Ala Shan platforms.

Tectonophysics↗

USGS2009

Crustal structure across the Three Gorges area of the Yangtze platform, central China, from seismic refraction/wide-angle reflection data

We present active-source seismic data recorded along a 300 km-long profile across the Three Gorges area of the western Yangtze platform, central China. From west to east, the profile crosses the Zigui basin, Huangling dome and Jianghan basin. The derived crustal P -wave velocity structure changes significantly across the Tongchenghe fault that lies at the transition from the Huangling dome to the Jianghan basin. West of the Tongchenghe fault, beneath the Zigui basin and the Huangling dome, we observe a ~ 42 km thick crust of relatively low average velocity (6.3–6.4 km/s). In contrast, east of the Tongchenghe fault, beneath the Jianghan basin, the crust is only 30 km thick and has a high average velocity (6.6–6.7 km/s). A west–east variation in crustal composition along the Tongchenghe fault is also inferred. West of the fault, P -wave velocities suggest a felsic composition with an intermediate layer at the base of the crust, whilst, east of the fault, felsic, intermediate, and mafic crustal layers are apparent. Our results suggest that the crust beneath the Jianghan basin has been thinned by rifting, accompanied by intrusion of the lower crust by mafic dikes and sills. The west-to-east division of the crust in the Three Gorges area coincides with first-order geophysical contrasts in gravity, topography, crustal and lithospheric thickness.

Tectonophysics↗