SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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

USGS2015

Direct measurement of asperity contact growth in quartz at hydrothermal conditions

Earthquake recurrence requires interseismic fault restrengthening which results from solid state deformation in room-temperature friction and indentation experiments. In contrast exhumed fault zones show solution-transport processes such as pressure solution and contact overgrowths influence fault zone properties . In the absence of fluid flow, overgrowths are driven by gradients in surface curvature where material is dissolved, diffuses, and precipitates at the contact without convergence normal to the contact. To determine the rate of overgrowth for quartz, we conducted single contact experiments in an externally heated pressure vessel. Convergence was continuously monitored using reflected-light interferometry through a long-working-distance microscope. Contact normal force was constant with an initial effective normal stress of 1.7 MPa, temperature was between 350 and 530{degree sign}C, and water pressure was constant at 150 MPa. Two control experiments were conducted: one dry at 425{degree sign}C and one bi-material (sapphire) at 425{degree sign}C and 150 MPa water pressure. No contact growth or convergence was observed in the controls. For wet single-phase contacts, growth was initially rapid and then decreased with time. No convergence was observed. Fluid inclusions indicate that the contact is not uniformly wetted. The contact is bounded by small regions of high aperture, reflecting local free-face dissolution as the source for the overgrowth. The apparent activation energy is ~125 kJ/mol. Extrapolation predicts rates of contact area increase orders of magnitude faster than in dry, room-temperature and hydrothermal friction experiments, suggesting that natural strength recovery near the base of the seismogenic zone could be dominated by contact overgrowth.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

The vertical structure of the circulation and dynamics in Hudson Shelf Valley

Hudson Shelf Valley is a 20–30 m deep, 5–10 km wide v-shaped submarine valley that extends across the Middle Atlantic Bight continental shelf. The valley provides a conduit for cross-shelf exchange via along-valley currents of 0.5 m s −1 or more. Current profile, pressure, and density observations collected during the winter of 1999–2000 are used to examine the vertical structure and dynamics of the flow. Near-bottom along-valley currents having times scales of a few days are driven by cross-shelf pressure gradients setup by wind stresses, with eastward (westward) winds driving onshore (offshore) flow within the valley. The along-valley momentum balance in the bottom boundary layer is predominantly between the pressure gradient and bottom stress because the valley bathymetry limits current veering. Above the bottom boundary layer, the flow veers toward an along-shelf (cross-valley) orientation and a geostrophic balance with some contribution from the wind stress (surface Ekman layer). The vertical structure and strength of the along-valley current depends on the magnitude and direction of the wind stress. During offshore flows driven by westward winds, the near-bottom stratification within the valley increases resulting in a thinner bottom boundary layer and weaker offshore currents. Conversely, during onshore flows driven by eastward winds the near-bottom stratification decreases resulting in a thicker bottom boundary layer and stronger onshore currents. Consequently, for wind stress magnitudes exceeding 0.1 N m −2 , onshore along-valley transport associated with eastward wind stress exceeds the offshore transport associated with westward wind stress of the same magnitude.

Hudson Shelf Valley↗

USGS2014

Formation of fine sediment deposit from a flash flood river in the Mediterranean Sea

We identify the mechanisms controlling fine deposits on the inner-shelf in front of the Besòs River, in the northwestern Mediterranean Sea. This river is characterized by a flash flood regime discharging large amounts of water (more than 20 times the mean water discharge) and sediment in very short periods lasting from hours to few days. Numerical model output was compared with bottom sediment observations and used to characterize the multiple spatial and temporal scales involved in offshore sediment deposit formation. A high-resolution (50 m grid size) coupled hydrodynamic-wave-sediment transport model was applied to the initial stages of the sediment dispersal after a storm-related flood event. After the flood, sediment accumulation was predominantly confined to an area near the coastline as a result of preferential deposition during the final stage of the storm. Subsequent reworking occurred due to wave-induced bottom shear stress that resuspended fine materials, with seaward flow exporting them toward the midshelf. Wave characteristics, sediment availability, and shelf circulation determined the transport after the reworking and the final sediment deposition location. One year simulations of the regional area revealed a prevalent southwestward average flow with increased intensity downstream. The circulation pattern was consistent with the observed fine deposit depocenter being shifted southward from the river mouth. At the southern edge, bathymetry controlled the fine deposition by inducing near-bottom flow convergence enhancing bottom shear stress. According to the short-term and long-term analyses, a seasonal pattern in the fine deposit formation is expected.

