Towards an improved understanding of the nitrate dynamics in lowland, permeable river-systems: Applications of INCA-N

The Integrated Catchment Model of Nitrogen (INCA-N) was applied to the Lambourn and Pang river-systems to integrate current process-knowledge and available-data to test two hypotheses and thereby determine the key factors and processes controlling the movement of nitrate at the catchment-scale in lowland, permeable river-systems: (i) that the in-stream nitrate concentrations were controlled by two end-members only: groundwater and soil-water, and (ii) that the groundwater was the key store of nitrate in these river-systems. Neither hypothesis was proved true or false. Due to equifinality in the model structure and parameters at least two alternative models provided viable explanations for the observed in-stream nitrate concentrations. One model demonstrated that the seasonal-pattern in the stream-water nitrate concentrations was controlled mainly by the mixing of ground- and soil-water inputs. An alternative model demonstrated that in-stream processes were important. It is hoped further measurements of nitrate concentrations made in the catchment soil- and ground-water and in-stream may constrain the model and help determine the correct structure, though other recent studies suggest that these data may serve only to highlight the heterogeneity of the system. Thus when making model-based assessments and forecasts it is recommend that all possible models are used, and the range of forecasts compared. In this study both models suggest that cereal production contributed approximately 50% the simulated in-stream nitrate toad in the two catchments, and the point-source contribution to the in-stream load was minimal. (c) 2006 Elsevier B.V. All rights reserved.

A numerical study of the probe method

The goal of this work is the numerical realization of the probe method suggested by Ikehata for the detection of an obstacle D in inverse scattering. The main idea of the method is to use probes in the form of point source (., z) with source point z to define an indicator function (I) over cap (z) which can be reconstructed from Cauchy data or far. eld data. The indicator function boolean AND (I) over cap (z) can be shown to blow off when the source point z tends to the boundary aD, and this behavior can be used to find D. To study the feasibility of the probe method we will use two equivalent formulations of the indicator function. We will carry out the numerical realization of the functional and show reconstructions of a sound-soft obstacle.

Probing the large-scale topology of the heliospheric magnetic field using Jovian electrons

Jupiter���s magnetosphere acts as a point source of near-relativistic electrons within the heliosphere. In this study, three solar cycles of Jovian electron data in near-Earth space are examined. Jovian electron intensity is found to peak for an ideal Parker spiral connection, but with considerable spread about this point. Assuming the peak in Jovian electron counts indicates the best magnetic connection to Jupiter, we ���nd a clear trend for fast and slow solar wind to be over- and under-wound with respect to the ideal Parker spiral, respectively. This is shown to be well explained in terms of solar wind stream interactions. Thus, modulation of Jovian electrons by corotating interaction regions (CIRs) may primarily be the result of changing magnetic connection, rather than CIRs acting as barriers to cross-���eld diffusion. By using Jovian electrons to remote sensing magnetic connectivity with Jupiter���s magnetosphere, we suggest that they provide a means to validate solar wind models between 1 and 5 AU, even when suitable in situ solar wind observations are not available. Furthermore, using Jovian electron observations as probes of heliospheric magnetic topology could provide insight into heliospheric magnetic ���eld braiding and turbulence, as well as any systematic under-winding of the heliospheric magnetic ���eld relative to the Parker spiral from footpoint motion of the magnetic ���eld.

Sewage effluent clean-up reduces phosphorus but not phytoplankton in lowland chalk stream (River Kennet, UK) impacted by water mixing from adjacent canal

Information is provided on phosphorus in the River Kennet and the adjacent Kennet and Avon Canal in southern England to assess their interactions and the changes following phosphorus reductions in sewage treatment work (STW) effluent inputs. A step reduction in soluble reactive phosphorus (SRP) concentration within the effluent (5 to 13 fold) was observed from several STWs discharging to the river in the mid-2000s. This translated to over halving of SRP concentrations within the lower Kennet. Lower Kennet SRP concentrations change from being highest under base-flow to highest under storm-flow conditions. This represented a major shift from direct effluent inputs to a within-catchment source dominated system characteristic of the upper part to the catchment. Average SRP concentrations in the lower Kennet reduced over time towards the target for good water quality. Critically, there was no corresponding reduction in chlorophyll-a concentration, the waters remaining eutrophic when set against standards for lakes. Following the up gradient input of the main water and SRP source (Wilton Water), SRP concentrations in the canal reduced down gradient to below detection limits at times near its junction with the Kennet downstream. However, chlorophyll concentrations in the canal were in an order of magnitude higher than in the river. This probably resulted from long water residence times and higher temperatures promoting progressive algal and suspended sediment generations that consumed SRP. The canal acted as a point source for sediment, algae and total phosphorus to the river especially during the summer months when boat traffic disturbed the canal's bottom sediments and the locks were being regularly opened. The short-term dynamics of this transfer was complex. For the canal and the supply source at Wilton Water, conditions remained hypertrophic when set against standards for lakes even when SRP concentrations were extremely low.

