Channel bed adjustments following major aggradation in a steep headwater setting: findings from Oyabu Creek, Kyushu, Japan

A typhoon in 1993 induced major aggradation along Oyabu Creek, a steep, gravel bed mountain stream in Kyushu, Japan. Processes of sediment reworking are inferred from a 7-year monitoring program that measured adjustments to channel cross-sections, the longitudinal profile, and the extent/distribution of bedrock outcrops along a 3-km study reach. Over time, the reach adopted a riffle and pool structure, with notable increase in the area of exposed bedrock on the bed. This adjustment process was characterised by progressive reduction in sediment storage change per unit flow. The relaxation pathway following disturbance induced by the typhoon was shaped by the magnitude and frequency of subsequent rainfall events, the capacity of these events to transport available sediments, and physical linkages between reaches. Adjacent subreaches demonstrated differing relaxation pathways in response to these influences, induced by spatial and temporal variability in threshold conditions along the channel. Longer-term evidence indicates that responses to major disturbance, such as the 1993 typhoon, occur as ‘cycles’ of around 20-year duration. A relaxation period of 7 years is required to attain a quasi-equilibrium bed configuration and rate of sediment flux. The timeframe of cycles is considered to reflect changes to hillslope–channel bed coupling, marking the period required to generate sufficient sediment stores to reactivate phases of aggradation and subsequent degradation.     

A multi-tiered earthquake hazard model for Australia

Earthquakes result from tectonic processes, and their distribution is strongly influenced by large-scale geology and the tectonic stress field. However, earthquake hazard estimates, particularly ground motion recurrence, have traditionally been computed using source models based primarily on instrumental and historical seismicity. In areas of low to moderate seismicity such as Australia, large earthquakes commonly occur in areas which have experienced little or no recent activity, making it difficult to develop source models based solely on seismicity.

The seismotectonic model developed for Australia that is presented here (AUS5) is based on geology, geophysics, tectonics and seismicity. The model was developed using a number of tiers of information, so that new information can easily be incorporated. The information used includes, but is not limited to, tectonic provinces, basins and ranges, gravity, magnetic, topography, and seismicity, all on a regional scale. On a local scale, for a site-specific earthquake hazard study, active faulting can be incorporated to provide fault source zones.

An earthquake hazard map showing peak ground acceleration with a 10% chance of exceedance in 50 years for southeastern Australia using the geologically defined seismotectonic model AUS5 is presented as an indication of how the model performs.     

Depth variation of seismic source scaling relations: implications for earthquake hazard in southeastern Australia

Advances in earthquake data acquisition and processing techniques have allowed for improved quantification of source parameters for local Australian earthquakes. Until recently, only hypocentral locations and local magnitudes (M_L) had been determined routinely, with little attention given to the inversion of additional source parameters. The present study uses these new source data (e.g. seismic moment, stress drop, source dimensions) to further extend our understanding of seismicity and the continental stress regime of the Australian landmass and its peripheral regions.

Earthquake activity within Australia is typically low, and the proportion of small to large events (i.e. the b value) is also low. It is observed that average stress drops for southeastern Australian earthquakes appear to increase with seismic moment to relatively high levels, up to approximately 10 MPa for M_L 5.0 earthquakes. This is thought to be indicative of high compressive crustal stress, coupled with strong rocks and fault asperities. Furthermore, the data indicates that shallow focus earthquakes (shallower than 6 km) appear to produce lower than average stress drops than deeper earthquakes (between 6 and 20 km) with similar moment.

Recurrence estimates were obtained for a discrete seismogenic zone in southeastern Australia. Decreasing b values with increasing focal depth for this zone indicate that larger earthquakes (with high stress drops) tend to occur deeper in the crust. This may offer an explanation for the apparent increase of stress drop with hypocentral depth. Consequently, earthquake hazard estimates that assume a uniform Gutenburg–Richter distribution with depth (i.e. constant b value) may be too conservative and therefore slightly overestimate seismic hazard for surface sites in southeastern Australia.     