Journal of Geophysical Research: Oceans↗

USGS2014

Small crater modification on Meridiani Planum and implications for erosion rates and climate change on Mars

A morphometric and morphologic catalog of ~100 small craters imaged by the Opportunity rover over the 33.5 km traverse between Eagle and Endeavour craters on Meridiani Planum shows craters in six stages of degradation that range from fresh and blocky to eroded and shallow depressions ringed by planed off rim blocks. The age of each morphologic class from <50–200 ka to ~20 Ma has been determined from the size-frequency distribution of craters in the catalog, the retention age of small craters on Meridiani Planum, and the age of the latest phase of ripple migration. The rate of degradation of the craters has been determined from crater depth, rim height, and ejecta removal over the class age. These rates show a rapid decrease from ~1 m/Myr for craters < 1 Ma to ~ <0.1 m/Myr for craters 10–20 Ma, which can be explained by topographic diffusion with modeled diffusivities of ~10 −6  m 2 /yr. In contrast to these relatively fast, short-term erosion rates, previously estimated average erosion rates on Mars over ~100 Myr and 3 Gyr timescales from the Amazonian and Hesperian are of order < 0.01 m/Myr, which is 3–4 orders of magnitude slower than typical terrestrial rates. Erosion rates during the Middle-Late Noachian averaged over ~250 Myr, and ~700 Myr intervals are around 1 m/Myr, comparable to slow terrestrial erosion rates calculated over similar timescales. This argues for a wet climate before ~3 Ga in which liquid water was the erosional agent, followed by a dry environment dominated by slow eolian erosion.

Journal of Geophysical Research E: Planets↗

USGS2015

Chemical variations in Yellowknife Bay formation sedimentary rocks analyzed by ChemCam on board the Curiosity rover on Mars

The Yellowknife Bay formation represents a ~5 m thick stratigraphic section of lithified fluvial and lacustrine sediments analyzed by the Curiosity rover in Gale crater, Mars. Previous works have mainly focused on the mudstones that were drilled by the rover at two locations. The present study focuses on the sedimentary rocks stratigraphically above the mudstones by studying their chemical variations in parallel with rock textures. Results show that differences in composition correlate with textures and both manifest subtle but significant variations through the stratigraphic column. Though the chemistry of the sediments does not vary much in the lower part of the stratigraphy, the variations in alkali elements indicate variations in the source material and/or physical sorting, as shown by the identification of alkali feldspars. The sandstones contain similar relative proportions of hydrogen to the mudstones below, suggesting the presence of hydrous minerals that may have contributed to their cementation. Slight variations in magnesium correlate with changes in textures suggesting that diagenesis through cementation and dissolution modified the initial rock composition and texture simultaneously. The upper part of the stratigraphy (~1 m thick) displays rocks with different compositions suggesting a strong change in the depositional system. The presence of float rocks with similar compositions found along the rover traverse suggests that some of these outcrops extend further away in the nearby hummocky plains.

Journal of Geophysical Research E: Planets↗

USGS2015

Mars Reconnaissance Orbiter and Opportunity observations of the Burns formation: crater hopping at Meridiani Planum

Compact Reconnaissance Imaging Spectrometer for Mars hyperspectral (1.0–2.65 µm) along-track oversampled observations covering Victoria, Santa Maria, Endeavour, and Ada craters were processed to 6 m/pixel and used in combination with Opportunity observations to detect and map hydrated Mg and Ca sulfate minerals in the Burns formation. The strongest spectral absorption features were found to be associated with outcrops that are relatively young and fresh (Ada) or preferentially scoured of dust, soil, and coatings by prevailing winds. At Victoria and Santa Maria, the scoured areas are on the southeastern rims and walls, opposite to the sides where wind-blown sands extend out of the craters. At Endeavour, the deepest absorptions are in Botany Bay, a subdued and buried rim segment that exhibits high thermal inertias, extensive outcrops, and is interpreted to be a region of enhanced wind scour extending up and out of the crater. Ada, Victoria, and Santa Maria outcrops expose the upper portion of the preserved Burns formation and show spectral evidence for the presence of kieserite. In contrast, gypsum is pervasive spectrally in the Botany Bay exposures. Gypsum, a relatively insoluble evaporative mineral, is interpreted to have formed close to the contact with the Noachian crust as rising groundwaters brought brines close to and onto the surface, either as a direct precipitate or during later diagenesis. The presence of kieserite at the top of the section is hypothesized to reflect precipitation from evaporatively concentrated brines or dehydration of polyhydrated sulfates.

