SEARCH · USGS Science

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.

SourceResearch connectionFrom yearThrough year

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

At least 577 records · Page 32Linked to original sources

USGS2006

Aftershock decay, productivity, and stress rates in Hawaii: Indicators of temperature and stress from magma sources

We examined dozens of aftershock sequences in Hawaii in terms of Gutenberg-Richter and modified Omori law parameters. We studied p, the rate of aftershock decay; A p , the aftershock productivity, defined as the observed divided by the expected number of aftershocks; and c, the time delay when aftershock rates begin to fall. We found that for earthquakes shallower than 20 km, p values >1.2 are near active magma centers. We associate this high decay rate with higher temperatures and faster stress relaxation near magma reservoirs. Deep earthquakes near Kilauea's inferred magma transport path show a range of p values, suggesting the absence of a large, deep magma reservoir. Aftershock productivity is >4.0 for flank earthquakes known to be triggered by intrusions but is normal (0.25 to 4.0) for isolated main shocks. We infer that continuing, post-main shock stress from the intrusion adds to the main shock's stress step and causes higher A p . High A p in other zones suggests less obvious intrusions and pulsing magma pressure near Kilauea's feeding conduit. We calculate stress rates and stress rate changes from pre-main shock and aftershock rates. Stress rate increased after many intrusions but decreased after large M7–8 earthquakes. Stress rates are highest in the seismically active volcano flanks and lowest in areas far from volcanic centers. We found sequences triggered by intrusions tend to have high A p , high (>0.10 day) c values, a stress rate increase, and sometimes a peak in aftershock rate hours after the main shock. We interpret these values as indicating continuing intrusive stress after the main shock.

Hawai'i↗

USGS2006

An evaluation of factors influencing pore pressure in accretionary complexes: Implications for taper angle and wedge mechanics

At many subduction zones, accretionary complexes form as sediment is off-scraped from the subducting plate. Mechanical models that treat accretionary complexes as critically tapered wedges of sediment demonstrate that pore pressure controls their taper angle by modifying basal and internal shear strength. Here, we combine a numerical model of groundwater flow with critical taper theory to quantify the effects of sediment and de??collement permeability, sediment thickness, sediment partitioning between accretion and underthrusting, and plate convergence rate on steady state pore pressure. Our results show that pore pressure in accretionary wedges can be viewed as a dynamically maintained response to factors which drive pore pressure (source terms) and those that limit flow (permeability and drainage path length). We find that sediment permeability and incoming sediment thickness are the most important factors, whereas fault permeability and the partitioning of sediment have a small effect. For our base case model scenario, as sediment permeability is increased, pore pressure decreases from near-lithostatic to hydrostatic values and allows stable taper angles to increase from ??? 2.5?? to 8??-12.5??. With increased sediment thickness in our models (from 100 to 8000 m), increased pore pressure drives a decrease in stable taper angle from 8.4??-12.5?? to <2.5-5??. In general, low-permeability and thick incoming sediment sustain high pore pressures consistent with shallowly tapered geometry, whereas high-permeability and thin incoming sediment should result in steep geometry. Our model results compare favorably with available data from active accretionary complexes. Active margins characterized by a significant proportion of fine-grained sediment within the incoming section, such as northern Antilles and eastern Nankai, exhibit thin taper angles, whereas those characterized by a higher proportion of sandy turbidites, such as Cascadia, Chile, and Mexico, have steep taper angles. Observations from active margins also indicate a strong trend of decreasing taper angle (from >15?? to <4??) with increased sediment thickness (from <1 to 7 km). One key implication is that hydrologic properties may strongly influence the strength of the crust in a wide range of geologic settings. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Dislocation pileup as a representation of strain accumulation on a strike-slip fault

