USGS Science⌕ Search

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.

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

At least 1,189 records · Page 66Linked to original sources
USGS2018

Benthic foraminifera from the Carnarvon Ramp reveal variability in Leeuwin Current activity (Western Australia) since the Pliocene

Benthic foraminiferal assemblages from a ~300 m deep core from an outer carbonate-ramp site off Western Australia (International Ocean Discovery Program Core U1460A) were examined to reconstruct the paleoceanographic evolution of the Carnarvon Ramp and the warm surficial Leeuwin Current (LC) for the last 3.54 Ma. Of the identified 179 benthic foraminiferal species, occurrences of the 15 most abundant taxa were assessed using Q-mode Cluster Analysis and Non-Metric Dimensional Scaling. Diversity, equitability, planktonic/benthic index, microhabitat preference, and sedimentary parameters such as lithology and sponge spicule content were analyzed to gather information about past intermediate- and surface-water circulation. Relative abundances of infaunal and epifaunal species were applied to indicate changes in organic-matter supply and oxygenation at the sea floor. Influence of upwelling was recognized by a high infaunal species ratio, with dominance by Uvigerina peregrina , Lagena annellatrachia and Trifarina bradyi . Epifaunal species such as Hanzawaia nipponica and Hyalinea florenceae gradually became more abundant around 1.14 Ma, indicating increased ventilation and establishment of the paleo-LC. A more substantial change was initiated by 0.91 Ma as marked by key species Spirorutilus carinatus and Rotorbinella sp. , together with increased faunal diversity, benthic foraminiferal accumulation rates, and evidence for suspension feeding sponges. With the LC flow suppressing upwelling, and better ventilated waters entering the shelf, the environment favored epifaunal agglutinates, rotalids, and miliolids, while buliminids decreased. Under high-flow conditions of the LC, sponge spicules and skeletal carbonate production reached an optimum at ~0.6 Ma before returning to modern conditions. Supported by these observations, we propose the following paleoceanographic evolution of the Carnarvon Ramp: During the late Pliocene to mid Pleistocene (3.54–0.91 Ma) conditions of deep-water upwelling from the Western Australian Current and Indian Ocean Gyre indicate the absence of the capping LC on the outer carbonate ramp. A transitional phase started in the mid Pleistocene (1.14–0.61 Ma). The paleo-LC triggered gradual oxygenation at the sediment-water interface, which coincided with an increase in carbonate sedimentation rates, and waning sea-surface productivity. During a third phase, mid Pleistocene to present (0.91–0 Ma), the LC’s intensity and flow rates peaked at ~0.6 Ma. Benthic foraminiferal accumulation rates reached a high, then decreased to present-day rates. For short periods, sea-surface productivity was moderately enhanced, likely due to fluctuating LC persistence or landward shift during glacial maxima.

Western Australia↗
USGS2009

The B7 family of immunoregulatory receptors: A comparative and evolutionary perspective

In mammals, T cell activation requires specific recognition of the peptide-MHC complex by the TcR and co-stimulatory signals. Important co-stimulatory receptors expressed by T cells are the molecules of the CD28 family, that regulate T cell activation, proliferation and tolerance. These receptors recognize B7s and B7-homologous (B7H) molecules that are typically expressed by the antigen presenting cells. In teleost fish, typical T cell responses have been described and the TcR, MHC and CD28/CTLA4 genes have been characterized. In contrast, the members of the B7 gene family have only been described in mammals and birds and have yet to be addressed in lower vertebrates. To learn more about the evolution of components guiding T cell activation in vertebrates, we performed a systematic genomic survey for the B7 co-stimulatory and co-inhibitory IgSF receptors in lower vertebrates with an emphasis on teleost fish. Our search identified fish sequences that are orthologous to B7, B7-H1/B7-DC, B7-H3 and B7-H4 as defined by sequence identity, phylogeny and combinations of short or long-range syntenic relationships. However, we were unable to identify clear orthologs for B7-H2 (CD275, ICOS ligand) in bony fish, which correlates with our prior inability to find ICOS in fish. Interestingly, our results indicate that teleost fish possess a single B7.1/B7.2 (CD80/86) molecule that likely interacts with CD28/CTLA4 as the ligand-binding regions seem to be conserved in both partners. Overall, our analyses implies that gene duplication (and loss) have shaped a molecular repertoire of B7-like molecules that was recruited for the refinement of T cell activation during the evolution of the vertebrates.