Journal of Geophysical Research E: Planets↗

USGS2015

Context of ancient aqueous environments on Mars from in situ geologic mapping at Endeavour Crater

Using the Mars Exploration Rover Opportunity , we have compiled one of the first field geologic maps on Mars while traversing the Noachian terrain along the rim of the 22 km diameter Endeavour Crater (Latitude −2°16′33″, Longitude −5°10′51″). In situ mapping of the petrographic, elemental, structural, and stratigraphic characteristics of outcrops and rocks distinguishes four mappable bedrock lithologic units. Three of these rock units predate the surrounding Burns formation sulfate-rich sandstones and one, the Matijevic Formation, represents conditions on early Mars predating the formation of Endeavour Crater. The stratigraphy assembled from these observations includes several geologic unconformities. The differences in lithologic units across these unconformities record changes in the character and intensity of the Martian aqueous environment over geologic time. Water circulated through fractures in the oldest rocks over periods long enough that texturally and elementally significant alteration occurred in fracture walls. These oldest pre-Endeavour rocks and their network of mineralized and altered fractures were preserved by burial beneath impact ejecta and were subsequently exhumed and exposed. The alteration along joints in the oldest rocks and the mineralized veins and concentrations of trace metals in overlying lithologic units is direct evidence that copious volumes of mineralized and/or hydrothermal fluids circulated through the early Martian crust. The wide range in intensity of structural and chemical modification from outcrop to outcrop along the crater rim shows that the ejecta of large (>8 km in diameter) impact craters is complex. These results imply that geologic complexity is to be anticipated in other areas of Mars where cratering has been a fundamental process in the local and regional geology and mineralogy.

Journal of Geophysical Research E: Planets↗

USGS2015

Evolution of Mars’ Northern Polar Seasonal CO 2 deposits: variations in surface brightness and bulk density

Small scale variations of seasonal ice are explored at different geomorphic units on the Northern Polar Seasonal Cap (NPSC). We use seasonal rock shadow measurements, combined with visible and thermal observations, to calculate density over time. The coupling of volume density and albedo allows us to determine the microphysical state of the seasonal CO 2 ice. We find two distinct endmembers across the NPSC: 1) Snow deposits may anneal to form an overlying slab layer that fractures. These low density deposits maintain relatively constant densities over springtime. 2) Porous slab deposits likely anneal rapidly in early spring and fracture in late spring. These high density deposits dramatically increase in density over time. The endmembers appear to be correlated with latitude.

Journal of Geophysical Research E: Planets↗

USGS2015

Atypical pit craters on Mars: new insights from THEMIS, CTX and HiRISE observations

More than 100 pit craters in the Tharsis region of Mars exhibit morphologies, diameters and thermal behaviors that diverge from the much larger bowl-shaped pit craters that occur in most regions across Mars. These Atypical Pit Craters (APCs) generally have sharp and distinct rims, vertical or overhanging walls that extend down to their floors, surface diameters of ~50-350 m, and high depth-to-diameter (d/D) ratios that are usually greater than 0.3 (which is an upper-range value for impacts and bowl-shaped pit craters), and can exceed values of 1.8. Observations by the Mars Odyssey THermal Emission Imaging System (THEMIS) show that APC floor temperatures are warmer at night, and fluctuate with much lower diurnal amplitudes than nearby surfaces or adjacent bowl-shaped pit craters. Kīlauea volcano, Hawai'i, hosts pit craters that formed through subsurface collapse into active volcanic dikes, resulting in pits that can appear morphologically analogous to either APCs or bowl-shaped pit craters. Partially-drained dikes are sometimes exposed within the lower walls and floors of these terrestrial APC analogs and can form extensive cave systems with unique microclimates. Similar caves in martian pit craters are of great interest for astrobiology. This study uses new observations by the Mars Reconnaissance Orbiter (MRO) High Resolution Imaging Science Experiment (HiRISE) and Context Camera (CTX) to refine previous work where seven APCs were described from lower-resolution THEMIS visible-wavelength (VIS) observations. Here, we identify locations of 115 APCs, map their distribution across the Tharsis region, characterize their internal morphologies with high-resolution observations, and discuss possible formation mechanisms.