The conventional model of strain accumulation on a vertical transform fault is a discrete screw dislocation in an elastic half-space with the Burgers vector of the dislocation increasing at the rate of relative plate motion. It would be more realistic to replace that discrete dislocation by a dislocation distribution, presumably a pileup in which the individual dislocations are in equilibrium. The length of the pileup depends upon the applied stress and the amount of slip that has occurred at depth. I argue here that the dislocation pileup (the transition on the fault from no slip to slip at the full plate rate) occupies a substantial portion of the lithosphere thickness. A discrete dislocation at an adjustable depth can reproduce the surface deformation profile predicted by a pileup so closely that it will be difficult to distinguish between the two models. The locking depth (dislocation depth) of that discrete dislocation approximation is substantially (???30%) larger than that (depth to top of the pileup) in the pileup model. Thus, in inverting surface deformation data using the discrete dislocation model, the locking depth in the model should not be interpreted as the true locking depth. Although dislocation pileup models should provide a good explanation of the surface deformation near the fault trace, that explanation may not be adequate at greater distances from the fault trace because approximating the expected horizontally distributed deformation at subcrustal depths by uniform slip concentrated on the fault is not justified.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

A new 1649-1884 catalog of destructive earthquakes near Tokyo and implications for the long-term seismic process

In order to assess the long-term character of seismicity near Tokyo, we construct an intensity-based catalog of damaging earthquakes that struck the greater Tokyo area between 1649 and 1884. Models for 15 historical earthquakes are developed using calibrated intensity attenuation relations that quantitatively convey uncertainties in event location and magnitude, as well as their covariance. The historical catalog is most likely complete for earthquakes M ??? 6.7; the largest earthquake in the catalog is the 1703 M ??? 8.2 Genroku event. Seismicity rates from 80 years of instrumental records, which include the 1923 M = 7.9 Kanto shock, as well as interevent times estimated from the past ???7000 years of paleoseismic data, are combined with the historical catalog to define a frequency-magnitude distribution for 4.5 ??? M ??? 8.2, which is well described by a truncated Gutenberg-Richter relation with a b value of 0.96 and a maximum magnitude of 8.4. Large uncertainties associated with the intensity-based catalog are propagated by a Monte Carlo simulation to estimations of the scalar moment rate. The resulting best estimate of moment rate during 1649-2003 is 1.35 ?? 1026 dyn cm yr-1 with considerable uncertainty at the 1??, level: (-0.11, + 0.20) ?? 1026 dyn cm yr-1. Comparison with geodetic models of the interseismic deformation indicates that the geodetic moment accumulation and likely moment release rate are roughly balanced over the catalog period. This balance suggests that the extended catalog is representative of long-term seismic processes near Tokyo and so can be used to assess earthquake probabilities. The resulting Poisson (or time-averaged) 30-year probability for M ??? 7.9 earthquakes is 7-11%.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Analysis of gravity and topography in the GLIMPSE study region: Isostatic compensation and uplift of the Sojourn and Hotu Matua Ridge systems

The Gravity Lieations Intraplate Melting Petrologic and Seismic Expedition (GLIMPSE) Experiment investigated the formation of a series of non-hot spot, intraplate volcanic ridges in the South Pacific and their relationship to cross-grain gravity lineaments detected by satellite altimetry. Using shipboard gravity measurements and a simple model of surface loading of a thin elastic plate, we estimate effective elastic thicknesses ranging from ???2 km beneath the Sojourn Ridge to a maximum of 10 km beneath the Southern Cross Seamount. These elastic thicknesses are lower than predicted for the 3-9 Ma seafloor on which the volcanoes lie, perhaps due to reheating and thinning of the plate during emplacement. Anomalously low apparent densities estimated for the Matua and Southern Cross seamounts 2050 and 2250 kg m-3, respectively, probably are artifacts caused by the assumption of only surface loading, ignoring the presence of subsurface loading in the form of underplated crust and/or low-density mantle. Using satellite free-air gravity and shipboard bathymetry, we calculate the age-detrended, residual mantle Bouguer anomaly (rMBA). The rMBA corrects the free-air anomaly for the direct effects of topography, including the thickening of the crust beneath the seamounts and volcanic ridges due to surface loading of the volcanic edifices. There are broad, negative rMBA anomalies along the Sojourn and Brown ridges and the Hotu Matua seamount chain that extend nearly to the East Pacific Rise. These negative rMBA anomalies connect to negative free-air anomalies in the western part of the study area that have been recognized previously as the beginnings of the cross-grain gravity lineaments. Subtracting the topographic effects of surface loading by the ridges and seamounts from the observed topography reveals that the ridges are built on broad bands of anomalously elevated seafloor. This swell topography and the negative rMBA anomalies contradict the predictions of lithospheric cracking models for the origin of gravity lineaments and associated volcanic ridges, favoring models with a dynamic mantle component such as small-scale convection or channelized asthenospheric return flow. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Distribution of recent volcanism and the morphology of seamounts and ridges in the GLIMPSE study area: Implications for the lithospheric cracking hypothesis for the origin of intraplate, non-hot spot volcanic chains

