A study of Holocene floodplain particle size characteristics with special reference to palaeochannel infills from the upper Severn basin,Wales, UK

Multiple sedimentary units from floodplain reaches at Welshpool on the upper River Severn and at the confluence of the Afon Tanat and Afon Vyrnwy (mid‐Wales, UK) were examined to ascertain if they have distinctive particle size characteristics. Changes in particle size characteristics and their possible relationship to known human and climatic impacts are also discussed. Ellipse plots of particle size characteristics from the River Severn floodplain at Welshpool show that coarse‐grained outwash deposits can be clearly discriminated from channel margin or palaeochannel sediments. In contrast, at the Afon Tanat–Vyrnwy study reach, this discrimination is not seen so clearly. The relationships between age and particle size characteristics from the most sampled sedimentary environment, palaeochannel infills, were also examined. The data from the River Severn floodplain at Welshpool show that palaeochannel sediments reveal a gradual but clear increase in particle size from the mid‐ to late Holocene towards the present day. Sediments deposited in the period 90–160 years BP are markedly coarser. It is suggested that these changes may be related to the combined effect of land‐use changes, metal mining impacts and changes in flood frequency and magnitude that occurred at this time within the upper Severn basin. In contrast, the particle size characteristics of post Late Devensian/Early Holocene units from Tanat–Vyrnwy palaeochannels were random with no discernible age–size patterns. It is suggested that the non‐systematic grain size distribution may be due to the steeper valley gradients of the Tanat–Vyrnwy system (and by inference higher stream powers) and its relatively narrow valley form enabling more effective coupling between coarser outwash deposits found on and at the edges of hillslopes and the valley floor. Although the two study reaches have undergone comparable environmental change during the Holocene and lie in the piedmont zone of their catchments, palaeochannel units of the same age possess distinctly different characteristics. Intrinsic reach‐scale geomorphic factors would appear to preclude the uniform application of particle size characteristics to determine alluvial response to environmental change. Consequently, care needs to be applied to the use of such data for environmental discrimination because the phenomenon of equifinality means that a specific set of sediment characteristics is not necessarily exclusive to specific fluvial environments in either space or time. Copyright © 2001 John Wiley & Sons, Ltd.     

Modelling alluvial fan morphology

A mathematical model which estimates the scale-independent sediment surface profile of alluvial fans has been developed. This model utilizes a diffusive sediment transport model and an unsteady, radial flow, conservation relationship. These equations are approximately solved assuming a quasi-steady-state closure with appropriate modelling assumptions for two end member fan types: (1) fans where most of the fan surface is depositionally active (denoted here as ‘homogeneous’) and (2) fans characterized by channelling and sediment sorting processes. The fundamental result for these two fan types is a dimensionless sediment profile relationship which approximates most fan surfaces. The model suggests that the overall dimensionless morphology of alluvial fans is governed more by fundamental diffusion principles in sediment deposition than by individual environmental or basin characteristics. Additionally, this work potentially can be extended to model temporal variation in fan development. Preliminary comparison with alluvial fan profiles is reasonable, indicating that this model provides useful qualitative and quantitative information relating to alluvial fan process and morphology. Copyright © 1999 John Wiley & Sons, Ltd.     

A pore-scale numerical model for flow through porous media

A pore-scale numerical model based on Smoothed Particle Hydrodynamics (SPH) is described for modelling fluid flow phenomena in porous media. Originally developed for astrophysics applications, SPH is extended to model incompressible flows of low Reynolds number as encountered in groundwater flow systems. In this paper, an overview of SPH is provided and the required modifications for modelling flow through porous media are described, including treatment of viscosity, equation of state, and no-slip boundary conditions. The performance of the model is demonstrated for two-dimensional flow through idealized porous media composed of spatially periodic square and hexagonal arrays of cylinders. The results are in close agreement with solutions obtained using the finite element method and published solutions in the literature. Copyright © 1999 John Wiley & Sons, Ltd.     

