SEARCH · Geology Reports

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

USGS2023

Intramolecular carbon isotope geochemistry of butane isomers from laboratory maturation and Monte-Carlo simulations of kerogen types I, II, and III

Position-specific (PS) carbon isotope compositions of light hydrocarbons such as propane and butane isomers (n-butane and i-butane) can provide a wealth of information on the history of natural gases in the subsurface reservoirs and other environments. For PS carbon isotope analysis of butane isomers, we have established a GC-pyrolysis-GC-isotope ratio mass spectrometry method with demonstrated accuracy. With this method, we analyzed PS δ 13 C values of butane isomers generated from the systematic laboratory pyrolysis experiments of three different kerogen types (I, II, and III) at temperatures of 310–430 °C with corresponding thermal maturity (Easy %R o ) ranging from 0.7 to 3.3. The observed evolution in the abundances of butane isomers can be interpreted and semi-quantitatively modeled based on the abundances of different C C bonds within the kerogens at low maturity and thermal degradation of butane isomers at high maturity. The δ 13 C values at the central sites of both nC 4 and iC 4 were heavier than those at the terminal positions, similar to our previous observations of propane. Their isotopic evolution with the maturity were controlled largely by kinetic isotope effects associated with breaking of different C C bonds during the generation and degradation of butane isomers. Kinetic Monte Carlo (kMC) simulations of n-butane generated from thermal cracking of model kerogens (I, II, and III) and an oil with a series of reactions (homolytic cleavage, β-scission, radical isomerization, H-abstraction, and termination by radical recombination) provided generally consistent results with the experimental observations, although the difference in PS δ 13 C values between the central and terminal positions are somewhat overestimated. On the other hand, the kMC simulation with homolytic cleavage and capping reactions alone produced significant deviations from the experimental results. Re-assessment of very limited data of PS δ 13 C values of natural butanes with our experimental and simulation results show that biodegradation significantly increased δ 13 C values at the central positions, not only of propane, but also of both butane isomers. This study lays a foundation and demonstrates the potential of PS isotope geochemistry of butane isomers to further improve our understanding of the sources, and geochemical and microbial processes of light hydrocarbons in the subsurface and other natural environments.

Geochimica et Cosmochimica Acta↗

USGS2013

Late Quaternary stream piracy and strath terrace formation along the Belle Fourche and lower Cheyenne Rivers, South Dakota and Wyoming

Stream piracy substantially affected the geomorphic evolution of the Missouri River watershed and drainages within, including the Little Missouri, Cheyenne, Belle Fourche, Bad, and White Rivers. The ancestral Cheyenne River eroded headward in an annular pattern around the eastern and southern Black Hills and pirated the headwaters of the ancestral Bad and White Rivers after ~ 660 ka. The headwaters of the ancestral Little Missouri River were pirated by the ancestral Belle Fourche River, a tributary to the Cheyenne River that currently drains much of the northern Black Hills. Optically stimulated luminescence (OSL) dating techniques were used to estimate the timing of this piracy event at ~ 22–21 ka. The geomorphic evolution of the Cheyenne and Belle Fourche Rivers is also expressed by regionally recognized strath terraces that include (from oldest to youngest) the Sturgis, Bear Butte, and Farmingdale terraces. Radiocarbon and OSL dates from fluvial deposits on these terraces indicate incision to the level of the Bear Butte terrace by ~ 63 ka, incision to the level of the Farmingdale terrace at ~ 40 ka, and incision to the level of the modern channel after ~ 12–9 ka. Similar dates of terrace incision have been reported for the Laramie and Wind River Ranges. Hypothesized causes of incision are the onset of colder climate during the middle Wisconsinan and the transition to the full-glacial climate of the late-Wisconsinan/Pinedale glaciation. Incision during the Holocene of the lower Cheyenne River is as much as ~ 80 m and is 3 to 4 times the magnitude of incision at ~ 63 ka and ~ 40 ka. The magnitude of incision during the Holocene might be due to a combined effect of three geomorphic processes acting in concert: glacial isostatic rebound in lower reaches (~ 40 m), a change from glacial to interglacial climate, and adjustments to increased watershed area resulting from piracy of the ancestral headwaters of the Little Missouri River.

