SEARCH · USGS Science

Results for “Evolution”

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

USGS2002

The Fish Canyon magma body, San Juan volcanic field, Colorado: Rejuvenation and eruption of an upper-crustal batholith

More than 5000 km 3 of nearly compositionally homogeneous crystal-rich dacite (∼68 wt % SiO 2 : ∼45% Pl + Kfs + Qtz + Hbl + Bt + Spn + Mag + Ilm + Ap + Zrn + Po) erupted from the Fish Canyon magma body during three phases: (1) the pre-caldera Pagosa Peak Dacite (an unusual poorly fragmented pyroclastic deposit, ∼200 km 3 ); (2) the syn-collapse Fish Canyon Tuff (one of the largest known ignimbrites, ∼5000 km 3 ); (3) the post-collapse Nutras Creek Dacite (a volumetrically minor lava). The late evolution of the Fish Canyon magma is characterized by rejuvenation of a near-solidus upper-crustal intrusive body (mainly crystal mush) of batholithic dimensions. The necessary thermal input was supplied by a shallow intrusion of more mafic magma represented at the surface by sparse andesitic enclaves in late-erupted Fish Canyon Tuff and by the post-caldera Huerto Andesite. The solidified margins of this intrusion are represented by holocrystalline xenoliths with Fish Canyon mineralogy and mineral chemistry and widely dispersed partially remelted polymineralic aggregates, but dehydration melting was not an important mechanism in the rejuvenation of the Fish Canyon magma. Underlying mafic magma may have evolved H 2 O–F–S–Cl-rich fluids that fluxed melting in the overlying crystal mush. Manifestations of the late up-temperature magma evolution are: (1) resorbed quartz, as well as feldspars displaying a wide spectrum of textures indicative of both resorption and growth, including Rapakivi textures and reverse growth zoning (An 27–28 to An 32–33 ) at the margins of many plagioclase phenocrysts; (2) high Sr, Ba, and Eu contents in the high-SiO 2 rhyolite matrix glass, which are inconsistent with extreme fractional crystallization of feldspar; (3) oscillatory and reverse growth zoning toward the margins of many euhedral hornblende phenocrysts (rimward increases from ∼5·5–6 to 7·7–8·5 wt % Al 2 O 3 ). Homogeneity in magma composition at the chamber-wide scale, contrasting with extreme textural and chemical complexities at the centimeter–millimeter scale, is consistent with a dynamic environment, wherein crystals with a variety of growth and resorption histories were juxtaposed shortly before eruption by convective currents.

Colorado↗

USGS2016

Probing the volcanic-plutonic connection and the genesis of crystal-rich rhyolite in a deeply dissected supervolcano in the Nevada Great Basin: Source of the late Eocene Caetano Tuff

