Search USGS⌕ Search

SEARCH · Search USGS

Results for “Evolution”

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 955 records · Page 53Linked to original sources
USGS2022

Variation in within-host replication kinetics among virus genotypes provides evidence of specialist and generalist infection strategies across three salmonid host species

Theory of the evolution of pathogen specialization suggests that a specialist pathogen gains high fitness in one host, but this comes with fitness loss in other hosts. By contrast, a generalist pathogen does not achieve high fitness in any host, but gains ecological fitness by exploiting different hosts, and has higher fitness than specialists in nonspecialized hosts. As a result, specialist pathogens are predicted to have greater variation in fitness across hosts, and generalists would have lower fitness variation across hosts. We test these hypotheses by measuring pathogen replicative fitness as within-host viral loads from the onset of infection to the beginning of virus clearance, using the rhabdovirus infectious hematopoietic necrosis virus (IHNV) in salmonid fish. Based on field prevalence and virulence studies, the IHNV subgroups UP, MD, and L are specialists, causing infection and mortality in sockeye salmon, steelhead, and Chinook salmon juveniles, respectively. The UC subgroup evolved naturally from a UP ancestor and is a generalist infecting all three host species but without causing severe disease. We show that the specialist subgroups had the highest peak and mean viral loads in the hosts in which they are specialized, and they had low viral loads in nonspecialized hosts, resulting in large variation in viral load across hosts. Viral kinetics show that the mechanisms of specialization involve the ability to both maximize early virus replication and avoid clearance at later times, with different mechanisms of specialization evident in different host–virus combinations. Additional nuances in the data included different fitness levels for nonspecialist interactions, reflecting different trade-offs for specialist viruses in other hosts. The generalist UC subgroup reached intermediate viral loads in all hosts and showed the smallest variation in fitness across hosts. The evolution of the UC generalist from an ancestral UP sockeye specialist was associated with fitness increases in steelhead and Chinook salmon, but only slight decreases in fitness in sockeye salmon, consistent with low- or no-cost generalism. Our results support major elements of the specialist–generalist theory, providing evidence of a specialist–generalist continuum in a vertebrate pathogen. These results also quantify within-host replicative fitness trade-offs resulting from the natural evolution of specialist and generalist virus lineages in multi-host ecosystems

Virus Evolution↗
USGS2012

Basins in ARC-continental collisions

Arc-continent collisions occur commonly in the plate-tectonic cycle and result in rapidly formed and rapidly collapsing orogens, often spanning just 5-15 My. Growth of continental masses through arc-continent collision is widely thought to be a major process governing the structural and geochemical evolution of the continental crust over geologic time. Collisions of intra-oceanic arcs with passive continental margins (a situation in which the arc, on the upper plate, faces the continent) involve a substantially different geometry than collisions of intra-oceanic arcs with active continental margins (a situation requiring more than one convergence zone and in which the arc, on the lower plate, backs into the continent), with variable preservation potential for basins in each case. Substantial differences also occur between trench and forearc evolution in tectonically erosive versus tectonically accreting margins, both before and after collision. We examine the evolution of trenches, trench-slope basins, forearc basins, intra-arc basins, and backarc basins during arc-continent collision. The preservation potential of trench-slope basins is low; in collision they are rapidly uplifted and eroded, and at erosive margins they are progressively destroyed by subduction erosion. Post-collisional preservation of trench sediment and trench-slope basins is biased toward margins that were tectonically accreting for a substantial length of time before collision. Forearc basins in erosive margins are usually floored by strong lithosphere and may survive collision with a passive margin, sometimes continuing sedimentation throughout collision and orogeny. The low flexural rigidity of intra-arc basins makes them deep and, if preserved, potentially long records of arc and collisional tectonism. Backarc basins, in contrast, are typically subducted and their sediment either lost or preserved only as fragments in melange sequences. A substantial proportion of the sediment derived from collisional orogenesis ends up in the foreland basin that forms as a result of collision, and may be preserved largely undeformed. Compared to continent-continent collisional foreland basins, arc-continent collisional foreland basins are short-lived and may undergo partial inversion after collision as a new, active continental margin forms outboard of the collision zone and the orogen whose load forms the basin collapses in extension.

Book chapter↗
USGS2017

