USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Oceans”

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

Effects of river regulation on aeolian landscapes, Colorado River, southwestern USA

Connectivity between fluvial and aeolian sedimentary systems plays an important role in the physical and biological environment of dryland regions. This study examines the coupling between fluvial sand deposits and aeolian dune fields in bedrock canyons of the arid to semiarid Colorado River corridor, southwestern USA. By quantifying significant differences between aeolian landscapes with and without modern fluvial sediment sources, this work demonstrates for the first time that the flow- and sediment-limiting effects of dam operations affect sedimentary processes and ecosystems in aeolian landscapes above the fluvial high water line. Dune fields decoupled from fluvial sand supply have more ground cover (biologic crust and vegetation) and less aeolian sand transport than do dune fields that remain coupled to modern fluvial sand supply. The proportion of active aeolian sand area also is substantially lower in a heavily regulated river reach (Marble–Grand Canyon, Arizona) than in a much less regulated reach with otherwise similar environmental conditions (Cataract Canyon, Utah). The interconnections shown here among river flow and sediment, aeolian sand transport, and biologic communities in aeolian dunes demonstrate a newly recognized means by which anthropogenic influence alters dryland environments. Because fluvial–aeolian coupling is common globally, it is likely that similar sediment-transport connectivity and interaction with upland ecosystems are important in other dryland regions to a greater degree than has been recognized previously.

Colorado↗

Characterizing the performance of ecosystem models across time scales: A spectral analysis of the North American Carbon Program site‐level synthesis

[1] Ecosystem models are important tools for diagnosing the carbon cycle and projecting its behavior across space and time. Despite the fact that ecosystems respond to drivers at multiple time scales, most assessments of model performance do not discriminate different time scales. Spectral methods, such as wavelet analyses, present an alternative approach that enables the identification of the dominant time scales contributing to model performance in the frequency domain. In this study we used wavelet analyses to synthesize the performance of 21 ecosystem models at 9 eddy covariance towers as part of the North American Carbon Program's site‐level intercomparison. This study expands upon previous single‐site and single‐model analyses to determine what patterns of model error are consistent across a diverse range of models and sites. To assess the significance of model error at different time scales, a novel Monte Carlo approach was developed to incorporate flux observation error. Failing to account for observation error leads to a misidentification of the time scales that dominate model error. These analyses show that model error (1) is largest at the annual and 20–120 day scales, (2) has a clear peak at the diurnal scale, and (3) shows large variability among models in the 2–20 day scales. Errors at the annual scale were consistent across time, diurnal errors were predominantly during the growing season, and intermediate‐scale errors were largely event driven. Breaking spectra into discrete temporal bands revealed a significant model‐by‐band effect but also a nonsignificant model‐by‐site effect, which together suggest that individual models show consistency in their error patterns. Differences among models were related to model time step, soil hydrology, and the representation of photosynthesis and phenology but not the soil carbon or nitrogen cycles. These factors had the greatest impact on diurnal errors, were less important at annual scales, and had the least impact at intermediate time scales.

Journal of Geophysical Research: Biogeosciences↗

Linkages between denitrification and dissolved organicmatter quality, Boulder Creek watershed, Colorado

Dissolved organic matter (DOM) fuels the majority of in-stream microbial processes, including the removal of nitrate via denitrification. However, little is known about how the chemical composition of DOM influences denitrification rates. Water and sediment samples were collected across an ecosystem gradient, spanning the alpine to plains, in central Colorado to determine whether the chemical composition of DOM was related to denitrification rates. Laboratory bioassays measured denitrification potentials using the acetylene block technique and carbon mineralization via aerobic bioassays, while organic matter characteristics were evaluated using spectroscopic and fractionation methods. Denitrification potentials under ambient and elevated nitrate concentrations were strongly correlated with aerobic respiration rates and the percent mineralized carbon, suggesting that information about the aerobic metabolism of a system can provide valuable insight regarding the ability of the system to additionally reduce nitrate. Multiple linear regressions (MLR) revealed that under elevated nitrate concentrations denitrification potentials were positively related to the presence of protein-like fluorophores and negatively related to more aromatic and oxidized fractions of the DOM pool. Using MLR, the chemical composition of DOM, carbon, and nitrate concentrations explained 70% and 78% of the observed variability in denitrification potential under elevated and ambient nitrate conditions, respectively. Thus, it seems likely that DOM optical properties could help to improve predictions of nitrate removal in the environment. Finally, fluorescence measurements revealed that bacteria used both protein and humic-like organic molecules during denitrification providing further evidence that larger, more aromatic molecules are not necessarily recalcitrant in the environment.

