We report on the concept of an innovative source to produce polarized proton/deuteron beams of a kinetic energy up to several GeV from a laser-driven plasma accelerator. Spin effects have been implemented into the particle-in-cell (PIC) simulation code VLPL (Virtual Laser Plasma Lab) to make theoretical predictions about the behavior of proton spins in laser-induced plasmas. Simulations of spin-polarized targets show that the polarization is conserved during the acceleration process. For the experimental realization, a polarized HCl gas-jet target is under construction using the fundamental wavelength of a Nd:YAG laser system to align the HCl bonds and simultaneously circularly polarized light of the fifth harmonic to photo-dissociate, yielding nuclear polarized H atoms. Subsequently, their degree of polarization is measured with a Lamb-shift polarimeter. The final experiments, aiming at the first observation of a polarized particle beam from laser-generated plasmas, will be carried out at the 10 PW laser system SULF at SIOM, Shanghai.

Key issues of glacier monitoring are changes in glacier geometry and glacier mass. As accurate direct measurements are costly and time-consuming, the use of various remote-sensing data for glacier monitoring is explored. One technology used and described here is airborne laser scanning. The method enables the derivation of high-quality digital elevation models (DEMs) with a vertical and horizontal accuracy in the sub-metre range. Between September 2001 and August 2002, three laser scanner data acquisition flights were carried out, covering the whole area of Engabreen, Norway, and corresponding well to the measurement dates for the mass-balance year 2001/02. The data quality of the DEMs is assessed (e.g. by comparing the values with a control area which has been surveyed independently or GPS ground profiles measured during the flights). For the whole glacier, surface elevation change and consequently volume change is calculated, quantified and compared with traditional mass-balance data for the same time interval. For the winter term, emergence/submergence velocity is determined from laser scanner data and snow-depth data and is compared with velocity measurements at stakes. The investigations reveal the high potential of airborne laser scanning for measuring the extent and the topography of glaciers as well as changes in geometry (Δarea, Δvolume).

A new continuum approach to snowdrift modelling is introduced. In addition, numerical studies are carried out to identify the influence of time-varying wind conditions on snowdrift simulations. We compare the snowdrift patterns at Grimming mountain, Austria, derived using a time-averaged wind field and a time-varying wind field obtained from the numerical weather prediction model INCA. The results show significant differences in the deposition patterns and snow depth even after a 6 hour drift period. Using time-averaged boundary conditions leads to an underprediction of the resulting snow depth caused by averaging the wind speed, which lets gusts of wind disappear while snow transport is a non-linear function of the wind speed. Using numerical weather prediction models for snowdrift simulation therefore provides enhanced knowledge of the snow depth for local avalanche warning services.

A new radioscopic imaging technique has been developed to measure firn density in unprecedented resolution and accuracy even when the porosity is low or the geometry of a core or piece of core is not perfect. The technique is based on an X-ray microfocus computer tomograph (ICE-CT) designed especially for ice-core applications. Applied on an archive piece of the Antarctic firn core B32 drilled in Dronning Maud Land in 1998, the obtained density profile shows a strong correlation with the calcium ion concentration as found previously in Greenland. Given the impurity–density relationship found previously in Greenland, our result suggests both improved accuracy of the new density measurements and an impurity–density relationship with a similar magnitude in Greenland to that on the Antarctic plateau. Our measurements provide first evidence that the impurity–density relationship is a universal feature of polar firn and that the calcium ion concentration can serve as a proxy to describe quantitatively the effect of the impurities on densification.

A new densification model, which simulates the effect of impurities on the densification of polar firn, is presented. The classical densification models of Herron and Langway (1980) and Pimienta and Barnola (Barnola and others, 1991) are modified by assuming that the activation energy for deformation is reduced by the impurities. Motivated by recent observations, the impurity effect is formulated on an empirical basis using the seasonally varying Ca2+ ion concentration. Excellent agreement between simulated and measured high-resolution density profiles confirms the new approach. The same parameterization applies for Greenland and Antarctica despite the one order of magnitude difference in impurity concentration. The new models allow us, for the first time, to simulate the density layering in firn down to the firn–ice transition. Our results emphasize the importance of impurities and density layering for the air entrapment and for dating gas records of deep ice cores, in particular for glacial climate conditions where the impurity concentrations are 10–100-fold higher than in modern firn.

Predictions of marine ice-sheet behaviour require models able to simulate grounding-line migration. We present results of an intercomparison experiment for plan-view marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no buttressing effects from lateral drag). Perturbation experiments specifying spatial variation in basal sliding parameters permitted the evolution of curved grounding lines, generating buttressing effects. The experiments showed regions of compression and extensional flow across the grounding line, thereby invalidating the boundary layer theory. Steady-state grounding-line positions were found to be dependent on the level of physical model approximation. Resolving grounding lines requires inclusion of membrane stresses, a sufficiently small grid size (<500 m), or subgrid interpolation of the grounding line. The latter still requires nominal grid sizes of <5 km. For larger grid spacings, appropriate parameterizations for ice flux may be imposed at the grounding line, but the short-time transient behaviour is then incorrect and different from models that do not incorporate grounding-line parameterizations. The numerical error associated with predicting grounding-line motion can be reduced significantly below the errors associated with parameter ignorance and uncertainties in future scenarios.

