SEARCH · Geology Reports

Results for “Earth Interactions”

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 451 records · Page 25Linked to original sources

USGS2013

Large wood budget and transport dynamics on a large river using radio telemetry

Despite the abundance of large wood (LW) river studies there is still a lack of understanding of LW transport dynamics on large low gradient rivers. This study used 290 radio frequency identification tagged (RFID) LW and 54 metal (aluminum) tagged LW, to quantify the percent of in-channel LW that moves per year and what variables play a role in LW transport dynamics. Aluminum tags were installed and monitored on LW in-transit during the rising limb of a flood, the mean distance traveled by those pieces during the week was 13.3 river kilometers (km) with a maximum distance of 72 km. RFID tagged LW moved a mean of 11.9 km/yr with a maximum observed at 101.1 km/yr. Approximately 41% of LW low on the bank moves per year. The high rate of transport and distance traveled is likely due to the lack of interaction between LW floating in the channel and the channel boundaries, caused primarily by the width of the channel relative to length of the LW. Approximately 80% of the RFID tags moved past a fixed reader during the highest 20% of river stage per year. LW transport and logjam dynamics are complicated at high flows as pieces form temporary jams that continually expand and contract. Unlike most other studies, key members that create a logjam were defined more by stability than jam size or channel/hydrologic conditions. Finally, using an existing geomorphic database for the river, and data from this study, we were able to develop a comprehensive LW budget showing that 5% of the in-channel LW population turns over each year (input from mass wasting and fluvial erosion equals burial, decomposition, and export out of system) and another 16% of the population moving within the system.

North Carolina;Virginia↗

USGS2022

Observation-constrained multicycle dynamic models of the southern San Andreas and the northern San Jacinto Faults: Addressing complexity in paleoearthquake extent and recurrence with realistic 2D fault geometry

Understanding mechanical conditions that lead to complexity in earthquakes is important to seismic hazard analysis. In this study, we simulate physics-based multicycle dynamic models of the San Andreas fault (Carrizo through San Bernardino sections) and the San Jacinto fault (Claremont and Clark strands). We focus on a complex fault geometry based on the Southern California Earthquake Center Community Fault Model and its effect over multiple earthquake cycles. Using geodetically derived strain rates, we validate the models against geologic slip rates and recurrence intervals at various paleoseismic sites. We find that the interactions among fault geometry, dynamic rupture and interseismic stress accumulation produce stress heterogeneities, leading to rupture segmentation and variability in earthquake recurrence. Our models produce earthquakes with rupture extents similar to a recent comprehensive paleoseismic catalog. The “earthquake gates” of the Big Bend and the Cajon Pass occasionally impede dynamic ruptures. The angle of compression, which is the subtraction of the maximum shear strain rate direction from the local fault strike, can better determine the likelihood of the impedance of restraining bends to dynamic ruptures. Because the Big Bend has an angle of compression of ∼20°, ruptures that traverse the Big Bend, like the 1857 Fort Tejon earthquake, are more frequent than expected based on empirical relations which predict the ∼40° restraining bend to terminate most ruptures. Our models indicate that large ruptures tend to initiate north of the Big Bend and propagate southwards, similar to the 1857 earthquake, providing critical information for ground shaking assessment in the region.

California↗

USGS2022

Multi-decadal simulation of marsh topography evolution under sea level rise and episodic sediment loads