Late Cenozoic faulting and large-magnitude extension in the Great Basin of the western USA has created locally deep windows into the upper crust, permitting direct study of volcanic and plutonic rocks within individual calderas. The Caetano caldera in north–central Nevada, formed during the mid-Tertiary ignimbrite flare-up, offers one of the best exposed and most complete records of caldera magmatism. Integrating whole-rock geochemistry, mineral chemistry, isotope geochemistry and geochronology with field studies and geologic mapping, we define the petrologic evolution of the magmatic system that sourced the >1100 km 3 Caetano Tuff. The intra-caldera Caetano Tuff is up to ∼5 km thick, composed of crystal-rich (30–45 vol. %), high-silica rhyolite, overlain by a smaller volume of comparably crystal-rich, low-silica rhyolite. It defies classification as either a monotonous intermediate or crystal-poor zoned rhyolite, as commonly ascribed to ignimbrite eruptions. Crystallization modeling based on the observed mineralogy and major and trace element geochemistry demonstrates that the compositional zonation can be explained by liquid–cumulate evolution in the Caetano Tuff magma chamber, with the more evolved lower Caetano Tuff consisting of extracted liquids that continued to crystallize and mix in the upper part of the chamber following segregation from a cumulate-rich, and more heterogeneous, source mush. The latter is represented in the caldera stratigraphy by the less evolved upper Caetano Tuff. Whole-rock major, trace and rare earth element geochemistry, modal mineralogy and mineral chemistry, O, Sr, Nd and Pb isotope geochemistry, sanidine Ar–Ar geochronology, and zircon U–Pb geochronology and trace element geochemistry provide robust evidence that the voluminous caldera intrusions (Carico Lake pluton and Redrock Canyon porphyry) are genetically equivalent to the least evolved Caetano Tuff and formed from magma that remained in the lower chamber after ignimbrite eruption and caldera collapse. Thus, the Caetano Tuff contradicts models for the mutually exclusive origins of voluminous volcanic and plutonic magmas in the upper crust. Crystal-scale O isotope data indicate that the Caetano Tuff is one of the most 18 O-enriched rhyolites in the Great Basin (δ 18 O magma = 10·2 ± 0·2‰), supporting anatexis of local metasedimentary basement crust. Metapelite xenoliths in the Carico Lake pluton and ubiquitous xenocrystic zircons in the Caetano Tuff provide constraints for the anatexis process; these data point to shallow (<15 km) dehydration melting of a protolith similar to the Proterozoic McCoy Creek Group siliciclastic sediments in eastern Nevada, projected beneath Caetano in fault-stacked shelf sediments that were thickened during Mesozoic crustal shortening. Mean zircon U–Pb ages for different stratigraphic levels of the intra-caldera Caetano Tuff are 34·2–34·5 Ma, 0·2–0·5 Myr older than the caldera sanidine 40 Ar/ 39 Ar age of 34·00 ± 0·03 Ma, documenting protracted duration of assembly and homogenization of isotopically diverse upper crustal melts, followed by crystallization and zonation to generate the Caetano Tuff magma chamber. Sanidine rims in the least evolved Caetano Tuff and in the Carico Lake pluton and Redrock Canyon porphyry have sharply zoned Ba domains that point to crystal growth during magmatic recharge events. The recharge magma is inferred to have been compositionally similar to the Caetano Tuff magma, with increased Ba resulting from remelting of Ba-rich sanidine cumulates. Mush reactivation to generate the Caetano Tuff eruption was sufficiently rapid to preserve compositional gradients in the intracaldera ignimbrite, calling into question models that predict homogeneity as a prerequisite for remobilizing crystal-rich ignimbrite magmas.

Journal of Petrology↗

USGS2014

A depth-averaged debris-flow model that includes the effects of evolving dilatancy. I. Physical basis

To simulate debris-flow behaviour from initiation to deposition, we derive a depth-averaged, two-phase model that combines concepts of critical-state soil mechanics, grain-flow mechanics and fluid mechanics. The model's balance equations describe coupled evolution of the solid volume fraction, m , basal pore-fluid pressure, flow thickness and two components of flow velocity. Basal friction is evaluated using a generalized Coulomb rule, and fluid motion is evaluated in a frame of reference that translates with the velocity of the granular phase, v s . Source terms in each of the depth-averaged balance equations account for the influence of the granular dilation rate, defined as the depth integral of ∇⋅ v s . Calculation of the dilation rate involves the effects of an elastic compressibility and an inelastic dilatancy angle proportional to m − m eq , where m eq is the value of m in equilibrium with the ambient stress state and flow rate. Normalization of the model equations shows that predicted debris-flow behaviour depends principally on the initial value of m − m eq and on the ratio of two fundamental timescales. One of these timescales governs downslope debris-flow motion, and the other governs pore-pressure relaxation that modifies Coulomb friction and regulates evolution of m . A companion paper presents a suite of model predictions and tests.

Proceedings of the Royal Society A↗

USGS2014

A depth-averaged debris-flow model that includes the effects of evolving dilatancy: II. Numerical predictions and experimental tests.