Molecular Immunology↗
USGS2012

Ocean-atmosphere dynamics during Hurricane Ida and Nor'Ida: An application of the coupled ocean-;atmosphere–wave–sediment transport (COAWST) modeling system

The coupled ocean–atmosphere–wave–sediment transport (COAWST) modeling system was used to investigate atmosphere–ocean–wave interactions in November 2009 during Hurricane Ida and its subsequent evolution to Nor'Ida, which was one of the most costly storm systems of the past two decades. One interesting aspect of this event is that it included two unique atmospheric extreme conditions, a hurricane and a nor'easter storm, which developed in regions with different oceanographic characteristics. Our modeled results were compared with several data sources, including GOES satellite infrared data, JASON-1 and JASON-2 altimeter data, CODAR measurements, and wave and tidal information from the National Data Buoy Center (NDBC) and the National Tidal Database. By performing a series of numerical runs, we were able to isolate the effect of the interaction terms between the atmosphere (modeled with Weather Research and Forecasting, the WRF model), the ocean (modeled with Regional Ocean Modeling System (ROMS)), and the wave propagation and generation model (modeled with Simulating Waves Nearshore (SWAN)). Special attention was given to the role of the ocean surface roughness. Three different ocean roughness closure models were analyzed: DGHQ (which is based on wave age), TY2001 (which is based on wave steepness), and OOST (which considers both the effects of wave age and steepness). Including the ocean roughness in the atmospheric module improved the wind intensity estimation and therefore also the wind waves, surface currents, and storm surge amplitude. For example, during the passage of Hurricane Ida through the Gulf of Mexico, the wind speeds were reduced due to wave-induced ocean roughness, resulting in better agreement with the measured winds. During Nor'Ida, including the wave-induced surface roughness changed the form and dimension of the main low pressure cell, affecting the intensity and direction of the winds. The combined wave age- and wave steepness-based parameterization (OOST) provided the best results for wind and wave growth prediction. However, the best agreement between the measured (CODAR) and computed surface currents and storm surge values was obtained with the wave steepness-based roughness parameterization (TY2001), although the differences obtained with respect to DGHQ were not significant. The influence of sea surface temperature (SST) fields on the atmospheric boundary layer dynamics was examined; in particular, we evaluated how the SST affects wind wave generation, surface currents and storm surges. The integrated hydrograph and integrated wave height, parameters that are highly correlated with the storm damage potential, were found to be highly sensitive to the ocean surface roughness parameterization.

Ocean Modelling↗
USGS2022

Smoothed Particle Hydrodynamics simulations of reef surf zone processes driven by plunging irregular waves

As waves interact with the slopes of coral reefs and other steep bathymetry profiles, plunging breaking usually occurs where the free surface overturns and violent water motion is triggered. Resolving these surf zone processes pose significant challenges for conventional mesh-based hydrodynamic models, due to the rapidly-deforming nature of the free surface and associated flows. Yet the accurate prediction of these surf zone hydrodynamics is critical for predicting a wide range of nearshore processes driven by wave breaking (e.g., wave dissipation and energy transfers; mean water levels and currents; and wave runup). In this study we assess the ability of the mesh-free, Lagrangian particle-based numerical modelling approach Smoothed Particle Hydrodynamics (SPH) based on DualSPHysics, to simulate the fine-scale hydrodynamic processes driven by irregular wave transformation over a fringing reef profile, by comparing results against detailed experimental observations from a physical modelling study. To greatly improve the computational efficiency, the SPH model was coupled to the mesh-based multi-layer nonhydrostatic wave-flow model SWASH. With this coupled approach, SWASH was used to efficiently simulate the evolution of non-breaking waves from the wavemaker up to the fore reef slope, with the SPH model then used to simulate the detailed hydrodynamic processes over the reef from just offshore of the breakpoint to the shoreline . The SPH model was able to accurately reproduce the complex free surface deformations during plunging breaking, the spectral evolution of waves across the reef flat (including nonlinear wave shape), the mean water levels and currents, and wave runup at the shoreline. Using the long duration simulations (>400 wave periods), the model was able to reproduce the full range of wave motions over the reef (from sea-swell to infragravity frequencies), including the increasing dominance of low frequency waves towards the shoreline and the large cross-reef standing wave motions excited by the reef geometry.