Switching on a two-dimensional time-harmonic scalar wave in the presence of a diffracting edge

This paper concerns the switching on of two-dimensional time-harmonic scalar waves. We ���rst review the switch-on problem for a point source in free space, then proceed to analyse the analogous problem for the diffraction of a plane wave by a half-line (the ���Sommerfeld problem���), determining in both cases the conditions under which the ���eld is well-approximated by the solution of the corresponding frequency domain problem. In both cases the rate of convergence to the frequency domain solution is found to be dependent on the strength of the singularity on the leading wavefront. In the case of plane wave diffraction at grazing incidence the frequency domain solution is immediately attained along the shadow boundary after the arrival of the leading wavefront. The case of non-grazing incidence is also considered.

Catchment phosphorous losses: an export coefficient modelling approach with scenario analysis for water management

This paper presents the development of an export coefficient model to characterise the rates and sources of P export from land to water in four reservoir systems located in a semi-arid rural region in southern of Portugal. The model was developed to enable effective management of these important water resource systems under the EU Water Framework Directive. This is the first time such an approach has been fully adapted for the semi-arid systems typical of Mediterranean Europe. The sources of P loading delivered to each reservoir from its catchment were determined and scenario analysis was undertaken to predict the likely impact of catchment management strategies on the scale of rate of P loading delivered to each water body from its catchment. The results indicate the importance of farming and sewage treatment works/collective septic tanks discharges as the main contributors to the total diffuse and point source P loading delivered to the reservoirs, respectively. A reduction in the total P loading for all study areas would require control of farming practices and more efficient removal of P from human wastes prior to discharge to surface waters. The scenario analysis indicates a strategy based solely on reducing the agricultural P surplus may result in only a slow improvement in water quality, which would be unlikely to support the generation of good ecological status in reservoirs. The model application indicates that a reduction of P-inputs to the reservoirs should first focus on reducing P loading from sewage effluent discharges through the introduction of tertiary treatment (P-stripping) in all major residential areas. The fully calibrated export coefficient modelling approach transferred well to semi-arid regions, with the only significant limitation being the availability of suitable input data to drive the model. Further studies using this approach in semi-arid catchments are now needed to increase the knowledge of nutrient export behaviours in semi-arid regions.

A time-domain probe method for three-dimensional rough surface reconstructions

The task of this paper is to develop a Time-Domain Probe Method for the reconstruction of impenetrable scatterers. The basic idea of the method is to use pulses in the time domain and the time-dependent response of the scatterer to reconstruct its location and shape. The method is based on the basic causality principle of timedependent scattering. The method is independent of the boundary condition and is applicable for limited aperture scattering data. In particular, we discuss the reconstruction of the shape of a rough surface in three dimensions from time-domain measurements of the scattered field. In practise, measurement data is collected where the incident field is given by a pulse. We formulate the time-domain fieeld reconstruction problem equivalently via frequency-domain integral equations or via a retarded boundary integral equation based on results of Bamberger, Ha-Duong, Lubich. In contrast to pure frequency domain methods here we use a time-domain characterization of the unknown shape for its reconstruction. Our paper will describe the Time-Domain Probe Method and relate it to previous frequency-domain approaches on sampling and probe methods by Colton, Kirsch, Ikehata, Potthast, Luke, Sylvester et al. The approach significantly extends recent work of Chandler-Wilde and Lines (2005) and Luke and Potthast (2006) on the timedomain point source method. We provide a complete convergence analysis for the method for the rough surface scattering case and provide numerical simulations and examples.

A study on orthogonality sampling.

The goal of this paper is to study and further develop the orthogonality sampling or stationary waves algorithm for the detection of the location and shape of objects from the far field pattern of scattered waves in electromagnetics or acoustics. Orthogonality sampling can be seen as a special beam forming algorithm with some links to the point source method and to the linear sampling method. The basic idea of orthogonality sampling is to sample the space under consideration by calculating scalar products of the measured far field pattern , with a test function for all y in a subset Q of the space , m = 2, 3. The way in which this is carried out is important to extract the information which the scattered fields contain. The theoretical foundation of orthogonality sampling is only partly resolved, and the goal of this work is to initiate further research by numerical demonstration of the high potential of the approach. We implement the method for a two-dimensional setting for the Helmholtz equation, which represents electromagnetic scattering when the setup is independent of the third coordinate. We show reconstructions of the location and shape of objects from measurements of the scattered field for one or several directions of incidence and one or many frequencies or wave numbers, respectively. In particular, we visualize the indicator function both with the Dirichlet and Neumann boundary condition and for complicated inhomogeneous media.