We evaluate a new depth-averaged mathematical model that is designed to simulate all stages of debris-flow motion, from initiation to deposition. A companion paper shows how the model’s five governing equations describe simultaneous evolution of flow thickness, solid volume fraction, basal pore-fluid pressure, and two components of flow momentum. Each equation contains a source term that represents the influence of state-dependent granular dilatancy. Here we recapitulate the equations and analyze their eigenstructure to show that they form a hyperbolic system with desirable stability properties. To solve the equations we use a shock-capturing numerical scheme with adaptive mesh refinement, implemented in an open-source software package we call D-Claw. As tests of D-Claw, we compare model output with results from two sets of large-scale debris-flow experiments. One set focuses on flow initiation from landslides triggered by rising pore-water pressures, and the other focuses on downstream flow dynamics, runout, and deposition. D-Claw performs well in predicting evolution of flow speeds, thicknesses, and basal pore-fluid pressures measured in each type of experiment. Computational results illustrate the critical role of dilatancy in linking coevolution of the solid volume fraction and pore-fluid pressure, which mediates basal Coulomb friction and thereby regulates debris-flow dynamics.

Proceedings of the Royal Society A↗

USGS2007

Chromosome painting among Proboscidea, Hyracoidea and Sirenia: Support for Paenungulata (Afrotheria, Mammalia) but not Tethytheria

Despite marked improvements in the interpretation of systematic relationships within Eutheria, particular nodes, including Paenungulata (Hyracoidea, Sirenia and Proboscidea), remain ambiguous. The combination of a rapid radiation, a deep divergence and an extensive morphological diversification has resulted in a limited phylogenetic signal confounding resolution within this clade both at the morphological and nucleotide levels. Cross-species chromosome painting was used to delineate regions of homology between Loxodonta africana (2n = 56), Procavia capensis (2n=54), Trichechus manatus latirostris (2n = 48) and an outgroup taxon, the aardvark (Orycteropus afer, 2n = 20). Changes specific to each lineage were identified and although the presence of a minimum of 11 synapomorphies confirmed the monophyly of Paenungulata, no change characterizing intrapaenungulate relationships was evident. The reconstruction of an ancestral paenungulate karyotype and the estimation of rates of chromosomal evolution indicate a reduced rate of genomic repatterning following the paenungulate radiation. In comparison to data available for other mammalian taxa, the paenungulate rate of chromosomal evolution is slow to moderate. As a consequence, the absence of a chromosomal character uniting two paenungulates (at the level of resolution characterized in this study) may be due to a reduced rate of chromosomal change relative to the length of time separating successive divergence events. ?? 2007 The Royal Society.

Proceedings of the Royal Society B: Biological Sci↗

USGS2019

Is there convergence of gut microbes in blood-feeding vertebrates?

Animal microbiomes play an important role in dietary adaptation, yet the extent to which microbiome changes exhibit parallel evolution is unclear. Of particular interest is an adaptation to extreme diets, such as blood, which poses special challenges in its content of proteins and lack of essential nutrients. In this study, we assessed taxonomic signatures (by 16S rRNA amplicon profiling) and potential functional signatures (inferred by Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt)) of haematophagy in birds and bats. Our goal was to test three alternative hypotheses: no convergence of microbiomes, convergence in taxonomy and convergence in function. We find a statistically significant effect of haematophagy in terms of microbial taxonomic convergence across the blood-feeding bats and birds, although this effect is small compared to the differences found between haematophagous and non-haematophagous species within the two host clades. We also find some evidence of convergence at the predicted functional level, although it is possible that the lack of metagenomic data and the poor representation of microbial lineages adapted to haematophagy in genome databases limit the power of this approach. The results provide a paradigm for exploring convergent microbiome evolution replicated with independent contrasts in different host lineages.

Philosophical Transactions of the Royal Society B:↗

USGS2021

The roles of antimicrobial resistance, phage diversity, isolation source, and selection in shaping the genomic architecture of Bacillus anthracis