Colorado↗

Peat accumulation in drained thermokarst lake basins in continuous, ice-rich permafrost, northern Seward Peninsula, Alaska

Thermokarst lakes and peat-accumulating drained lake basins cover a substantial portion of Arctic lowland landscapes, yet the role of thermokarst lake drainage and ensuing peat formation in landscape-scale carbon (C) budgets remains understudied. Here we use measurements of terrestrial peat thickness, bulk density, organic matter content, and basal radiocarbon age from permafrost cores, soil pits, and exposures in vegetated, drained lake basins to characterize regional lake drainage chronology, C accumulation rates, and the role of thermokarst-lake cycling in carbon dynamics throughout the Holocene on the northern Seward Peninsula, Alaska. Most detectable lake drainage events occurred within the last 4,000 years with the highest drainage frequency during the medieval climate anomaly. Peat accumulation rates were highest in young (50–500 years) drained lake basins (35.2 g C m −2 yr −1 ) and decreased exponentially with time since drainage to 9 g C m −2 yr −1 in the oldest basins. Spatial analyses of terrestrial peat depth, basal peat radiocarbon ages, basin geomorphology, and satellite-derived land surface properties (Normalized Difference Vegetation Index (NDVI); Minimum Noise Fraction (MNF)) from Landsat satellite data revealed significant relationships between peat thickness and mean basin NDVI or MNF. By upscaling observed relationships, we infer that drained thermokarst lake basins, covering 391 km 2 (76%) of the 515 km 2 study region, store 6.4–6.6 Tg organic C in drained lake basin terrestrial peat. Peat accumulation in drained lake basins likely serves to offset greenhouse gas release from thermokarst-impacted landscapes and should be incorporated in landscape-scale C budgets.

Alaska↗

Persistence of soil organic matter in eroding versus depositional landform positions

Soil organic matter (SOM) processes in dynamic landscapes are strongly influenced by soil erosion and sedimentation. We determined the contribution of physical isolation of organic matter (OM) inside aggregates, chemical interaction of OM with soil minerals, and molecular structure of SOM in controlling storage and persistence of SOM in different types of eroding and depositional landform positions. By combining density fractionation with elemental and spectroscopic analyses, we showed that SOM in depositional settings is less transformed and better preserved than SOM in eroding landform positions. However, which environmental factors exert primary control on storage and persistence of SOM depended on the nature of the landform position considered. In an annual grassland watershed, protection of SOM by physical isolation inside aggregates and chemical association of organic matter (complexation) with soil minerals, as assessed by correlation with radiocarbon concentration, were more effective in the poorly drained, lowest-lying depositional landform positions, compared to well-drained landform positions in the upper parts of the watershed. Results of this study demonstrated that processes of soil erosion and deposition are important mechanisms of long-term OM stabilization.

Journal of Geophysical Research: Biogeosciences↗

Landscape controls on total and methyl Hg in the Upper Hudson River basin, New York, USA

Approaches are needed to better predict spatial variation in riverine Hg concentrations across heterogeneous landscapes that include mountains, wetlands, and open waters. We applied multivariate linear regression to determine the landscape factors and chemical variables that best account for the spatial variation of total Hg (THg) and methyl Hg (MeHg) concentrations in 27 sub-basins across the 493 km 2 upper Hudson River basin in the Adirondack Mountains of New York. THg concentrations varied by sixfold, and those of MeHg by 40-fold in synoptic samples collected at low-to-moderate flow, during spring and summer of 2006 and 2008. Bivariate linear regression relations of THg and MeHg concentrations with either percent wetland area or DOC concentrations were significant but could account for only about 1/3 of the variation in these Hg forms in summer. In contrast, multivariate linear regression relations that included metrics of (1) hydrogeomorphology, (2) riparian/wetland area, and (3) open water, explained about 66% to >90% of spatial variation in each Hg form in spring and summer samples. These metrics reflect the influence of basin morphometry and riparian soils on Hg source and transport, and the role of open water as a Hg sink. Multivariate models based solely on these landscape metrics generally accounted for as much or more of the variation in Hg concentrations than models based on chemical and physical metrics, and show great promise for identifying waters with expected high Hg concentrations in the Adirondack region and similar glaciated riverine ecosystems.