Ocean Modelling↗
USGS2019

Crustal architecture beneath the southern Midcontinent (USA) and controls on Mesoproterozoic iron-oxide mineralization from 3D geophysical models

Several types of critical mineral-bearing ore deposits in the southern Midcontinent region of the U.S. are hosted in Mesoproterozoic igneous rocks largely concealed beneath Paleozoic cover. Discerning the architecture of igneous intrusions and volcanic centers in the crust is fundamental to understanding the geologic evolution of this vast region and its mineral resources. To advance the understanding of the geologic framework beneath the Southeast Missouri Iron Metallogenic Province, we invert continental-scale magnetic and gravity anomaly data to three-dimensional (3D) physical property models. The regional models image altered and mineralized igneous rocks near the Precambrian basement surface and underlying intrusive complexes that extend down to the Moho. At shallow crustal levels, our models confirm that iron oxide-apatite ± rare earth element (IOA±REE) deposits and iron oxide-copper-gold ± cobalt (IOCG) deposits occur within or near the edges of large low density/low susceptibility early Mesoproterozoic (ca 1.4 Ga) silicic calderas and (ca 1.3 Ga) granitic plutons. Previous isotopic and geochemical studies conclude that the iron deposits and their volcanic host rocks originated from mantle-derived and crustal melts that erupted during regional extension. Extension was associated with thermal event(s) that produced the large-scale silicic magmatism related to the ca 1.45 Ga Eastern Granite Rhyolite Province (EGRP) and the 1.35 Ga Southern Granite Rhyolite Province (SGRP). We postulate that early in the evolution of the EGRP, several trans-crustal magmatic plumbing systems developed that are evident in the 3D models. The Southeast Missouri Metallogenic Province is underlain by one such magmatic system that is expressed as a northwest-trending ~ 50 km-wide by 200 km-long elongate track of high susceptibility at deep crustal levels. The high susceptibility corridor splays upward through the crust to the Precambrian surface where the iron deposits are the epigenetic manifestation of this magmatic event. Our findings confirm that the iron deposits, with no distinct connection at the surface, are connected to one large magmatic system at depth. We propose that other similar susceptibility tracks, which are present along the top of the mantle, mark additional feeder zones that allowed magma to ascend to the main eruptive centers that produced the Granite Rhyolite Provinces. The early Mesoproterozoic extensional tectonic framework established crustal-scale pathways that controlled the distribution of subsequent magmatic activity, including the ca 1.4 Ga calderas and underlying intrusions, ca 1.3 Ga silicic plutons and Phanerozoic alkaline intrusions. If these interpretations are correct, our study has identified large areas that are prospective for critical mineral-bearing ore deposits and, importantly, suggests that the Mesoproterozoic architecture may have influenced subsequent magmatism and hydrothermal activity in the southern Midcontinent of the U.S.

Ore Geology Reviews↗
USGS2017

Sedimentology, sequence-stratigraphy, and geochemical variations in the Mesoproterozoic Nonesuch Formation, northern Wisconsin, USA