CS2: A PIECEWISE-LINEAR MODEL FOR LARGE STRAIN CONSOLIDATION

This paper presents a piecewise-linear finite-difference model for one-dimensional large strain consolidation called CS2. CS2 is developed using a fixed Eulerian co-ordinate system and constitutive relationships which are defined by discrete data points. The model is dimensionless such that solutions are independent of the initial height of the compressible layer and the absolute magnitude of the hydraulic conductivity of the soil. The capability of CS2 is illustrated using four example problems involving small strain, large strain, self-weight, and non-linear constitutive relationships. In each case, the performance of the model is comparable to other available analytical and numerical solutions. Using CS2, correction factors are developed for the conventional Terzaghi theory which account for the effect of vertical strain on computed values by elapsed time and maximum excess pore pressure during consolidation. © 1997 John Wiley & Sons, Ltd.     

Flow separation on Zongo Glacier,Cordillera Real,Bolivia

Meltwaters collected from the proglacial stream escaping from Zongo Glacier (2·1 km²), Bolivia (16°S), have been monitored in order to analyse the internal drainage system of an Andean glacier. Electrical conductivity has been measured sporadically between February 1995 and March 1996, during 16 one-day field surveys, under various meteorological conditions in summer and winter. The mixing-model technique based on the electrical conductivity is used for a quantitative separation of discharge which is derived from continuous water level registration. Tracer experiments (mainly uranine dye and NaCl salt) have been carried out from March to June 1997 to obtain information about the internal drainage system. In the tropical Andes, accumulation only occurs in austral summer, whereas ablation occurs throughout the year and is higher during the accumulation season, between November and March. The assumptions involved in the use of mixing models for analysis of glacial drainage structure are applicable for tropical glaciers because glacial conduits do not suffer complete closure, and are permanently supplied by meltwaters, even in wintertime. Two components of discharge are separated: an englacial flow originating from surface meltwater which is routed without chemical enrichment, and offering low electrical conductivity; and a subglacial one routed in contact with bedrock or sediments showing high ionic concentrations. Electrical conductivity of meltwater varies diurnally, inversely to discharge fluctuations. According to this behaviour, total discharge is mainly formed by the englacial component. The drainage structures for englacial and subglacial flow have to be widely interconnected, as indicated by diurnal variations of the subglacial discharge. Comparison of hydrograph separation based on conductivity and on ¹⁸O isotope confirms that the subglacial flow is influenced by surface melting. A hydrograph separation of the subglacial flow is proposed, between a diurnal variable component, composed of water coming from the englacial network, and a base flow, which may vary seasonally. The dye tracing experiments confirm the drainage complexity of Zongo Glacier and demonstrate the interest of identifying three main drainage components. © 1998 John Wiley & Sons, Ltd.     

A fine balance: Accommodation dominated control of contemporaneous cool-carbonate shelf-edge clinoforms and tropical reef-margin trajectories,North Carnarvon Basin,Northwestern Australia

The concurrent development of a cool-carbonate Miocene clinoform system and the tropical reef which developed on its shelf in the North Carnarvon Basin is studied. The study, based on seismic interpretation and geometrical analysis, seeks to investigate how the architecture of the clinoforms develops in relation to the advance of the reef-margin, providing a proxy for discussing contemporaneous shoreline versus shelf-edge development. The progradation of the reef and shelf-edge often display a closely mirrored development, although the reef twice advances an order of two to three times the concurrent advance of the shelf-edge. The forced regression of the second advance, as compared to the normal regression during the first, is observed in proportionally higher input of sediment towards advance of the shelf-edge and toe, along with a gentler slope. The inability of the shelf-edge to keep pace with the reef-margin (and by proxy the shoreline) during lower accommodation/sedimentation is a result of the increased volume of sediment required to match reef-margin advance beyond the shelf-edge. Increased accommodation/sedimentation ratios promote higher trajectories where the volumes on shelf and slope are more balanced and the development more closely matched. The observed matched development of reef and shelf-edge during both limited and increased slope sedimentation, suggest that accommodation is the dominant control on the location and trajectory of both ‘shoreline’ and shelf-edge, and that excess sediment is deposited along the slope.     