Understanding lake and catchment history as a tool for integrated lake management

Sustainable lake management for nutrient-enriched lakes must be underpinned by an understanding of both the functioning of the lake, and the origins of changes in nutrient loading from the catchment. To date, limnologists have tended to focus on studying the impact of nutrient enrichment on the lake biota, and the dynamics of nutrient cycling between the water column, biota and sediments within the lake. Relatively less attention has been paid to understanding the specific origins of nutrient loading from the catchment and nutrient transport pathways linking the lake to its catchment. As such, when devising catchment management strategies to reduce nutrient loading on enriched lakes, assumptions have been made regarding the relative significance of non-point versus point sources in the catchment. These are not always supported by research conducted on catchment nutrient dynamics in other fields of freshwater science. Studies on nutrient enrichment in lakes need to take account of the history of catchment use and management specific to each lake in order to devise targeted and sustainable management strategies to reduce nutrient loading to enriched lakes. Here a modelling approach which allows quantification of the relative contribution of nutrients from each specific point and non-point catchment source over the course of catchment history is presented. The approach has been applied to three contrasting catchments in the U.K. for the period 1931 to present. These are the catchment of Slapton Ley in south Devon, the River Esk in Cumbria and the Deben Estuary in Suffolk. Each catchment showed marked variations in the nature and intensity of land use and management. The model output quantifies the relative importance of point source versus non-point livestock and land use sources in each of the catchments, and demonstrates the necessity for an understanding of site-specific catchment history in devising suitable management strategies for the reduction of nutrient loading on enriched lakes.

Steady state and dynamic modelling of nitrogen in the River Kennet: impacts of land use change since the 1930s

Steady state and dynamic models have been developed and applied to the River Kennet system. Annual nitrogen exports from the land surface to the river have been estimated based on land use from the 1930s and the 1990s. Long term modelled trends indicate that there has been a large increase in nitrogen transport into the river system driven by increased fertiliser application associated with increased cereal production, increased population and increased livestock levels. The dynamic model INCA Integrated Nitrogen in Catchments. has been applied to simulate the day-to-day transport of N from the terrestrial ecosystem to the riverine environment. This process-based model generates spatial and temporal data and reproduces the observed instream concentrations. Applying the model to current land use and 1930s land use indicates that there has been a major shift in the short term dynamics since the 1930s, with increased river and groundwater concentrations caused by both non-point source pollution from agriculture and point source discharges. ���

Nitrogen speciation and phosphorus fractionation dynamics in a lowland chalk catchment

A detailed analysis of temporal and spatial trends in nitrogen (N) speciation and phosphorus (P) fractionation in the Wylye, a lowland Chalk sub-catchment of the Hampshire Avon, UK is presented, identifying the sources contributing to nutrient enrichment, and temporal variability in the fractionation of nutrients in transit from headwaters to lower reaches of the river. Samples were collected weekly from ten monitoring stations with daily sampling at three further sites over one year, and monthly inorganic N and total reactive P (TRP) concentrations at three of the ten weekly monitoring stations over a ten year period are also presented. The data indicate significant daily and seasonal variation in nutrient fractionation in the water column, resulting from plant uptake of dissolved organic and inorganic nutrient fractions in the summer months, increased delivery of both N and P from diffuse sources in the autumn to winter period and during high flow events, and lack of dilution of point source discharges to the Wylye from septic tank, small package Sewage Treatment Works (STW) and urban Waste Water Treatment Works (WwTW) during the summer low flow period. Weekly data show that contributing source areas vary along the river with headwater N and P strongly influenced by diffuse inorganic N and particulate P fluxes, and SRP and organic-rich point source contributions from STW and WwTW having a greater influence in the lower reaches. Long-term data show a decrease in TRP concentrations at all three monitoring stations, with the most pronounced decrease occurring downstream from Warminster WwTW, following the introduction of P stripping at the works in 2001. Inorganic N demonstrates no statistically significant change over the ten year period of record in the rural headwaters, but an increase in the lower reaches downstream from the WwTW which may be due to urban expansion in the lower catchment.

The interactive responses of water quality and hydrology to changes in multiple stressors, and implications for the long-term effective management of phosphorus