We use core descriptions and portable X-ray fluorescence analyses to identify lithofacies and stratigraphic surfaces for the Mesoproterozoic Nonesuch Formation within the Ashland syncline, Wisconsin. We group lithofacies into facies associations and construct a sequence stratigraphic framework based on lithofacies stacking and stratigraphic surfaces. The fluvial-alluvial facies association (upper Copper Harbor Conglomerate) is overlain across a transgressive surface by the fluctuating-profundal facies association (lower Nonesuch Formation). The fluctuating-profundal facies association comprises a retrogradational sequence set overlain across a maximum flooding surface by an aggradational-progradational sequence set comprising fluctuating-profundal, fluvial-lacustrine, and fluvial-alluvial facies associations (middle Nonesuch through lower Freda Formations). Lithogeochemistry supports sedimentologic and stratigraphic interpretations. Fe/S molar ratios reflect the oxidation state of the lithofacies; values are most depleted above the maximum flooding surface where lithofacies are chemically reduced and are greatest in the chemically oxidized lithofacies. Si/Al and Zr/Al molar ratios reflect the relative abundance of detrital heavy minerals vs. clay minerals; greater values correlate with larger grain size. Vertical facies association stacking records depositional environments that evolved from fluvial and alluvial, to balanced-fill lake, to overfilled lake, and returning to fluvial and alluvial. Elsewhere in the basin, where accommodation was greatest, some volume of fluvial-lacustrine facies is likely present below the transgressive stratigraphic surface. This succession of continental and lake-basin types indicates a predominant tectonic driver of basin evolution. Lithofacies distribution and geochemistry indicate deposition within an asymmetric half-graben bounded on the east by a west-dipping growth fault. While facies assemblages are lacustrine and continental, periodic marine incursions are probable, especially across maximum transgressive surfaces. We demonstrate a sequence-stratigraphic approach may be applied to fine-grained Precambrian sediments using traditional rock description and supporting lithogeochemistry. Identification of a characteristic lithofacies succession in Mesoproterozoic sediments demonstrates fundamental controls commonly interpreted for Phanerozoic lake systems may be extended into the Precambrian. These controls result in a predictable association of lithofacies, with distinct physical, biological, and geochemical properties. This has regional significance for carbon sequestration and the distribution of mineral and hydrocarbon resources and broader significance for addressing Mesoproterozoic paleogeographic reconstructions and questions related to the evolution of terrestrial life.

Wisconsin↗
USGS2019

Middle Pleistocene formation of the Rio Grande Gorge, San Luis Valley, south-central Colorado and north-central New Mexico, USA: Process, timing, and downstream implications

The Rio Grande is the fourth longest river in North America extending over 3,000 km from the Rocky Mountains to the Gulf of Mexico. The Pleistocene evolution of this river from individual subbasins into a coalesced fluvial system has been long debated. Herein, we constrain the middle Pleistocene evolution of the northernmost and largest Rio Grande basin, the San Luis basin, and the timing of incision of the Rio Grande Gorge, based on new geologic mapping, 3 He surface exposure dating, and U-series dating on pedogenic carbonates. 3 He dating of shoreline deposits of closed basin Lake Alamosa and fluvially scoured bedrock along the Rio Grande gorge rim indicate the San Luis basin was connected with southerly Rio Grande basins and progressively incised since ∼400 ka. Waning tectonic activity, coupled with the hydrologic response to intensified middle Pleistocene glaciation, drove southward spillover of Lake Alamosa ∼400 ka, expanding the river's drainage basin by nearly 24,000 km 2 , and adding recharge areas in the high-altitude, glaciated San Juan, Tusas, and Sangre de Cristo Mountains. Geologic barriers to streamflow on the Taos Plateau were breached and basin tributary drainages were integrated and inset into their present canyons ∼250-200 ka, with approximately 120 m of incision during ∼400-200 ka. Previous work demonstrates that southerly streamflow through basins south of the SLB terminated into a large, rift-related bolson complex until the middle Pleistocene, and most axial Rio Grande river incision has occurred since 630 ka. We propose that middle Pleistocene integration of the SLB was the major geomorphic event that eroded Plio-Pleistocene tectonic and volcanic barriers and re-established flow of the Rio Grande from the Rocky Mountains to the Gulf of Mexico.

Colorado, New Mexico↗
USGS1999

Multivariate analyses of Erzgebirge granite and rhyolite composition: Implications for classification of granites and their genetic relations