Lithospheric cracking by remotely applied stresses or thermoelastic stresses has been suggested to be the mechanism responsible for the formation of intraplate volcanic ridges in the Pacific that clearly do not form above fixed hot spots. As part of the Gravity Lineations Intraplate Melting Petrology and Seismic Expedition (GLIMPSE) project designed to investigate the origin of these features, we have mapped two volcanic chains that are actively forming to the west of the East Pacific Rise using multibeam echo sounding and side-scan sonar. Side-scan sonar reveals the distribution of rough seafloor corresponding to recent, unsedimented lava flows. In the Hotu Matua volcanic complex, recent flows and volcanic edifices are distributed over a region 450 km long and up to 65 km wide, with an apparent, irregular age progression from older flows in the west to younger in the east. The 550-km-long Southern Cross Seamount/Sojourn Ridge/Brown Ridge chain appears to have been recently active only at its eastern end near the East Pacific Rise. A third region of recent flows is found 120 km north of Southern Cross Seamount in seafloor approximately 9 Myr old. No indication of lithospheric extension in the form of faulting or graben formation paralleling the trend of the volcanic chains is found in the vicinity of recent flows or anywhere else in the study area. Thermoelastic cracking could be a factor in the formation of a few small, very narrow volcanic ridges, but most of the volcanic activity is broadly distributed in wide swaths with no indication of formation along narrow cracks. The Sojourn and Brown chains appear to begin as distributed zones of small seamounts that later develop into segmented ridges, perhaps under the influence of membrane stresses from self-loading. We suggest that the linear volcanic chains are created by moving melting anomalies in the asthenosphere and that lithospheric cracking plays at most a secondary role. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Regional P wave velocity structure of the Northern Cascadia Subduction Zone

This paper presents the first regional three-dimensional, P wave velocity model for the Northern Cascadia Subduction. Zone (SW British Columbia and NW Washington State) constructed through tomographic inversion of first-arrival traveltime data from active source experiments together with earthquake traveltime data recorded at permanent stations. The velocity model images the structure of the subducting Juan de Fuca plate, megathrust, and the fore-arc crust and upper mantle. Beneath southern Vancouver Island the megathrust above the Juan de Fuca plate is characterized by a broad zone (25-35 km depth) having relatively low velocities of 6.4-6.6 km/s. This relative low velocity zone coincides with the location of most of the episodic tremors recently mapped beneath Vancouver Island, and its low velocity may also partially reflect the presence of trapped fluids and sheared lower crustal rocks. The rocks of the Olympic Subduction Complex are inferred to deform aseismically as evidenced by the lack of earthquakes withi the low-velocity rocks. The fore-arc upper mantle beneath the Strait of Georgia and Puget Sound is characterized by velocities of 7.2-7.6 km/s. Such low velocities represent regional serpentinization of the upper fore-arc mantle and provide evidence for slab dewatering and densification. Tertiary sedimentary basins in the Strait of Georgia and Puget Lowland imaged by the velocity model lie above the inferred region of slab dewatering and densification and may therefore partly result from a higher rate of slab sinking. In contrast, sedimentary basins in the Strait of Juan de Fuca lie in a synclinal depression in the Crescent Terrane. The correlation of in-slab earthquake hypocenters M>4 with P wave velocities greater than 7.8 km/s at the hypocenters suggests that they originate near the oceanic Moho of the subducting Juan de Fuca plate. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Damped regional-scale stress inversions: Methodology and examples for southern California and the Coalinga aftershock sequence