Bacillus anthracis, the causative agent of anthrax disease, is a worldwide threat to livestock, wildlife and public health. While analyses of genetic data from across the globe have increased our understanding of this bacterium’s population genomic structure, the influence of selective pressures on this successful pathogen is not well understood. In this study, we investigate the effects of antimicrobial resistance, phage diversity, geography and isolation source in shaping population genomic structure. We also identify a suite of candidate genes potentially under selection, driving patterns of diversity across 356 globally extant B. anthracis genomes. We report ten antimicrobial resistance genes and 11 different prophage sequences, resulting in the first large-scale documentation of these genetic anomalies for this pathogen. Results of random forest classification suggest genomic structure may be driven by a combination of antimicrobial resistance, geography and isolation source, specific to the population cluster examined. We found strong evidence that a recombination event linked to a gene involved in protein synthesis may be responsible for phenotypic differences between comparatively disparate populations. We also offer a list of genes for further examination of B. anthracis evolution, based on high-impact single nucleotide polymorphisms (SNPs) and clustered mutations. The information presented here sheds new light on the factors driving genomic structure in this notorious pathogen and may act as a road map for future studies aimed at understanding functional differences in terms of B. anthracis biogeography, virulence and evolution.

Microbial Genomics↗

USGS2025

Diversity of polyomaviruses and papillomaviruses in penguins from eastern and western Antarctica

Polyomaviruses and papillomaviruses are icosahedral viruses with small circular dsDNA genomes. Limited information on their diversity and evolution in avian hosts is available, with even less known regarding Antarctic penguins. Prior to this study, only one polyomavirus and two papillomaviruses had been identified in Adélie penguins ( Pygoscelis adeliae ). To expand our knowledge of these viruses in Antarctic penguins, we collected faecal and cloacal swab samples from 246 Adélie penguins over 3 breeding seasons (2021–2024) and 10 emperor penguins ( Aptenodytes forsteri ) during the 2023–2024 season on Ross Island (Ross Sea). Additionally, we sampled 66 Adélie, 40 chinstrap ( Pygoscelis antarcticus ) and 71 gentoo ( Pygoscelis papua ) penguins during the 2022–2023 season across various sites on the Antarctic Peninsula. All samples were screened for papillomaviruses and polyomaviruses. We identified 31 polyomaviruses in Adélie, gentoo and chinstrap penguins and 4 papillomaviruses in Adélie penguins sampled in both eastern and western Antarctica. The 31 penguin polyomaviruses belong to a single species but form four distinct variants that are host species specific with strong geographic clustering. The four papillomaviruses represent three different types, of which two are new types from Adélie penguins sampled on Yalour Island in the West Antarctic Peninsula. Co-occurrence of two polyomavirus variants was identified in two individual gentoo penguins. Both of these variants appear to be circulating in gentoo penguins at Cierva Cove, Hope Bay in Trinity Peninsula along the Antarctic Peninsula, and at Hannah Point on Livingstone Island and Stranger Point on King George Island in the South Shetland Islands. Here, we expand the known diversity, host and geographical ranges of penguin polyomaviruses and, together with a previously identified polyomavirus on Ross Island from 2012 to 2013, show that they form five distinct lineages. The four papillomaviruses identified in this study, together with two previously identified from Ross Island in 2012 and 2013 breeding seasons, show substantial diversity reflecting four papillomavirus types across three viral species and two distinct genera. Continued surveillance and viral genomic analysis across a larger geographical framework will help understand the evolution, transmission and incidence rates of these viruses.

Microbial Genomics↗

USGS2003

Phylogeography of infectious haematopoietic necrosis virus in North America

Infectious hematopoietic necrosis virus (IHNV) is a rhabdoviral pathogen that infects wild and cultured salmonid fish throughout the Pacific Northwest of North America. IHNV causes severe epidemics in young fish and can cause disease or occur asymptomatically in adults. In a broad survey of 323 IHNV field isolates, sequence analysis of a 303 nucleotide variable region within the glycoprotein gene revealed a maximum nucleotide diversity of 8·6 %, indicating low genetic diversity overall for this virus. Phylogenetic analysis revealed three major virus genogroups, designated U, M and L, which varied in topography and geographical range. Intragenogroup genetic diversity measures indicated that the M genogroup had three- to fourfold more diversity than the other genogroups and suggested relatively rapid evolution of the M genogroup and stasis within the U genogroup. We speculate that factors influencing IHNV evolution may have included ocean migration ranges of their salmonid host populations and anthropogenic effects associated with fish culture.