High-precision major, minor and trace element analyses for 44 elements have been made of 329 Late Variscan granitic and rhyolitic rocks from the Erzgebirge metallogenic province of Germany. The intrusive histories of some of these granites are not completely understood and exposures of rock are not adequate to resolve relationships between what apparently are different plutons. Therefore, it is necessary to turn to chemical analyses to decipher the evolution of the plutons and their relationships. A new classification of Erzgebirge plutons into five major groups of granites, based on petrologic interpretations of geochemical and mineralogical relationships (low-F biotite granites; low-F two-mica granites; high-F, high-P2O5 Li-mica granites; high-F, low-P2O5 Li-mica granites; high-F, low-P2O5 biotite granites) was tested by multivariate techniques. Canonical analyses of major elements, minor elements, trace elements and ratio variables all distinguish the groups with differing amounts of success. Univariate ANOVA's, in combination with forward-stepwise and backward-elimination canonical analyses, were used to select ten variables which were most effective in distinguishing groups. In a biplot, groups form distinct clusters roughly arranged along a quadratic path. Within groups, individual plutons tend to be arranged in patterns possibly reflecting granitic evolution. Canonical functions were used to classify samples of rhyolites of unknown association into the five groups. Another canonical analysis was based on ten elements traditionally used in petrology and which were important in the new classification of granites. Their biplot pattern is similar to that from statistically chosen variables but less effective at distinguishing the five groups of granites. This study shows that multivariate statistical techniques can provide significant insight into problems of granitic petrogenesis and may be superior to conventional procedures for petrological interpretation.

Computers & Geosciences↗
USGS2007

Chapter 7 Magmatic-hydrothermal fluid interaction and mineralization in alkali-syenite nodules from the Breccia Museo pyroclastic deposit, Naples, Italy

The Breccia Museo, a pyroclastic flow that crops out in the Campi Flegrei volcanic complex (Naples, Italy), contains alkali-syenite (trachyte) nodules with enrichment in Cl and incompatible elements (e.g., U, Zr, Th, and rare-earth elements). Zircon was dated at ≈52 ka, by U-Th isotope systematics using a SHRIMP. Scanning electron microscope and electron microprobe analysis of the constituent phases have documented the mineralogical and textural evolution of the nodules of feldspar and mafic accumulations on the magma chamber margins. Detailed electron microprobe data are given for alkali and plagioclase feldspar, salite to ferrosalite clinopyroxene, pargasite, ferrogargasite, magnesio-hastingsite hornblende amphibole, biotite mica, Cl-rich scapolite, and a member (probable davyne-type) of the cancrinite group. Detailed whole rock, major and minor element data are also presented for selected nodules. A wide variety of common and uncommon accessory minerals were identified such as zircon, baddeleyite, zirconolite, pollucite, sodalite, titanite, monazite, cheralite, apatite, titanomagnetite and its alteration products, scheelite, ferberite, uraninite/thorianite, uranpyrochlore, thorite, pyrite, chalcopyrite, and galena. Scanning electron microscope analysis of opened fluid inclusions identified halite, sylvite, anhydrite, tungstates, carbonates, silicates, sulfides, and phosphates; most are probably daughter minerals. Microthermometric determinations on secondary fluid inclusions hosted by alkali feldspar define a temperature regime dominated by hypersaline aqueous fluids. Fluid-inclusion temperature data and mineral-pair geothermometers for coexisting feldspars and hornblende and plagioclase were used to construct a pressure-temperature scenario for the development and evolution of the nodules. We have compared the environment of porphyry copper formation and the petrogenetic environment constructed for the studied nodules. The suite of ore minerals observed in the nodules supports a potential for mineralization, which is similar to that observed in the alkaline volcanic systems of southern Italy (Pantelleria, Pontine Archipelago, Mt. Somma-Vesuvius).

Campi Flegrei↗
USGS2011

Climate, Fire and Geology in the Convergence of Mediterranean-type Climate Ecosystems

Integrating Climate, Fire and Geology in a Fire-prone World Fire challenges the long-standing hegemony of ecology, biogeography and paleoecology that climate and soils are sufficient to explain the origin and distribution of plant species. In a world where half of the land surface is fire-prone (Krawchuk et al . 2009), understanding the past and predicting the future requires a close integration of climate, fire and geology. The dogma that fire is an anthropogenic phenomenon of little use in understanding paleoecology (Axelrod 1980, 1989), or merely incidental to vegetation development (Hopper 2009), is rapidly being replaced with a better understanding of paleofire's impact on land plant evolution (Scott 2000; Pausas & Keeley 2009). Attempts to model future global vegetation patterns have been demonstrated to be inadequate without including both natural and anthropogenic fire regimes (Bond et al . 2005). Bond and Keeley (2005) outlined the conundrum posed by alternative explanations for the present distribution of vegetation and assembly of communities. Classical explanations have invoked resource-based mechanisms that are driven by climate and soils. There are ecosystems where resource-based mechanisms may be sufficient, but on many seasonally dry landscapes ecosystem processes such as fire play a major role in the organization and evolution of vegetation.

