Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research - Planets”

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 55 records · Page 3Linked to original sources

Ice core evidence of rapid air temperature increases since 1960 in alpine areas of the Wind River Range, Wyoming, United States

Site-specific transfer functions relating delta oxygen 18 (δ 18 O) values in snow to the average air temperature ( T A ) during storms on Upper Fremont Glacier (UFG) were used in conjunction with δ 18 O records from UFG ice cores to reconstruct long-term trends in air temperature from alpine areas in the Wind River Range, Wyoming. Transfer functions were determined by using data collected from four seasonal snowpacks (1989–1990, 1997–1998, 1998–1999, and 1999–2000). The timing and amount of each storm was determined from an automated snowpack telemetry (SNOTEL) site, 22 km northeast of UFG, and ∼1060 m in elevation below UFG. Statistically significant and positive correlations between δ 18 O values in the snow and T A were consistently found in three of the four seasonal snowpacks. The snowpack with the poor correlation was deposited in 1997–1998 during the 1997–1998 El Niño Southern Oscillation (ENSO). An ultrasonic snow-depth sensor installed on UFG provided valuable insights into site-specific storms and postdepositional processes that occur on UFG. The timing of storms recorded at the UFG and Cold Springs SNOTEL sites were similar; however, selected storms did not correlate. Snow from storms occurring after mid-October and followed by high winds was most susceptible to redeposition of snow. This removal of lower temperature snowfall could potentially bias the δ 18 O values preserved in ice core records to environmental conditions reflecting higher air temperatures and lower wind speeds. Transfer functions derived from seasonal snow cover on UFG were used to reconstruct T A values from δ 18 O values determined from two ice cores collected from UFG. Reconstructed air temperatures from the ice core data indicate an increase in T A of ∼3.5°C from the mid-1960s to the early 1990s in the alpine areas of northwestern Wyoming. Reconstructed T A from the ice core records between the end of the Little Ice Age (LIA), mid-1800s, and the early 1990s indicate a T A increase of ∼5°C. The historically reconstructed T A values from the UFG were significantly higher than the global average observed during the 20th Century but were in agreement with T A increases observed at selected, high-altitude and high-latitude sites in other parts of the world. Additional research is required to determine if part of the observed trend toward heavier δ 18 O values in ice from the UFG since the LIA (and increased T A ) is due to an increased proportion of snowfall from southerly storm tracks and moisture sources, as seems to have been the situation in 1997–1998.

Wyoming↗

Sediment and nutrient delivery from thermokarst features in the foothills of the North Slope, Alaska: Potential impacts on headwater stream ecosystems

Permafrost is a defining characteristic of the Arctic environment. However, climate warming is thawing permafrost in many areas leading to failures in soil structure called thermokarst. An extensive survey of a 600 km2 area in and around the Toolik Lake Natural Research Area (TLNRA) revealed at least 34 thermokarst features, two thirds of which were new since ???1980 when a high resolution aerial survey of the area was done. Most of these thermokarst features were associated with headwater streams or lakes. We have measured significantly increased sediment and nutrient loading from thermokarst features to streams in two well-studied locations near the TLNRA. One small thermokarst gully that formed in 2003 on the Toolik River in a 0.9 km2 subcatchment delivered more sediment to the river than is normally delivered in 18 years from 132 km2 in the adjacent upper Kuparuk River basin (a long-term monitoring reference site). Ammonium, nitrate, and phosphate concentrations downstream from a thermokarst feature on Imnavait Creek increased significantly compared to upstream reference concentrations and the increased concentrations persisted over the period of sampling (1999-2005). The downstream concentrations were similar to those we have used in a long-term experimental manipulation of the Kuparuk River and that have significantly altered the structure and function of that river. A subsampling of other thermokarst features from the extensive regional survey showed that concentrations of ammonium, nitrate, and phosphate were always higher downstream of the thermokarst features. Our previous research has shown that even minor increases in nutrient loading stimulate primary and secondary production. However, increased sediment loading could interfere with benthic communities and change the responses to increased nutrient delivery. Although the terrestrial area impacted by thermokarsts is limited, the aquatic habitat altered by these failures can be extensive. If warming in the Arctic foothills accelerates thermokarst formation, there may be substantial and wide-spread impacts on arctic stream ecosystems that are currently poorly understood. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research G: Biogeosciences↗

