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.     

11 million years of Oligocene geomagnetic field behaviour

An 11 million year long record of the Oligocene geomagnetic field has been obtained from pelagic sediments of DSDP Hole 522 An average sample spacing of 4 cm yielded approximately one specimen per 4 to 8 kyr. The rock magnetics are remarkabh consistent across the entire interval. Previous work demonstrated a magnetic mineralogy dominated by magnetically stable magnetite. The natural remanent magnetism (NRM) carries an Oligocene polarity timescale that is in excellent agreement with the Oligocene reversal record as determined from marine magnetic anomalies (MMAs), including many of the so-called 'crypto-chrons'. Normalized NRM intensities from the undisturbed portions of the record yield a time series of variations with features consistent with a number of other palaeointensity time series derived from both sedimentary and lava sequences. These features include consistent, major decreases in palaeointensity (DIPs) at reversal boundaries, and occasional DIPs between reversal boundaries that could correspond to lineated 'tiny wiggles' in the MMA records. The data set suggests that the overall field strength was 40 per cent higher in the first half of the Oligocene when the average reversal frequency was 1.6 Myr^-1 than in the second half when the reversal frequency was 4 Myr^-1. There is also a weak dependence of average field strength on length of polarity interval. Finally, in the three cores suited to spectral analysis (of coherent polarity and relative intensity independent of lithological contamination), there is a persistent ca. 30–50ka periodicity in the variations of the relative intensity, suggesting that the geomagnetic field 'pulses' at about this frequency, not only during the Brunhes (as demonstrated by Tauxe & Shackleton 1994), but in the Oligocene as well.     

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.