Book chapter↗
USGS2011

Relationships between silicic plutonism and volcanism: Geochemical evidence

Field associations (voluminous ash flow deposits, rhyolitic stocks and dykes, ring complexes), evidence of repeated influxes of mafic magma, and thermal constraints indicate that many high-level silicic plutons (magma chambers) acted as open systems for considerable parts of their history. The long thermal lifetime, as well as other evidence from the volcanic record, suggests that some such systems reached a quasi-steady state in which magma input was balanced by magma output for times longer than those required for crystallisation. Reconstruction of the evolution of large, long-lived caldera-forming systems, such as that of the Jemez Mountains, New Mexico, indicates that many chambers have lost a highly fractionated silicic cap, in some cases cyclically. Crystallised plutons may contain no obvious record of this evolutionary phase. Geochemical data from silicic ash flow deposits can be used to reconstruct the volcanic stage of pluton development. Many silicic systems, especially of alkaline affinity, apparently pass from a stage in which melt evolution is dominated by crystal-liquid processes to one in which other processes may also contribute to differentiation. Apparently, the transition is most readily achieved in volatile-rich, alkaline silicic systems emplaced in complex, ancient sialic crust of the cratons. Once established, the preservation of highly fractionated caps on magma chambers requires a balance between thermal input and cooling-induced crystallisation. If heat enters the system too quickly, the cap may get stirred into the dominant magma volume by convection. If heat input is too slow, the magma body will crystallise inward from the margins, and the plutonic-consolidation stage will begin. © 1988, Royal Society of Edinburgh. All rights reserved.

Transactions of the Royal Society of Edinburgh, Ea↗
USGS2000

On rate-state and Coulomb failure models

We examine the predictions of Coulomb failure stress and rate-state frictional models. We study the change in failure time (clock advance) Δt due to stress step perturbations (i.e., coseismic static stress increases) added to "background" stressing at a constant rate (i.e., tectonic loading) at time t 0 . The predictability of Δt implies a predictable change in seismicity rate r(t)/r 0 , testable using earthquake catalogs, where r 0 is the constant rate resulting from tectonic stressing. Models of r(t)/r 0 , consistent with general properties of aftershock sequences, must predict an Omori law seismicity decay rate, a sequence duration that is less than a few percent of the mainshock cycle time and a return directly to the background rate. A Coulomb model requires that a fault remains locked during loading, that failure occur instantaneously, and that Δt is independent of t 0 . These characteristics imply an instantaneous infinite seismicity rate increase of zero duration. Numerical calculations of r(t)/r 0 for different state evolution laws show that aftershocks occur on faults extremely close to failure at the mainshock origin time, that these faults must be "Coulomb-like," and that the slip evolution law can be precluded. Real aftershock population characteristics also may constrain rate-state constitutive parameters; a may be lower than laboratory values, the stiffness may be high, and/or normal stress may be lower than lithostatic. We also compare Coulomb and rate-state models theoretically. Rate-state model fault behavior becomes more Coulomb-like as constitutive parameter a decreases relative to parameter b. This is because the slip initially decelerates, representing an initial healing of fault contacts. The deceleration is more pronounced for smaller a, more closely simulating a locked fault. Even when the rate-state Δt has Coulomb characteristics, its magnitude may differ by some constant dependent on b. In this case, a rate-state model behaves like a modified Coulomb failure model in which the failure stress threshold is lowered due to weakening, increasing the clock advance. The deviation from a non-Coulomb response also depends on the loading rate, elastic stiffness, initial conditions, and assumptions about how state evolves.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Importance of a stochastic distribution of floods and erosion thresholds in the bedrock river incision problem