Montana, Nebraska, North Dakota, South Dakota, Wyo↗

USGS2017

Excursions in fluvial (dis)continuity

Lurking below the twin concepts of connectivity and disconnectivity are their first, and in some ways, richer cousins: continuity and discontinuity. In this paper we explore how continuity and discontinuity represent fundamental and complementary perspectives in fluvial geomorphology, and how these perspectives inform and underlie our conceptions of connectivity in landscapes and rivers. We examine the historical roots of continuum and discontinuum thinking, and how much of our understanding of geomorphology rests on contrasting views of continuity and discontinuity. By continuum thinking we refer to a conception of geomorphic processes as well as geomorphic features that are expressed along continuous gradients without abrupt changes, transitions, or thresholds. Balance of forces, graded streams, and hydraulic geometry are all examples of this perspective. The continuum view has played a prominent role in diverse disciplinary fields, including ecology, paleontology, and evolutionary biology, in large part because it allows us to treat complex phenomena as orderly progressions and invoke or assume equilibrium processes that introduce order and prediction into our sciences. In contrast the discontinuous view is a distinct though complementary conceptual framework that incorporates non-uniform, non-progressive, and non-equilibrium thinking into understanding geomorphic processes and landscapes. We distinguish and discuss examples of three different ways in which discontinuous thinking can be expressed: 1) discontinuous spatial arrangements or singular events; 2) specific process domains generally associated with thresholds, either intrinsic or extrinsic; and 3) physical dynamics or changes in state, again often threshold-linked. In moving beyond the continuous perspective, a fertile set of ideas comes into focus: thresholds, non-equilibrium states, heterogeneity, catastrophe. The range of phenomena that is thereby opened up to scientific exploration similarly expands: punctuated episodes of cutting and filling, discretization of landscapes into hierarchies of structure and control, the work of extreme events. Orderly and progressive evolution towards a steady or ideal state is replaced by chaotic episodes of disturbance and recovery. Recent developments in the field of geomorphology suggest that we may be on the cusp of a new paradigm that recognizes that both continuous and discontinuous processes and mechanisms play a role in fluvial processes and landscape evolution with neither holding sway over the other and both needed to see rivers as they are.

Geomorphology↗

USGS2020

Modeling the spatial dynamics of marsh ponds in New England salt marshes

Ponds are common features on salt marshes, yet it is unclear how they affect large-scale marsh evolution. We developed a spatially explicit model that combines cellular automata for pond formation, expansion, and drainage, and partial differential equations for elevation dynamics. We use the mesotidal Barnstable marsh (MA, USA) as a case study, for which we measured pond expansion rate by remote sensing analysis over a 41-year time span. We estimated pond formation rate by comparing observed and modeled pond size distribution, and predicted pond deepening by comparing modeled and measured pond depth. The Barnstable marsh is currently in the pond recovery regime, i.e.,every pond revegetates and recovers the necessary elevation to support plant growth after re-connecting to the channel network. This pond dynamic creates an equivalent (i.e.,spatially and temporally averaged over the whole marsh) 0.5–2 mm/yr elevation loss that needs to be supplemented by excess vertical accretion. We explore how the pond regime would change with decreased sediment supply and increased relative sea-level rise (RSLR) rate, focusing on the case in which the vegetated marsh keeps pace with RSLR. When the RSLR rate remains below the minimum unvegetated deposition rate, the pond dynamics is nearly unaltered and ponds always occupy ~10% of the marsh area. However, when RSLR rate exceeds this threshold, the ponds in the marsh interior – which receive the least amount of suspended sediment – do not recover after drainage. These ponds transition to mudflats and permanently occupy up to 30% of the marsh area depending on RSLR rate. For marshes with a small tidal range, such as the microtidal Sage Lot Pond marsh on the opposite side of the peninsula from Barnstable marsh, high RSLR rates could bring every portion of the marsh into the pond runaway regime, with the whole marsh eventually converting into mudflats. In this regime, the existing marsh would disappear within centuries to millennia depending on the RSLR rate. Because of the spatial and temporal components of marsh evolution, a single RSLR threshold value applied across the entire marsh landscape provides a limited description of the marsh vulnerability to RSLR.

Massachusetts↗

