Physico-chemical controls on phosphorus cycling in two lowland streams. Part 1 - the water column

This paper investigates phosphorus (P) transport and transformation dynamics in two contrasting sub-catchments of the River Kennel, England. Samples were collected daily under baseflow and hourly under stormflow conditions using autosamplers for 2 years and analysed for a range of determinands (full P fractionation, suspended sediment (SS), cations, pH, alkalinity, temperature and oxygen). Concentrations of SRP, SUP, PP and SS were higher in the flashy River Enborne (means of 0.186, 0.071, 0.101 and 34 mg l(-1), respectively) than the groundwater-fed River Lambourn (0.079, 0.057, 0.028 and 9 mg l(-1), respectively). A seasonal trend in the daily P dataset was evident, with lower concentrations during intermediate flows and the spring (caused by a dilution effect and macrophyte uptake) than during baseflow conditions. However, in the hourly P dataset, highest concentrations were observed during storm events in the autumn and winter (reflecting higher scour with increased capacity to entrain particles). Storm events were more significant in contributing to the total P load in the River Enborne than the River Lambourn, especially during August to October, when dry antecedent conditions were observed in the catchment. Re-suspension of P-rich sediment that accumulated within the channel during summer low flows might account for these observations. It is suggested that a P-calcite co-precipitation mechanism was operating during summer in the River Lambourn, while adsorption by metal oxyhydroxide groups was an important mechanism controlling P fractionation in the River Enborne. The influence of flow conditions and channel storage/release mechanisms on P dynamics in these two lowland rivers is assessed. (C) 2004 Elsevier B.V. All rights reserved.

Hiatal surfaces from the Miocene Globigerina Limestone Formation of Malta: Biostratigraphy, sedimentology, trace fossils and early diagenesis

The Miocene Globigerina Limestone of the Maltese islands contains widespread omission surfaces with very different characteristics and origins. The terminal Lower Globigerina Limestone hardground (TLGLHg) formed during a period of falling sea level. Coccolith assemblages suggest shallowness. Sedimentary structures and trace fossil assemblages, indicate increasing frequency of storm events and erosional episodes, towards the surface. Calcite cementation which took place around Thalassinoides burrows and formed irregular nodules was followed by dissolution of aragonite. It is suggested that lithification was linked to microbial reactions involving organic matter. In contrast two later surfaces, the terminal Middle Globigerina Limestone omissionground (TMGLOg), which marks the Lower to Middle Miocene boundary, and the Fomm-ir-Rih local hardground (FiRLHg) both contain early diagenetic dolomite. Lithification took place in two phases. The dolomite is interpreted to have formed beneath the sea floor: it was subsequently exhumed and partially corroded as the precipitation of calcitic and phosphatic cements took place around burrows open to the circulation of sea water. (C) 2008 Elsevier B.V. All rights reserved.

Does reactive surface area depend on grain size? Results from pH 3, 25 circle C far-from-equilibrium flow-through dissolution experiments on anorthite and biotite