The thermal history of the intergalactic medium

At redshifts z ≳2, most of the baryons reside in the smooth intergalactic medium which is responsible for the low column density Ly α forest. This photoheated gas follows a tight temperature–density relation which introduces a cut-off in the distribution of widths of the Ly α absorption lines ( b -parameters) as a function of column density. We have measured this cut-off in a sample of nine high-resolution, high signal-to-noise ratio quasar spectra and determined the thermal evolution of the intergalactic medium in the redshift range 2.0–4.5. At a redshift z ∼3, the temperature at the mean density shows a peak and the gas becomes nearly isothermal. We interpret this as evidence for the reionization of He  ii .     

Assessment of observational networks with the Representer Matrix Spectra method—application to a 3D coastal model of the Bay of Biscay

The development of coastal ocean modeling in the recent years has allowed an improved representation of the associated complex physics. Such models have become more realistic, to the point that they can now be used to design observation networks in coastal areas, with the idea that a “good” network is a network that controls model state error. To test this ability without performing data assimilation, we set up a technique called Representer Matrix Spectra (RMS) technique that combines the model state and observation error covariance matrices into a single scaled representer matrix. Examination of the spectrum and the eigenvectors of that matrix informs us on which model state error modes a network can detect and constrain amidst the observation error background. We applied our technique to a 3D coastal model in the Bay of Biscay, with a focus on mesoscale activity, and tested the performance of various altimetry networks and an in situ array deployment strategy. It appears that a single nadir altimeter is not efficient enough at capturing coastal mesoscale physics, while a wide swath altimeter would do a much better job. Testing various local in situ array configurations confirms that adding a current meter to a vertical temperature measurement array improves the detection of secondary variability modes, while shifting the array higher on the shelf break would obviously enhance the model constraint along the coast. The RMS technique is easily set up and used as a “black box,” but the utility of its results is maximized by previous knowledge of model state error physics. The technique provides both quantitative (eigenvalues) and qualitative (eigenvectors) tools to study and compare various network options. The qualitative approach is essential to discard possibly inconsistent modes.     

Sensitivity-guided reduction of parametric dimensionality for multi-objective calibration of watershed models

Problem complexity for watershed model calibration is heavily dependent on the number of parameters that can be identified during model calibration. This study investigates the use of global sensitivity analysis as a screening tool to reduce the parametric dimensionality of multi-objective hydrological model calibration problems while maximizing the information extracted from hydrological response data. This study shows that by expanding calibration problem formulations beyond traditional, statistical error metrics to also include metrics that capture indices or signatures of hydrological function, it is possible to reduce the complexity of calibration while maintaining high quality model predictions. The sensitivity-guided calibration is demonstrated using the Sacramento Soil Moisture Accounting (SAC-SMA) conceptual rainfall–runoff model of moderate complexity (i.e., up to 14 freely varying parameters). Using both statistical and hydrological metrics, optimization results demonstrate that parameters controlling at least 20% of the model output variance (through individual effects and interactions) should be included in the calibration process. This threshold generally yields 30–40% reductions in the number of SAC-SMA parameters requiring calibration – setting the others to a priori values – while maintaining high quality predictions. Two parameters are recommended to be calibrated in all cases (percent impervious area and lower zone tension water storage), three parameters are needed in drier watersheds (additional impervious area, riparian zone vegetation, and percent of percolation going to tension storage), and the lower zone parameters are crucial unless the watershed is very dry. Overall, this study demonstrates that a coupled, multi-objective sensitivity and calibration analysis better captures differences between watersheds during model calibration and serves to maximize the value of available watershed response time series. These contributions are particularly important given the ongoing development of more complex integrated models, which will require new tools to address the growing discrepancy between the information content of hydrological data and the number of model parameters that have to be estimated.