We present a new focal mechanism stress inversion technique to produce regional-scale models of stress orientation containing the minimum complexity necessary to fit the data. Current practice is to divide a region into small subareas and to independently fit a stress tensor to the focal mechanisms of each subarea. This procedure may lead to apparent spatial variability that is actually an artifact of overfitting noisy data or nonuniquely fitting data that does not completely constrain the stress tensor. To remove these artifacts while retaining any stress variations that are strongly required by the data, we devise a damped inversion method to simultaneously invert for stress in all subareas while minimizing the difference in stress between adjacent subareas. This method is conceptually similar to other geophysical inverse techniques that incorporate damping, such as seismic tomography. In checkerboard tests, the damped inversion removes the stress rotation artifacts exhibited by an undamped inversion, while resolving sharper true stress rotations than a simple smoothed model or a moving-window inversion. We show an example of a spatially damped stress field for southern California. The methodology can also be used to study temporal stress changes, and an example for the Coalinga, California, aftershock sequence is shown. We recommend use of the damped inversion technique for any study examining spatial or temporal variations in the stress field.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

Imaging the transition from Aleutian subduction to Yakutat collision in central Alaska, with local earthquakes and active source data

In southern and central Alaska the subduction and active volcanism of the Aleutian subduction zone give way to a broad plate boundary zone with mountain building and strike-slip faulting, where the Yakutat terrane joins the subducting Pacific plate. The interplay of these tectonic elements can be best understood by considering the entire region in three dimensions. We image three-dimensional seismic velocity using abundant local earthquakes, supplemented by active source data. Crustal low-velocity correlates with basins. The Denali fault zone is a dominant feature with a change in crustal thickness across the fault. A relatively high-velocity subducted slab and a low-velocity mantle wedge are observed, and high Vp/Vs beneath the active volcanic systems, which indicates focusing of partial melt. North of Cook Inlet, the subducted Yakutat slab is characterized by a thick low-velocity, high-Vp/Vs, crust. High-velocity material above the Yakutat slab may represent a residual older slab, which inhibits vertical flow of Yakutat subduction fluids. Alternate lateral flow allows Yakutat subduction fluids to contribute to Cook Inlet volcanism and the Wrangell volcanic field. The apparent northeast edge of the subducted Yakutat slab is southwest of the Wrangell volcanics, which have adakitic composition consistent with melting of this Yakutat slab edge. In the mantle, the Yakutat slab is subducting with the Pacific plate, while at shallower depths the Yakutat slab overthrusts the shallow Pacific plate along the Transition fault. This region of crustal doubling within the shallow slab is associated with extremely strong plate coupling and the primary asperity of the Mw 9.2 great 1964 earthquake. Copyright 2006 by the American Geophysical Union.

Alaska↗

USGS2007

Microplate model for the present-day deformation of Tibet

[1] Site velocities from 349 Global Positioning System (GPS) stations are used to construct an 11-element quasi-rigid block model of the Tibetan Plateau and its surroundings. Rigid rotations of five major blocks are well determined, and average translation velocities of six smaller blocks can be constrained. Where data are well distributed the velocity field can be explained well by rigid block motion and fault slip across block boundaries. Residual misfits average 1.6 mm/yr compared to typical one standard deviation velocity uncertainties of 1.3 mm/yr. Any residual internal straining of the blocks is small and heterogeneous. However, residual substructure might well represent currently unresolved motions of smaller blocks. Although any smaller blocks must move at nearly the same rate as the larger blocks within which they lie, undetected relative motions between them could be significant, particularly where there are gaps in GPS coverage. Predicted relative motions between major blocks agree with the observed sense of slip and along-strike partitioning of motion across major faults. However, predicted slip rates across Tibet's major strike-slip faults are low, only 5–12 mm/yr, a factor of 2–3 smaller than most rates estimated from fault offset features dated by radiometric methods as ∼2000 to ∼100,000 year old. Previous work has suggested that both GPS data and low fault slip rates are incompatible with rigid block motions of Tibet. The results reported here overcome these objections.