The effect of lithology on valley width, terrace distribution, and coarse sediment provenance in a tectonically stable catchment with flat-lying stratigraphy

How rock resistance or erodibility affects fluvial landforms and processes is an outstanding question in geomorphology that has recently garnered attention owing to the recognition that the erosion rates of bedrock channels largely set the pace of landscape evolution. In this work, we evaluate valley width, terrace distribution, and sediment provenance in terms of reach scale variation in lithology in the study reach and discuss the implications for landscape evolution in a catchment with relatively flat2 lying stratigraphy and very little uplift. A reach of the 21 Buffalo National River in Arkansas was partitioned into lithologic reaches and the mechanical and chemical resistance of the main lithologies making up the catchment was measured. Valley width and the spatial distribution of terraces were compared among the different lithologic reaches. The surface grain size and provenance of coarse (2-90 mm) sediment of both modern gravel bars and older terrace deposits were measured and defined. The results demonstrate a strong impact of lithology upon valley width, terrace distribution, and coarse sediment provenance and therefore, upon landscape evolution processes. Channel down-cutting through different lithologies creates variable patterns of resistance across catchments and continents. Particularly in post-tectonic and nontectonic landscapes, the variation in resistance that arises from the exhumation of different rocks in channel longitudinal profiles can impact local base levels, initiating responses that can be propagated through channel networks. The rate at which that response is transmitted through channels is potentially amplified and/or mitigated by differences between the resistance of channel beds and sediment loads. In the study 36 reach, variation in lithologic resistance influences the prevalence of lateral and vertical 37 processes, thus producing a spatial pattern of terraces that reflects rock type rather than 38 climate, regional base level change, or hydrologic variability.

Arkansas, Missouri↗
USGS2019

Advances in computational morphodynamics using the International River Interface Cooperative (iRIC) software

Results from computational morphodynamics modeling of coupled flow-bed-sediment systems are described for ten applications as a review of recent advances in the field. Each of these applications is drawn from solvers included in the public-domain International River Interface Cooperative (iRIC) software package. For mesoscale river features such as bars, predictions of alternate and higher mode river bars are shown for flows with equilibrium sediment supply and for a single case of oversupplied sediment. For microscale bed features such as bedforms, computational results are shown for the development and evolution of two-dimensional bedforms using a simple closure-based two-dimensional model, for two- and three-dimensional ripples and dunes using a three-dimensional large-eddy simulation flow model coupled to a physics-based particle transport model, and for the development of bed streaks using a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver with a simple sediment-transport treatment. Finally, macroscale or channel evolution treatments are used to examine the temporal development of meandering channels, a failure model for cantilevered banks, the effect of bank vegetation on channel width, the development of channel networks in tidal systems, and the evolution of bedrock channels. In all examples, computational morphodynamics results from iRIC solvers are shown to compare well to observations of natural bed morphology. For each of the three scales investigated here, brief suggestions for future work and potential research directions are offered.

Earth Surface Processes and Landforms↗
USGS2021

Geometry of obstacle marks at instream boulders-Integration of laboratory investigations and field observations

Obstacle marks are instream bedforms, typically composed of an upstream frontal scour hole and a downstream sediment accumulation in the vicinity of an obstacle. Local scouring at infrastructure (e.g. bridge piers) is a well‐studied phenomenon in hydraulic engineering, while less attention is given to the time‐dependent evolution of frontal scour holes at instream boulders and their geometric relations (depth to width, and length ratio). Furthermore, a comparison between laboratory studies and field observations is rare. Therefore, the morphodynamic importance of such scour features to fluvial sediment transport and morphological change is largely unknown. In this study, obstacle marks at boulder‐like obstructions were physically modelled in 30 unscaled process‐focused flume experiments (runtime per experiment ≥ 5760 min) at a range of flows (subcritical, clear‐water conditions, emergent and submerged water levels) and boundary conditions designed to represent the field setting (i.e. obstacle tilting, and limited thickness of the alluvial layer). Additionally, geometries of scour holes at 90 in‐situ boulders (diameter ≥ 1 m) located in a 50‐km segment of the Colorado River in Marble Canyon (AZ) were measured from a 1 m‐resolution digital elevation model. Flume experiments reveal similar evolution of local scouring, irrespective of hydraulic conditions, controlled by the scour incision, whereas the thickness of the alluvial layer and obstacle tilting into the evolving frontal scour hole limit incision. Three temporal evolution phases—(1) rapid incision, (2) decreasing incision, and (3) scour widening—are identified based on statistical analysis of spatiotemporal bed elevation time series. A quantitative model is presented that mechanistically predicts enlargement in local scour length and width based on (1) scour depth, (2) the inclination of scour slopes, and (3) the planform area of the frontal scour hole bottom. The comparison of field observations and laboratory results demonstrates scale invariance of geometry, which implies similitude of processes and form rather than equifinality.