New York↗

Carbon and sediment accumulation in the Everglades (USA) during the past 4000 years: rates, drivers, and sources of error

Tropical and sub-tropical wetlands are considered to be globally important sources for greenhouse gases but their capacity to store carbon is presumably limited by warm soil temperatures and high rates of decomposition. Unfortunately, these assumptions can be difficult to test across long timescales because the chronology, cumulative mass, and completeness of a sedimentary profile are often difficult to establish. We therefore made a detailed analysis of a core from the principal drainage outlet of the Everglades of South Florida, to assess these problems and determine the factors that could govern carbon accumulation in this large sub-tropical wetland. Accelerator mass spectroscopy dating provided direct evidence for both hard-water and open-system sources of dating errors, whereas cumulative mass varied depending upon the type of method used. Radiocarbon dates of gastropod shells, nevertheless, seemed to provide a reliable chronology for this core once the hard-water error was quantified and subtracted. Long-term accumulation rates were then calculated to be 12.1 g m -2 yr -1 for carbon, which is less than half the average rate reported for northern and tropical peatlands. Moreover, accumulation rates remained slow and relatively steady for both organic and inorganic strata, and the slow rate of sediment accretion ( 0.2 mm yr -1 ) tracked the correspondingly slow rise in sea level (0.35 mm yr -1 ) reported for South Florida over the past 4000 years. These results suggest that sea level and the local geologic setting may impose long-term constraints on rates of sediment and carbon accumulation in the Everglades and other wetlands.

Florida↗

Spatiotemporal analysis of black spruce forest soils and implications for the fate of C

Post-fire storage of carbon (C) in organic-soil horizons was measured in one Canadian and three Alaskan chronosequences in black spruce forests, together spanning stand ages of nearly 200 yrs. We used a simple mass balance model to derive estimates of inputs, losses, and accumulation rates of C on timescales of years to centuries. The model performed well for the surface and total organic soil layers and presented questions for resolving the dynamics of deeper organic soils. C accumulation in all study areas is on the order of 20–40 gC/m 2 / yr for stand ages up to ∼200 yrs. Much larger fluxes, both positive and negative, are detected using incremental changes in soil C stocks and by other studies using eddy covariance methods for CO 2 . This difference suggests that over the course of stand replacement, about 80% of all net primary production (NPP) is returned to the atmosphere within a fire cycle, while about 20% of NPP enters the organic soil layers and becomes available for stabilization or loss via decomposition, leaching, or combustion. Shifts toward more frequent and more severe burning and degradation of deep organic horizons would likely result in an acceleration of the carbon cycle, with greater CO 2 emissions from these systems overall.

Journal of Geophysical Research↗

Reoccupation of floodplains by rivers and its relation to the age structure of floodplain vegetation

River channel dynamics over many decades provide a physical control on the age structure of floodplain vegetation as a river occupies and abandons locations. Floodplain reoccupation by a river, in particular, determines the interval of time during which vegetation can establish and mature. A general framework for analyzing floodplain reoccupation and a time series model are developed and applied to five alluvial rivers in the United States. Channel dynamics in these rivers demonstrate time-scale dependence with short-term oscillation in active channel area in response to floods and subsequent vegetation growth and progressive lateral movement that accounts for much of the cumulative area occupied by the rivers over decades. Rivers preferentially reoccupy locations recently abandoned causing a decreasing probability of reoccupation with time since abandonment. For a typical case, a river is 10 times more likely to reoccupy an area it abandoned in the past decade than it is to reoccupy an area it abandoned 30 yrs ago. The decreasing probability of reoccupation over time is consistent with observations of persistent stands of late seral stage floodplain forest. A power function provides a robust approach for estimating the cumulative area occupied by a river and the age structure of riparian forests resulting from a specific historical sequence of streamflow in comparison to either linear or exponential alternatives.

Journal of Geophysical Research↗

Dissolved organic carbon and chromophoric dissolved organic matter properties of rivers in the USA