Journal of Geophysical Research↗

USGS2015

Examining spectral variations in localized lunar dark mantle deposits

The localized lunar dark mantle deposits (DMDs) in Alphonsus, J. Herschel, and Oppenheimer craters were analyzed using visible-near-infrared spectroscopy data from the Moon Mineralogy Mapper. Spectra of these localized DMDs were analyzed for compositional and mineralogical variations within the deposits and were compared with nearby mare basalt units. Spectra of the three localized DMDs exhibited mafic absorption features indicating iron-rich compositions, although the DMDs were spectrally distinct from nearby mare basalts. All of the DMDs contained spectral signatures of glassy materials, suggesting the presence of volcanic glass in varying concentrations across the individual deposits. In addition, the albedo and spectral signatures were variable within the Alphonsus and Oppenheimer crater DMDs, suggesting variable deposit thickness and/or variations in the amount of mixing with the local substrate. Two previously unidentified localized DMDs were discovered to the northeast of Oppenheimer crater. The identification of high concentrations of volcanic glass in multiple localized DMDs in different locations suggests that the distribution of volcanic glass across the lunar surface is much more widespread than has been previously documented. The presence of volcanic glass implies an explosive, vulcanian eruption style for localized DMDs, as this allows volcanic glass to rapidly quench, inhibiting crystallization, compared to the larger hawaiian-style eruptions typical of regional DMD emplacement where black beads indicate a higher degree of crystallization. Improved understanding of the local and global distributions of volcanic glass in lunar DMDs will further constrain lunar degassing and compositional evolution throughout lunar volcanic history.

Journal of Geophysical Research E: Planets↗

USGS2015

Limited role for thermal erosion by turbulent lava in proximal Athabasca Valles, Mars

The Athabasca Valles flood lava is among the most recent (<50 Ma) and best preserved effusive lava flows on Mars and was probably emplaced turbulently. The Williams et al. (2005) model of thermal erosion by lava has been applied to what we term “proximal Athabasca,” the 75 km long upstream portion of Athabasca Valles. For emplacement volumes of 5000 and 7500 km 3 and average flow thicknesses of 20 and 30 m, the duration of the eruption varies between ~11 and ~37 days. The erosion of the lava flow substrate is investigated for three eruption temperatures (1270°C, 1260°C, and 1250°C), and volatile contents equivalent to 0–65 vol % bubbles. The largest erosion depths of ~3.8–7.5 m are at the lava source, for 20 m thick and bubble-free flows that erupted at their liquidus temperature (1270°C). A substrate containing 25 vol % ice leads to maximum erosion. A lava temperature 20°C below liquidus reduces erosion depths by a factor of ~2.2. If flow viscosity increases with increasing bubble content in the lava, the presence of 30–50 vol % bubbles leads to erosion depths lower than those relative to bubble-free lava by a factor of ~2.4. The presence of 25 vol % ice in the substrate increases erosion depths by a factor of 1.3. Nevertheless, modeled erosion depths, consistent with the emplacement volume and flow duration constraints, are far less than the depth of the channel (~35–100 m). We conclude that thermal erosion does not appear to have had a major role in excavating Athabasca Valles.

Journal of Geophysical Research E: Planets↗

USGS2014

Residual shear strength variability as a primary control on movement of landslides reactivated by earthquake-induced ground motion: Implications for coastal Oregon, U.S.

Most large seismogenic landslides are reactivations of preexisting landslides with basal shear zones in the residual strength condition. Residual shear strength often varies during rapid displacement, but the response of residual shear zones to seismic loading is largely unknown. We used a ring shear apparatus to perform simulated seismic loading tests, constant displacement rate tests, and tests during which shear stress was gradually varied on specimens from two landslides to improve understanding of coseismic landslide reactivation and to identify shear strength models valid for slow gravitational failure through rapid coseismic failure. The landslides we studied represent many along the Oregon, U.S., coast. Seismic loading tests resulted in (1) catastrophic failure involving unbounded displacement when stresses represented those for the existing landslides and (2) limited to unbounded displacement when stresses represented those for hypothetical dormant landslides, suggesting that coseismic landslide reactivation may be significant during future great earthquakes occurring near the Oregon Coast. Constant displacement rate tests indicated that shear strength decreased exponentially during the first few decimeters of displacement but increased logarithmically with increasing displacement rate when sheared at 0.001 cm s −1 or greater. Dynamic shear resistance estimated from shear strength models correlated well with stresses observed during seismic loading tests, indicating that displacement rate and amount primarily controlled failure characteristics. We developed a stress-based approach to estimate coseismic landslide displacement that utilizes the variable shear strength model. The approach produced results that compared favorably to observations made during seismic loading tests, indicating its utility for application to landslides.