Understanding the amplitudes of noise correlation measurements

Cross correlation of ambient seismic noise is known to result in time series from which station-station travel-time measurements can be made. Part of the reason that these cross-correlation travel-time measurements are reliable is that there exists a theoretical framework that quantifies how these travel times depend on the features of the ambient noise. However, corresponding theoretical results do not currently exist to describe how the amplitudes of the cross correlation depend on such features. For example, currently it is not possible to take a given distribution of noise sources and calculate the cross correlation amplitudes one would expect from such a distribution. Here, we provide a ray-theoretical framework for calculating cross correlations. This framework differs from previous work in that it explicitly accounts for attenuation as well as the spatial distribution of sources and therefore can address the issue of quantifying amplitudes in noise correlation measurements. After introducing the general framework, we apply it to two specific problems. First, we show that we can quantify the amplitudes of coherency measurements, and find that the decay of coherency with station-station spacing depends crucially on the distribution of noise sources. We suggest that researchers interested in performing attenuation measurements from noise coherency should first determine how the dominant sources of noise are distributed. Second, we show that we can quantify the signal-to-noise ratio of noise correlations more precisely than previous work, and that these signal-to-noise ratios can be estimated for given situations prior to the deployment of seismometers. It is expected that there are applications of the theoretical framework beyond the two specific cases considered, but these applications await future work.

Journal of Geophysical Research↗

Surface temperature patterns in complex terrain: Daily variations and long-term change in the central Sierra Nevada, California

A realistic description of how temperatures vary with elevation is crucial for ecosystem studies and for models of basin-scale snowmelt and spring streamflow. This paper explores surface temperature variability using temperature data from an array of 37 sensors, called the Yosemite network, which traverses both slopes of the Sierra Nevada in the vicinity of Yosemite National Park, California. These data indicate that a simple lapse rate is often a poor description of the spatial temperature structure. Rather, the spatial pattern of temperature over the Yosemite network varies considerably with synoptic conditions. Empirical orthogonal functions (EOFs) were used to identify the dominant spatial temperature patterns and how they vary in time. Temporal variations of these surface temperature patterns were correlated with large-scale weather conditions, as described by National Centers for Environmental Prediction-National Center for Atmospheric Research Reanalysis data. Regression equations were used to downscale larger-scale weather parameters, such as Reanalysis winds and pressure, to the surface temperature structure over the Yosemite network. These relationships demonstrate that strong westerly winds are associated with relatively warmer temperatures on the east slope and cooler temperatures on the west slope of the Sierra, and weaker westerly winds are associated with the opposite pattern. Reanalysis data from 1948 to 2005 indicate weakening westerlies over this time period, a trend leading to relatively cooler temperatures on the east slope over decadal timescale's. This trend also appears in long-term observations and demonstrates the need to consider topographic effects when examining long-term changes in mountain regions. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

Introduction to special section on the California-Arizona crustal transect: Part II

In May 1988, a multidisciplinary conference was held in Flagstaff, Arizona, to synthesize recent research on the crustal structural and evolution of a transect across the southern Cordillera. This transect extends from the southern Colorado Plateau southwestward across the Arizona Transition Zone, through the Basin and Range province in the Mojave and Sonoran deserts, to the San Andreas fault system, and thence westward across the Imperial Valley and Peninsular Ranges and into the southern California Borderland. Most discussion concerned the eastern part of the transect, from the Colorado Plateau to the San Andreas fault [ Sass et al ., 1988; Howard et al ., 1990]. For this meeting, the name CACTIS (California‐Arizona Crustal Transect‐Interim Synthesis) was coined by J. Sass.

Journal of Geophysical Research↗

Integrating perspectives to understand lake ice dynamics in a changing world