New observations on the geomorphology and origins of Mountain Lake,Virginia

Mountain Lake is the only natural lake of significance in the unglaciated southern Appalachian Highlands. It is located near the summit of Salt Pond Mountain, Giles County Virginia, at an elevation of 1177 m. It is underlain by Ordovician and Silurian non‐calcareous shale and sandstone of the Martinsburg, Juniata and Clinch formations. Historical (250 years) and sediment (6000 years) records suggest that the size of the lake has varied periodically. In the 1930s lake origin was proposed as due to valley damming by a lateral landslide (Hutchinson and Pickford, 1932) or damming by scree (Sharp, 1933). A later theory modified the landslide hypothesis to the primarily vertical collapse of a canyon feature in the Clinch (Parker et al., 1975). Fracture trace analysis now reveals a regional lineation feature associated with the lake. This feature is present surficially both downgradient from the lake to the northwest, and upgradient to the southeast. Sonar bathymetry and diver reconnaissance show it expressed as a (relatively sediment‐free) narrow open crevice in the deepest (33 m) portion of the lake, probably a fault. Hydrologic observation and resistivity suggest preferential water movement along this fracture, as well as leakage directly from the lake. The present study suggests conduit erosion within this feature and periodic vertical downsettling of overlying Clinch material as the primary mechanism of lake origin and water‐level fluctuations through time. Copyright © 2001 John Wiley & Sons, Ltd.     

Processing And Interpretation Of Electromagnetic Induction Array Data

The simultaneous nature of array data can be exploitedin electromagnetic induction studiesfor three general purposes.First, one or more reference sites can be usedto reduce bias, improve signal-to-noise ratios, and provide bettercontrol over source complications and coherent noisein estimates of MT impedances and otherEM transfer functions (TFs).Although a single good reference site can dramatically improveTF estimates, improvements due to multiple sites are often rather modest, because local noise is the limiting factor.Secondly, arrays allow for estimation of inter-station transfer functions, and maps of anomalous horizontal field variations. Relatively straightforward modifications to inversion codes wouldallow quantitative interpretation of these additionalconstraints on resistivity variations.Finally, with arrays it is possible to estimate the response of theEarth to a richer spectrum of external source excitations. In particular,the natural extension to the usual uniform source assumption implicit inthe MT method allows for three curl-free magnetic gradient sources.Quantitative interpretation of the response of a three-dimensionalEarth to these sources could provide additionalconstraints on large scale variations in crustal and uppermantle resistivity, and might help to overcome problemsdue to aliasing of near surface distortion in widely spaced MT data.     

Pool-fills: a window to palaeoflood history and response in bedrock-confined rivers

Channel‐scale sedimentary units associated with bedrock‐controlled riffle‐pool morphology are examined in detail along Sandy Creek gorge, an ephemeral stream in arid south‐eastern central Australia. Pool‐fills comprise cut‐and‐fill assemblages of poorly sorted sediments ranging in texture from muds to boulders. Five unit types are defined based on particle size, sedimentary structures, geometry and bounding surface character: (1) coarse‐grained bar platform; (2) fine‐grained bar supraplatform; (3) fine‐grained pool‐fill; (4) fine‐grained bench; and (5) modern pool‐fill. The last coarse‐grained unit currently lining the pools suggests an altered sedimentation style over the post‐settlement period (post‐ad 1860s). Situated at bedrock valley constrictions, pool‐fills are compared with other sedimentary units associated with recirculating currents: eddy bars and slackwater deposits. But only the fine‐grained bench units reflect eddy recirculation; the pool‐fills are principally forced‐bars associated with bedrock‐controlled or ‘forced’ riffle‐pool morphology. A late Holocene palaeoflood history is proposed based on radiocarbon ages from the pool‐fills: multiple phases of cut‐and‐fill activity were preceded by a superflood 3400–1900 years ago that eroded the pool‐fills to bedrock. The resilience of the pool‐fills was illustrated by the passage of a 1‐in‐100‐year flood in 1992, which caused only minor erosion. The presence of pool‐fills may provide a window to past phases of river activity that cannot be extracted from either historical records/observations or palaeoflood slackwater sediment analyses. The formation and sedimentary preservation potential of these landforms reflect a combination of hydraulic and structural influences, but the occurrence of high‐magnitude floods exerts the dominant control.     