Coastal marsh within Mediterranean climate zones is exposed to episodic watershed runoff and sediment loads that occur during storm events. Simulating future marsh accretion under sea level rise calls for attention to: (a) physical processes acting over the time scale of storm events and (b) biophysical processes acting over time scales longer than storm events. Using the upper Newport Bay in Southern California as a case study, we examine the influence of event-scale processes on simulated change in marsh topography by comparing: (a) a biophysical model that integrates with an annual time step and neglects event-scale processes (BP-Annual), (b) a physical model that resolves event-scale processes but neglects biophysical interactions (P-Event), and (c) a biophysical model that resolves event-scale physical processes and biophysical processes at annual and longer time scales (BP-Event). A calibrated BP-Event model shows that large (>20-year return period) episodic storm events are major drivers of marsh accretion, depositing up to 30 cm of sediment in one event. Greater deposition is predicted near fluvial sources and tidal channels and less on marshes further from fluvial sources and tidal channels. In contrast, the BP-Annual model poorly resolves spatial structure in marsh accretion as a consequence of neglecting event-scale processes. Furthermore, the P-Event model significantly overestimates marsh accretion as a consequence of neglecting marsh surface compaction driven by annual scale biophysical processes. Differences between BP-Event and BP-Annual models translate up to 20 cm per century in marsh surface elevation.

California↗

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↗

USGS2021

The making of the NEAM Tsunami Hazard Model 2018 (NEAMTHM18)