USGS2011

Titan's cloud seasonal activity from winter to spring with Cassini/VIMS

Since Saturn orbital insertion in July 2004, the Cassini orbiter has been observing Titan throughout most of the northern winter season (October 2002–August 2009) and the beginning of spring, allowing a detailed monitoring of Titan’s cloud coverage at high spatial resolution with close flybys on a monthly basis. This study reports on the analysis of all the near-infrared images of Titan’s clouds acquired by the Visual and Infrared Mapping Spectrometer (VIMS) during 67 targeted flybys of Titan between July 2004 and April 2010. The VIMS observations show numerous sporadic clouds at southern high and mid-latitudes, rare clouds in the equatorial region, and reveal a long-lived cloud cap above the north pole, ubiquitous poleward of 60°N. These observations allow us to follow the evolution of the cloud coverage during almost a 6-year period including the equinox, and greatly help to further constrain global circulation models (GCMs). After 4 years of regular outbursts observed by Cassini between 2004 and 2008, southern polar cloud activity started declining, and completely ceased 1 year before spring equinox. The extensive cloud system over the north pole, stable between 2004 and 2008, progressively fractionated and vanished as Titan entered into northern spring. At southern mid-latitudes, clouds were continuously observed throughout the VIMS observing period, even after equinox, in a latitude band between 30°S and 60°S. During the whole period of observation, only a dozen clouds were observed closer to the equator, though they were slightly more frequent as equinox approached. We also investigated the distribution of clouds with longitude. We found that southern polar clouds, before disappearing in mid-2008, were systematically concentrated in the leading hemisphere of Titan, in particular above and to the east of Ontario Lacus, the largest reservoir of hydrocarbons in the area. Clouds are also non-homogeneously distributed with longitude at southern mid-latitudes. The n = 2-mode wave pattern of the distribution, observed since 2003 by Earth-based telescopes and confirmed by our Cassini observations, may be attributed to Saturn’s tides. Although the latitudinal distribution of clouds is now relatively well reproduced and understood by the GCMs, the non-homogeneous longitudinal distributions and the evolution of the cloud coverage with seasons still need investigation. If the observation of a few single clouds at the tropics and at northern mid-latitudes late in winter and at the start of spring cannot be further interpreted for the moment, the obvious shutdown of the cloud activity at Titan’s poles provides clear signs of the onset of the general circulation turnover that is expected to accompany the beginning of Titan’s northern spring. According to our GCM, the persistence of clouds at certain latitudes rather suggests a ‘sudden’ shift in near future of the meteorology into the more illuminated hemisphere. Finally, the observed seasonal change in cloud activity occurred with a significant time lag that is not predicted by our model. This may be due to an overall methane humidity at Titan’s surface higher than previously expected.

Icarus↗

USGS2012

Changes in sources and storage in a karst aquifer during a transition from drought to wet conditions

Understanding the sources and processes that control groundwater compositions and the timing and magnitude of groundwater vulnerability to potential surface-water contamination under varying meteorologic conditions is critical to informing groundwater protection policies and practices. This is especially true in karst terrains, where infiltrating surface water can rapidly affect groundwater quality. We analyzed the evolution of groundwater compositions (major ions and Sr isotopes) during the transition from extreme drought to wetconditions, and used inverse geochemical modeling (PHREEQC) to constrain controls on groundwater compositions during this evolution. Spring water and groundwater from two wells dominantly receiving diffuse and conduit flow (termed diffuse site and conduit site, respectively) in the Barton Springs segment of the Edwards aquifer (central Texas, USA) and surface water from losing streams that recharge the aquifer were sampled every 3–4 weeks during November 2008–March 2010. During this period, water compositions at the spring and conduit sites changed rapidly but there was no change at the diffuse site, illustrating the dual nature (i.e., diffuse vs. conduit) of flow in this karst system. Geochemical modeling demonstrated that, within a month of the onset of wetconditions, the majority of spring water and groundwater at the conduit site was composed of surface water, providing quantitative information on the timing and magnitude of the vulnerability of groundwater to potential surface-water contamination. The temporal pattern of increasing spring discharge and changing pattern of covariation between spring discharge and surface-water (steam) recharge indicates that that there were two modes of aquifer response—one with a small amount of storage and a second that accommodates more storage.