Ice cover plays a critical role in physical, biogeochemical, and ecological processes in lakes. Despite its importance, winter limnology remains relatively understudied. Here, we provide a primer on the predominant drivers of freshwater lake ice cover and the current methodologies used to study lake ice, including in situ and remote sensing observations, physical based models, and experiments. We highlight opportunities for future research by integrating these four disciplines to address key knowledge gaps in our understanding of lake ice dynamics in changing winters. Advances in technology, data integration, and interdisciplinary collaboration will allow the field to move toward developing global forecasts of lake ice cover for small to large lakes across broad spatial and temporal scales, quantifying ice quality and ice thickness, moving from binary to continuous ice records, and determining how winter ice conditions and quality impact ecosystem processes in lakes over winter. Ultimately, integrating disciplines will improve our ability to understand the impacts of changing winters on lake ice.

Journal of Geophysical Research: Biogeosciences↗

Coseismic deformation of the Wrights tunnel during the 1906 San Francisco earthquake: A key to understanding 1906 fault slip and 1989 surface ruptures in the southern Santa Cruz Mountains, California

The Wrights tunnel is an abandoned railroad tunnel that crosses the San Andreas fault in the southern Santa Cruz Mountains in the vicinity of the 1989 Loma Prieta earthquake. The tunnel was damaged and deformed during the 1906 San Francisco earthquake and a plot showing postearthquake measurements made in the tunnel is given by Lawson [1908]. The amount of offset shown on this plot (1.5 m) has been used in several studies as being representative of the amount of fault offset along this segment of the San Andreas fault in 1906. Our historical research shows that different observers reported different amounts of fault offset in the tunnel and that the 1.5 m given on the plot is not a surveyed measurement. In addition, the plot of the tunnel has been interpreted in several previous studies as evidence of a broad (1.5 km) zone of faulting beneath Summit Ridge. Our analysis shows that this plot need not indicate a broad zone of deformation. Our historical research and modeling of the tunnel measurements indicate that faulting was confined to a zone less than 400 m wide and that 60–85% of the coseismic slip occurred across a single fault plane. There is no evidence for offset across a second shear zone beneath Summit Ridge in 1906. This implies that surface fractures reported on Summit Ridge in 1906 were not associated with significant deformation of the tunnel, implying that they were shallow, surficial features. By analogy, the very similar fractures that occurred on Summit Ridge in 1989 were also probably the result of shallow gravitational, rather than deep-seated tectonic, processes. Our modeling also indicates that total coseismic, near-surface slip across the San Andreas fault zone in the Wrights tunnel in 1906 was at least 1.7–1.8 m.

Journal of Geophysical Research B: Solid Earth↗

Arctic transitions in the Land - Atmosphere System (ATLAS): Background, objectives, results, and future directions

This paper briefly reviews the background, objectives, and results of the Arctic Transitions in the Land-Atmosphere System (ATLAS) Project to date and provides thoughts on future directions. The key goal of the ATLAS Project is to improve understanding of controls over spatial and temporal variability of terrestrial processes in the Arctic that have potential consequences for the climate system, i.e., processes that affect the exchange of water and energy with the atmosphere, the exchange of radiatively active gases with the atmosphere, and the delivery of freshwater to the Arctic Ocean. Three important conclusions have emerged from research associated with the ATLAS Project. First, associated with the observation that the Alaskan Arctic has warmed significantly in the last 30 years, permafrost is warming, shrubs are expanding, and there has been a temporary release of carbon dioxide from tundra soils. Second, the winter is a more important period of biological activity than previously appreciated. Biotic processes, including shrub expansion and decomposition, affect snow structure and accumulation and affect the annual carbon budget of tundra ecosystems. Third, observed vegetation changes can have a significant positive feedback to regional warming. These vegetation effects are, however, less strong than those exerted by land-ocean heating contrasts and the topographic constraints on air mass movements. The papers of this special section provide additional insights related to these conclusions and to the overall goal of ATLAS.

Journal of Geophysical Research D: Atmospheres↗

A detached eddy simulation model for the study of lateral separation zones along a large canyon-bound river