Approaches to surface complexation modeling of Uranium(VI) adsorption on aquifer sediments

Uranium(VI) adsorption onto aquifer sediments was studied in batch experiments as a function of pH and U(VI) and dissolved carbonate concentrations in artificial groundwater solutions. The sediments were collected from an alluvial aquifer at a location upgradient of contamination from a former uranium mill operation at Naturita, Colorado (USA). The ranges of aqueous chemical conditions used in the U(VI) adsorption experiments (pH 6.9 to 7.9; U(VI) concentration 2.5 · 10⁻⁸ to 1 · 10⁻⁵ M; partial pressure of carbon dioxide gas 0.05 to 6.8%) were based on the spatial variation in chemical conditions observed in 1999-2000 in the Naturita alluvial aquifer. The major minerals in the sediments were quartz, feldspars, and calcite, with minor amounts of magnetite and clay minerals. Quartz grains commonly exhibited coatings that were greater than 10 nm in thickness and composed of an illite-smectite clay with occluded ferrihydrite and goethite nanoparticles. Chemical extractions of quartz grains removed from the sediments were used to estimate the masses of iron and aluminum present in the coatings. Various surface complexation modeling approaches were compared in terms of the ability to describe the U(VI) experimental data and the data requirements for model application to the sediments. Published models for U(VI) adsorption on reference minerals were applied to predict U(VI) adsorption based on assumptions about the sediment surface composition and physical properties (e.g., surface area and electrical double layer). Predictions from these models were highly variable, with results overpredicting or underpredicting the experimental data, depending on the assumptions used to apply the model. Although the models for reference minerals are supported by detailed experimental studies (and in ideal cases, surface spectroscopy), the results suggest that errors are caused in applying the models directly to the sediments by uncertain knowledge of: 1) the proportion and types of surface functional groups available for adsorption in the surface coatings; 2) the electric field at the mineral-water interface; and 3) surface reactions of major ions in the aqueous phase, such as Ca²⁺, Mg²⁺, HCO₃⁻, SO₄²⁻, H₄SiO₄, and organic acids. In contrast, a semi-empirical surface complexation modeling approach can be used to describe the U(VI) experimental data more precisely as a function of aqueous chemical conditions. This approach is useful as a tool to describe the variation in U(VI) retardation as a function of chemical conditions in field-scale reactive transport simulations, and the approach can be used at other field sites. However, the semi-empirical approach is limited by the site-specific nature of the model parameters.     

Evaluation of a pumping test of the Snake River Plain aquifer using axial-flow numerical modeling

The Snake River Plain aquifer in southeast Idaho is hosted in a thick sequence of layered basalts and interbedded sediments. The degree to which the layering impedes vertical flow has not been well understood, yet is a feature that may exert a substantial control on the movement of contaminants. An axial-flow numerical model, RADFLOW, was calibrated to pumping test data collected by a straddle-packer system deployed at 23 depth intervals in four observation wells to evaluate conceptual models and estimate properties of the Snake River Plain aquifer at the Idaho National Engineering and Environmental Laboratory. A delayed water-table response observed in intervals beneath a sediment interbed was best reproduced with a three-layer simulation. The results demonstrate the hydraulic significance of this interbed as a semi-confining layer. Vertical hydraulic conductivity of the sediment interbed was estimated to be about three orders of magnitude less than vertical hydraulic conductivity of the lower basalt and upper basalt units. The numerical model was capable of representing aquifer conceptual models that could not be represented with any single analytical technique. The model proved to be a useful tool for evaluating alternative conceptual models and estimating aquifer properties in this application. Electronic Publication