The NEAM Tsunami Hazard Model 2018 (NEAMTHM18) is a probabilistic hazard model for tsunamis generated by earthquakes. It covers the coastlines of the North-eastern Atlantic, the Mediterranean, and connected seas (NEAM). NEAMTHM18 was designed as a three-phase project. The first two phases were dedicated to the model development and hazard calculations, following a formalized decision-making process based on a multiple-expert protocol. The third phase was dedicated to documentation and dissemination. The hazard assessment workflow was structured in Steps and Levels. There are four Steps: Step-1) probabilistic earthquake model; Step-2) tsunami generation and modeling in deep water; Step-3) shoaling and inundation; Step-4) hazard aggregation and uncertainty quantification. Each Step includes a different number of Levels. Level-0 always describes the input data; the other Levels describe the intermediate results needed to proceed from one Step to another. Alternative datasets and models were considered in the implementation. The epistemic hazard uncertainty was quantified through an ensemble modeling technique accounting for alternative models’ weights and yielding a distribution of hazard curves represented by the mean and various percentiles. Hazard curves were calculated at 2,343 Points of Interest (POI) distributed at an average spacing of ∼20 km. Precalculated probability maps for five maximum inundation heights (MIH) and hazard intensity maps for five average return periods (ARP) were produced from hazard curves. In the entire NEAM Region, MIHs of several meters are rare but not impossible. Considering a 2% probability of exceedance in 50 years (ARP≈2,475 years), the POIs with MIH >5 m are fewer than 1% and are all in the Mediterranean on Libya, Egypt, Cyprus, and Greece coasts. In the North-East Atlantic, POIs with MIH >3 m are on the coasts of Mauritania and Gulf of Cadiz. Overall, 30% of the POIs have MIH >1 m. NEAMTHM18 results and documentation are available through the TSUMAPS-NEAM project website (http://www.tsumaps-neam.eu/), featuring an interactive web mapper. Although the NEAMTHM18 cannot substitute in-depth analyses at local scales, it represents the first action to start local and more detailed hazard and risk assessments and contributes to designing evacuation maps for tsunami early warning.

Frontiers in Earth Science↗

USGS2025

A time-space model of graphite mineral systems

Increasing demand for graphite in energy storage systems warrants review of graphite ore genesis in a mineral systems framework. Orogenic graphite encompasses the metamorphic and orogenic mineral systems that produce flake graphite and hydrothermal vein (lump and chip) graphite deposits, respectively. A common feature of orogenic graphite deposits is an association with upper amphibolite- to granulite-facies metasedimentary rocks in continent-continent or continent-island arc collisional orogens. Orogenic flake graphite deposits form primarily through graphitization of organic carbon during regional metamorphism, but strain localization and partial melting of pelitic protoliths are likely important processes for graphite grade and quality enrichment. Orogenic vein graphite deposits precipitate from hydrothermal fluids, possibly derived from metamorphism or anatexis at depth. Decarbonation reactions in mixed calcareous-carbonaceous metasedimentary sequences are the most likely carbon sources for the veins. In contrast, intrusion-related graphite includes magmatic-hydrothermal and metamorphic mineral systems that form primarily in continental arc settings via the interaction of magmas with carbonaceous sedimentary packages. Magmatic-hydrothermal flake graphite deposits are hosted in plutonic and volcanic rocks, and result from the exsolution of CO 2 -CH 4 -rich fluids from melts contaminated by such packages. Contact metamorphism of carbonaceous sedimentary rocks by plutons produces some microcrystalline (amorphous) graphite deposits, including many in China. Compilation of geologic data from known graphite deposits globally suggests that pulses of carbon deposition in the Paleoproterozoic, Mesoproterozoic, and Neoproterozoic provided source carbon material. Subsequent supercontinent orogenesis at ca. 2,100 to 1,700 Ma (Columbia), ca. 1,300 to 1,000 Ma (Rodinia), and ca. 650 to 500 Ma (Gondwana) resulted in the genesis of orogenic flake and vein graphite deposits, where favorable geologic components overlapped with organic ± carbonate carbon-rich strata. Cryogenian deposition of graphite protoliths and Cryogenian – Cambrian metamorphic mineralization account for nearly 75% of all known resources globally and coincide with profound carbon isotope excursions and climate variability, implying a link with the global carbon budget. Comparatively few graphite deposits are associated with Pangea-forming orogens, attributed to less exhumation and/or denudation. High-temperature metasedimentary belts containing organic carbon-rich protoliths are most favorable for hosting orogenic flake graphite deposits, whereas sequences that also contain carbonate protoliths are favorable for orogenic graphite veins. Continent-scale orogenic belts may host both deposit types along with vanadium deposits. Use of a time-space mineral systems framework for graphite deposits can improve exploration models needed to ensure future supply of this critical mineral and provide insights into Earth’s long-term carbon cycle.

Mineralium Deposita↗

USGS2024

Predicting characteristic length scales of barrier island segmentation in microtidal environments

Segmented barrier islands can be found in regions with small tidal ranges. In contrast to tidally dominated barriers, where inlet dynamics are thought to control island length scales, the controls on barrier island length scales in wave-dominated environments have not been quantified. These microtidal barriers typically have a curved shoreline, suggesting the influence of wave-driven alongshore sediment transport. Microtidal barriers are also typically hydrodynamically isolated from one another, as weak tidal flows limit interactions between adjoining barriers. To better understand the controls on and scales of barrier segmentation in the relative absence of tides, here we develop a theoretical framework to estimate the alongshore length scales at which a barrier will either breach or heal following a disturbance in the barrier morphology. The non-dimensional framework compares the timescales of overwash (advective) and alongshore sediment transport (diffusive) processes along barrier island chains. We then apply this framework to modern barrier islands in the microtidal Gulf of Mexico using wave hindcast data and the lengths, widths, heights, and lagoon depths measured from remotely sensed geospatial data and topobathymetric data. We find that most of these barriers are currently longer than their critical length scale, often as a result of coastal restoration efforts. Our critical length scale analysis suggests that most of the Gulf of Mexico barriers are vulnerable to segmentation despite coastal restoration efforts intended to protect fisheries and the mainland coasts.

Florida, Louisiana, Mississippi↗

USGS2009

Changing sources of strontium to soils and ecosystems across the Hawaiian Islands

Strontium isotope ratios assist ecosystem scientists in constraining the sources of alkaline earth elements, but their interpretation can be difficult because of complexities in mineral weathering and in the geographical and environmental controls on elemental additions and losses. Hawaii is a "natural laboratory" where a number of important biogeochemical variables have either limited ranges or vary in systematic ways, providing a unique opportunity to understand the impact of time, climate, and atmospheric inputs on the evolution of base cation sources to ecosystems. There are three major sources of strontium (Sr) to these ecosystems, each with distinct isotopic compositions: basalt lava, Asian dust, and rainfall. We present Sr isotope and concentration data on both bulk soil digests and NH 4 Ac extracts from soil profiles covering a wide range of environments and substrate ages. Bulk soil material from dry climates and/or young substrate ages with > 80????g g - 1 Sr retain basalt-like Sr isotopic signatures, whereas those with Sr concentrations < 80????g g - 1 can have isotope signatures that range from basalt-like values to the more radiogenic values associated with continental dust. Although both dust accumulation and lava weathering are time- and rainfall-dependent, the overall concentration of Sr drops with increasing leaching even as quartz and mica derived from continental dust sources increase to > 40% by mass. At elevated dust levels, lava-derived Sr is low and dust-derived Sr is the dominant control of 87 Sr/ 86 Sr in bulk soils; however, 87 Sr/ 86 Sr of NH 4 Ac-extractable Sr largely reflects atmospheric deposition of marine aerosol in these situations. Overall, whole-soil Sr isotope values are controlled by complex interactions between Sr provided by lava weathering but partially lost by leaching, and Sr provided by dust but held in more resistant minerals. The isotopic composition of NH 4 Ac-extractable Sr and of the biota is controlled by lava weathering and rainfall contribution of Sr with only minor contributions from radiogenic dust sources. ?? 2009 Elsevier B.V.

Chemical Geology↗

USGS2022

PHREEQ-N-AMDTreat+REYs water-quality modeling tools to evaluate acid mine drainage treatment strategies for recovery of rare-earth elements

The PHREEQ-N-AMDTreat+REYs water-quality modeling tools have the fundamental capability to simulate aqueous chemical reactions and predict the formation of metal-rich solids during the treatment of acid mine drainage (AMD). These new user-friendly, publicly available tools were expanded from the PHREEQ-N-AMDTreat tools to include the precipitation of rare-earth elements plus yttrium (REYs) and the adsorption of REYs onto hydrous Fe, Al, and Mn oxides. The tool set consists of a caustic titration model that indicates equilibrium surface and aqueous speciation of REYs as functions of pH and caustic agent, and a kinetics+adsorption model that simulates progressive changes in pH, major ions, and REYs in water and solids during sequential steps through passive and/or active treatment. Each model has a user interface (UI) that facilitates the input of water-quality data and adjustment to geochemical or treatment system variables; for example, retention time and aeration rate are adjustable parameters in the kinetics model. On-screen graphs display results of changes in metals and associated solute concentrations as functions of pH or retention time; details are summarized in output tables. A goal of such modeling is to identify strategies that could produce a concentrated REYs extract from AMD or mine waste leachate. For example, if REYs could be concentrated after first removing substantial Fe and Al, the final REYs-bearing phase(s) could be more efficiently processed for REYs recovery and, therefore, may represent a more valuable commodity. Preliminary modeling supports the hypothesis that Fe and Al can be removed at pH < 5.5 using conventional sequential oxidation and neutralization treatment processes without removing REYs, and that further increasing pH can promote the adsorption of REYs by hydrous Mn oxides. Alternatively, chemicals such as oxalate or phosphate may be added to precipitate REYs compounds following initial steps to decrease Fe and Al concentrations. The aqueous geochemical model framework is comprehensive and permits evaluation of effects from interactive chemical and physical variables. Field studies that demonstrate REYs attenuation from AMD and corresponding solid-phase formation during specific treatment steps plus laboratory studies of aqueous/solid interactions are helpful to corroborate, refine, and constrain modelin parameters.

Conference Paper↗

USGS1976

Tectonic implications of space-time patterns of Cenozoic magmatism in the western United States

Locations of 2,100 radiometrically dated igneous rocks were plotted on a series of 20 maps, each representing an interval within the period 80 m.y. B.P. to present. Derivative maps showing the distributions in space and time of dated granitic intrusive rocks, silicic lavas and domes, ash-flow tuffs, andesitic-dacitic rocks, and basalts depict well the two main petrogenetic assemblages noted previously by others: (1) mainly intermediate andesitic-dacitic suites, including associated granitic intrusive rocks, silicic extrusive rocks, and minor basaltic lavas, are interpreted as reflecting plate interactions related to subduction along the continental margin; and (2) bimodal suites, dominantly basaltic but with minor silicic extrusive rocks, are interpreted as reflecting extensional tectonics. Space-time distribution of the two assemblages suggests that magmatic arcs extended continously parallel to the continental margin from Canada to Mexico in latest Mesozoic and in Oligocene times. An early Cenozoic null in magmatism in the Great Basin may delineate the region where subduction was arrested temporarily by development of the proto-San Andreas fault as a transform in coastal California or, alternatively, may reflect complex subsurface configurations of subducted plates. The late Cenozoic transition from subduction-related magmatism to extention-related basaltic volcanism in the southern Cordillera occurred at different times in different areas in harmony with current concepts about the migration of the Mendocino triple junction as the modern San Andreas transform fault was formed. The plots also reveal the existence of several discrete magmatic loci where igneous activity of various kinds was characteristically more intense and long-lived than elsewhere.

Earth and Planetary Science Letters↗

USGS2019

Introduction to special issue on gas hydrate in porous media: Linking laboratory and field‐scale phenomena

The proliferation of drilling expeditions focused on characterizing natural gas hydrate as a potential energy resource has spawned widespread interest in gas hydrate reservoir properties and associated porous media phenomena. Between 2017 and 2019, a Special Section of this journal compiled contributed papers elucidating interactions between gas hydrate and sediment based on laboratory, numerical modeling, and field studies. Motivated mostly by field observations in the northern Gulf of Mexico and offshore Japan, several papers focus on the mechanisms for gas hydrate formation and accumulation, particularly with vapor phase gas, not dissolved gas, as the precursor to hydrate. These studies rely on numerical modeling or laboratory experiments using sediment packs or benchtop micromodels. A second focus of the Special Section is the role of fines in inhibiting production of gas from methane hydrate, controlling the distribution of hydrate at a pore scale, and influencing the bulk behavior of seafloor sediments. Other papers fill knowledge gaps related to the physical properties of hydrate-bearing sediments and advance new approaches in coupled thermal-mechanical modeling of these sediments during hydrate dissociation. Finally, one study addresses the long-standing question about the fate of methane hydrate at the molecular level when CO2 is injected into natural reservoirs under hydrate-forming conditions.

Journal of Geophysical Research B: Solid Earth↗

USGS2022

Stress heterogeneity as a driver of aseismic slip during the 2011 Prague, Oklahoma aftershock sequence

The interaction of aseismic and seismic slip before and after an earthquake is fundamental for both earthquake nucleation and postseismic stress relaxation. However, it can be difficult to determine where and when aseismic slip occurs within the seismogenic zone because geodetic techniques are limited to detecting moderate to large slip amplitudes or long duration small slip amplitudes. Here, we use repeating earthquakes (earthquakes that re-rupture the same fault patch) as a proxy for aseismic slip during the 2011 Prague, Oklahoma earthquake sequence. We find that aseismic slip in the Prague earthquake sequence occurs both within the granitic basement and the overlying sedimentary rocks. The repeating earthquakes show that patches of aseismic slip are mostly located at fault intersections. These fault intersections hosted possible mainshock slip, abundant aftershocks, and afterslip. We estimate that ∼40% of the aftershocks are driven by afterslip. We interpret that aseismic slip occurs at fault intersections where stress heterogeneity creates patches of lower stress that are stable within a nonsteady state, rate-state framework.

Oklahoma↗

USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Journal of Geophysical Research B: Solid Earth↗

USGS2020

USGS “Did You Feel It?” — Science and lessons from twenty years of citizen science-based macroseismology

The U.S. Geological Survey (USGS) “Did You Feel It?” (DYFI) system is an automatic method for rapidly collecting macroseismic intensity data from Internet users’ shaking and damage reports and for generating intensity maps immediately following felt earthquakes. DYFI has been in operation for nearly two decades (1999-2019) in the United States, and for nearly 15 years globally. During that period the amount of data collected is astounding: Over 5 million individual DYFI intensity reports—spanning all magnitude and distance ranges—have been amassed and archived. DYFI allows for macroseismic data collection at rates and quantities never before imagined, and thus high-quality macroseismic intensity maps can be made almost immediately, and with more complete coverage at higher resolution than in the past. DYFI also allows for valuable positive interactions of the citizenry with a Federal science agency. In essence, the widespread adoption of DYFI—along with ShakeMap—has facilitated the general acceptance of the very concept of shaking intensity, fundamentally improving our agency’s ability to communicate both hazard and risk to the population. DYFI effectively confirms the importance of reporting and inculcating the public’s understanding of intensity—in addition to magnitude—for a proper perspective of earthquake risk-related decision-making. Furthermore, the vast amount of DYFI data allow for data-rich analyses of otherwise intractable seismological, sociological, and earthquake impact studies that range from quantifying the shaking due to induced earthquakes, to human response and risk perception, to relating recorded shaking metrics to macroseismic effects, and the attenuation of intensity with magnitude and distance. Naturally, web-based data collection also poses challenges. After two decades of experience acquiring data with the DYFI system, we address some of these hallenges by documenting refinements to our algorithmic and operational procedures that have evolved over that time. Lastly, we outline new opportune research and development directions for our DYFI approach to citizen seismology.

Frontiers in Earth Science↗

USGS2026

Diverse novel and avian-associated viruses in the ileal viromes of northern mockingbird (Mimus polyglottos)

Viruses are the most abundant and diverse organisms on Earth, though only a small portion cause disease. Understanding viral diversity is key to understanding and predicting pathogen emergence and zoonotic spillover. Here, we use meta-transcriptomic sequencing to examine the viral communities in the ileum of 25 Northern Mockingbirds ( Mimus polyglottos ) from various locations across Texas. We assembled high-quality genomes of 43 viral species (40 species identified to 13 families, one to kingdom, and two to realm), 38 of which were novel. They tentatively represent avian- (n = 3), arthropod- (n = 21), plant- (n = 5) and fungi- (n = 4) associated, or other (n = 10) viruses. The arthropod-associated Dicistroviridae family was the most dominant, comprising known and potentially new species. Of potential epidemiological importance were three novel and avian-associated viruses: members of the families Hepeviridae and Picornaviridae , and a new Matryoshka RNA virus. The Matryoshka RNA virus 8 (MaRNAV-8) is sister to other Matryoshka RNA viruses, and its co-occurrence with haemosporida further supports the nested virus-parasite-vector-vertebrate host relationship of this group of viruses, with potential implications for parasite evolution, fitness and load and vector competence. The Picornaviridae virus is a member of an avian hepatovirus clade, found nested within a clade containing both the mammalian pathogens Hepatovirus A – I and the avian Tremovirus pathogens, suggestive of a newly discovered pathogen of Northern Mockingbird. Although the recovered Hepeviridae virus is of unknown pathology, its family members include the Hepatitis E viruses. With the great diversity and novelty described from ileal viromes, discriminating potential pathogens and commensal microbiota from viruses associated with food items remains challenging. A deeper understanding of virus transmission and the risk of potential zoonosis can be enhanced by tracking viruses through the food web and via inter-specific and predator-prey interactions, particular in areas subject to land-use change, where human-wildlife interactions are increased and the risks from emerging pathogens of veterinary and medical importance are more pronounced.

Texas↗

USGS2012

Lead and strontium isotopes and related trace elements as genetic tracers in the Upper Cenozoic rhyolite-basalt association of the Yellowstone Plateau volcanic field

Supported by various field geologic and petrologic data, the contents of Pb, U, Th, Rb, and Sr and the isotopic compositions of Pb and Sr for upper Cenozoic volcanic rocks of the Yellowstone Plateau volcanic field are consistent with the hypothesis of derivation of the basaltic and rhyolitic magmas by partial melting of distinct source regions in the upper mantle and lower crust, respectively. All the basalt samples analyzed but one have systematically lower values of 207 Pb/ 204 Pb and 87 Sr/ 86 Sr than the rhyolites. The values of 206 Pb/ 204 Pb are smaller, and 87 Sr/ 86 Sr are mostly larger than known values in oceanic basalts. In all but one case, the values of 207 Pb/ 204 Pb are higher than expected from an extrapolation of known values in oceanic basalts to less radiogenic values of 206 Pb/ 204 Pb. Because there are no xenoliths, phenocrysts are only moderate to sparse in abundance, REE patterns are low and flat at the radiogenic end of lead isotopic compositions, several values of Rb/Sr are low, and 80% of the basalt samples form a well-developed secondary isochron separate from the rhyolites, we favor an interpretation for basalt genesis wherein isotopic signatures of most mafic magmas were attained in a continental ‘keel’ of mantlelike character about 2.6 b.y. old or somewhat older attached to the crust, and these signatures were unaltered by magma passage through the crust. At the very least, the current data continue to cast serious doubt as to the inevitability of crustal contamination for basaltic magma intruding the continental environment and postulate that much can be learned about the mantle under continents through the study of continental basalts. One basalt unit with an unusually low value of 207 Pb/ 204 Pb and an 87 Ar/ 86 Ar less than 0.704 may represent subcontinental ‘keel’-derived magma that rose unaltered to the surface. Our data also are not consistent with formation of this rhyolite-basalt association primarily by such processes as crystal fractionation, separation of immiscible silicate liquids from a common parental magma, or fractional melting of a homogeneous source. Rather as a conceptual model, we envision large mafic intrusions to have been injected into the lower crust resulting in rhyolite generation through partial anatexis of the adjacent wall rocks which probably had a 206 Pb/ 204 Pb < 17 and 87 Sr/ 86 Sr > 0.709; a model that has much in common with that proposed by Holmes (1931). All the other hypotheses listed have the necessary added complication that either the basalt or the rhyolite or both become contaminated after the two magma types separated, have problems accounting for the lack of igneous rocks of intermediate compositions or production of such large volumes of rhyolitic material (∼5000 km 3 ), and fail to explain why rhyolitic magma is not a more common occurrence in the ocean basin. We appeal to bouyancy of rhyolites to generate a barrier for basalt magma migration and account for the great preponderance of rhyolite relative to basalt at the surface. Furthermore, the complex isotopic picture in the rhyolites indicates that many of these magmas interacted with the upper crustal geologic units that they traversed. The interactions involved diverse processes, probably including reacton with hydrothermal fluids or hydrothermally altered rocks at high levels as well as by contamination with Phanerozoic sedimentary and Precambrian crystalline rocks at deeper levels. At the very least, we feel our study adds a cautionary note to the currently increasingly popular hypothesis that differentiation of basalt or gabbro magmas to rhyolite or granite (as distinct from tonalite or dacite) is a common occurrence and is therefore an important continential building process. Models for formation of rhyolite and granite predominantly by reworking of crust (anatexis) must still be considered. The primitive Archean mantle of the region was characterized by higher Rb/Sr, U/Pb, and Th/U values than are typical of modern suboceanic mantle. The mantle residuum within the continental subcrustal lithosperic ‘keel’ that resulted from the Archean crustal differentiation event probably was depleted in Rb/Sr and U/Pb, and the crust was correspondingly enriched in these ratios. The crust probably was further differentiated by an Archean high-grade metamorphism, during or after the primary event, into a granulitic lower crust depleted in U/Pb and Rb/Sr and a lower-grade upper crust enriched in these ratios.

Journal of Geophysical Research Solid Earth↗

USGS2021

Warming and microbial uptake influence the fate of added soil carbon across a Hawai'ian weathering gradient

Tropical forest soils contain some of the largest carbon (C) stocks on Earth, yet the effects of warming on the fate of fresh C entering tropical soils are still poorly understood. This research sought to understand how the fate of fresh C entering soils is influenced by warming, soil weathering status, and C chemistry. We hypothesized that compounds that are quickly incorporated into microbial biomass (i.e., greater C use efficiency [CUE]) subsequently have longer-term (255 days) retention in soil. We also hypothesized that relatively weathered soils with greater sorptive capacity also retain more fresh C in the short and longer-terms, and that C in these soils is more resistant to weathering loss compared with less weathered soils. We tested these hypotheses by adding two 13 C-labeled compounds (glucose and glycine) to three tropical forest soils from a weathering gradient in Hawai'i, and then incubating soils at ambient (16 °C), +5 °C, and +10 °C for 255 days. We found that 255-day 13 C retention in mineral soil across sites and temperatures was best predicted by two factors: initial retention of 13 C in mineral soil and initial microbial 13 CUE (Adjusted R 2 = 0.78). Carbon compound type influenced 13 C initial retention, with greater glucose- 13 C retention versus glycine- 13 C retention in mineral soils and microbial biomass, corresponding to greater glucose- 13 C retention in soil at 255 days. Warming had a negative longer-term effect on the retention of 13 C only in the least-weathered soil, supporting our hypothesis. These results show that initial retention of fresh C in soils via mineral sorption and microbial uptake is a strong predictor of longer-term retention, indicating that immediate C losses are a major hurdle for soil C storage. Also, retention of fresh C appears most sensitive to warming in less-weathered tropical soils, supporting the idea that mineral sorption may provide some protections against warming. Understanding the interaction between soil sorptive properties and warming for C cycling could improve predictions of forest-climate feedbacks for tropical regions.

Soil Biology and Biochemistry↗

USGS2016

Joint analysis of geodetic and earthquake fault-plane solution data to constrain magmatic sources: A case study from Kīlauea Volcano

A joint analysis of geodetic and seismic datasets from Kīlauea Volcano during a period of magmatic unrest in 2006 demonstrates the effectiveness of this combination for testing and constraining models of magma dynamics for a complex, multi-source system. At the end of 2003, Kīlauea's summit began a four-year-long period of inflation due to a surge in magma supply to the volcano. In 2006, for the first time since 1982, Kīlauea's Southwest Rift Zone (SWRZ) also experienced inflation. To investigate the characteristics of active magma sources and the nature of their interactions with faults in the SWRZ during 2006, we integrate, through Coulomb stress modeling, contemporary geodetic data from InSAR and GPS with a new catalogue of double-couple fault-plane solutions for volcano-tectonic earthquakes. We define two periods of inflation during 2006 based on the rate of deformation measured in daily GPS data, spanning February to 15 March 2006 (Period 1) and 16 March to 30 September 2006 (Period 2). InSAR data for these two periods are inverted to determine the position, change in size, and shape of inflation sources in each period. Our new models are consistent with microseismic activity from each period. They suggest that, during Period 1, deformation in the SWRZ can be explained by pressurization of magma in a spherical reservoir beneath the south caldera, and that, during Period 2, magma was also aseismically intruded farther to the southwest into the SWRZ along a sub-horizontal plane. Our Coulomb stress analysis shows that the microseismicity recorded in the SWRZ is induced by overpressurization of the south caldera reservoir, and not by magma intrusion into the SWRZ. This study highlights the importance of a joint analysis of independent geophysical datasets to fully constrain the nature of magma accumulation.

Hawaii↗