Journal of Geophysical Research B: Solid Earth↗

USGS2006

February 2003 marine atmospheric conditions and the bora over the northern Adriatic

A winter oceanographic field experiment provided an opportunity to examine the atmospheric marine conditions over the northern Adriatic. Mean February winds are from a northeasterly direction over most of the Adriatic and a more northerly direction along the western coast. Wind speeds are fastest in jets over the NE coast during bora events and weakest in the mid-northwestern Adriatic. Diurnal air temperature cycles are smallest on the NE coast and largest in the midwestern Adriatic. The maximum sea-air difference is +10??C on the eastern coast and near zero on the midwestern Adriatic. Boras are northeasterly (from) wind events that sweep off Croatia and Slovenia, bringing slightly colder and drier air over the northern Adriatic. The main bora season is December to March. Winter 2002-2003 was normal for bora events. Synoptic-scale temporal variations are correlated over the northern Adriatic. Fastest Bora winds and highest wind stress over the northern Adriatic is concentrated in four topographically controlled jets. The strongest is the Senj Jet, while the Trieste Jet extends across the entire northern Adriatic. Between each two jets is a weak wind zone. The greatest mean net heat loss is in bora jets in the NE Adriatic, where it was -438 W m-2 and is weakest in the midwestern northern Adriatic, where it was near zero. Wind stress is concentrated over the NE half of Adriatic in four bora jets, while wind stress is weak in the NW Adriatic. There is significant variation in wind stress mean and standard deviation structure over the northern Adriatic with each bora event. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗

USGS2006

Surface drifter derived circulation in the northern and middle Adriatic Sea: Response to wind regime and season

More than 120 satellite-tracked drifters were deployed in the northern and middle Adriatic (NMA) Sea between September 2002 and November 2003, with the purpose of studying the surface circulation at mesoscale to seasonal scale in relation to wind forcing, river runoff, and bottom topography. Pseudo-Eulerian and Lagrangian statistics were calculated from the low-pass-filtered drifter velocity data between September 2002 and December 2003. The structure of the mean circulation is determined with unprecedented high horizontal resolution by the new data. In particular, mean currents, velocity variance, and kinetic energy levels are shown to be maximal in the Western Adriatic Current (WAC). Separating data into seasons, we found that the mean kinetic energy is maximal in fall, with high values also in winter, while it is significantly weaker in summer. High-resolution Local Area Model Italy winds were used to relate the drifter velocities to the wind fields. The surface currents appear to be significantly influenced by the winds. The mean flow during the northeasterly bora regime shows an intensification of the across-basin recirculating currents. In addition, the WAC is strongly intensified both in intensity and in its offshore lateral extension. In the southeasterly sirocco regime, northward flow without recirculation dominates in the eastern half of the basin, while during northwesterly maestro the WAC is enhanced. Separating the data into low and high Po River discharge rates for low-wind conditions shows that the WAC and the velocity fluctuations in front of the Po delta are stronger for high Po River runoff. Lagrangian covariance, diffusivity, and integral time and space scales are larger in the along-basin direction and are maximal in the southern portion of the WAC. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗

USGS2006

Globally synchronous ice core volcanic tracers and abrupt cooling during the last glacial period

We perform a Monte Carlo pattern recognition analysis of the coincidence between three regional volcanic histories from ice coring of Greenland and Antarctica over the period 2 to 45 ka, using SO4 anomalies in Greenland and East Antarctica determined by continuous core chemistry, together with West Antarctic volcanic ash layers determined by remote optical borehole logging and core assays. We find that the Antarctic record of volcanism correlates with Glacial abrupt climate change at a 95% to >99.8% (???3??) significance level and that volcanic depositions at the three locations match at levels exceeding 3??, likely indicating that many common horizons represent single eruptive events which dispersed material world wide. These globally coincident volcanics were associated with abrupt cooling, often simultaneous with onsets or sudden intensifications of millennial cold periods. The striking agreement between sites implies that the consistency of current timescales obtained by isotopic and glaciological dating methods is better than estimated. Copyright 2006 by the American Geogphysical Union.

