Atmospheric deposition of heavy metals to North European marginal seas: Scenarios and trend for lead

The contribution of the atmosphere to the input of heavy metals to marine environments of Northern Europe is determined by applying an EMEP-type Lagrangian trajectory model. The results show that the model is capable of simulating long-range transport of heavy metals from European anthropogenic sources to the North Sea and the Baltic Sea including emission reduction scenarios for lead. Model calculations for lead indicate maximum annual inputs from the atmosphere of more than 2,400 tonnes for the North Sea and 1,300 tonnes for the baltic sea. It was calculated that in the 1985–1990 time period, mainly as a consequence of use of unleaded gasoline, the trend show a pronounced decrease of the deposition fluxes of lead from the atmosphere.     

Stratigraphy and facies distribution of the Utsira formation and the Pliocene sequences in the northern North Sea

The Upper Cenozoic deposits of the northern North Sea have been analysed in order to establish a regional and detailed stratigraphy. The Utsira Formation is subdivided into four log-units and mapped, and two main depocentres are outlined. The lower part of the Utsira Formation consists of thick marine, mounded sand bodies, interpreted as overall stacked lowstand fan deposits, while the upper part of the formation consists of more clayey-silty intervals, indicating increased relative sea level. The succeeding progradational Pliocene deposits are subdivided into 13 high-frequency depositional sequences and are mapped. The sequences are grouped into four composite sequences. Each of the four Pliocene composite sequences is composed of one or two rather locally distributed, prograding sequences (lowstand sequence set), succeeded by one or two more widely distributed aggrading-prograding sequences (transgressive-highstand sequence set). Boundaries between the composite sequences are recorded as marked changes in distribution of depocentres and sequence architecture. The regional uplift of Scandinavia is believed to be the main control on sediment input, feeding the succeeding general prograding Pliocene sequences. Oscillations of the eustatic sea level punctuated the tectonically controlled progradation and affected variations in the accommodation space, and thus created the high-frequency sequences.     

Correction of aerosol influence in Landsat 5 Thematic Mapper data

The influence of aerosols in the visible and near infrared part of the electromagnetic spectrum was studied by simulations of Landsat 5 Thematic Mapper measurements. The radiative transfer model used is based on the matrix-operator-method and was applied to different surface types, represented by specific spectral albedo values. On the basis of a single scattering approach for atmospheric correction, an algorithm was developed to correct for the influences of aerosols, air molecules and athmospheric trace gases on Thematic Mapper measurements above land surfaces using additional measurements above nearby located ocean surfaces to estimate the optical properties. The optical thickness of a cloud-free atmosphere has therefore been varied in the model for different aerosol types and surface reflectances.     

Applying SAR-IRSL methodology for dating fine-grained sediments from Lake El’gygytgyn, north-eastern Siberia

Lake El’gygytgyn is situated in a 3.6 Ma old impact crater in north-eastern Siberia and probably represents one of the most complete archives of Arctic climate change. Investigated here is the potential of infra-red stimulated luminescence (IRSL) using the single-aliquot regenerative-dose (SAR) approach for dating sediments from this lake. Independent age control is available from a published age model of a parallel core that is based on tuning sediment proxies with regional insolation and the results of previous multiple aliquot IRSL dating. Although the site is located within volcanic bedrock, anomalous fading seems to have little effect on the calculated ages. The modelled water content for the entire time of burial is seen as the most prominent uncertainty at this particular site. Despite these potential error sources, SAR-IRSL ages are in acceptable agreement with the given timeframe and clearly point to the possibility to establish independent chronologies at this site up to at least 400,000 years.     

Resetting of sediments mobilised by the LGM ice-sheet in southern Norway

Former geological field investigations in the Rondane area, east-central southern Norway, have proposed that the maximum Fennoscandian ice-sheet coverage occurred during the Late Weichselian Glacial Maximum (LGM, ca. 20 ka) and that subsequent glaciofluvial sediments were first deposited in the early Holocene (after 10 ka). However, recent field investigations with ages from three internally consistent quartz optically stimulated luminescence (OSL) age series show an apparent deglaciation of northern Rondane in the period 20.0–13.8 ka. We examine here the possibility that these ages are too old because the sediment was not completely zeroed prior to deposition. Our investigations of incomplete bleaching use modern analogues, small aliquots, and single grains of quartz. First, the symmetric shape of small aliquot equivalent dose distributions suggests that the sediment was probably well bleached at deposition. This is supported by 5 modern analogue equivalent doses (D_e) of 0.6 Gy, 1.5% of the typical D_e from the deglaciation sediments. Finally, from single grain studies on three samples, we conclude that there is no evidence for poor bleaching in these samples; thus the weighted mean gives the best estimates of D_e, and these are completely consistent with both large and small aliquot estimates for these samples. These comparisons between large aliquots, modern analogues, small aliquots and single-grain analyses help to validate the OSL ages and confirm the complete resetting of these sediments prior to deposition.     

Intercomparison of Two Hydroxyl Radical Measurement Techniques at the Atmosphere Simulation Chamber SAPHIR