Texas↗

USGS2009

A lacustrine record from Lop Nur, Xinjiang, China: Implications for paleoclimate change during Late Pleistocene

Climate variability during the Late Pleistocene is studied from the proxies in core CK-2 drilled from the Luobei Depression (91??03???E, 40??47???N), Lop Nur in the eastern Tarim Basin, Xinjiang, China. Geophysical and geochemical properties, including magnetic susceptibility, granularity, chroma, carbonate content, loss on ignition and trace elements, have been determined to reconstruct the environmental evolution of the area during 32-9 ka BP. The chronology is established by uranium-thorium disequilibrium dating techniques. Our data suggest four paleoclimate stages, indicating glacial variations between cold-humid and warm-arid environments. A period of extreme humidity occurred during 31,900-19,200 yr BP is attributed the last glacial maximum (LGM). The period was followed by a warm-arid episode during 19,200-13,500 yr BP. Then a cold-humid interval during 13,500-12,700 yr BP may correspond to another cooling phases at high latitudes of the Northern Hemisphere. The last stage from 12,700 to 9000 yr BP has a trend that the climate turned warm and arid. The Lop Nur region is characterized by particularly humid stadials and arid interstadials. The climate variability in Lop Nur was constrained by global climate change because it is correlated with Dansgaard-Oeschger and Heinrich events, which were observed at the northern high latitudes. The synchroneity of the palaeoclimatic events suggested that cold air activity at the northern high latitudes was the most important factor that influenced the climate evolution in the Lop Nur region. A probable mechanism that involves the migration of westerly winds is proposed to interpret this synchroneity. ?? 2008 Elsevier Ltd.

Journal of Asian Earth Sciences↗

USGS2010

The thermal signature of volcanic eruptions on Io and Earth,

We investigate a spectrum-based technique to identify the style of active volcanic eruptions on Jupiter's moon Io. Thermal remote sensing of Io has had to rely primarily on low-spatial-resolution data, similar to low-spatial-resolution satellite data applied to detecting and charting the temporal evolution of terrestrial hot spots. These terrestrial analyses use data from sensors designed to monitor the weather and sea surface temperature. On Io, such low-spatial-resolution data are used to classify eruption styles (modes of emplacement) by means of several criteria related to the temporal evolution of the infrared spectrum associated with the eruptive activity at each hot spot, which we term “thermal signature.” We find that the ratio of the emission at 2 and 5 µm, and how this ratio changes with time, is often diagnostic of eruption style, even in low-spatial-resolution data. Tests using thermal data for terrestrial “ground truth” cases show that our classification system is valid on Earth. The results of our analysis can be used to aid in the design of future space-based instruments that can be used for volcano monitoring on Io, as well as Earth.

Journal of Volcanology and Geothermal Research↗

USGS2014

Brittle deformation and slope failure at the North Menan Butte tuff cone, Eastern Snake River Plain, Idaho

The manifestation of brittle deformation within inactive slumps along the North Menan Butte, a basaltic tuff cone in the Eastern Snake River Plain, is investigated through field and laboratory studies. Microstructural observations indicate that brittle strain is localized along deformation bands, a class of structural discontinuity that is predominant within moderate to high-porosity, clastic sedimentary rocks. Various subtypes of deformation bands are recognized in the study area based on the sense of strain they accommodate. These include dilation bands (no shear displacement), dilational shear bands, compactional shear bands and simple shear bands (no volume change). Measurements of the host rock permeability between the deformation bands indicate that the amount of brittle strain distributed throughout this part of the rock is negligible, and thus deformation bands are the primary means by which brittle strain is manifest within this tuff. Structural discontinuities that are similar in appearance to deformation bands are observed in other basaltic tuffs. Therefore deformation bands may represent a common structural feature of basaltic tuffs that have been widely misclassified as fractures. Slumping and collapse along the flanks of active volcanoes strongly influence their eruptive behavior and structural evolution. Therefore characterizing the process of deformation band and fault growth within basaltic tuff is key to achieving a more complete understanding of the evolution of basaltic volcanoes and their associated hazards.

Idaho↗

USGS2022

