A GIS based approach to back trajectory analysis for the source apportionment of aerosol constituents and its first application

A new method has been developed to combine back trajectory statistics with a detailed land cover analysis. It provides numeric proxies for the residence times of sampled air masses above certain land cover classes (marine, natural vegetation, agricultural lands, urban areas, and bare areas), as well as further meteorological parameters (mean trajectory length, solar radiation along trajectory, and local height of the boundary mixing layer). The method has been implemented into a GIS-enabled database system to allow for an efficient processing of large datasets with low computational demands. A principal component analysis was performed on a dataset including the modelled residence times, the modelled meteorological parameters, some measured meteorological parameters (wind speed and temperature), and the concentrations of 10 particle constituents (inorganic ions and organic and elemental carbon) in 5 particle size ranges for 29 winter- and summertime samples at an urban background site in Leipzig, Germany. Six principal components could be extracted which together explained about 80% of the total variance in the dataset. The factors could be attributed to the influence of meteorology to continental background pollution, secondary formation processes in polluted air masses, wood burning, aged sea-salt, local traffic, and long-range transported crustal material. The modelled residence times and the meteorological parameters were generally consistent with the existing knowledge of specific particle sources and thereby facilitated and strengthened the interpretation of the factors. Moreover, they allowed for a clear distinction between continental background pollution and secondary formation processes, which has not been possible in previous source apportionment studies. The results demonstrate that the combined usage of back trajectory, land cover, and meteorological data by the presented method yields valuable additional information on the history of sampled air masses, which can improve the quality of source apportionment of atmospheric aerosol constituents.     

Particle tracking in the vicinity of Helgoland,North Sea: a model comparison

Station Helgoland Roads in the south-eastern North Sea (German Bight) hosts one of the richest long-term time series of marine observations. Hydrodynamic transport simulations can help understand variability in the local data brought about by intermittent changes of water masses. The objective of our study is to estimate to which extent the outcome of such transport simulations depends on the choice of a specific hydrodynamic model. Our basic experiment consists of 3,377 Lagrangian simulations in time-reversed mode initialized every 7 h within the period Feb 2002–Oct 2004. Fifty-day backward simulations were performed based on hourly current fields from four different hydrodynamic models that are all well established but differ with regard to spatial resolution, dimensionality (2D or 3D), the origin of atmospheric forcing data, treatment of boundary conditions, presence or absence of baroclinic terms, and the numerical scheme. The particle-tracking algorithm is 2D; fields from 3D models were averaged vertically. Drift simulations were evaluated quantitatively in terms of the fraction of released particles that crossed each cell of a network of receptor regions centred at the island of Helgoland. We found substantial systematic differences between drift simulations based on each of the four hydrodynamic models. Sensitivity studies with regard to spatial resolution and the effects of baroclinic processes suggest that differences in model output cannot unambiguously be assigned to certain model properties or restrictions. Therefore, multi-model simulations are needed for a proper identification of uncertainties in long-term Lagrangian drift simulations.     

Seasonal to decadal forcing of high water level percentiles in the German Bight throughout the last century

For the purpose of coastal planning and management, especially under changing climatic conditions, enhanced knowledge about the evolution of extreme sea levels in the past, present, and future is required. This paper presents statistical analyses of high seasonal water level percentiles of 13 tide gauges in the German Bight, spanning over a period of up to 109 years throughout the twentieth and twenty-first centuries. Seasonal and annual high percentile time series of water levels were investigated in comparison to the mean sea level (MSL) for changes on seasonal, inter-annual, and decadal timescales. While throughout the first half of the twentieth century extreme water levels generally followed changes in MSL, during the second half of the century, linear extreme sea level trends exceeded those in MSL in the order of 9–64 cm per century. The largest, although insignificant, contribution to the magnitude of these trends occurs in the winter season (January to March), while smaller but, due to the generally lower atmospheric variability, significant changes are observed during spring (April to June). The observed multi-decadal trends are generally in good agreement with multi-decadal trends in the corresponding percentiles of local zonal surface winds. Only small parts of the trends remain unexplained. It is suggested that these remaining trends result from modifications in the local tidal regime. For the aspects of coastal planning, the findings clarify that in the German Bight, in addition to changes in MSL, potential changes in storminess and in the tidal regime significantly contribute to the development of extreme water levels. Since these factors have influenced the characteristic of extremes throughout the recent past, they also have to be taken into account when estimating design water levels for, e.g., dikes (in a warming climate) under changing greenhouse gas emissions.     

Interaction of waves,currents and tides,and wave-energy impact on the beach area of Sylt Island

Erosion due to waves is an important and actual problem for most coastal areas of the North Sea. The objective of this study was to estimate the impact of wave action on the coastline of Sylt Island. From a 2-year time series (November 1999 to October 2001) of hydrological and wave parameters generated with a coupled wave–current modelling system, a period comprising storm ‘Anatol’ (3–4 December 1999) is used to investigate the effects of waves on currents and water levels and the input of wave energy into the coastline. The wave-induced stress causes an increase of the current velocity of 1 m/s over sand and an additional drift along the coast of about 20 cm/s. This produces a water level increase of more than 20 cm in parts of the tidal basin. The model system also calculates the wave energy input into the coastline. Scenario runs for December 1999 with a water level increase of 50 cm and wind velocity increased by 10% show that the input of the wave energy into the west coast of Sylt Island increases by 30% compared to present conditions. With regard to the forecasted near-future (Woth et al., Ocean Dyn 56:3–15, 2006) increase of strong storm surges, the scenario results indicate an increased risk of coastal erosion in the surf zone of Sylt Island.     

Ten Years of Development of Equipment for Measurement of Atmospheric Radioactive Xenon for the Verification of the CTBT