Laboratory determined mineral weathering rates need to be normalised to allow their extrapolation to natural systems. The principle normalisation terms used in the literature are mass, and geometric- and BET specific surface area (SSA). The purpose of this study was to determine how dissolution rates normalised to these terms vary with grain size. Different size fractions of anorthite and biotite ranging from 180-150 to 20-10 mu m were dissolved in pH 3, HCl at 25 degrees C in flow through reactors under far from equilibrium conditions. Steady state dissolution rates after 5376 h (anorthite) and 4992 h (biotite) were calculated from Si concentrations and were normalised to initial- and final- mass and geometric-, geometric edge- (biotite), and BET SSA. For anorthite, rates normalised to initial- and final-BET SSA ranged from 0.33 to 2.77 X 10(-10) mol(feldspar) m(-2) s(-1), rates normalised to initial- and final-geometric SSA ranged from 5.74 to 8.88 X 10(-10) mol(feldspar) m(-2) s(-1) and rates normalised to initial- and final-mass ranged from 0.11 to 1.65 mol(feldspar) g(-1) s(-1). For biotite, rates normalised to initial- and final-BET SSA ranged from 1.02 to 2.03 X 10(-12) mol(biotite) m(-2) s(-1), rates normalised to initial- and final-geometric SSA ranged from 3.26 to 16.21 X 10(-12) mol(biotite) m(-2) s(-1), rates normalised to initial- and final-geometric edge SSA ranged from 59.46 to 111.32 x 10(-12) mol(biotite) m(-2) s(-1) and rates normalised to initial- and final-mass ranged from 0.81 to 6.93 X 10(-12) mol(biotite) g(-1) s(-1). For all normalising terms rates varied significantly (p <= 0.05) with grain size. The normalising terms which gave least variation in dissolution rate between grain sizes for anorthite were initial BET SSA and initial- and final-geometric SSA. This is consistent with: (1) dissolution being dominated by the slower dissolving but area dominant non-etched surfaces of the grains and, (2) the walls of etch pits and other dissolution features being relatively unreactive. These steady state normalised dissolution rates are likely to be constant with time. Normalisation to final BET SSA did not give constant ratios across grain size due to a non-uniform distribution of dissolution features. After dissolution coarser grains had a greater density of dissolution features with BET-measurable but unreactive wall surface area than the finer grains. The normalising term which gave the least variation in dissolution rates between grain sizes for biotite was initial BET SSA. Initial- and final-geometric edge SSA and final BET SSA gave the next least varied rates. The basal surfaces dissolved sufficiently rapidly to influence bulk dissolution rate and prevent geometric edge SSA normalised dissolution rates showing the least variation. Simple modelling indicated that biotite grain edges dissolved 71-132 times faster than basal surfaces. In this experiment, initial BET SSA best integrated the different areas and reactivities of the edge and basal surfaces of biotite. Steady state dissolution rates are likely to vary with time as dissolution alters the ratio of edge to basal surface area. Therefore they would be more properly termed pseudo-steady state rates, only appearing constant because the time period over which they were measured (1512 h) was less than the time period over wich they would change significantly. (c) 2006 Elsevier Inc. All rights reserved.

Salt marsh morphodynamics: An investigation of tidal flows and marsh channel equilibrium

Although numerous field studies have evaluated flow and transport processes in salt marsh channels, the overall role of channels in delivering and removing material from salt marsh platforms is still poorly characterised. In this paper, we consider this issue based on a numerical hydrodynamic model for a prototype marsh system and on a field survey of the cross-sectional geometry of a marsh channel network. Results of the numerical simulations indicate that the channel transfers approximately three times the volume of water that would be estimated from mass balance considerations alone. Marsh platform roughness exerts a significant influence on the partitioning of discharge between the channel and the marsh platform edge, alters flow patterns on the marsh platform due to its effects on channel-to-platform transfer and also controls the timing of peak discharge relative to marsh-edge overtopping. Although peak channel discharges and velocities are associated with the flood tide and marsh inundation, a larger volume of water is transferred by the channel during ebb flows, a portion of which transfer takes place after the tidal height is below the marsh platform. Detailed surveys of the marsh channels crossing a series of transects at Upper Stiffkey Marsh, north Norfolk, England, show that the total channel cross-sectional area increases linearly with catchment area in the inner part of the marsh, which is consistent with the increase in shoreward tidal prism removed by the channels. Toward the marsh edge, however, a deficit in the total cross-sectional area develops, suggesting that discharge partitioning between the marsh channels and the marsh platform edge may also be expressed in the morphology of marsh channel systems.

Characterization of mineral surfaces using FIB and TEM: a case study of naturally weathered alkali feldspars