Earth Surface Processes and Landforms↗
USGS2000

Europa's Crust and Ocean: Origin, Composition, and the Prospects for Life

We have considered a wide array of scenarios for Europa's chemical evolution in an attempt to explain the presence of ice and hydrated materials on its surface and to understand the physical and chemical nature of any ocean that may lie below. We postulate that, following formation of the jovian system, the europan evolutionary sequence has as its major links: (a) initial carbonaceous chondrite rock, (b) global primordial aqueous differentiation and formation of an impure primordial hydrous crust, (c) brine evolution and intracrustal differentiation, (d) degassing of Europa's mantle and gas venting, (e) hydrothermal processes, and (f) chemical surface alteration. Our models were developed in the context of constraints provided by Galileo imaging, near infrared reflectance spectroscopy, and gravity and magnetometer data. Low-temperature aqueous differentiation from a carbonaceous CI or CM chondrite precursor, without further chemical processing, would result in a crust/ocean enriched in magnesium sulfate and sodium sulfate, consistent with Galileo spectroscopy. Within the bounds of this simple model, a wide range of possible layered structures may result; the final state depends on the details of intracrustal differentiation. Devolatilization of the rocky mantle and hydrothermal brine reactions could have produced very different ocean/crust compositions, e.g., an ocean/crust of sodium carbonate or sulfuric acid, or a crust containing abundant clathrate hydrates. Realistic chemical-physical evolution scenarios differ greatly in detailed predictions, but they generally call for a highly impure and chemically layered crust. Some of these models could lead also to lateral chemical heterogeneities by diapiric upwellings and/or cryovolcanism. We describe some plausible geological consequences of the physical-chemical structures predicted from these scenarios. These predicted consequences and observed aspects of Europa's geology may serve as a basis for further analys is and discrimination among several alternative scenarios. Most chemical pathways could support viable ecosystems based on analogy with the metabolic and physiological versatility of terrestrial microorganisms. ?? 2000 Academic Press.

Icarus↗
USGS2001

Major, trace element and isotope geochemistry (Sr-Nd-Pb) of interplinian magmas from Mt. Somma-Vesuvius (Southern Italy)