At the atmosphere simulation chamber SAPHIR in Jülich both Laser-Induced Fluorescence Spectroscopy (LIF) and Long-Path Differential Optical Laser Absorption Spectroscopy (DOAS) are operational for the detection of OH radicals at tropospheric levels. The two different spectroscopic techniques were compared within the controlled environment of SAPHIR based on all simultaneous measurements acquired in 2003 (13 days). Hydroxyl radicals were scavenged by added CO during four of these days in order to experimentally check the calculated precisions at the detection limit. LIF measurements have a higher precision (σ= 0.88×10⁶ cm^–3) and better time resolution (Δt = 60 s), but the DOAS method (σ= 1.24×10⁶ cm^–3, Δt = 135 s) is regarded as primary standard for comparisons because of its good accuracy. A high correlation coefficient of r = 0.95 was found for the whole data set highlighting the advantage of using a simulation chamber. The data set consists of two groups. The first one includes 3 days, where the LIF measurements yield (1 – 2) ×10⁶ cm^–3 higher OH concentrations than observed by the DOAS instrument. The experimental conditions during these days are characterized by increased NO_x concentration and a small dynamic range in OH. Excellent agreement is found within the other group of 6 days. The regression to the combined data of this large group yields unity slope without a significant offset.     

Temporal-moment matching for truncated breakthrough curves for step or step-pulse injection

The method of temporal moments is an efficient approach for analyzing breakthrough curves (BTCs). By matching the moments of the BTCs computed through parametric transfer-function models or one-dimensional transport models to those of the data, one can estimate the parameters characterizing the transfer function or apparent transport parameters. The classical method of moments presumes infinite duration. However, the measurement of BTCs is usually terminated prematurely, before the concentration has reached zero. Unless this truncation of the BTCs has been taken into account, the estimates of the parameters may be in error. Truncated measured BTCs are sometimes extrapolated assuming exponential decay. In this study, we use the concept of moments of the truncated impulse–response function [Jawitz JW. Moments of truncated continuous univariate distributions. Adv Water Res 2004;27:269–81] in the analysis of truncated BTCs corresponding to the commonly encountered step and step-pulse injection modes. The method is straightforward, based on the relation, which we derive, between truncated moments of the impulse–response function and the measured BTC. It is practical to apply and does not require the extrapolation of the measured BTC. The method is also accurate. In a numerical study we discuss how short a step-pulse injection may be so that we can approximate it as instantaneous. Finally, we apply the method to the analysis of a field-scale tracer test.     

Routine analysis of sediment pore water of high ionic strength

Pore water of sediments plays an important role in aquatic systems as mediator and as the reactive zone between the sediment and surface water. Sediment pore waters with high ionic strength from acidified pit lakes were investigated to obtain information about the influence on the lake water quality. The analysis of soluble reactive phosphorus, nitrate, ammonium, silica, dissolved organic carbon, ferrous iron, sulfate, chloride, sodium, potassium, calcium, magnesium, manganese, and total dissolved iron was predetermined as the dataset required for evaluation of water quality. The data collection procedure was optimised by designing a methodology for stabilisation, dilution of pore water samples and adaptation of analytical methods. The developed methodology was evaluated with respect to the effort required in the laboratory under routine conditions. In the first round of analyses, 72% of 638 individual analyses from a random selection of pore water profiles were found to be within the calibration ranges. Remedial actions to handle the remaining 28% of invalid analytical results are exemplified. Differences between comparative analyses of some ions by continuous flow analysis, ion chromatography, and atomic emission spectroscopy were evaluated. The majority of results measured by ion chromatography differed on a highly significant level from results measured by atomic emission spectroscopy. Possible reasons, originating from the extreme sample matrix, are discussed. Finally, the designed methodology and the results of the method comparison are used to recommend the selection of analytical methods under specific conditions.     

Aluminous granulites from the Pulur complex, NE Turkey: a case of partial melting, efficient melt extraction and crystallisation

In the Pulur complex, NE Turkey, a heterogeneous rock sequence ranging from quartz-rich mesocratic gneisses to silica- and alkali-deficient, Fe-, Mg- and Al-rich melanocratic rocks is characterized by granulite-facies assemblages involving garnet, cordierite, sillimanite, ilmenite, ±spinel, ±plagioclase, ±quartz, ±biotite, ±corundum, rutile and monazite. Textural evidence for partial melting in the aluminous granulites, particularly leucosomes, is largely absent or strongly obliterated by a late-stage hydrothermal overprint. However, inclusion relations, high peak P–T conditions, the refractory modes, bulk and biotite compositions of the melanocratic rocks strongly support a model of partial melting. The melt was almost completely removed from the melanocratic rocks and crystallised within the adjacent mesocratic gneisses which are silica-rich, bear evidence of former feldspar and show a large range in major element concentrations as well as a negative correlation of most elements with SiO₂. Peak conditions are estimated to be ≥800 °C and 0.7–0.8 GPa. Subsequent near-isothermal decompression to 0.4–0.5 GPa at 800–730 °C is suggested by the formation of cordierite coronas and cordierite–spinel symplectites around garnet and in the matrix. Sm–Nd, Rb–Sr and ⁴⁰Ar/³⁹Ar isotope data indicate peak conditions at 330 Ma and cooling below 300 °C at 310 Ma.     

Efficient Computation of Linearized Cross-Covariance and Auto-Covariance Matrices of Interdependent Quantities

In many geostatistical applications, spatially discretized unknowns are conditioned on observations that depend on the unknowns in a form that can be linearized. Conditioning takes several matrix–matrix multiplications to compute the cross-covariance matrix of the unknowns and the observations and the auto-covariance matrix of the observations. For large numbers n of discrete values of the unknown, the storage and computational costs for evaluating these matrices, proportional to n ², become strictly inhibiting. In this paper, we summarize and extend a collection of highly efficient spectral methods to compute these matrices, based on circulant embedding and the fast Fourier transform (FFT). These methods are applicable whenever the unknowns are a stationary random variable discretized on a regular equispaced grid, imposing an exploitable structure onto the auto-covariance matrix of the unknowns. Computational costs are reduced from O(n ²) to O(nlog₂ n) and storage requirements are reduced from O(n ²) to O(n).