Olivine and glass chemistry record cycles of plumbing system recovery after summit collapse events at Kīlauea Volcano, Hawai‘i

The eruptive activity of Kīlauea Volcano (Hawai‘i) in the past 2500 years has alternated between centuries-long periods dominated either by explosive or effusive eruptions. The onset of explosive periods appears to be marked by caldera collapse events at the volcano's summit accompanied by draining of Kīlauea's magmatic plumbing system. Here we leverage >1800 olivine forsterite (Fo) contents, >900 glass MgO contents, and estimated magma supply rates from the past six centuries to describe the relationships between summit collapse and the composition of erupted material. On a first order basis, the major element chemistry of the centuries-long eruptive periods largely originates from fundamental differences between fractional crystallization of shallowly stored magmas during high-supply effusive-dominated periods versus little evolution of mafic recharge magmas during low-supply explosive-dominated periods. The modern effusive period (1820s-present) is dominated by relatively evolved olivine forsterite contents (Fo81–82) for Kīlauea, which is interpreted to reflect a buffered crustal reservoir system in which shallow storage and fractional crystallization control the composition of magmas. In contrast, olivine crystals from the explosive Keanakāko‘i Tephra (1500 - early 1800s C.E.) are dominated by higher olivine forsterite contents (Fo89) which are interpreted to reflect more primitive compositions, are correlated with glass MgO compositions extending to high values (e.g.,11.0 wt%), and show signs of magma mixing (zoned olivine, bimodal Fo populations). These signatures reflect a disrupted reservoir system in which high-MgO recharge melts mix with melts left over from draining of the shallow (<5 km) magma plumbing. Superimposed on these explosive-effusive periods are three decades- to centuries long periods of progressively evolving olivine and glass compositions. Eruptions that occur after caldera collapse in ~1500C.E. and smaller scale crater collapse events in 1790 (inferred) and 1924 have heterogeneous olivine populations dominated by ≥Fo88 and typically high MgO glasses. These compositions reflect inefficient mixing of stored and primitive recharge magmas after the disruption of the shallow plumbing system. After these collapses, olivine Fo and glass MgO subsequently evolve to <Fo82 and <7.0 wt% compositions, reflecting the recovery of the crustal plumbing system to an end-member system state characterized by efficient mixing of recharge and stored magmas that serve to buffer the shallow magma reservoirs. These evolved signatures suggest that a mature and buffered reservoir system may be a key condition for significant disruptions of volcanic plumbing systems. Plumbing system recovery is slower following large-scale caldera collapse (hundreds of years) compared to recovery following smaller crater collapse (tens of years), which may be modulated by differences in magma supply rates. Following the 2018 crater collapse olivine populations have high-Fo but glasses are low MgO, suggesting that this collapse might have disrupted shallow magma pathways but not strongly impacted the reservoir(s). Ultimately, olivine and glass major element chemistry record the impacts of caldera and smaller but significant summit crater collapse events at Kīlauea and could be used to provide a framework for better characterizing long-term volcano evolution in Hawai‘i and shield volcanoes elsewhere.

Hawaii↗

USGS2019

Morphodynamics of a field of crescent-shaped rippled scour depressions: Northern Monterey Bay, CA

Despite the prevalence of rippled scour depression (RSD) on the world's continental shelves and their importance as nursery habitats for many commercially-important species, the processes responsible for their formation and geomorphic evolution are still not well understood. Most studies that focused on RSD evolution have been based on data acquired over multiple years to decades, and often during calmer summer months. Here, multiple geophysical seafloor mapping surveys of a field of RSD's off Santa Cruz, California, USA, were conducted over the course of two winters with distinctly different oceanographic forcing; these results were compared to time series oceanographic measurements made within the survey area, and put into the context of long-term trends from survey data collected during previous decades. Although the migration of these features across the seafloor occurs on decadal time scales, significant change was detected within a given year, depending on the wave climate and fluvial sediment input. RSDs shrank up to 16% or grew up to 28% in a span of a few months. Such change was predominantly the result of redistribution of fine-grained sediment input to the system by local rivers. Migration of the seafloor features appears to be the result of accumulation of fine-grained sediment at the boundary between the RSD and the surrounding fine-grained cap.

California↗

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↗