Atmospheric measurement of radioactive xenon isotopes (radioxenon) plays a key role in remote monitoring of nuclear explosions, since it has a high capability to capture radioactive debris for a wide range of explosion scenarios. It is therefore a powerful tool in providing evidence for nuclear testing, and is one of the key components of the verification regime of the Comprehensive Nuclear-Test-Ban Treaty (CTBT). The reliability of this method is largely based on a well-developed measurement technology. In the 1990s, with the prospect of the build-up of a monitoring network for the CTBT, new development of radioxenon equipment started. This article summarizes the physical and technical principles upon which the radioxenon technology is based and the advances the technology has undergone during the last 10 years. In contrast to previously used equipment, which was manually operated, the new generation of radioxenon monitoring equipment is designed for automated and continuous operation in remote field locations. Also the analytical capabilities of the equipment were strongly enhanced. Minimum detectable concentrations of the recently developed systems are well below 1 mBq/m³ for the key nuclide ¹³³Xe for sampling periods between 8 and 24 h. All the systems described here are also able to separately measure with low detection limits the radioxenon isotopes ^131mXe, ^133mXe and ¹³⁵Xe, which are also relevant for the detection of nuclear tests. The equipment has been extensively tested during recent years by operation in a laboratory environment and in field locations, by performing comparison measurements with laboratory type equipment and by parallel operation. These tests demonstrate that the equipment has reached a sufficiently high technical standard for deployment in the global CTBT verification regime.     

Laboratory seismic measurements: an aid in the interpretation of seismic field data

This contribution reports data and evidence from laboratory measurements of compressional and shear-wave velocities and mineralogical, microstructural and chemical analyses in crustal rocks that may provide an aid to understanding and interpreting geophysical data. It is shown that the correlation of laboratory seismic data with stratigraphic field data can provide the necessary data base for constructing a seismic model of the Earth's crust in a particular geological area.     

Syngenetic, thin graphite-rich horizons in lower crustal rocks from the Serre San Bruno, Calabria (Italy), and implications for the nature of high-conducting deep crustal layers

In order to better understand the nature of deep crustal high electrical conductivity, we studied the electrical properties of a tilted section of a former lower continental crust exposed in the Calabrian arc of the Alpine-Apennine mountain system. Geoelectric field measurements and impedance measurements on rock samples showed that these high-grade metamorphic rocks are generally highly resistive as expected for crystalline, electrolytically conducting rocks of low porosity. This holds for graphite-free metabasites as well as for metapelites which generally contain accessory, up to 3% biogenic graphite in the form of isolated grains. Clearly as an exception, a group of thin stratiform black horizons with thicknesses of 1-15 cm within the metapelitic series was detected by means of self-potential (SP) measurements. Rock samples from these horizons exhibit high, quasi-metallic bulk conductivities of up to 50 S/m (0.02 Ωm) in agreeement with up to 20% syngenetic graphite, forming a network of interconnected streaks or crack fillings. The high amount of carbon most probably originates from organic matter of C_org-rich black shales. Relative enrichment of the low mobility graphitic matter compared to the carbon content of the assumed protoliths may have been due to pressure solution and partial melting during prograde metamorphism, without major contribution of a fluid phase, resulting in isolated graphite flakes. Although enriched, graphite in this form has little effect on electrical conductivity. For the Calabrian black horizons, microscopic analyses make conceivable that, in a further decisive step, isolated graphite grains were mechanically smeared to continuous pathways during uplift by shearing, producing hereby the observed graphitic network which is needed to generate high conductivity. As C_org-rich black shales are common members of sedimentary sequences throughout the earth’s history, good conductors of this type may be expected in the continental crust at any depth depending on tectonic and metamorphic history, with the exception of magmatic protoliths. Regarding the extremely high conductivity of the meta-black shale samples containing syngenetic sheared graphite, a total thickness of a few meters of such rocks is sufficient to explain magnetotelluric high conductivity anomalies in the deep crust.     

Validating a turbulence closure against estuarine microstructure measurements

A simple k–ε turbulence closure is introduced which has no stability functions but instead a Richardson number-dependent turbulent Prandtl number. Its free parameters are determined in a comparison with microstructure observations from a stratified and sheared tidal estuary and laboratory measurements. The closure is able to simulate observed turbulent dissipation rates (ε) and turbulent length scales (l_th) in regions of strong mean shear and small gradient Richardson number (R_g) to within factors of 2–3. It fails in regions of small shear and large R_g, presumably because of the dominance of internal wave-driven mixing. Additional simulations with a k–ε closure with stability functions taken from Canuto et al. [Canuto, V.M., Howard, A., Cheng, Y., Dubovikov, M.S., 2001. Ocean turbulence I: one-point closure model. Momentum and heat vertical diffusivities. J. Phys. Oceanogr. 31, 1413–1426] and with the closure of Baumert and Peters [Baumert, H., Peters, H., 2004. Turbulence closure, steady state, and collapse into waves. J. Phys. Oceanogr. 34, 505–512] show poor performance. Establishing a valid 1:1 comparison of simulated and observed ε and l_th requires nudging the model velocity and density toward observed values because free model integrations quickly diverge from the observations. Steady state gradient Richardson numbers are constrained to a range of 0.18–0.25, while flux Richardson numbers are constrained to the range of 0.1–0.22. The closure output is rather insensitive to such parameter variations. The simulations are sensitive, however, to the treatment of the observed velocity and density used to nudge the model. Good closure performance requires averaging the measured tidal flow over about an hour, a time scale for which conventional numerical models of estuarine circulations should be able to match observed shears. In the closure simulations the TKE balance stays close to a production–dissipation balance. The time rate of change and vertical diffusion of TKE are small, of the same order of magnitude, and vary in magnitude relative to each other systematically across the water column.