Oregon↗

USGS2014

Characterizing riverbed sediment using high-frequency acoustics 1: spectral properties of scattering

Bed-sediment classification using high-frequency hydro-acoustic instruments is challenging when sediments are spatially heterogeneous, which is often the case in rivers. The use of acoustic backscatter to classify sediments is an attractive alternative to analysis of topography because it is potentially sensitive to grain-scale roughness. Here, a new method is presented which uses high-frequency acoustic backscatter from multibeam sonar to classify heterogeneous riverbed sediments by type (sand, gravel,rock) continuously in space and at small spatial resolution. In this, the first of a pair of papers that examine the scattering signatures from a heterogeneous riverbed, methods are presented to construct spatially explicit maps of spectral properties from geo-referenced point clouds of geometrically and radiometrically corrected echoes. Backscatter power spectra are computed to produce scale and amplitude metrics that collectively characterize the length scales of stochastic measures of riverbed scattering, termed ‘stochastic geometries’. Backscatter aggregated over small spatial scales have spectra that obey a power-law. This apparently self-affine behavior could instead arise from morphological- and grain-scale roughnesses over multiple overlapping scales, or riverbed scattering being transitional between Rayleigh and geometric regimes. Relationships exist between stochastic geometries of backscatter and areas of rough and smooth sediments. However, no one parameter can uniquely characterize a particular substrate, nor definitively separate the relative contributions of roughness and acoustic impedance (hardness). Combinations of spectral quantities do, however, have the potential to delineate riverbed sediment patchiness, in a data-driven approach comparing backscatter with bed-sediment observations (which is the subject of part two of this manuscript).

Journal of Geophysical Research F: Earth Surface↗

USGS2014

Characterizing riverbed sediment using high-frequency acoustics 2: scattering signatures of Colorado River bed sediment in Marble and Grand Canyons

In this, the second of a pair of papers on the statistical signatures of riverbed sediment in high-frequency acoustic backscatter, spatially explicit maps of the stochastic geometries (length- and amplitude-scales) of backscatter are related to patches of riverbed surfaces composed of known sediment types, as determined by geo-referenced underwater video observations. Statistics of backscatter magnitudes alone are found to be poor discriminators between sediment types. However, the variance of the power spectrum, and the intercept and slope from a power-law spectral form (termed the spectral strength and exponent, respectively) successfully discriminate between sediment types. A decision-tree approach was able to classify spatially heterogeneous patches of homogeneous sands, gravels (and sand-gravel mixtures), and cobbles/boulders with 95, 88, and 91% accuracy, respectively. Application to sites outside the calibration, and surveys made at calibration sites at different times, were plausible based on observations from underwater video. Analysis of decision trees built with different training data sets suggested that the spectral exponent was consistently the most important variable in the classification. In the absence of theory concerning how spatially variable sediment surfaces scatter high-frequency sound, the primary advantage of this data-driven approach to classify bed sediment over alternatives is that spectral methods have well understood properties and make no assumptions about the distributional form of the fluctuating component of backscatter over small spatial scales.

Arizona↗

USGS2015

Glaciological and marine geological controls on terminus dynamics of Hubbard Glacier, southeast Alaska

Hubbard Glacier, located in southeast Alaska, is the world's largest non-polar tidewater glacier. It has been steadily advancing since it was first mapped in 1895; occasionally, the advance creates an ice or sediment dam that blocks a tributary fjord (Russell Fiord). The sustained advance raises the probability of long-term closure in the near-future, which will strongly impact the ecosystem of Russell Fiord and the nearby community of Yakutat. Here, we examine a 43-year record of flow speeds and terminus position to understand the large-scale dynamics of Hubbard Glacier. Our long-term record shows that the rate of terminus advance has increased slightly since 1895, with the exception of a slowed advance between approximately 1972 and 1984. The short-lived closure events in 1986 and 2002 were not initiated by perturbations in ice velocity or environmental forcings, but were likely due to fluctuations in sedimentation patterns at the terminus. This study points to the significance of a coupled system where short-term velocity fluctuations and morainal shoal development control tidewater glacier terminus position.