Lateral flow separation occurs in rivers where banks exhibit strong curvature. In canyon-boundrivers, lateral recirculation zones are the principal storage of fine-sediment deposits. A parallelized,three-dimensional, turbulence-resolving model was developed to study the flow structures along lateralseparation zones located in two pools along the Colorado River in Marble Canyon. The model employs thedetached eddy simulation (DES) technique, which resolves turbulence structures larger than the grid spacingin the interior of the flow. The DES-3D model is validated using Acoustic Doppler Current Profiler flowmeasurements taken during the 2008 controlled flood release from Glen Canyon Dam. A point-to-pointvalidation using a number of skill metrics, often employed in hydrological research, is proposed here forfluvial modeling. The validation results show predictive capabilities of the DES model. The model reproducesthe pattern and magnitude of the velocity in the lateral recirculation zone, including the size and position ofthe primary and secondary eddy cells, and return current. The lateral recirculation zone is open, havingcontinuous import of fluid upstream of the point of reattachment and export by the recirculation returncurrent downstream of the point of separation. Differences in magnitude and direction of near-bed andnear-surface velocity vectors are found, resulting in an inward vertical spiral. Interaction between therecirculation return current and the main flow is dynamic, with large temporal changes in flow direction andmagnitude. Turbulence structures with a predominately vertical axis of vorticity are observed in the shearlayer becoming three-dimensional without preferred orientation downstream.

Arizona↗

Chronological refinement of an ice core record at Upper Fremont Glacier in south central North America

The potential to use ice cores from alpine glaciers in the midlatitudes to reconstruct paleoclimatic records has not been widely recognized. Although excellent paleoclimatic records exist for the polar regions, paleoclimatic ice core records are not common from midlatitude locations. An ice core removed from the Upper Fremont Glacier in Wyoming provides evidence for abrupt climate change during the mid-1800s. Volcanic events (Krakatau and Tambora) identified from electrical conductivity measurements (ECM) and isotopic and chemical data from the Upper Fremont Glacier were reexamined to confirm and refine previous chronological estimates of the ice core. At a depth of 152 m the refined age-depth profile shows good agreement (1736 ± 10 A.D.) with the 14C age date (1729 ± 95 A.D.). The δ 18 O profile of the Upper Fremont Glacier (UFG) ice core indicates a change in climate known as the Little Ice Age (LIA). However, the sampling interval for δ 18 O is sufficiently large (20 cm) such that it is difficult to pinpoint the LIA termination on the basis of δ 18 O data alone. Other research has shown that changes in the δ18O variance are generally coincident with changes in ECM variance. The ECM data set contains over 125,000 data points at a resolution of 1 data point per millimeter of ice core. A 999-point running average of the ECM data set and results from f tests indicates that the variance of the ECM data decreases significantly at about 108 m. At this depth, the age-depth profile predicts an age of 1845 A.D. Results indicate the termination of the LIA was abrupt with a major climatic shift to warmer temperatures around 1845 A.D. and continuing to present day. Prediction limits (error bars) calculated for the profile ages are ±10 years (90% confidence level). Thus a conservative estimate for the time taken to complete the LIA climatic shift to present-day climate is about 10 years, suggesting the LIA termination in alpine regions of central North America may have occurred on a relatively short (decadal) timescale. Copyright 2000 by the American Geophysical Union.

Wyoming↗

Volcano monitoring using the Global Positioning System: Filtering strategies

Permanent Global Positioning System (GPS) networks are routinely used for producing improved orbits and monitoring secular tectonic deformation. For these applications, data are transferred to an analysis center each day and routinely processed in 24-hour segments. To use GPS for monitoring volcanic events, which may last only a few hours, real-time or near real-time data processing and subdaily position estimates are valuable. Strategies have been researched for obtaining station coordinates every 15 min using a Kalman filter; these strategies have been tested on data collected by a GPS network on Kilauea Volcano. Data from this network are tracked continuously, recorded every 30 s, and telemetered hourly to the Hawaiian Volcano Observatory. A white noise model is heavily impacted by data outages and poor satellite geometry, but a properly constrained random walk model fits the data well. Using a borehole tiltmeter at Kilauea's summit as ground-truth, solutions using different random walk constraints were compared. This study indicates that signals on the order of 5 mm/h are resolvable using a random walk standard deviation of 0.45 cm/√h. Values lower than this suppress small signals, and values greater than this have significantly higher noise at periods of 1–6 hours.

Hawai'i↗

Archie’s saturation exponent for natural gas hydrate in coarse-grained reservoirs