Soluble reactive phosphorus (SRP) plays a key role in eutrophication, a global problem decreasing habitat quality and in-stream biodiversity. Mitigation strategies are required to prevent SRP fluxes from exceeding critical levels, and must be robust in the face of potential changes in climate, land use and a myriad of other influences. To establish the longevity of these strategies it is therefore crucial to consider the sensitivity of catchments to multiple future stressors. This study evaluates how the water quality and hydrology of a major river system in the UK (the River Thames) respond to alterations in climate, land use and water resource allocations, and investigates how these changes impact the relative performance of management strategies over an 80-year period. In the River Thames, the relative contributions of SRP from diffuse and point sources vary seasonally. Diffuse sources of SRP from agriculture dominate during periods of high runoff, and point sources during low flow periods. SRP concentrations rose under any future scenario which either increased a) surface runoff or b) the area of cultivated land. Under these conditions, SRP was sourced from agriculture, and the most effective single mitigation measures were those which addressed diffuse SRP sources. Conversely, where future scenarios reduced flow e.g. during winters of reservoir construction, the significance of point source inputs increased, and mitigation measures addressing these issues became more effective. In catchments with multiple point and diffuse sources of SRP, an all-encompassing effective mitigation approach is difficult to achieve with a single strategy. In order to attain maximum efficiency, multiple strategies might therefore be employed at different times and locations, to target the variable nature of dominant SRP sources and pathways.

Efficient calculation of the green function for acoustic propagation above a homogeneous impedance plane

This paper is concerned with the problem of propagation from a monofrequency coherent line source above a plane of homogeneous surface impedance. The solution of this problem occurs in the kernel of certain boundary integral equation formulations of acoustic propagation above an impedance boundary, and the discussion of the paper is motivated by this application. The paper starts by deriving representations, as Laplace-type integrals, of the solution and its first partial derivatives. The evaluation of these integral representations by Gauss-Laguerre quadrature is discussed, and theoretical bounds on the truncation error are obtained. Specific approximations are proposed which are shown to be accurate except in the very near field, for all angles of incidence and a wide range of values of surface impedance. The paper finishes with derivations of partial results and analogous Laplace-type integral representations for the case of a point source.

Preferential dust sources: a geomorphological classification designed and for use in global aerosol models

We present a simple theoretical land-surface classification that can be used to determine the location and temporal behavior of preferential sources of terrestrial dust emissions. The classification also provides information about the likely nature of the sediments, their erodibility and the likelihood that they will generate emissions under given conditions. The scheme is based on the dual notions of geomorphic type and connectivity between geomorphic units. We demonstrate that the scheme can be used to map potential modern-day dust sources in the Chihuahuan Desert, the Lake Eyre Basin and the Taklamakan. Through comparison with observed dust emissions, we show that the scheme provides a reasonable prediction of areas of emission in the Chihuahuan Desert and in the Lake Eyre Basin. The classification is also applied to point source data from the Western Sahara to enable comparison of the relative importance of different land surfaces for dust emissions. We indicate how the scheme could be used to provide an improved characterization of preferential dust sources in global dust-cycle models.

A self-similar expansion model for use in solar wind transient propogation studies

Since the advent of wide-angle imaging of the inner heliosphere, a plethora of techniques have been developed to investigate the three-dimensional structure and kinematics of solar wind transients, such as coronal mass ejections, from their signatures in single- and multi-spacecraft imaging observations. These techniques, which range from the highly complex and computationally intensive to methods based on simple curve ���tting, all have their inherent advantages and limitations. In the analysis of single-spacecraft imaging observations, much use has been made of the ���xed �� ���tting (FPF) and harmonic mean ���tting (HMF) techniques, in which the solar wind transient is considered to be a radially propagating point source (���xed ��, FP, model) and a radially expanding circle anchored at Sun centre (harmonic mean, HM, model), respectively. Initially, we compare the radial speeds and propagation directions derived from application of the FPF and HMF techniques to a large set of STEREO/Heliospheric Imager (HI) observations. As the geometries on which these two techniques are founded constitute extreme descriptions of solar wind transients in terms of their extent along the line of sight, we describe a single-spacecraft ���tting technique based on a more generalized model for which the FP and HM geometries form the limiting cases. In addition to providing estimates of a transient���s speed and propagation direction, the self-similar expansion ���tting (SSEF) technique provides, in theory, the capability to estimate the transient���s angular extent in the plane orthogonal to the ���eld of view. Using the HI observations, and also by performing a Monte Carlo simulation, we assess the potential of the SSEF technique.