Rayleigh–Bénard convection in rotating spherical shells can be considered as a simplified analogue of many astrophysical and geophysical fluid flows. Here, we use three-dimensional direct numerical simulations to study this physical process. We construct a dataset of more than 200 numerical models that cover a broad parameter range with Ekman numbers spanning $3\times 10^{-7}\leqslant E\leqslant 10^{-1}$ , Rayleigh numbers within the range $10^{3}<Ra<2\times 10^{10}$ and a Prandtl number of unity. The radius ratio $r_{i}/r_{o}$ is 0.6 in all cases and gravity is assumed to be proportional to $1/r^{2}$ . We investigate the scaling behaviours of both local (length scales, boundary layers) and global (Nusselt and Reynolds numbers) properties across various physical regimes from onset of rotating convection to weakly rotating convection. Close to critical, the convective flow is dominated by a triple force balance between viscosity, Coriolis force and buoyancy. For larger supercriticalities, a small subset of our numerical data approach the asymptotic diffusivity-free scaling of rotating convection $Nu\sim Ra^{3/2}E^{2}$ in a narrow fraction of the parameter space delimited by $6\,Ra_{c}\leqslant Ra\leqslant 0.4\,E^{-8/5}$ . Using a decomposition of the viscous dissipation rate into bulk and boundary layer contributions, we establish a theoretical scaling of the flow velocity that accurately describes the numerical data. In rapidly rotating turbulent convection, the fluid bulk is controlled by a triple force balance between Coriolis, inertia and buoyancy, while the remaining fraction of the dissipation can be attributed to the viscous friction in the Ekman layers. Beyond $Ra\simeq E^{-8/5}$ , the rotational constraint on the convective flow is gradually lost and the flow properties continuously vary to match the regime changes between rotation-dominated and non-rotating convection. We show that the quantity $RaE^{12/7}$ provides an accurate transition parameter to separate rotating and non-rotating convection.

The 14C Quality Assurance Programme coordinated by the IAEA (Rozanski et al. 1992) prepared a set of five new intercomparison materials, including 40–50 ka old subfossil wood excavated from New Zealand peat bogs (IAEA C-4 standard). Statistical analysis of 79 14C measurements made on the wood indicated considerable variation in the results, with a marked skewness toward more modern values. The wide range of results and the possibility of inhomogeneity within the standard prompted the recovery and analysis of replacement material. The new subfossil wood sample is kauri (Agathis australis), at least 50 ka old, excavated from a swamp in Northland. It is in the form of a single plank, 6 m long, weighing 80 kg. It will be forwarded to the IAEA in Vienna for milling and distribution. Subsamples were obtained from both ends of the plank and analyzed by six laboratories. We present here the results of these analyses and compare them with the previous IAEA intercalibration results for the C-4 standard.

New radiocarbon calibration curves, IntCal04 and Marine04, have been constructed and internationally ratified to replace the terrestrial and marine components of IntCal98. The new calibration data sets extend an additional 2000 yr, from 0–26 cal kyr BP (Before Present, 0 cal BP = AD 1950), and provide much higher resolution, greater precision, and more detailed structure than IntCal98. For the Marine04 curve, dendrochronologically-dated tree-ring samples, converted with a box diffusion model to marine mixed-layer ages, cover the period from 0–10.5 cal kyr BP. Beyond 10.5 cal kyr BP, high-resolution marine data become available from foraminifera in varved sediments and U/Th-dated corals. The marine records are corrected with site-specific 14C reservoir age information to provide a single global marine mixed-layer calibration from 10.5–26.0 cal kyr BP. A substantial enhancement relative to IntCal98 is the introduction of a random walk model, which takes into account the uncertainty in both the calendar age and the 14C age to calculate the underlying calibration curve (Buck and Blackwell, this issue). The marine data sets and calibration curve for marine samples from the surface mixed layer (Marine04) are discussed here. The tree-ring data sets, sources of uncertainty, and regional offsets are presented in detail in a companion paper by Reimer et al. (this issue).

The first meeting of the IntCal04 working group took place at Queen's University Belfast from April 15 to 17, 2002. The participants are listed as co-authors of this report. The meeting considered criteria for the acceptance of data into the next official calibration dataset, the importance of including reliable estimates of uncertainty in both the radiocarbon ages and the cal ages, and potential methods for combining datasets. This preliminary report summarizes the criteria that were discussed, but does not yet give specific recommendations for inclusion or exclusion of individual datasets.