Major, trace element and isotopic (Sr, Nd, Pb) data are reported for representative samples of interplinian (Protohistoric, Ancient Historic and Medieval Formations) activity of Mt. Somma-Vesuvius volcano during the last 3500 years. Tephra and lavas exhibit significant major, trace element and isotopic variations. Integration of these data with those obtained by previous studies on the older Somma suites and on the latest activity, allows to better trace a complete petrological and geochemical evolution of the Mt. Somma-Vesuvius magmatism. Three main groups of rocks are recognized. A first group is older than 12.000 yrs, and includes effusive-explosive activity of Mt. Somma. The second group (8000-2700 yrs B.P.) includes the products emitted by the Ottaviano (8000 yrs. B.P.) and Avellino (3550 yrs B.P.) plinian eruptions and the interplinian activity associated with the Protohistoric Formation. Ancient Historic Formation (79-472 A.D.), Medieval Formation (472-1139 A.D.) and Recent interplinian activity (1631-1944 A.D.) belong to the third group of activity (79-1944 A.D.). The three groups of rocks display distinct positive trends of alkalis vs. silica, which become increasingly steeper with age. In the first group there is an increase in silica and alkalis with time, whereas an opposite tendency is observed in the two younger groups. Systematic variations are also evident among the incompatible (Pb, Zr, Hf, Ta, Th, U, Nb, Rb, Cs, Ba) and compatible elements (Sr, Co, Cr). REE document variable degrees of fractionation, with recent activity displaying higher La/Yb ratios than Medieval and Ancient Historic products with the same degree of evolution. N-MORB normalized multi-element diagrams for interplinian rocks show enrichment in Rb, Th, Nb, Zr and Sm (> *10 N-MORB). Sr isotope ratios are variable, with Protohistoric rocks displaying 87Sr/86Sr= 0.70711-0.70810, Ancient Historic 87Sr/86Sr=0.70665-0.70729, and Medieval 87Sr/86Sr=0.70685-0.70803. Neodymium isotopic compositions in the interplinian rocks show a tendency to become slightly more radiogenic with age, from the Protohistoric (143Nd/144Nd=0.51240-0.51247) to Ancient Historic (143Nd/144Nd=0.51245-0.51251). Medieval interplinian activity (143Nd/144Nd: 0.51250-0.51241) lacks meaningful internal trends. All the interplinian rocks have virtually homogeneous compositions of 207Pb/204Pb and 208Pb/204Pb in acid-leached residues (207Pb/204Pb ???15.633 to 15.687, 208Pb/204Pb ???38.947 to 39.181). Values of 206Pb/204Pb are very distinctive, however, and discriminate among the three interplinian cycles of activity (Protohistoric: 18.929-18.971, Ancient Historic: 19.018-19.088, Medieval: 18.964-19.053). Compositional trends of major, trace element and isotopic compositions clearly demonstrate strong temporal variations of the magma types feeding the Somma-Vesuvius activity. These different trends are unlikely to be related only to low pressure evolutionary processes, and reveal variations of parental melt composition. Geochemical data suggest a three component mixing scheme for the interplinian activity. These involve HIMU-type and DMM-type mantle and Calabrian-type lower crust. Interaction between these components has taken place in the source; however, additional quantitative constraints must be acquired in order to better discriminate between magma characteristics inherited from the sources and those acquired during shallow level evolution.

Mineralogy and Petrology↗
USGS1995

Beach-steps: An evolutionary perspective

Field observation of contrasting beach-step behavior at Canaveral National Seashore on two subsequent days when incident-wave conditions in the inner surf zone were similar prompted this re-examination of our conceptual and quantitative understanding of beach steps. These lower-foreshore features are more complex than previously assumed, evolving through erosional as well as accretional phases, and displaying equifinality in geometric form but not necessarily internal sedimentary structure. Past and recent evidence is reviewed that links beach steps to incident waves at the surging-plunging transition and to the action of a backwash vortex. Tides and low-frequency waves likely play no direct role in beach-step initiation, although their presence can have pronounced influences on modulating nearshore hydrodynamics, and thus, on beach-step maintenance and evolution. A generalized, conceptual model capturing these aspects of beach-step dynamics is presented. Beach-step initiation proceeds via step “carving”, “excavation”, or “building” depending on the erosional-accretional character of the beach-foreshore system. Subsequent evolution of the step form may take one of several alternative morphodynamic pathways including stepface “retreat”, step “drag down”, or step “infilling/elimination” depending on tidal stage/range or wave set-up/setdown. Additional data on equilibrium beach-step forms and associated morphodynamic and hydrodynamic conditions in the field are necessary before quantitative models of beach-step existence and evolution can be formulated with realistic results.

Marine Geology↗
USGS2022

Maturation study of vitrinite in carbonaceous shales and coals: Insights from hydrous pyrolysis

