A record of accelerated erosion in the recent sediments of Blelham Tarn in the English Lake district

Two frozen cores from Blelham Tarn were subsampled and measured using mineral magnetic, loss-on-ignition (LOI), radiometric, granulometric and diatom analyses. A detailed chronology was established using varves, radioisotopes and diatoms. This has enabled an accurately dated reconstruction of sedimentation over the past forty years. Despite a large increase in lake productivity, evidence suggests that the observed exponential increase in sedimentation rates can be attributed to erosion within the catchment. The predominant sediment source has been identified as surface soil. A comparison between the trend of accelerated sedimentation and the record of increased sheep stocking density for the area within which the most of the catchment lies, as well as observations of contemporary surface processes within the catchment, both suggest that much of the recent erosion is a direct response to increased pressure from sheep grazing.     

Flexible spectral methods for the generation of random fields with power-law semivariograms

Random field generators serve as a tool to model heterogeneous media for applications in hydrocarbon recovery and groundwater flow. Random fields with a power-law variogram structure, also termed fractional Brownian motion (fBm) fields, are of interest to study scale dependent heterogeneity effects on one-phase and two-phase flow. We show that such fields generated by the spectral method and the Inverse Fast Fourier Transform (IFFT) have an incorrect variogram structure and variance. To illustrate this we derive the prefactor of the fBm spectral density function, which is required to generate the fBm fields. We propose a new method to generate fBm fields that introduces weighting functions into the spectral method. It leads to a flexible and efficient algorithm. The flexibility permits an optimal choice of summation points (that is points in frequency space at which the weighting function is calculated) specific for the autocovariance structure of the field. As an illustration of the method, comparisons between estimated and expected statistics of fields with an exponential variogram and of fBm fields are presented. For power-law semivariograms, the proposed spectral method with a cylindrical distribution of the summation points gives optimal results.     

The identifiability of parameters in a water quality model of the Biebrza River, Poland

The identifiability of model parameters of a steady state water quality model of the Biebrza River and the resulting variation in model results was examined by applying the Monte Carlo method which combines calibration, identifiability analysis, uncertainty analysis, and sensitivity analysis. The water quality model simulates the steady state concentration profiles of chloride, phosphate, ammonium, and nitrate as a function of distance along a river. The water quality model with the best combination of parameter values simulates the observed concentrations very well. However, the range of possible modelled concentrations obtained for other more or less equally eligible combinations of parameter values is rather wide. This range in model outcomes reflects possible errors in the model parameters. Discrepancies between the range in model outcomes and the validation data set are only caused by errors in model structure, or (measurement) errors in boundary conditions or input variables. In this sense the validation procedure is a test of model capability, where the effects of calibration errors are filtered out. It is concluded that, despite some slight deviations between model outcome and observations, the model is successful in simulating the spatial pattern of nutrient concentrations in the Biebrza River.     

Influence of the Load Wavelength on the Permeability of a Viscosity Interface in the Mantle

Assuming a radially stratified Newtonian mantle in a steady-state approximation, we demonstrate that the permeability of a viscosity interface at 660-km depth strongly depends on the wavelength of buoyancy forces driving the flow. The flow induced by long-wavelength loads penetrates through the boundary freely even if the viscosity increases by two orders. In contrast, the boundary is practically impermeable for short-wavelength loads located in the upper mantle. Thus, a stepwise increase of viscosity is a significant obstacle for small descending features in the upper mantle, but huge upper mantle downwellings, or upwellings formed in the-lower mantle can overcome it easily. This indicates that certain care is necessary in interpreting the seismic structure of the mantle by means of flow models. The global tomographic image includes only the first few degrees of the harmonic series and, consequently, its interpretation in terms of a present-day flow field results in a predominantly whole-mantle circulation even for extreme viscosity contrasts.     

Detecting lateral resistivity inhomogeneities with the Schlumberger array