Accurately quantifying the amount of naturally occurring gas hydrate in marine and permafrost environments is important for assessing its resource potential and understanding the role of gas hydrate in the global carbon cycle. Electrical resistivity well logs are often used to calculate gas hydrate saturations, S h , using Archie's equation. Archie's equation, in turn, relies on an empirical saturation parameter, n . Though n = 1.9 has been measured for ice‐bearing sands and is widely used within the hydrate community, it is highly questionable if this n value is appropriate for hydrate‐bearing sands. In this work, we calibrate n for hydrate‐bearing sands from the Canadian permafrost gas hydrate research well, Mallik 5L‐38, by establishing an independent downhole S h profile based on compressional‐wave velocity log data. Using the independently determined S h profile and colocated electrical resistivity and bulk density logs, Archie's saturation equation is solved for n, and uncertainty is tracked throughout the iterative process. In addition to the Mallik 5L‐38 well, we also apply this method to two marine, coarse‐grained reservoirs from the northern Gulf of Mexico Gas Hydrate Joint Industry Project: Walker Ridge 313‐H and Green Canyon 955‐H. All locations yield similar results, each suggesting n ≈ 2.5 ± 0.5. Thus, for the coarse‐grained hydrate bearing ( S h > 0.4) of greatest interest as potential energy resources, we suggest that n = 2.5 ± 0.5 should be applied in Archie's equation for either marine or permafrost gas hydrate settings if independent estimates of n are not available.

Journal of Geophysical Research B: Solid Earth↗

Testing the accuracy of a 1-D volcanic plume model in estimating mass eruption rate

During volcanic eruptions, empirical relationships are used to estimate mass eruption rate from plume height. Although simple, such relationships can be inaccurate and can underestimate rates in windy conditions. One-dimensional plume models can incorporate atmospheric conditions and give potentially more accurate estimates. Here I present a 1-D model for plumes in crosswind and simulate 25 historical eruptions where plume height H obs was well observed and mass eruption rate M obs could be calculated from mapped deposit mass and observed duration. The simulations considered wind, temperature, and phase changes of water. Atmospheric conditions were obtained from the National Center for Atmospheric Research Reanalysis 2.5° model. Simulations calculate the minimum, maximum, and average values ( M min , M max , and M avg ) that fit the plume height. Eruption rates were also estimated from the empirical formula M empir = 140 H obs 4.14 ( M empir is in kilogram per second, H obs is in kilometer). For these eruptions, the standard error of the residual in log space is about 0.53 for M avg and 0.50 for M empir . Thus, for this data set, the model is slightly less accurate at predicting M obs than the empirical curve. The inability of this model to improve eruption rate estimates may lie in the limited accuracy of even well-observed plume heights, inaccurate model formulation, or the fact that most eruptions examined were not highly influenced by wind. For the low, wind-blown plume of 14–18 April 2010 at Eyjafjallajökull, where an accurate plume height time series is available, modeled rates do agree better with M obs than M empir .

Journal of Geophysical Research D: Atmospheres↗

Association between atmospheric circulation patterns and firn-ice core records from the Inilchek glacierized area, central Tien Shan, Asia

Glacioclimatological research in the central Tien Shan was performed in the summers of 1998 and 1999 on the South Inilchek Glacier at 5100-5460 m. A 14.36 m firn-ice core and snow samples were collected and used for stratigraphic, isotopic, and chemical analyses. The firn-ice core and snow records were related to snow pit measurements at an event scale and to meteorological data and synoptic indices of atmospheric circulation at annual and seasonal scales. Linear relationships between the seasonal air temperature and seasonal isotopic composition in accumulated precipitation were established. Changes in the ??18O air temperature relationship, in major ion concentration and in the ratios between chemical species, were used to identify different sources of moisture and investigate changes in atmospheric circulation patterns. Precipitation over the central Tien Shan is characterized by the lowest ionic content among the Tien Shan glaciers and indicates its mainly marine origin. In seasons of minimum precipitation, autumn and winter, water vapor was derived from the and and semiarid regions in central Eurasia and contributed annual maximal solute content to snow accumulation in Tien Shan. The lowest content of major ions was observed in spring and summer layers, which represent maximum seasonal accumulation when moisture originates over the Atlantic Ocean and Mediterranean and Black Seas. Copyright 2004 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