Using a focused ion beam (FIB) instrument, electron-transparent samples (termed foils) have been cut from the naturally weathered surfaces of perthitic alkali feldspars recovered from soils overlying the Shap granite, northwest England. Characterization of these foils by transmission electron microscopy (TEM) has enabled determination of the crystallinity and chemical composition of near-surface regions of the feldspar and an assessment of the influence of intragranular microtextures on the microtopography of grain surfaces and development of etch pits. Damage accompanying implantation of the 30 kV Ga+ ions used for imaging and deposition of protective platinum prior to ion milling creates amorphous layers beneath outer grain surfaces, but can be overcome by coating grains with > 85 nm of gold before FIB work. The sidewalls of the foil and feldspar surrounding original voids are also partially amorphized during later stages of ion milling. No evidence was found for the presence of amorphous or crystalline weathering products or amorphous "leached layers" immediately beneath outer grain surfaces. The absence of a leached layer indicates that chemical weathering of feldspar in the Shap soils is stoichiometric, or if non-stoichiometric, either the layer is too thin to resolve by the TEM techniques used (i.e., <=similar to 2.5 nm) or an insufficient proportion of ions have been leached from near-surface regions so that feldspar crystallinity is maintained. No evidence was found for any difference in the mechanisms of weathering where a microbial filament rests on the feldspar surface. Sub-micrometer-sized steps on the grain surface have formed where subgrains and exsolution lamellae have influenced the propagation of fractures during physical weathering, whereas finer scale corrugations form due to compositional or strain-related differences in dissolution rates of albite platelets and enclosing tweed orthoclase. With progressive weathering, etch pits that initiated at the grain surface extend into grain interiors as etch tubes by exploiting preexisting networks of nanopores that formed during the igneous history of the grain. The combination of FIB and TEM techniques is an especially powerful way of exploring mechanisms of weathering within the "internal zone" beneath outer grain surfaces, but results must be interpreted with caution owing to the ease with which artifacts can be created by the high-energy ion and electron beams used in the preparation and characterization of the foils.

Laurentide ice streaming on the Canadian Shield: A conflict with the soft-bedded ice stream paradigm?

Elucidating the controls on the location and vigor of ice streams is crucial to understanding the processes that lead to fast disintegration of ice flows and ice sheets. In the former North American Laurentide ice sheet, ice stream occurrence appears to have been governed by topographic troughs or areas of soft-sediment geology. This paper reports robust evidence of a major paleo-ice stream over the northwestern Canadian Shield, an area previously assumed to be incompatible with fast ice flow because of the low relief and relatively hard bedrock. A coherent pattern of subglacial bedforms (drumlins and megascalle glacial lineations) demarcates the ice stream flow set, which exhibits a convergent onset zone, a narrow main trunk with abrupt lateral margins, and a lobate terminus. Variations in bedform elongation ratio within the flow set match theoretical expectations of ice velocity. In the center of the ice stream, extremely parallel megascalle glacial lineations tens of kilometers long with elongation ratios in excess of 40:1 attest to a single episode of rapid ice flow. We conclude that while bed properties are likely to be influential in determining the occurrence and vigor of ice streams, contrary to established views, widespread soft-bed geology is not an essential requirement for those ice streams without topographic control. We speculate that the ice stream acted as a release valve on ice-sheet mass balance and was initiated by the presence of a proglacial lake that destabilized the ice-sheet margin and propagated fast ice flow through a series of thermomechanical feedbacks involving ice flow and temperature.

Impacts of climate change on in-stream nitrogen in a lowland chalk stream: An appraisal of adaptation strategies