Journal of Geophysical Research D: Atmospheres↗

USGS2006

Atmospheric mercury emissions from substrates and fumaroles associated with three hydrothermal systems in the western United States

This paper quantifies atmospheric mercury (Hg) emissions from substrates and fumaroles associated with three hydrothermal systems: Lassen Volcanic Center, California (LVC); Yellowstone Caldera, Wyoming (YC); and Dixie Valley, Nevada (DV). Substrate Hg fluxes were measured using field chamber methods at thermal and nonthermal sites. The highest Hg fluxes (up to 541 ng m −2 h −1 ) were measured at thermal active areas. Fluxes from altered and unaltered nonthermal sites were <10 ng m −2 h −1 and were comparable to those measured in natural low-Hg background regions for YC and DV, and at LVC they were slightly elevated. Similarly, reactive gaseous mercury concentrations were higher in thermal active areas. Using a Geographic Information System framework for scaling, estimated area-average Hg emissions from substrates were 0.9–3.8 ng m −2 h −1 at LVC, 0.8–2.8 ng m −2 h −1 at YC, and 0.4–0.5 ng m −2 h −1 at DV. At LVC, nonthermal areas were the primary substrate source of atmospheric Hg (>98%). At YC, substrate Hg emissions were dominated (50 to 90%) by acidically altered thermal areas. Substrate emissions at DV were low and primarily from nonthermal areas (66% to 75%). Fumarole emissions at LVC (91–146 kg yr −1 ) and YC (0.18–1.6 kg yr −1 for Mud Volcano) were estimated by applying Hg:H 2 O and Hg:CO 2 ratios in hydrothermal gas samples to H 2 O and CO 2 emissions. Applying total area-average emissions from substrates and thermal features at LVC, YC, and DV to similar systems across the conterminous United States, yearly atmospheric Hg emissions from active hydrothermal systems are projected to be 1.3–2.1 Mg.

Journal of Geophysical Research D: Atmospheres↗

USGS2006

Cloud-to-ground lightning and surface rainfall in warm-season Florida thunderstorms

Relationships between cloud-to-ground (CG) lightning and surface rainfall have been examined in nine isolated, warm-season thunderstorms on the east coast of central Florida. CG flashes and the associated rain volumes were measured as a function of time in storm-centered reference frames that followed each storm over a network of rain gauges. Values of the storm-average rain volume per CG flash ranged from 0.70 ?? 104 to 6.4 ?? 104 m3/CG flash, with a mean (and standard deviation) of 2.6 ?? 104 ?? 2.1 ?? 104 m3/CG flash. Values of the rain volume concurrent with CG flashes ranged from 0.11 ?? 104 to 4.9 ?? 104 m3/CG flash with a mean of 2.1 ?? 104 ?? 2.0 ?? 104 m3/CG flash. The lag-time between the peak CG flash rate and the peak rainfall rate (using 5 min bins), and the results of a lag correlation analysis, show that surface rainfall tends to follow the lightning (positive lag) by up to 20 min in six storms. In one storm the rainfall preceded the lightning by 5 min, and two storms had nonsignificant lags. Values of the lagged rain volume concurrent with CG flashes ranged from 0.43 ?? 104 to 4.9 ?? 104 m3/CG flash, and the mean was 1.9 ?? 104 ?? 1.7 ?? 104 m3/CG flash. For the five storms that produced 12 or more flashes and had significant lags, a plot of the optimum lag time versus the total number of CG flashes shows a linear trend (R2 = 0.56). The number of storms is limited, but the lag results do indicate that large storms tend to have longer lags. A linear fit to the lagged rain volume vs. the number of concurrent CG flashes has a slope of 1.9 ?? 104 m3/CG flash (R2 = 0.83). We conclude that warm-season Florida thunderstorms produce a roughly constant rain volume per CG flash and that CG lightning can be used to estimate the location and intensity of convective rainfall in that weather regime. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗

USGS2005

Evolving force balance at Columbia Glacier, Alaska, during its rapid retreat