Fluvial erosion of bedrock occurs during occasional flood events when boundary shear stress exceeds a critical threshold to initiate incision. Therefore efforts to model the evolution of topography over long timescales should include an erosion threshold and should be driven by a stochastic distribution of erosive events. However, most bedrock incision models ignore the threshold as a second‐order detail. In addition, climate is poorly represented in most landscape evolution models, so the quantitative relationship between erosion rate and measurable climatic variables has been elusive. Here we show that the presence of an erosion threshold, when combined with a well‐constrained, probabilistic model of storm and flood occurrence, has first‐order implications for the dynamics of river incision in tectonically active areas. First, we make a direct calculation of the critical shear stress required to pluck bedrock blocks for a field site in New York. Second, we apply a recently proposed stochastic, threshold, bedrock incision model to a series of streams in California, with known tectonic and climatic forcing. Previous work in the area has identified a weak relationship between channel gradient or relief and rock uplift rate that is not easily explained by simpler detachment‐limited models. The results with the stochastic threshold model show that even low erosion thresholds, which are exceeded in steep channels during high‐frequency flood events, fundamentally affect the predicted relationship between gradient and uplift rate in steady state rivers, in a manner consistent with the observed topography. This correspondence between theory and data is, however, nonunique; models in which a thin alluvial cover may act to inhibit channel incision in the low uplift rate zone also provide plausible explanations for the observed topography. Third, we explore the broader implications of the stochastic threshold model to the development of fluvial topography in active tectonic settings. We suggest that continued field applications of geomorphic models, including physically meaningful thresholds and stochastic climate distributions, are required to advance our knowledge of interactions among surficial, climatic, and crustal processes.

California↗
USGS2007

Bed forms created by simulated waves and currents in a large flume

The morphology and evolution of bed forms created by combinations of waves and currents were investigated using an oscillating plate in a 4-m-wide flume. Current speed ranged from 0 to 30 cm/s, maximum oscillatory velocity ranged from 20 to 48 cm/s, oscillation period was 8 s (except for one run with 12 s period), and the median grain size was 0.27 mm. The angle between oscillations and current was 90°, 60°, or 45°. At the end of each run the sand bed was photographed and ripple dimensions were measured. Ripple wavelength was also determined from sonar images collected throughout the runs. Increasing the ratio of current to wave (i.e., oscillatory) velocity decreased ripple height and wavelength, in part because of the increased fluid excursion during the wave period. Increasing the ratio of current to waves, or decreasing the angle between current and waves, increased the three-dimensionality of bed forms. During the runs, ripple wavelength increased by a factor of about 2. The average number of wave periods for evolution of ripple wavelength to 90% of its final value was 184 for two-dimensional ripples starting from a flat bed. Bed form orientations at the end of each run were compared to four potential controlling factors: the directions of waves, current, maximum instantaneous bed shear stress, and maximum gross bed form normal transport (MGBNT). The directions of waves and of MGBNT were equally good predictors of bed form orientations, and were significantly better than the other two factors.

Journal of Geophysical Research C: Oceans↗
USGS2008

Probabilistic seismic hazard in the San Francisco Bay area based on a simplified viscoelastic cycle model of fault interactions

We construct a viscoelastic cycle model of plate boundary deformation that includes the effect of time-dependent interseismic strain accumulation, coseismic strain release, and viscoelastic relaxation of the substrate beneath the seismogenic crust. For a given fault system, time-averaged stress changes at any point (not on a fault) are constrained to zero; that is, kinematic consistency is enforced for the fault system. The dates of last rupture, mean recurrence times, and the slip distributions of the (assumed) repeating ruptures are key inputs into the viscoelastic cycle model. This simple formulation allows construction of stress evolution at all points in the plate boundary zone for purposes of probabilistic seismic hazard analysis (PSHA). Stress evolution is combined with a Coulomb failure stress threshold at representative points on the fault segments to estimate the times of their respective future ruptures. In our PSHA we consider uncertainties in a four-dimensional parameter space: the rupture peridocities, slip distributions, time of last earthquake (for prehistoric ruptures) and Coulomb failure stress thresholds. We apply this methodology to the San Francisco Bay region using a recently determined fault chronology of area faults. Assuming single-segment rupture scenarios, we find that fature rupture probabilities of area faults in the coming decades are the highest for the southern Hayward, Rodgers Creek, and northern Calaveras faults. This conclusion is qualitatively similar to that of Working Group on California Earthquake Probabilities, but the probabilities derived here are significantly higher. Given that fault rupture probabilities are highly model-dependent, no single model should be used to assess to time-dependent rupture probabilities. We suggest that several models, including the present one, be used in a comprehensive PSHA methodology, as was done by Working Group on California Earthquake Probabilities.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Diurnal and vertical variability of the sensible heat and carbon dioxide budgets in the atmospheric surface layer