Journal of General Virology↗

USGS2007

In situ time-series measurements of subseafloor sediment properties

The capabilities and diversity of subsurface sediment sensors lags significantly from what is available for the water column, thereby limiting progress in understanding time-dependent seabed exchange and high-frequency acoustics. To help redress this imbalance, a new instrument, the autonomous sediment profiler (ASP), is described herein. ASP consists of a four-electrode, Wenner-type resistivity probe and a thermistor that log data at 0.1-cm vertical intervals over a 58-cm vertical profile. To avoid resampling the same spot on the seafloor, the probes are moved horizontally within a 20 times 100-cm -2 area in one of three preselected patterns. Memory and power capacities permit sampling at hourly intervals for up to 3-mo duration. The system was tested in a laboratory tank and shown to be able to resolve high-frequency sediment consolidation, as well as changes in sediment roughness. In a field test off the southern coast of France, the system collected resistivity and temperature data at hourly intervals for 16 d. Coupled with environmental data collected on waves, currents, and suspended sediment, the ASP is shown to be useful for understanding temporal evolution of subsurface sediment porosity, although no large depositional or erosional events occurred during the deployment. Following a rapid decrease in bottom-water temperature, the evolution of the subsurface temperature field was consistent with the 1-D thermal diffusion equation coupled with advection in the upper 3-4 cm. Collectively, the laboratory and field tests yielded promising results on time-dependent seabed change.

IEEE Journal of Oceanic Engineering↗

USGS2023

Landscape transcriptomics as a tool for addressing global change effects across diverse species

Landscape transcriptomics is an emerging field studying how genome-wide expression patterns reflect dynamic landscape-scale environmental drivers, including habitat, weather, climate, and contaminants, and the subsequent effects on organismal function. This field is benefitting from advancing and increasingly accessible molecular technologies, which in turn are allowing the necessary characterization of transcriptomes from wild individuals distributed across natural landscapes. This research is especially important given the rapid pace of anthropogenic environmental change and potential impacts that span levels of biological organization. We discuss three major themes in landscape transcriptomic research: connecting transcriptome variation across landscapes to environmental variation, generating and testing hypotheses about the mechanisms and evolution of transcriptomic responses to the environment, and applying this knowledge to species conservation and management. We discuss challenges associated with this approach and suggest potential solutions. We conclude that landscape transcriptomics has great promise for addressing fundamental questions in organismal biology, ecology, and evolution, while providing tools needed for conservation and management of species.

Molecular Ecology Resources↗

USGS2024

The underlying causes of differential migration: Assumptions, hypotheses, and predictions

Mechanisms governing the migratory decisions of birds have long fascinated ecologists and sparked considerable debate. Identifying factors responsible for variation in migration distance, also known as differential migration, has been a popular approach to understanding the mechanisms underlying migratory behaviour more generally. However, research progress has been slowed by the continued testing of overlapping, non-mechanistic, and circular predictions among a small set of historically entrenched hypotheses. We highlight the body size hypothesis and suggest that the predictions commonly tested have impeded progress because body size relationships with migration distance are predictions made by several distinct hypotheses with contrasting mechanisms. The cost of migration itself has not been adequately accounted for in most hypotheses, and we propose two flight efficiency hypotheses with time- and energy-minimizing mechanisms that allow individuals to mitigate the risks inherent to longer migrations. We also advance two conceptual versions of the social dominance hypothesis based on two distinct underlying mechanisms related to distance minimization and food maximization that will help clarify the role of competition in driving migratory decisions. Overall, we describe and refine 12 mechanistic hypotheses proposed to explain differential migration (along with several other special-case hypotheses), seven of which have underlying mechanisms related to food limitation as past research has identified this to be an important driver of differential migration. We also thoroughly reviewed 145 publications to assess the amount of support for 10 critical assumptions underlying alternative hypotheses for differential migration in birds. Our review reveals that surprisingly few studies explicitly evaluate assumptions within a differential migration context. Generating and testing strong predictions and critical assumptions underlying mechanisms of alternative hypotheses will improve our ability to differentiate among these explanations of differential migration. Additionally, future intraspecific progress will be greatest if investigators continue to focus on mechanisms underlying variation in migration distance within rather than among demographic classes, as previous research has found differing mechanisms to be responsible for differential migration among demographic classes. Interspecifically, a thorough comparative analysis that seeks to explain variation in migration distance among species would broaden both our understanding of the mechanisms regulating current differential migration patterns and those that led to the evolution of migration more generally. Collectively, we provide a framework that, together with advances in animal-borne tracking and other technology, can be used to advance our understanding of the causes of differential migration distance, and migratory decisions more generally.