The presence of vitrinite in sedimentary rocks of post-Silurian age allows its reflectance to be used to estimate the thermal maturation of organic matter in petroleum systems. Increasing reflectance of vitrinite, which is primarily driven by aromaticity, depends primarily on the time and temperature attributes of its evolutionary pathway. This study evaluated carbonaceous shales proximal to coal measures and coal samples via isothermal hydrous pyrolysis (HP) to compare differences in the maturation pathways of vitrinite. Sample residues were analysed via vitrinite reflectance (VRo), geochemical screening tests (organic carbon and programmed temperature pyrolysis), and infrared spectroscopy. The study included samples from Indian and North American basins, to observe differences in vitrinite evolution with respect to enclosing mineral matrix, starting degree of aromaticity, organic matter types, stratigraphic age, and depositional environment. The organic content of HP residues shows an intuitive response to the thermal stress of HP, e.g., a general depletion of total organic carbon (TOC) content, pyrolyzate (S2), and hydrogen index with increasing HP temperature. Infrared proxies including C-factor and CH2/CH3 generally decrease with increasing thermal maturity indicating loss of O via CO2 generation and the thermal cracking of aliphatic organic matter. Tmax, production index (PI), and VRo show intuitive increasing values with respect to HP temperature. The least mature sample (0.48 ± 0.05% VRo) generally experienced the maximum change in these parameters during maturation, whereas the most mature sample (0.99 ± 0.06% VRo) generally showed the least change. This observation is consistent with higher kinetic barriers to reaction in more aromatic vitrinite which contains higher bond dissociation energies. Devolatilization of vitrinite during HP causes formation of gas evacuation vacuoles and contraction cracks in the vitrinite grains of both coal and carbonaceous shale. Similarities in vitrinite response to HP between coal and carbonaceous shale suggest that thermal evolution of the vitrinite maceral is principally controlled by inherent rate-limiting kinetic parameters related to its molecular structure. Whereas, the stratigraphic age, sedimentary environment, surrounding organic matter, lithology, and catalysis by mineral composition have less effect. To further improve our understanding of vitrinite aromatization and kinetic parameters, future studies of vitrinite reflectance thermal evolution with temperature should include coal and carbonaceous shale from the same stratigraphic section and extant woody tissue from modern vascular plants.

International Journal of Coal Geology↗
USGS2019

Morphodynamic modelling of the wilderness breach, Fire Island, New York. Part I: Model set-up and validation

On October 29, 2012, storm surge and large waves produced by Hurricane 13 Sandy resulted in the formation of a breach in eastern Fire Island, NY. The goals of this study 14 are to gain a better understanding of the physical processes that govern breach behavior and 15 to assess whether process-based models can be used to forecast the evolution of future 16 breaches. The Wilderness Breach grew rapidly in size during the first winter following 17 formation. Growth of the breach was accompanied by the formation of a complex of flood 18 shoals inside Great South Bay, a primary channel that flowed through the eastern part of the 19 flood shoals, and an ebb shoal on the ocean side of the breach. From the summer of 2013 20 through late 2015, the breach continued to change and evolve, albeit at a much slower pace 21 than in the first year after formation. A hybrid combination of Delft3D and XBeach models is 22 used to hindcast the morphodynamic evolution of the Wilderness Breach over the first three 23 years after formation. The formation of the breach during Hurricane Sandy is not part of the 24 simulations. Model simulations are initiated with a post-storm topography in which the 25 breach is already present. The models are capable of hindcasting the main morphodynamic 26 changes of the Wilderness Breach. The spatial patterns, as well as the bulk statistics, such as 27 2 breach geometry and sediment volume changes, are reasonably 28 well reproduced by the model. 29 The model sheds light on previously unknown processes of breach evolution, especially 30 regarding sediment transport and flow regimes within the breach complex.

New York↗
USGS2022

Trace elements in olivine fingerprint the source of 2018 magmas and shed light on explosive-effusive eruption cycles at Kīlauea Volcano

Understanding magma genesis and the evolution of intensive parameters (temperature, pressure, composition, degree of melting) in the mantle source of highly active volcanic systems is crucial for interpreting magma supply changes over time and recognizing cyclic behavior to anticipate future volcanic behavior. Major and trace elements in olivine are commonly used to study variations in mantle lithologies and melting conditions (e.g., temperature, pressure, oxygen fugacity) affecting the mantle over time. Here, we track the temporal evolution of primary melts through the most recent cycle of explosive and effusive eruptions at Kīlauea (Hawai‘i), which spans the last ∼500 years. We report major and trace elements in olivine from the last explosive period (∼1500 – early 1820’s Keanakāko‘i Tephra) and the most recent decade of the current effusive period (2018 LERZ, 2015–2018 Pu‘u‘ō‘ō, 2008–2018 lava lake and 2020 eruption in Halema‘uma‘u). Scandium concentrations in olivine allow characterizing changes in mantle source between 1500 and 2018, and suggest that the recent (2015–2018) magma feeding the Pu‘u‘ō‘ō cone did not significantly interact with the magma that erupted in the LERZ in 2018. The evolution of olivine and melt compositions over the past 500 years is not easily reconcilable with variations in mantle potential temperature, pressure of mantle melt pooling and storage, or oxygen fugacity. Instead, Sc, Mn, and Co concentrations and Ni/Mg ratio in high forsterite (Fo >87) olivine advocate for an increase in the proportion of clinopyroxene in the mantle source associated with a slightly higher degree of partial melting from 1500 to 2018. Changes in primitive melt compositions and degrees of mantle melting may well modulate magma supply to the crust and formation-replenishment of steady or ephemeral summit reservoirs, and thereby control transitions between explosive and effusive periods at Kīlauea. Analyzing trace elements in olivine at Kīlauea and elsewhere could therefore provide important clues on subtle changes occurring at the mantle level that might herald changes in volcanic behavior.