Change in atmospheric mineral aerosols in response to climate: Last glacial period, preindustrial, modern, and doubled carbon dioxide climates

Desert dust simulations generated by the National Center for Atmospheric Research's Community Climate System Model for the current climate are shown to be consistent with present day satellite and deposition data. The response of the dust cycle to last glacial maximum, preindustrial, modern, and doubled-carbon dioxide climates is analyzed. Only natural (non-land use related) dust sources are included in this simulation. Similar to some previous studies, dust production mainly responds to changes in the source areas from vegetation changes, not from winds or soil moisture changes alone. This model simulates a +92%, +33%, and -60% change in dust loading for the last glacial maximum, preindustrial, and doubled-carbon dioxide climate, respectively, when impacts of carbon dioxide fertilization on vegetation are included in the model. Terrestrial sediment records from the last glacial maximum compiled here indicate a large underestimate of deposition in continental regions, probably due to the lack of simulation of glaciogenic dust sources. In order to include the glaciogenic dust sources as a first approximation, we designate the location of these sources, and infer the size of the sources using an inversion method that best matches the available data. The inclusion of these inferred glaciogenic dust sources increases our dust flux in the last glacial maximum from 2.1 to 3.3 times current deposition. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

Bora event variability and the role of air-sea feedback

A two-way interacting high resolution numerical simulation of the Adriatic Sea using the Navy Coastal Ocean Model (NCOM) and Coupled Ocean/ Atmosphere Mesoscale Prediction System (COAMPS??) was conducted to improve forecast momentum and heat flux fields, and to evaluate surface flux field differences for two consecutive bora events during February 2003. (COAMPS?? is a registered trademark of the Naval Research Laboratory.) The strength, mean positions and extensions of the bora jets, and the atmospheric conditions driving them varied considerably between the two events. Bora 1 had 62% stronger heat flux and 51% larger momentum flux than bora 2. The latter displayed much greater diurnal variability characterized by inertial oscillations and the early morning strengthening of a west Adriatic barrier jet, beneath which a stronger west Adriatic ocean current developed. Elsewhere, surface ocean current differences between the two events were directly related to differences in wind stress curl generated by the position and strength of the individual bora jets. The mean heat flux bias was reduced by 72%, and heat flux RMSE reduced by 30% on average at four instrumented over-water sites in the two-way coupled simulation relative to the uncoupled control. Largest reductions in wind stress were found in the bora jets, while the biggest reductions in heat flux were found along the north and west coasts of the Adriatic. In bora 2, SST gradients impacted the wind stress curl along the north and west coasts, and in bora 1 wind stress curl was sensitive to the Istrian front position and strength. The two-way coupled simulation produced diminished surface current speeds of ???12% over the northern Adriatic during both bora compared with a one-way coupled simulation. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗

GCIP water and energy budget synthesis (WEBS)

As part of the World Climate Research Program's (WCRPs) Global Energy and Water-Cycle Experiment (GEWEX) Continental-scale International Project (GCIP), a preliminary water and energy budget synthesis (WEBS) was developed for the period 1996-1999 fromthe "best available" observations and models. Besides this summary paper, a companion CD-ROM with more extensive discussion, figures, tables, and raw data is available to the interested researcher from the GEWEX project office, the GAPP project office, or the first author. An updated online version of the CD-ROM is also available at http://ecpc.ucsd.edu/gcip/webs.htm/. Observations cannot adequately characterize or "close" budgets since too many fundamental processes are missing. Models that properly represent the many complicated atmospheric and near-surface interactions are also required. This preliminary synthesis therefore included a representative global general circulation model, regional climate model, and a macroscale hydrologic model as well as a global reanalysis and a regional analysis. By the qualitative agreement among the models and available observations, it did appear that we now qualitatively understand water and energy budgets of the Mississippi River Basin. However, there is still much quantitative uncertainty. In that regard, there did appear to be a clear advantage to using a regional analysis over a global analysis or a regional simulation over a global simulation to describe the Mississippi River Basin water and energy budgets. There also appeared to be some advantage to using a macroscale hydrologic model for at least the surface water budgets. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