Alaska↗

USGS2015

Radar attenuation and temperature within the Greenland Ice Sheet

The flow of ice is temperature-dependent, but direct measurements of englacial temperature are sparse. The dielectric attenuation of radio waves through ice is also temperature-dependent, and radar sounding of ice sheets is sensitive to this attenuation. Here we estimate depth-averaged radar-attenuation rates within the Greenland Ice Sheet from airborne radar-sounding data and its associated radiostratigraphy. Using existing empirical relationships between temperature, chemistry, and radar attenuation, we then infer the depth-averaged englacial temperature. The dated radiostratigraphy permits a correction for the confounding effect of spatially varying ice chemistry. Where radar transects intersect boreholes, radar-inferred temperature is consistently higher than that measured directly. We attribute this discrepancy to the poorly recognized frequency dependence of the radar-attenuation rate and correct for this effect empirically, resulting in a robust relationship between radar-inferred and borehole-measured depth-averaged temperature. Radar-inferred englacial temperature is often lower than modern surface temperature and that of a steady state ice-sheet model, particularly in southern Greenland. This pattern suggests that past changes in surface boundary conditions (temperature and accumulation rate) affect the ice sheet's present temperature structure over a much larger area than previously recognized. This radar-inferred temperature structure provides a new constraint for thermomechanical models of the Greenland Ice Sheet.

Journal of Geophysical Research F: Earth Surface↗

USGS2015

Coevolution of bed surface patchiness and channel morphology: 1. Mechanisms of forced patch formation

Riverbeds frequently display a spatial structure where the sediment mixture composing the channel bed has been sorted into discrete patches of similar grain size. Even though patches are a fundamental feature in gravel bed rivers, we have little understanding of how patches form, evolve, and interact. Here we present a two-dimensional morphodynamic model that is used to examine in greater detail the mechanisms responsible for the development of forced bed surface patches and the coevolution of bed morphology and bed surface patchiness. The model computes the depth-averaged channel hydrodynamics, mixed-grain-size sediment transport, and bed evolution by coupling the river morphodynamic model Flow and Sediment Transport with Morphological Evolution of Channels (FaSTMECH) with a transport relation for gravel mixtures and the mixed-grain-size Exner equation using the active layer assumption. To test the model, we use it to simulate a flume experiment in which the bed developed a sequence of alternate bars and temporally and spatially persistent forced patches with a general pattern of coarse bar tops and fine pools. Cross-stream sediment flux causes sediment to be exported off of bars and imported into pools at a rate that balances downstream gradients in the streamwise sediment transport rate, allowing quasi-steady bar-pool topography to persist. The relative importance of lateral gravitational forces on the cross-stream component of sediment transport is a primary control on the amplitude of the bars. Because boundary shear stress declines as flow shoals over the bars, the lateral sediment transport is increasingly size selective and leads to the development of coarse bar tops and fine pools.

Journal of Geophysical Research F: Earth Surface↗

USGS2015

Coevolution of bed surface patchiness and channel morphology: 2. Numerical experiments

In gravel bed rivers, bed topography and the bed surface grain size distribution evolve simultaneously, but it is not clear how feedbacks between topography and grain sorting affect channel morphology. In this, the second of a pair of papers examining interactions between bed topography and bed surface sorting in gravel bed rivers, we use a two-dimensional morphodynamic model to perform numerical experiments designed to explore the coevolution of both free and forced bars and bed surface patches. Model runs were carried out on a computational grid simulating a 200 m long, 2.75 m wide, straight, rectangular channel, with an initially flat bed at a slope of 0.0137. Over five numerical experiments, we varied (a) whether an obstruction was present, (b) whether the sediment was a gravel mixture or a single size, and (c) whether the bed surface grain size feeds back on the hydraulic roughness field. Experiments with channel obstructions developed a train of alternate bars that became stationary and were connected to the obstruction. Freely migrating alternate bars formed in the experiments without channel obstructions. Simulations incorporating roughness feedbacks between the bed surface and flow field produced flatter, broader, and longer bars than simulations using constant roughness or uniform sediment. Our findings suggest that patches are not simply a by-product of bed topography, but they interact with the evolving bed and influence morphologic evolution.

Journal of Geophysical Research F: Earth Surface↗