The impacts of climate change on nitrogen (N) in a lowland chalk stream are investigated using a dynamic modelling approach. The INCA-N model is used to simulate transient daily hydrology and water quality in the River Kennet using temperature and precipitation scenarios downscaled from the General Circulation Model (GCM) output for the period 1961-2100. The three GCMs (CGCM2, CSIRO and HadCM3) yield very different river flow regimes with the latter projecting significant periods of drought in the second half of the 21st century. Stream-water N concentrations increase over time as higher temperatures enhance N release from the soil, and lower river flows reduce the dilution capacity of the river. Particular problems are shown to occur following severe droughts when N mineralization is high and the subsequent breaking of the drought releases high nitrate loads into the river system. Possible strategies for reducing climate-driven N loads are explored using INCA-N. The measures include land use change or fertiliser reduction, reduction in atmospheric nitrate and ammonium deposition, and the introduction of water meadows or connected wetlands adjacent to the river. The most effective strategy is to change land use or reduce fertiliser use, followed by water meadow creation, and atmospheric pollution controls. Finally, a combined approach involving all three strategies is investigated and shown to reduce in-stream nitrate concentrations to those pre-1950s even under climate change. (c) 2006 Elsevier B.V. All rights reserved.

Integrated modelling of climate change impacts on water resources and quality in a lowland catchment: River Kennet, UK

An integrated approach to climate change impact assessment is explored by linking established models of regional climate (SDSM), water resources (CATCHMOD) and water quality (INCA) within a single framework. A case study of the River Kennet illustrates how the system can be used to investigate aspects of climate change uncertainty, deployable water resources, and water quality dynamics in upper and lower reaches of the drainage network. The results confirm the large uncertainty in climate change scenarios and freshwater impacts due to the choice of general circulation model (GCM). This uncertainty is shown to be greatest during summer months as evidenced by large variations between GCM-derived projections of future tow river flows, deployable yield from groundwater, severity of nutrient flushing episodes, and Long-term trends in surface water quality. Other impacts arising from agricultural land-use reform or delivery of EU Water Framework Directive objectives under climate change could be evaluated using the same framework. (c) 2006 Elsevier B.V. All rights reserved.

Solvent Effects on the Formation of Nanoparticles and Multilayered Coatings Based on Hydrogen-Bonded Interpolymer Complexes of Poly(acrylic acid) with Homo- and Copolymers of N-Vinyl Pyrrolidone

The formation of hydrogen-bonded interpolymer complexes between poly(acrylic acid) and poly(N-vinyl pyrrolidone) as well as amphiphilic copolymers of N-vinyl pyrrolidone with vinyl propyl ether has been studied in aqueous and organic solutions. It was demonstrated that introduction of vinyl propyl ether units into the macromolecules of the nonionic polymer enhances their ability to form complexes in aqueous solutions due to more significant contribution of hydrophobic effects. The complexation was found to be a multistage process that involves the formation of primary polycomplex particles, which further aggregate to form spherical nanoparticles. Depending on the environmental factors (pH, solvent nature), these nanoparticles may either form stable colloidal solutions or undergo further aggregation, resulting in precipitation of interpolymer complexes. In organic solvents, the intensity of complex formation increases in the following order: methanol < ethanol < isopropanol < dioxane. The multilayered coatings were developed using layer-by-layer deposition of interpolymer complexes on glass surfaces. It was demonstrated that the solvent nature affects the efficiency of coating deposition.

The three-dimensional vortical nature of atmospheric and oceanic turbulent flows

Using a novel numerical method at unprecedented resolution, we demonstrate that structures of small to intermediate scale in rotating, stratified flows are intrinsically three-dimensional. Such flows are characterized by vortices (spinning volumes of fluid), regions of large vorticity gradients, and filamentary structures at all scales. It is found that such structures have predominantly three-dimensional dynamics below a horizontal scale LLR, where LR is the so-called Rossby radius of deformation, equal to the characteristic vertical scale of the fluid H divided by the ratio of the rotational and buoyancy frequencies f/N. The breakdown of two-dimensional dynamics at these scales is attributed to the so-called "tall-column instability" [D. G. Dritschel and M. de la Torre Juárez, J. Fluid. Mech. 328, 129 (1996)], which is active on columnar vortices that are tall after scaling by f/N, or, equivalently, that are narrow compared with LR. Moreover, this instability eventually leads to a simple relationship between typical vertical and horizontal scales: for each vertical wave number (apart from the vertically averaged, barotropic component of the flow) the average horizontal wave number is equal to f/N times the vertical wave number. The practical implication is that three-dimensional modeling is essential to capture the behavior of rotating, stratified fluids. Two-dimensional models are not valid for scales below LR. ©1999 American Institute of Physics.