Dissolved organic carbon (DOC) concentration and chromophoric dissolved organic matter (CDOM) parameters were measured over a range of discharge in 30 U.S. rivers, covering a diverse assortment of fluvial ecosystems in terms of watershed size and landscape drained. Relationships between CDOM absorption at a range of wavelengths (a 254 , a 350 , a 440 ) and DOC in the 30 watersheds were found to correlate strongly and positively for the majority of U.S. rivers. However, four rivers (Colorado, Colombia, Rio Grande and St. Lawrence) exhibited statistically weak relationships between CDOM absorption and DOC. These four rivers are atypical, as they either drain from the Great Lakes or experience significant impoundment of water within their watersheds, and they exhibited values for dissolved organic matter (DOM) parameters indicative of autochthonous or anthropogenic sources or photochemically degraded allochthonous DOM and thus a decoupling between CDOM and DOC. CDOM quality parameters in the 30 rivers were found to be strongly correlated to DOM compositional metrics derived via XAD fractionation, highlighting the potential for examining DOM biochemical quality from CDOM measurements. This study establishes the ability to derive DOC concentration from CDOM absorption for the majority of U.S. rivers, describes characteristics of riverine systems where such an approach is not valid, and emphasizes the possibility of examining DOM composition and thus biogeochemical function via CDOM parameters. Therefore, the usefulness of CDOM measurements, both laboratory-based analyses and in situ instrumentation, for improving spatial and temporal resolution of DOC fluxes and DOM dynamics in future studies is considerable in a range of biogeochemical studies.

Journal of Geophysical Research G: Biogeosciences↗

Uncovering the nonadiabatic response of geosynchronous electrons to geomagnetic disturbance

We describe an energy spectrum method for scaling electron integral flux, which is measured at a constant energy, to phase space density at a constant value of the first adiabatic invariant which removes much of the variation due to reversible adiabatic effects. Applying this method to nearly a solar cycle (1995 - 2006) of geosynchronous electron integral flux (E>2.0MeV) from the GOES satellites, we see that much of the diurnal variation in electron phase space density at constant energy can be removed by the transformation to phase space density at constant μ (4000 MeV/G). This allows us a clearer picture of underlying non-adiabatic electron population changes due to geomagnetic activity. Using scaled phase space density, we calculate the percentage of geomagnetic storms resulting in an increase, decrease or no change in geosynchronous electrons as 38%, 7%, and 55%, respectively. We also show examples of changes in the electron population that may be different than the unscaled fluxes alone suggest. These examples include sudden electron enhancements during storms which appear during the peak of negative Dst for μ-scaled phase space density, contrary to the slow increase seen during the recovery phase for unscaled phase space density for the same event.

Journal of Geophysical Research↗

Problem of the Love‐Gannon relation between the asymmetric disturbance field and Dst

Love and Gannon (2009) discovered that statistically, over a fifty year period the difference in the dawn and dusk disturbance‐field H component at low latitudes (hourly averaged) is linearly proportional to Dst. If the difference is designated by δ DD in units of nT/R E , then the Love‐Gannon (L‐G) relation is δ DD = − 0.2 Dst. At any time departures from the relation can be large. Nonetheless, the relation is evident for all values of Dst and persists throughout magnetic storms, both the main phase and the recovery phase. The Love‐Gannon discovery presents a problem to current understanding of the relation between the causes of δ DD and Dst because the dawn dusk asymmetry in the disturbance field is presumably governed by a long‐established magnetosphere‐ionosphere coupling theory which predicts a characteristic time scale (the shielding time) of less than an hour whereas the characteristic time scale for Dst (the ring current decay time) is more like ten hours. Thus, without forcing both time scales toward each other to the limits of their ranges, a linear proportionality between δ DD and Dst cannot be derived from the current understanding of the causes of the asymmetry and the ring current. This conclusion is the paper's main contribution. In addition, we attempt to get around the conflict of time scales by looking at other possibilities for generating δ DD that depend directly on the ring current. The most promising of these is the possibility that the ring current decay mechanism creates a quasi‐permanent, local‐time modification of the ring current compared to what it would be in the absence of the decay mechanism and that this modification causes a field‐aligned current that closes through the ionosphere and generates the asymmetry δ DD . This idea has the virtue of coupling the asymmetry directly to the ring current and of accounting for the persistence of the L‐G proportionality through the recovery phase of magnetic storms.

Journal of Geophysical Research A: Space Physics↗

Laboratory triggering of stick-slip events by oscillatory loading in the presence of pore fluid with implications for physics of tectonic tremor