Hawaii↗
USGS2026

Lunar analog study using portable gamma-ray neutron detector: Radiochemical mapping of silicic-to-basaltic volcanic terrains in the San Francisco Volcanic Field, Arizona

Studying lunar silicic volcanism provides key insights into the Moon’s crustal evolution, magmatic processes, and volcanic history. The Lunar-VISE (Lunar Vulkan Imaging and Spectroscopy Explorer) mission will investigate the Gruithuisen domes, a unique lunar region hypothesized to have formed through silicic volcanism. Using instruments on a Firefly Aerospace lander and a Honeybee Robotics rover, Lunar-VISE will analyze mineralogy, geochemistry, and surface properties to determine the origin and evolution of the domes, with a gamma-ray and neutron spectrometer (LV-GRNS) among its payload instruments. In preparation for this mission, we conducted preliminary fieldwork using a handheld gamma-ray neutron detector with NaI(Tl) and Cs LiYCl :Ce (CLYC) scintillators. This study focuses on various rhyolitic and basaltic volcanic centers in the San Francisco Volcanic Field (SFVF) near Flagstaff, Arizona—specifically Sugarloaf Peak, Bonito Lava Flow, and Robinson Mountain, a region containing several well-characterized lunar analog sites. The SFVF was selected for this study due to its broad compositional diversity spanning basaltic to rhyolitic compositions, its well-preserved volcanic morphologies, and its extensive use in previous NASA field campaigns and astronaut training exercises, making it an ideal terrestrial laboratory for testing planetary exploration techniques. We measured natural radioactivity, specifically from potassium (K), thorium (Th), and uranium (U) and other elements in their decay chains, which serve as diagnostic tracers of magmatic differentiation and crustal evolution processes, to assess geochemical variability across compositionally diverse terrains. The detector’s sensitivity was assessed across varying concentration levels. Regionally averaged concentrations of radioisotopes were determined by gamma-ray spectroscopy at selected sites. Our measurements reveal significant compositional variations between sites, with Sugarloaf Peak (rhyolitic, silicic dome) exhibiting the highest average radioisotope concentrations (K: 3.29 ± 0.24 wt%, U: 14.35 ± 1.81 ppm, Th: 27.14 ± 2.43 ppm), while Bonito Lava Flow (K: 1.13 ± 0.08 wt%, U: 3.86 ± 0.52 ppm, Th: 7.52 ± 1.02 ppm), and Robinson Mountain (K: 1.46 ± 0.16 wt%, U: 5.59 ± 1.06 ppm, Th: 11.75 ± 1.98 ppm) show lower concentrations aligned with basaltic compositions. Gamma-ray fluxes were elevated by a factor of approximately three to four at Sugarloaf Peak relative to nearby basaltic terrains, consistent with expected geochemical differentiation patterns. Furthermore, at meter scales, proximity to geological features significantly affects measurements. Th concentrations adjacent to a cliff face at Sugarloaf’s base were 28% higher than values measured 3-5 meters away from the same feature, demonstrating localized compositional heterogeneity. Analysis of station-by-station measurements reveals that K, U, and Th concentrations show an elevation trend at Robinson Mountain and generally higher concentrations near Sugarloaf Peak’s summit, suggesting progressive magmatic differentiation and/or the presence of more evolved lithologies at higher elevations. By mapping these elements on Earth using GRNS, we aim to optimize measurement techniques for Lunar-VISE and similar rover-borne missions by demonstrating the applicability of portable GRNS instruments for planetary surface exploration in an analog environment.

Arizona↗
USGS2007