Acoustic scattering by mildly rough unbounded surfaces in three dimensions

For a nonlocally perturbed half- space we consider the scattering of time-harmonic acoustic waves. A second kind boundary integral equation formulation is proposed for the sound-soft case, based on a standard ansatz as a combined single-and double-layer potential but replacing the usual fundamental solution of the Helmholtz equation with an appropriate half- space Green's function. Due to the unboundedness of the surface, the integral operators are noncompact. In contrast to the two-dimensional case, the integral operators are also strongly singular, due to the slow decay at infinity of the fundamental solution of the three-dimensional Helmholtz equation. In the case when the surface is sufficiently smooth ( Lyapunov) we show that the integral operators are nevertheless bounded as operators on L-2(Gamma) and on L-2(Gamma G) boolean AND BC(Gamma) and that the operators depend continuously in norm on the wave number and on G. We further show that for mild roughness, i.e., a surface G which does not differ too much from a plane, the boundary integral equation is uniquely solvable in the space L-2(Gamma) boolean AND BC(Gamma) and the scattering problem has a unique solution which satisfies a limiting absorption principle in the case of real wave number.

On the initiation of surface waves by turbulent shear flow

An analytical model is developed for the initial stage of surface wave generation at an air-water interface by a turbulent shear flow in either the air or in the water. The model treats the problem of wave growth departing from a flat interface and is relevant for small waves whose forcing is dominated by turbulent pressure fluctuations. The wave growth is predicted using the linearised and inviscid equations of motion, essentially following Phillips [Phillips, O.M., 1957. On the generation of waves by turbulent wind. J. Fluid Mech. 2, 417-445], but the pressure fluctuations that generate the waves are treated as unsteady and related to the turbulent velocity field using the rapid-distortion treatment of Durbin [Durbin, P.A., 1978. Rapid distortion theory of turbulent flows. PhD thesis, University of Cambridge]. This model, which assumes a constant mean shear rate F, can be viewed as the simplest representation of an oceanic or atmospheric boundary layer. For turbulent flows in the air and in the water producing pressure fluctuations of similar magnitude, the waves generated by turbulence in the water are found to be considerably steeper than those generated by turbulence in the air. For resonant waves, this is shown to be due to the shorter decorrelation time of turbulent pressure in the air (estimated as proportional to 1/Gamma), because of the higher shear rate existing in the air flow, and due to the smaller length scale of the turbulence in the water. Non-resonant waves generated by turbulence in the water, although being somewhat gentler, are still steeper than resonant waves generated by turbulence in the air. Hence, it is suggested that turbulence in the water may have a more important role than previously thought in the initiation of the surface waves that are subsequently amplified by feedback instability mechanisms.

Conservative upwind difference schemes for the Euler equations for real gas flows in a duct

In a recent paper [P. Glaister, Conservative upwind difference schemes for compressible flows in a Duct, Comput. Math. Appl. 56 (2008) 1787–1796] numerical schemes based on a conservative linearisation are presented for the Euler equations governing compressible flows of an ideal gas in a duct of variable cross-section, and in [P. Glaister, Conservative upwind difference schemes for compressible flows of a real gas, Comput. Math. Appl. 48 (2004) 469–480] schemes based on this philosophy are presented for real gas flows with slab symmetry. In this paper we seek to extend these ideas to encompass compressible flows of real gases in a duct. This will incorporate the handling of additional terms arising out of the variable geometry and the non-ideal nature of the gas.