Form and growth of bars in a wandering gravel‐bed river

The changing form of developing alluvial river bars has rarely been studied in the field, especially in the context of the fixed, compound, mainly alternate gravel bars that are the major morphological feature of the wandering style. Century scale patterns of three‐dimensional growth and development, and the consequent scaling relations of such bars, are examined along the gravel‐bed reach of lower Fraser River, British Columbia, Canada. A retrospective view based on maps and aerial photographs obtained through the twentieth century shows that individual bars have a life history of about 100 years, except in certain, protected positions. A newly formed gravel bar quickly assumes its ultimate thickness and relatively quickly approaches its equilibrium length. Growth continues mainly by lateral accretion of unit bars, consistent with the lateral style of instability of the river. Bar growth is therefore allometric. Mature bars approach equilibrium dimensions and volume that scale with the overall size of the channel. Accordingly, the bars conform with several published criteria for the ultimate dimensions of alternate barforms. Sand bars, observed farther downstream, have notably different morphology. Fraser River presents a typical wandering channel planform, exhibiting elements of both meandered and low‐order braided channels. Hydraulic criteria to which the Fraser bars conform illustrate why this planform develops and persists. The modest rate of bed material transfer along the channel – typical of the wandering type – determines a century‐length time scale for bar development. This time scale is consistent with estimates that have been made for change of the macroform elements that determine the overall geometry of alluvial channels. Copyright © 2009 John Wiley & Sons, Ltd.     

Geophysical Finite-Element Simulation Tool (GeoFEST): Algorithms and Validation for Quasistatic Regional Faulted Crust Problems

GeoFEST (Geophysical Finite Element Simulation Tool) is a two- and three-dimensional finite element software package for the modeling of solid stress and strain in geophysical and other continuum domain applications. It is one of the featured high-performance applications of the NASA QuakeSim project. The program is targeted to be compiled and run on UNIX systems, and is running on diverse systems including sequential and message-passing parallel systems. Solution to the elliptical partial differential equations is obtained by finite element basis sampling, resulting in a sparse linear system primarily solved by conjugate gradient iteration to a tolerance level; on sequential systems a Crout factorization for the direct inversion of the linear system is also supported. The physics models supported include isotropic linear elasticity and both Newtonian and power-law viscoelasticity, via implicit quasi-static time stepping. In addition to triangular, quadrilateral, tetrahedral and hexahedral continuum elements, GeoFEST supports split-node faulting, body forces, and surface tractions. This software and related mesh refinement strategies have been validated on a variety of test cases with rigorous comparison to analytical solutions. These include a box-shaped domain with imposed motion on one surface, a pair of strike slip faults in stepover arrangement, and two community-agreed benchmark cases: a strike slip fault in an enclosing box, and a quarter-domain circular fault problem. Scientific applications of the code include the modeling of static and transient co- and post-seismic earth deformation, Earth response to glacial, atmospheric and hydrological loading, and other scenarios involving the bulk deformation of geologic media.