Generation and propagation of nonlinear internal waves in Massachusetts Bay

During the summer, nonlinear internal waves (NLIWs) are commonly observed propagating in Massachusetts Bay. The topography of the area is unique in the sense that the generation area (over Stellwagen Bank) is only 25 km away from the shoaling area, and thus it represents an excellent natural laboratory to study the life cycle of NLIWs. To assist in the interpretation of the data collected during the 1998 Massachusetts Bay Internal Wave Experiment (MBIWE98), a fully nonlinear and nonhydrostatic model covering the generation/shoaling region was developed, to investigate the response of the system to the range of background and driving conditions observed. Simplified models were also used to elucidate the role of nonlinearity and dispersion in shaping the NLIW field. This paper concentrates on the generation process and the subsequent evolution in the basin. The model was found to reproduce well the range of propagation characteristics observed (arrival time, propagation speed, amplitude), and provided a coherent framework to interpret the observations. Comparison with a fully nonlinear hydrostatic model shows that during the generation and initial evolution of the waves as they move away from Stellwagen Bank, dispersive effects play a negligible role. Thus the problem can be well understood considering the geometry of the characteristics along which the Riemann invariants of the hydrostatic problem propagate. Dispersion plays a role only during the evolution of the undular bore in the middle of Stellwagen Basin. The consequences for modeling NLIWs within hydrostatic models are briefly discussed at the end.

Massachusetts↗
USGS2008

Seismicity and deformation induced by magma accumulation at three basaltic volcanoes

We analyzed the evolution of volcano‐tectonic (VT) seismicity and deformation at three basaltic volcanoes (Kilauea, Mauna Loa, Piton de la Fournaise) during phases of magma accumulation. We observed that the VT earthquake activity displays an accelerating evolution at the three studied volcanoes during the time of magma accumulation. At the same times, deformation rates recorded at the summit of Kilauea and Mauna Loa volcanoes were not accelerating but rather tend to decay. To interpret these observations, we propose a physical model describing the evolution of pressure produced by the accumulation of magma into a reservoir. This variation of pressure is then used to force a simple model of damage, where damage episodes are equivalent to earthquakes. This model leads to an exponential increase of the VT activity and to an exponential decay of the deformation rate during accumulation phases. Seismicity and deformation data are well fitted by such an exponential model. The time constant, deduced from the exponential increase of the seismicity, is in agreement with the time constant predicted by the model of magma accumulation. This VT activity can thus be a direct indication of the accumulation of magma at depth, and therefore can be seen as a long‐term precursory phenomenon, at least for the three studied basaltic volcanoes. Unfortunately, it does not allow the prediction of the onset of future eruptions, as no diverging point (i.e., critical time) is present in the model.

Journal of Geophysical Research B: Solid Earth↗
USGS2019

Fault initiation in serpentinite

Serpentinite fault rheology is fundamental to tectonic and earthquake processes, yet links between deformation textures and strength evolution during fault initiation are poorly constrained. Here I present field and petrographic microstructural observations of unsheared and sheared serpentinite that demonstrate a progression of fault development. I compliment observations with a clast size distribution analysis to investigate the evolution of fault rigidity, and a numerical model to query the stress distribution of a clast-in-matrix geometry. Unsheared microstructures reveal well-aligned, elongate serpentine in the matrix and short, randomly-oriented serpentine in clasts. Sheared matrix displays cataclastic textures, discrete brittle surfaces, dissolution bands and ductile textures defined by anastomosing networks of well-aligned, fine-grained serpentine. During fault initiation, matrix serpentine anisotropy promotes slip on basal planes or fiber aggregates, and clast-on-clast interactions drive a high bulk viscosity prone to stick-slip behavior. As deformation progresses clast fracturing is focused at clast tips and smaller clasts are preferentially removed by dissolution-precipitation processes, increasing the relative abundance of matrix. Strain is continually focused in the matrix and as clast content reduces the bulk viscosity drops. This study reveals that viscosity contrasts formed by primary serpentinization textures are a major driver for the development and strength evolution of faults. On a continental-scale, similar processes may govern earthquake distributions, fault growth and segmentation patterns on young, serpentinite-hosted faults.

California↗
USGS2019

Azimuthal seismic anisotropy of 70 Ma Pacific‐plate upper mantle