The physical mechanism by which the low-frequency earthquakes (LFEs) that make up portions of tectonic (also called non-volcanic) tremor are created is poorly understood. In many areas of the world, tectonic tremor and LFEs appear to be strongly tidally modulated, whereas ordinary earthquakes are not. Anomalous seismic wave speeds, interpreted as high pore fluid pressure, have been observed in regions that generate tremor. Here we build upon previous laboratory studies that investigated the response of stick-slip on artificial faults to oscillatory, tide-like loading. These previous experiments were carried out using room-dry samples of Westerly granite, at one effective stress. Here we augment these results with new experiments on Westerly granite, with the addition of varying effective stress using pore fluid at two pressures. We find that raising pore pressure, thereby lowering effective stress can significantly increase the degree of correlation of stick-slip to oscillatory loading. We also find other pore fluid effects that become important at higher frequencies, when the period of oscillation is comparable to the diffusion time of pore fluid into the fault. These results help constrain the conditions at depth that give rise to tidally modulated LFEs, providing confirmation of the effective pressure law for triggering and insights into why tremor is tidally modulated while earthquakes are at best only weakly modulated.

Journal of Geophysical Research B: Solid Earth↗

Heat flow in vapor dominated areas of the Yellowstone Plateau volcanic field: implications for the thermal budget of the Yellowstone Caldera

Characterizing the vigor of magmatic activity in Yellowstone requires knowledge of the mechanisms and rates of heat transport between magma and the ground surface. We present results from a heat flow study in two vapor dominated, acid-sulfate thermal areas in the Yellowstone Caldera, the 0.11 km 2 Obsidian Pool Thermal Area (OPTA) and the 0.25 km 2 Solfatara Plateau Thermal Area (SPTA). Conductive heat flux through a low permeability layer capping large vapor reservoirs is calculated from soil temperature measurements at >600 locations and from laboratory measurements of soil properties. The conductive heat output is 3.6 ± 0.4 MW and 7.5 ± 0.4 MW from the OPTA and the SPTA, respectively. The advective heat output from soils is 1.3 ± 0.3 MW and 1.2 ± 0.3 MW from the OPTA and the SPTA, respectively and the heat output from thermal pools in the OPTA is 6.8 ± 1.4 MW. These estimates result in a total heat output of 11.8 ± 1.4 MW and 8.8 ± 0.4 MW from OPTA and SPTA, respectively. Focused zones of high heat flux in both thermal areas are roughly aligned with regional faults suggesting that faults in both areas serve as conduits for the rising acid vapor. Extrapolation of the average heat flux from the OPTA (103 ± 2 W·m −2 ) and SPTA (35 ± 3 W·m −2 ) to the ~35 km 2 of vapor dominated areas in Yellowstone yields 3.6 and 1.2 GW, respectively, which is less than the total heat output transported by steam from the Yellowstone Caldera as estimated by the chloride inventory method (4.0 to 8.0 GW).

Journal of Geophysical Research B: Solid Earth↗

Externally triggered renewed bubble nucleation in basaltic magma: the 12 October 2008 eruption at Halema‘uma‘u Overlook vent, Kīlauea, Hawai‘i, USA

From October 2008 until present, dozens of small impulsive explosive eruptions occurred from the Overlook vent on the southeast side of Halema&lsquo;uma&lsquo;u Crater, at Kīlauea volcano, USA. These eruptions were triggered by rockfalls from the walls of the volcanic vent and conduit onto the top of the lava column. Here we use microtextural observations and data from clasts erupted during the well-characterized 12 October 2008 explosive eruption at Halema&lsquo;uma&lsquo;u to extend existing models of eruption triggering. We present a potential mechanism for this eruption by combining microtextural observations with existing geophysical and visual data sets. We measure the size and number density of bubbles preserved in juvenile ejecta using 2D images and X-ray microtomography. Our data suggest that accumulations of large bubbles with diameters of >50 &mu; m to at least millimeters existed at shallow levels within the conduit prior to the 12 October 2008 explosion. Furthermore, a high number density of small bubbles <50 &mu; m is measured in the clasts, implying very rapid nucleation of bubbles. Visual observations, combined with preexisting geophysical data, suggest that the impact of rockfalls onto the magma free surface induces pressure changes over short timescales that (1) nucleated new additional bubbles in the shallow conduit leading to high number densities of small bubbles and (2) expanded the preexisting bubbles driving upward acceleration. The trigger of eruption and bubble nucleation is thus external to the degassing system.