Biological Reviews↗

USGS2018

Phylogeny and species traits predict bird detectability

Avian acoustic communication has resulted from evolutionary pressures and ecological constraints. We therefore expect that auditory detectability in birds might be predictable by species traits and phylogenetic relatedness. We evaluated the relationship between phylogeny, species traits, and field‐based estimates of the two processes that determine species detectability (singing rate and detection distance) for 141 bird species breeding in boreal North America. We used phylogenetic mixed models and cross‐validation to compare the relative merits of using trait data only, phylogeny only, or the combination of both to predict detectability. We found a strong phylogenetic signal in both singing rates and detection distances; however the strength of phylogenetic effects was less than expected under Brownian motion evolution. The evolution of behavioural traits that determine singing rates was found to be more labile, leaving more room for species to evolve independently, whereas detection distance was mostly determined by anatomy (i.e. body size) and thus the laws of physics. Our findings can help in disentangling how complex ecological and evolutionary mechanisms have shaped different aspects of detectability in boreal birds. Such information can greatly inform single‐ and multi‐species models but more work is required to better understand how to best correct possible biases in phylogenetic diversity and other community metrics.

Ecography↗

USGS2014

Hybridization of an invasive shrub affects tolerance and resistance to defoliation by a biological control agent

Evolution has contributed to the successful invasion of exotic plant species in their introduced ranges, but how evolution affects particular control strategies is still under evaluation. For instance, classical biological control, a common strategy involving the utilization of highly specific natural enemies to control exotic pests, may be negatively affected by host hybridization because of shifts in plant traits, such as root allocation or chemical constituents. We investigated introgression between two parent species of the invasive shrub tamarisk ( Tamarix spp.) in the western United States, and how differences in plant traits affect interactions with a biological control agent. Introgression varied strongly with latitude of origin and was highly correlated with plant performance. Increased levels of T. ramosissima introgression resulted in both higher investment in roots and tolerance to defoliation and less resistance to insect attack. Because tamarisk hybridization occurs predictably on the western U.S. landscape, managers may be able to exploit this information to maximize control efforts. Genetic differentiation in plant traits in this system underpins the importance of plant hybridization and may explain why some biological control releases are more successful than others.

Colorado↗

USGS2017

Incorporating evolutionary insights to improve ecotoxicology for freshwater species

Ecotoxicological studies have provided extensive insights into the lethal and sublethal effects of environmental contaminants. These insights are critical for environmental regulatory frameworks, which rely on knowledge of toxicity for developing policies to manage contaminants. While varied approaches have been applied to ecotoxicological questions, perspectives related to the evolutionary history of focal species or populations have received little consideration. Here, we evaluate chloride toxicity from the perspectives of both macroevolution and contemporary evolution. First, by mapping chloride toxicity values derived from the literature onto a phylogeny of macroinvertebrates, fish, and amphibians, we tested whether macroevolutionary relationships across species and taxa are predictive of chloride tolerance. Next, we conducted chloride exposure tests for two amphibian species to assess whether potential contemporary evolutionary change associated with environmental chloride contamination influences chloride tolerance across local populations. We show that explicitly evaluating both macroevolution and contemporary evolution can provide important and even qualitatively different insights from those obtained via traditional ecotoxicological studies. While macroevolutionary perspectives can help forecast toxicological end points for species with untested sensitivities, contemporary evolutionary perspectives demonstrate the need to consider the environmental context of exposed populations when measuring toxicity. Accounting for divergence among populations of interest can provide more accurate and relevant information related to the sensitivity of populations that may be evolving in response to selection from contaminant exposure. Our data show that approaches accounting for and specifically examining variation among natural populations should become standard practice in ecotoxicology.