Plate formation and evolution processes are predicted to generate upper mantle seismic anisotropy and negative vertical velocity gradients in oceanic lithosphere. However, predictions for upper mantle seismic velocity structure do not fully agree with the results of seismic experiments. The strength of anisotropy observed in the upper mantle varies widely. Further, many refraction studies observe a fast direction of anisotropy rotated several degrees with respect to the paleospreading direction, suggesting that upper mantle anisotropy records processes other than 2D corner flow and plate‐driven shear near mid‐ocean ridges. We measure 6.0 ± 0.3% anisotropy at the Moho in 70 Ma lithosphere in the central Pacific with a fast direction parallel to paleospreading, consistent with mineral alignment by 2D mantle flow near a mid‐ocean ridge. We also find an increase in the strength of anisotropy with depth, with vertical velocity gradients estimated at 0.02 km/s/km in the fast direction and 0 km/s/km in the slow direction. The increase in anisotropy with depth can be explained by mechanisms for producing anisotropy other than intrinsic effects from mineral fabric, such as aligned cracks or other structures. This measurement of seismic anisotropy and gradients reflects the effects of both plate formation and evolution processes on seismic velocity structure in mature oceanic lithosphere, and can serve as a reference for future studies to investigate the processes involved in lithospheric formation and evolution.

Journal of Geophysical Research B: Solid Earth↗
USGS2019

Insights into pāhoehoe lava emplacement using visible and thermal structure-from-motion photogrammetry

We present the evolution over 3 months of a 2016–2017 pāhoehoe flow at Kīlauea as it changed from a narrow sheet flow into a compound lava field fed by a stable system of tubes. The portion of the flow located on Kīlauea's coastal plain was characterized using helicopter‐based visible and thermal structure‐from‐motion photogrammetry to construct a series of georeferenced digital surface models and thermal maps on eight different days. Results reveal key influences on the emplacement and evolution of such long‐lived pāhoehoe flows. This region of the flow grew by ~12 × 10 6 m 3 with a near‐constant time‐average discharge rate of 1.2–2.7 m 3 /s. The development of two tube systems is captured and shows an initial nascent tube enhanced by a narrow topographic confinement, which later inflated and created a topographic inversion that modulated the emplacement of a second flow lobe with its own tube system. The analysis of breakouts at various stages of the field's life suggests that the evolution of the thermal and morphological properties of the flow surface reflect its maturity. Thermal properties of breakouts were used to expand the empirical relationship of breakout cooling to longer timescales. This study contributes to the long‐term development and validation of more accurate predictive models for pāhoehoe, required during the management of long‐lasting lava flow crises in Hawai'i and elsewhere.

Hawaii↗
USGS2021

Fluid-earthquake and earthquake-earthquake interactions in southern Kansas, USA

An increase in injection activity associated with energy production in southern Kansas starting in 2013 has been linked to the occurrence of more than 130,000 earthquakes (M −1.5 to 4.9) between 2014 and 2017. Studies suggest that the dramatic increase in seismicity rate is related to wastewater injection into the highly permeable Arbuckle formation. Most of the seismicity is located in the underlying crystalline basement, for which hydrological properties and specific fault geometries are unknown. Additionally, some earthquake clusters occurred relatively far (tens of kilometers) from the main injection wells. Therefore, the effect of pore pressure diffusion may be insufficient to explain the relation between the volume of injected fluids and the spatiotemporal evolution of seismicity. Combining physical models (static stress and poroelasticity) and a statistical cluster analysis applied to a high-resolution relocated catalog, we analyze the evolution of seismicity in southern Kansas. We find that pore pressure changes (Δ p ) and Coulomb stress changes (ΔCFS) due to fluid diffusion smaller than 0.1 MPa are enough to initiate seismic sequences, which then evolve depending on their distance from the major injection wells. However, we find that earthquake sequences have different seismogenic responses to Δ p and ΔCFS in terms of triggering threshold. In regions located close to disposal wells (Harper area) our cluster analysis suggests that both earthquake interactions and fluid diffusion control the evolution of seismicity. On the other hand, at greater distances (Milan area), where clustering behavior suggests greater earthquake interactions, we find that coseismic ΔCFS are larger than Δ p .

Kansas, Oklahoma↗