A method of detecting lateral resistivity inhomogeneities with a multi-electrode cable for the Schlumberger array is presented. Using such a cable system, two dipole-dipole soundings are gathered in addition to the Schlumberger sounding at each point. Offset differences calculated from the dipole-dipole data give qualitative information about lateral resistivity inhomogeneities. Results from 1D modelling can give quantitative information about the 2D resistivity distribution since the dipole-dipole soundings are made to either side of the Schlumberger sounding point. Examples from two different locations in Norway are shown. At Reinøya, northern Norway, a dipping layer, confirmed by refraction seismic data, was identified. In a sedimentary basin with a contamination plume at Haslemoen, southern Norway, the method has revealed lateral variations in resistivity. In both cases, the Schlumberger soundings could be fitted with a 1D resistivity model.     

Space-time mapping of soil salinity using probabilistic bayesian maximum entropy

The mapping of saline soils is the first task before any reclamation effort. Reclamation is based on the knowledge of soil salinity in space and how it evolves with time. Soil salinity is traditionally determined by soil sampling and laboratory analysis. Recently, it became possible to complement these hard data with soft secondary data made available using field sensors like electrode probes. In this study, we had two data sets. The first includes measurements of field salinity (EC_a) at 413 locations and 19 time instants. The second, which is a subset of the first (13 to 20 locations), contains, in addition to EC_a, salinity determined in the laboratory (EC_2.5). Based on a procedure of cross-validation, we compared the prediction performance in the space-time domain of 3 methods: kriging using either only hard data (HK) or hard and mid interval soft data (HMIK), and Bayesian maximum entropy (BME) using probabilistic soft data. We found that BME was less biased, more accurate and giving estimates, which were better correlated with the observed values than the two kriging techniques. In addition, BME allowed one to delineate with better detail saline from non-saline areas.     

Gravity reduction with three-dimensional atmospheric pressure data for precise ground gravity measurements

The redistribution of air masses induces gravity variations (atmospheric pressure effect) up to about 20 μgal. These variations are disturbing signals in gravity records and they must be removed very carefully for detecting weak gravity signals. In the past, different methods have been developed for modelling of the atmospheric pressure effect. These methods use local or two-dimensional (2D) surface atmospheric pressure data and a standard height-dependent air density distribution. The atmospheric pressure effect is consisting of the elastic deformation and attraction term. The deformation term can be well modelled with 2D surface atmospheric pressure data, for instance with the Green's function method. For modelling of the attraction term, three-dimensional (3D) data are required. Results with 2D data are insufficient.From European Centre for Medium-Range Weather Forecasts (ECMWF) 3D atmospheric pressure data are now available. The ECMWF data used here are characterised by a spacing of Δ and Δλ = 0.5°, 60 pressure levels up to a height of 60 km and an interval of 6 h. These data are used for modelling of the atmospheric attraction term. Two attraction models have been developed based on the point mass attraction of air segments and the gravity potential of the air masses. The modelling shows a surface pressure-independent part of gravity variations induced by mass redistribution of the atmosphere in the order of some μgal. This part can only be determined by using 3D atmospheric pressure data. It has been calculated for the Vienna Superconducting Gravimeter site.From this follows that the gravity reduction can be improved by applying the 3D atmospheric attraction model for analysing long-periodic tidal waves including the polar tide. The same improvement is expected for reduction of long-term absolute gravity measurements or comparison of gravity measurements at different seasonal times. By using 3D atmospheric pressure data, the gravity correction can be improved up to some μgal.     

Managed water levels and the expansion of emergent vegetation along a lakeshore

Water-level fluctuations may be used to promote the expansion of emergent vegetation along lakeshores. We present the case of the lake Volkerak-Zoommeer in the Netherlands, a fresh-water lake created in 1987 after the enclosure of an estuary. Using an experimental area in which the water level could be manipulated, it was shown that partial summer drawdown of the shoreline created suitable conditions for germination and growth of tall emergent species (in particular Phragmites australis). Plant survival and growth depended on subsequent water-level fluctuations and grazing by waterbirds. Based on the experiment and empirical data, a model was developed to predict the effects of the water-level regime on potential reed bed development. The model was applied for four hydrological scenarios that have been considered for the water-level management of the lake.