Changes in driving and resistive stresses play an essential role in governing the buoyancy forces that are important controls on the speed and irreversibility of tidewater glacier retreats. We describe changes in geometry, velocity, and strain rate and present a top-down force balance analysis performed over the lower reach of Columbia Glacier. Our analysis uses new measurements and estimates of basal topography and photogrammetric surface velocity measurements made between 1977 and 2001, while assuming depth-independent strain. Sensitivity tests show that the method is robust and insensitive to small changes in the calculation parameters. Spatial distributions of ice speed show little correspondence with driving stress. Instead, spatial patterns of ice speed exhibit a nonlinear correspondence with basal drag. Primary resistance to flow comes from basal drag, but lateral drag becomes increasingly more important throughout the retreat, which may account for observed increases in speed. Maximum basal drag is always located in a prominent constriction located ~12 km upstream from the preretreat terminus. Once the terminus retreated into deep water off the terminal moraine marking the modern maximum extent, the upstream location of this maximum basal drag helped to promote thinning and decrease effective pressure in the lower region by limiting replenishing ice flow from upstream. An increase in both ice velocity and calving resulted, initiating what appears to be an irreversible retreat. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

USGS2005

Modeling effects of bank friction and woody bank vegetation on channel flow and boundary shear stress in the Rio Puerco, New Mexico

[1] We have applied a physically based model for steady, horizontally uniform flow to calculate reach-averaged velocity and boundary shear-stress distributions in a natural stream with woody vegetation on the channel banks. The model calculates explicitly the form drag on woody plant stems and includes the effects of vegetation on the boundary shear stress, velocity, and turbulence fields. Average channel shapes, bed gradients, and shrub characteristics were determined for four long, internally similar reaches covering an 81 -km segment of the lower Rio Puerco, in north-central New Mexico. Hydraulic geometries for each of three flow events were determined from reach-average channel shapes, slopes, and measured silt-line heights above the bed. Model results show that friction on the lateral boundaries reduced the boundary shear stress in the center of the Rio Puerco channel by as much as 20%. In reaches with moderate to dense bank shrubs, perimeter-averaged boundary shear stresses were reduced by almost 40% during near bankfull flows. Model-calculated discharges for all three flow events indicate there was a loss of about 40% of the upstream discharge through the 81-km river segment. Infiltration rates estimated from discharge losses during the quasi-steady peak flows and areas of the bed and banks suggest rates of flow loss into the silty sand of the upper banks are about 8 times greater than rates of loss into the clayey silt-covered bed and lower banks.

Journal of Geophysical Research F: Earth Surface↗

USGS2006

Acoustic bed velocity and bed load dynamics in a large sand bed river

Development of a practical technology for rapid quantification of bed load transport in large rivers would represent a revolutionary advance for sediment monitoring and the investigation of fluvial dynamics. Measurement of bed load motion with acoustic Doppler current profiles (ADCPs) has emerged as a promising approach for evaluating bed load transport. However, a better understanding of how ADCP data relate to conditions near the stream bed is necessary to make the method practical for quantitative applications. In this paper, we discuss the response of ADCP bed velocity measurements, defined as the near-bed sediment velocity detected by the instrument's bottom-tracking feature, to changing sediment-transporting conditions in the lower Missouri River. Bed velocity represents a weighted average of backscatter from moving bed load particles and spectral reflections from the immobile bed. The ratio of bed velocity to mean bed load particle velocity depends on the concentration of the particles moving in the bed load layer, the bed load layer thickness, and the backscatter strength from a unit area of moving particles relative to the echo strength from a unit area of unobstructed bed. A model based on existing bed load transport theory predicted measured bed velocities from hydraulic and grain size measurements with reasonable success. Bed velocities become more variable and increase more rapidly with shear stress when the transport stage, defined as the ratio of skin friction to the critical shear stress for particle entrainment, exceeds a threshold of about 17. This transition in bed velocity response appears to be associated with the appearance of longer, flatter bed forms at high transport stages.

Journal of Geophysical Research F: Earth Surface↗