Hawaii↗

Hydrate morphology: Physical properties of sands with patchy hydrate saturation

The physical properties of gas hydrate-bearing sediments depend on the volume fraction and spatial distribution of the hydrate phase. The host sediment grain size and the state of effective stress determine the hydrate morphology in sediments; this information can be used to significantly constrain estimates of the physical properties of hydrate-bearing sediments, including the coarse-grained sands subjected to high effective stress that are of interest as potential energy resources. Reported data and physical analyses suggest hydrate-bearing sands contain a heterogeneous, patchy hydrate distribution, whereby zones with 100% pore-space hydrate saturation are embedded in hydrate-free sand. Accounting for patchy rather than homogeneous hydrate distribution yields more tightly constrained estimates of physical properties in hydrate-bearing sands and captures observed physical-property dependencies on hydrate saturation. For example, numerical modeling results of sands with patchy saturation agree with experimental observation, showing a transition in stiffness starting near the series bound at low hydrate saturations but moving toward the parallel bound at high hydrate saturations. The hydrate-patch size itself impacts the physical properties of hydrate-bearing sediments; for example, at constant hydrate saturation, we find that conductivity (electrical, hydraulic and thermal) increases as the number of hydrate-saturated patches increases. This increase reflects the larger number of conductive flow paths that exist in specimens with many small hydrate-saturated patches in comparison to specimens in which a few large hydrate saturated patches can block flow over a significant cross-section of the specimen.

Journal of Geophysical Research B: Solid Earth↗

Clustering of GPS velocities in the Mojave Block, southeastern California

We find subdivisions within the Mojave Block using cluster analysis to identify groupings in the velocities observed at GPS stations there. The clusters are represented on a fault map by symbols located at the positions of the GPS stations, each symbol representing the cluster to which the velocity of that GPS station belongs. Fault systems that separate the clusters are readily identified on such a map. The most significant representation as judged by the gap test involves 4 clusters within the Mojave Block. The fault systems bounding the clusters from east to west are 1) the faults defining the eastern boundary of the Northeast Mojave Domain extended southward to connect to the Hector Mine rupture, 2) the Calico-Paradise fault system, 3) the Landers-Blackwater fault system, and 4) the Helendale-Lockhart fault system. This division of the Mojave Block is very similar to that proposed by Meade and Hager. However, no cluster boundary coincides with the Garlock Fault, the northern boundary of the Mojave Block. Rather, the clusters appear to continue without interruption from the Mojave Block north into the southern Walker Lane Belt, similar to the continuity across the Garlock Fault of the shear zone along the Blackwater-Little Lake fault system observed by Peltzer et al. Mapped traces of individual faults in the Mojave Block terminate within the block and do not continue across the Garlock Fault [Dokka and Travis, ].

California↗

Bathymetric controls on sediment transport in the Hudson River estuary: Lateral asymmetry and frontal trapping

Analyses of field observations and numerical model results have identified that sediment transport in the Hudson River estuary is laterally segregated between channel and shoals, features frontal trapping at multiple locations along the estuary, and varies significantly over the spring-neap tidal cycle. Lateral gradients in depth, and therefore baroclinic pressure gradient and stratification, control the lateral distribution of sediment transport. Within the saline estuary, sediment fluxes are strongly landward in the channel and seaward on the shoals. At multiple locations, bottom salinity fronts form at bathymetric transitions in width or depth. Sediment convergences near the fronts create local maxima in suspended-sediment concentration and deposition, providing a general mechanism for creation of secondary estuarine turbidity maxima at bathymetric transitions. The lateral bathymetry also affects the spring-neap cycle of sediment suspension and deposition. In regions with broad, shallow shoals, the shoals are erosional and the channel is depositional during neap tides, with the opposite pattern during spring tides. Narrower, deeper shoals are depositional during neaps and erosional during springs. In each case, the lateral transfer is from regions of higher to lower bed stress, and depends on the elevation of the pycnocline relative to the bed. Collectively, the results indicate that lateral and along-channel gradients in bathymetry and thus stratification, bed stress, and sediment flux lead to an unsteady, heterogeneous distribution of sediment transport and trapping along the estuary rather than trapping solely at a turbidity maximum at the limit of the salinity intrusion.

Journal of Geophysical Research↗