Evolutionary Applications↗

USGS2015

PHT3D-UZF: A reactive transport model for variably-saturated porous media

A modified version of the MODFLOW/MT3DMS-based reactive transport model PHT3D was developed to extend current reactive transport capabilities to the variably-saturated component of the subsurface system and incorporate diffusive reactive transport of gaseous species. Referred to as PHT3D-UZF, this code incorporates flux terms calculated by MODFLOW's unsaturated-zone flow (UZF1) package. A volume-averaged approach similar to the method used in UZF-MT3DMS was adopted. The PHREEQC-based computation of chemical processes within PHT3D-UZF in combination with the analytical solution method of UZF1 allows for comprehensive reactive transport investigations (i.e., biogeochemical transformations) that jointly involve saturated and unsaturated zone processes. Intended for regional-scale applications, UZF1 simulates downward-only flux within the unsaturated zone. The model was tested by comparing simulation results with those of existing numerical models. The comparison was performed for several benchmark problems that cover a range of important hydrological and reactive transport processes. A 2D simulation scenario was defined to illustrate the geochemical evolution following dewatering in a sandy acid sulfate soil environment. Other potential applications include the simulation of biogeochemical processes in variably-saturated systems that track the transport and fate of agricultural pollutants, nutrients, natural and xenobiotic organic compounds and micropollutants such as pharmaceuticals, as well as the evolution of isotope patterns.

Ground Water↗

USGS2007

Sediment compaction rates and subsidence in deltaic plains: Numerical constraints and stratigraphic influences

Natural sediment compaction in deltaic plains influences subsidence rates and the evolution of deltaic morphology. Determining compaction rates requires detailed knowledge of subsurface geotechnical properties and depositional history, neither of which is often readily available. To overcome this lack of knowledge, we numerically forward model the incremental sedimentation and compaction of stochastically generated stratigraphies with geotechnical properties typical of modern depositional environments in the Mississippi River delta plain. Using a Monte Carlo approach, the range of probable compaction rates for stratigraphies with compacted thicknesses <150 m and accumulation times <20 kyr. varies, but maximum values rarely exceed a few mm yr −1 . The fastest compacting stratigraphies are composed primarily of peat and bar sand, whereas the slowest compacting stratigraphies are composed of prodelta mud and natural levee deposits. These results suggest that compaction rates can significantly influence vertical and lateral stratigraphic trends during deltaic evolution.

Basin Research↗

USGS2004

A viscoelastic damage model with applications to stable and unstable fracturing

A viscoelastic damage rheology model is presented that provides a generalization of Maxwell viscoelasticity to a non-linear continuum mechanics framework incorporating material degradation and recovery, transition from stable to unstable fracturing and gradual accumulation of non-reversible deformation. The model is a further development of the damage rheology framework of Lyakhovsky et al. for evolving effective elasticity. The framework provides a quantitative treatment for macroscopic effects of evolving distributed cracking with local density represented by an intensive state variable. The formulation, based on thermodynamic principles, leads to a system of kinetic equations for the evolution of damage. An effective viscosity inversely proportional to the rate of damage increase is introduced to account for gradual accumulation of irreversible deformation due to dissipative processes. A power-law relation between the damage variable and elastic moduli leads to a non-linear coupling between the rate of damage evolution and the damage variable itself. This allows the model to reproduce a transition from stable to unstable fracturing of brittle rocks and the Kaiser effect. 3-D numerical simulations based on the model formulation for homogeneous and heterogeneous materials account for the main features of rock behaviour under large strain. The model coefficients are constrained, using triaxial laboratory experiments with low-porosity Westerly granite and high-porosity Berea sandstone samples.

Geophysical Journal International↗