The diurnal and vertical variability of heat and carbon dioxide (CO2) in the atmospheric surface layer are studied by analyzing measurements from a 213 in tower in Cabauw (Netherlands). Observations of thermodynamic variables and CO2 mixing ratio as well as vertical profiles of the turbulent fluxes are used to retrieve the contribution of the budget terms in the scalar conservation equation. On the basis of the daytime evolution of turbulent fluxes, we calculate the budget terms by assuming that turbulent fluxes follow a linear profile with height. This assumption is carefully tested and the deviation ftom linearity is quantified. The budget calculation allows us to assess the importance of advection of heat and CO2 during day hours for three selected days. It is found that, under nonadvective conditions, the diurnal variability of temperature and CO2 is well reproduced from the flux divergence measurements. Consequently, the vertical transport due to the turbulent flux plays a major role in the daytime evolution of both scalars and the advection is a relatively small contribution. During the analyzed days with a strong contribution of advection of either heat or carbon dioxide, the flux divergence is still an important contribution to the budget. For heat, the quantification of the advection contribution is in close agreement with results from a numerical model. For carbon dioxide, we qualitatively corroborate the results with a Lagrangian transport model. Our estimation of advection is compared with, traditional estimations based on the Net Ecosystem-atmosphere Exchange (NEE). Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
USGS2009

Numerical simulation of dune-flat bed transition and stage‐discharge relationship with hysteresis effect

This work presents recent advances on morphodynamic modeling of bed forms under unsteady discharge. This paper includes further development of a morphodynamic model proposed earlier by Giri and Shimizu (2006a). This model reproduces the temporal development of river dunes and accurately replicates the physical properties associated with bed form evolution. Model results appear to provide accurate predictions of bed form geometry and form drag over bed forms for arbitrary steady flows. However, accurate predictions of temporal changes of form drag are key to the prediction of stage‐discharge relation during flood events. Herein, the model capability is extended to replicate the dune–flat bed transition, and in turn, the variation of form drag produced by the temporal growth or decay of bed forms under unsteady flow conditions. Some numerical experiments are performed to analyze hysteresis of the stage‐discharge relationship caused by the transition between dune and flat bed regimes during rising and falling stages of varying flows. The numerical model successfully simulates dune–flat bed transition and the associated hysteresis of the stage‐discharge relationship; this is in good agreement with physical observations but has been treated in the past only using empirical methods. A hypothetical relationship for a sediment parameter (the mean step length) is proposed to a first level of approximation that enables reproduction of the dune–flat bed transition. The proposed numerical model demonstrates its ability to address an important practical problem associated with bed form evolution and flow resistance in varying flows.

Water Resources Research↗
USGS2019

Mechanisms of methane hydrate formation in geological systems

Natural gas hydrates are ice-like mixtures of gas (mostly methane) and water that are widely found in sediments along the world’s continental margins and within and beneath permafrost in a near-surface depth interval where the pressure is sufficiently high and temperature sufficiently low for gas hydrate to be stable. Beneath this interval, gas hydrate is not stable and free gas may be present. This paper reviews the multiple quantitative models that have proposed to describe the genesis of gas hydrate in geological systems. We emphasize the importance of coupling multi-phase flow (vapor and liquid) and multicomponent reactive transport with geological history to describe the dynamical processes of gas hydrate formation and evolution in geological systems. By understanding the generation and evolution of gas hydrate through time, we will better understand their role in the carbon cycle, their potential to contribute to climate change and geohazards, and how to design optimal strategies for the environmentally safe production of gas from hydrate reservoirs.

Reviews of Geophysics↗