Selected topics in arctic atmosphere and climate

This paper summarizes the main elements of four IPY projects that examine the Arctic Atmosphere. All four projects focus on present conditions with a view to anticipating possible climate change. All four investigate the Arctic atmosphere, ocean, ice, and land interfacial surfaces. One project uses computer models to simulate the dynamics of the Arctic atmosphere, storms, and their interactions with the ocean and ice interface. Another project uses statistical methods to infer transports of pollutants as simulated in large-scale global atmospheric and oceanic models verifying results with available observations. A third project focuses on measurements of pollutants at the ice-ocean?Catmosphere interface, with reference to model estimates. The fourth project is concerned with multiple, high accuracy measurements at Eureka in the Canadian Archipelago. While these projects are distinctly different, led by different teams and interdisciplinary collaborators, with different technical approaches and methodologies, and differing objectives, they all strive to understand the processes of the Arctic atmosphere and climate, and to lay the basis for projections of future changes. Key findings include: ? Decreased sea ice leads to more intense storms, higher winds, reduced surface albedo, increased surface air temperature, and enhanced vertical mixing in the upper ocean. ? Arctic warming may affect toxic chemicals by remobilizing persistent organic pollutants and augmenting mercury deposition/retention in the environment. ? Changes in sea ice can dramatically change processes in and at the ice surface related to ozone, mercury and bromine oxide and related chemical/physical properties. ? Structure and properties of the Arctic atmospheric??troposphere to stratosphere??and tracking of transport of pollution and smoke plumes from mid-latitudes to the poles.     

The natural radiochemical and growth records in contemporary hermatypic corals from the Atlantic and Caribbean

The skeletal banding observed in vertical cross-sections of two species of coral,Montastrea annularis from Jamaica andDiploria labrynthiformis from Bermuda, has been confirmed to be annual by ingrowth of²²⁸Th towards radioactive equilibrium with²²⁸Ra taken up from seawater. In addition the Bermuda coral is shown to have incorporated²¹⁰Pb at essentially constant specific activity over at least the last thirty years.The individual annual patterns of growth banding in a suite of specimens ofD. Labrynthiformis from Bermuda are shown to be correlative. The mass of CaCO₃ accumulated each year is positively correlated with the band width for the corresponding year.     

Oxidation and Origin of Organic Matter in Surficial Eastern Mediterranean Hemipelagic Sediments

Aerobic mineralisation of C_org in surface sedimentsof the deep (>2000 m water depth) eastern Mediterranean Sea has been quantified by analysis of detailedbox core C_org concentration versus depth profiles and the modelling environment for early diageneticproblems MEDIA. The reactive fraction comprises 60–80% of the total C_org reachingthe sediments and is largely oxidised within the surficial 10 cm. A non-reactive C _orgfraction (G_NR) dominates at depths >10 cm, and makes up20–40% of the total C _org flux to the sediments. First-order rateconstants for decomposition of the reactive fraction calculated from theC _org profiles range from 5.4 × 10^-3 to8.0 × 10^-3 y^-1 to 8.0 × 10^-3 y^-1. Total mineralization rates in thesurface sediment are between 1.7 and 2.6 mol C cm^-2 y^-1 and thus are typical for oligotrophic, deep-seaenvironments. The low fluxes and rapid remineralisation of C _org are accompanied by²¹⁰Pb_excess surface mixed layers which are only 2 cm deep, among the thinnest reported for oxygenated marine sediments.Model results indicate a mismatch between the C _org profiles and O₂ microprofileswhich were measured onboard ship. This can be attributed to a combination of decompression artefactsaffecting onboard measurement of the O₂ profiles or the leakage ofoxygen into the core during handling on deck. Furthermore, the used D_b values, based on ²¹⁰Pb, may not befully appropriate; calculations with higher D_b values improve the O₂ fits. The surficial sediment¹³C _org values of -22 become less negative with increasing depth and decreasing C _orgconcentrations. The major ¹³C change occurs in the top 3 to 4 cm and coincides with the interval weremost of the organic carbon oxidation takes place. This indicates that the reactive fractionof organic matter, commonly assumed to be marine, has a more negative ¹³C _orgthan the refractory fraction, usually held to be terrestrial. Palaeoproductivity estimates calculated from thesediment data by means of literature algorithms yield low surface productivities(12–88 gC m^-2 y^-1), which are in good agreement with field measurements of primary productivity in otherstudies. Such values are, however, significantly lower than those indicated by recent productivitymaps of the area derived from satellite imagery (>100 gC m^-2 y^-1).     

Complex Rays and Wave Packets for Decaying Signals in Inhomogeneous,Anisotropic and Anelastic Media

Diffraction and anelasticity problems involving decaying, “evanescent” or “inhomogeneous” waves can be studied and modelled using the notion of “complex rays”. The wavefront or “eikonal” equation for such waves is in general complex and leads to rays in complex position-slowness space. Initial conditions must be specified in that domain: for example, even for a wave originating in a perfectly elastic region, the ray to a real receiver in a neighbouring anelastic region generally departs from a complex point on the initial-values surface. Complex ray theory is the formal extension of the usual Hamilton equations to complex domains. Liouville's phase-space-incompressibility theorem and Fermat's stationary-time principle are formally unchanged. However, an infinity of paths exists between two fixed points in complex space all of which give the same final slowness, travel time, amplitude, etc. This does not contradict the fact that for a given receiver position there is a unique point on the initial-values surface from which this infinite complex ray family emanates.In perfectly elastic media complex rays are associated with, for example, evanescent waves in the shadow of a caustic. More generally, caustics in anelastic media may lie just outside the real coordinate subspace and one must trace complex rays around the complex caustic in order to obtain accurate waveforms nearby or the turning waves at greater distances into the lit region. The complex extension of the Maslov method for computing such waveforms is described. It uses the complex extension of the Legendre transformation and the extra freedom of complex rays makes pseudocaustics avoidable. There is no need to introduce a Maslov/KMAH index to account for caustics in the geometrical ray approximation, the complex amplitude being generally continuous. Other singular ray problems, such as the strong coupling around acoustic axes in anisotropic media, may also be addressed using complex rays.Complex rays are insightful and practical for simple models (e.g. homogeneous layers). For more complicated numerical work, though, it would be desirable to confine attention to real position coordinates. Furthermore, anelasticity implies dispersion so that complex rays are generally frequency dependent. The concept of group velocity as the velocity of a spatial or temporal maximum of a narrow-band wave packet does lead to real ray/Hamilton equations. However, envelope-maximum tracking does not itself yield enough information to compute synthetic seismogramsFor anelasticity which is weak in certain precise senses, one can set up a theory of real, dispersive wave-packet tracking suitable for synthetic seismogram calculations in linearly visco-elastic media. The seismologically-accepiable constant-Q rheology of Liu et al. (1976), for example, satisfies the requirements of this wave-packet theory, which is adapted from electromagnetics and presented as a reasonable physical and mathematical basis for ray modelling in inhomogeneous, anisotropic, anelastic media. Dispersion means that one may need to do more work than for elastic media. However, one can envisage perturbation analyses based on the ray theory presented here, as well as extensions like Maslov's which are based on the Hamiltonian properties.     

Records of radionuclide deposition in two salt marshes in the United Kingdom with contrasting redox and accumulation conditions

Single cores from two salt marshes in the United Kingdom located near different nuclear facilities were investigated to compare chronostratigraphic estimates derived from the natural radionuclide ²¹⁰Pb_excess with estimates from the known times of introduction of artificial radionuclides to the environment. Both cores selected had clear visual indications of redox zonation, and evidence for diagenetic redox perturbation of the radionuclide records was also sought. In the core from Beaulieu Marsh on the south coast of England, the redox zonation was revealed by the profiles of the redox-sensitive elements Mn + I, Fe + P + As, and S, and the records of nuclear discharges were entirely contained within oxidized conditions in the upper 40 cm. The constant flux/constant sedimentation ²¹⁰Pb_excess accumulation estimate was 76% of that derived from the 1963 fallout ¹³⁷Cs level (0.35 vs. 0.46 g cm⁻² yr⁻¹ dry mass), but the constant flux ²¹⁰Pb_excess method indicated that accumulation rates were lower at Beaulieu before ∼1950. On any timescale, ¹³⁷Cs appears earlier in the sediment record than its introduction to the environment, but although downward diffusion of ¹³⁷Cs relative to ²⁴¹Am has clearly occurred, the ¹³⁷Cs peak still appeared in place and there was negligible penetration of ¹³⁷Cs into underlying reduced conditions. This core also contained a peak of the ⁶⁰Co discharges from either or both the Winfrith and La Hague nuclear plants that peaked in 1980 and 1984, respectively. The sediments in the core from Wyre Marsh on the eastern coast of the Irish Sea had accumulated more rapidly than those at Beaulieu, and in this case the redox zonation could be established only from Mn and S profiles. Here, the constant initial activity ²¹⁰Pb_excess accumulation rate estimate was 125% of that derived from the ¹³⁷Cs peak correlated with the 1975 Sellafield discharge maximum (0.79 vs. 0.64 g cm⁻² yr⁻¹). Sellafield discharge ¹³⁷Cs swamps fallout or Chernobyl ¹³⁷Cs signals in this core, but the ¹³⁷Cs and ²⁴¹Am sediment records are well separated and remain consistent with the slightly different discharge patterns over time. This is so despite the fact that the maximum activity levels of both isotopes are now located well into reducing conditions out of which Mn must have migrated. The ²¹⁰Pb profile appeared similarly unaffected by the oxidized/reduced boundary in this case. This core was too short to define the limits of any downward ¹³⁷Cs migration. ²¹⁰Pb_excess accumulation rate estimates for salt marshes should be viewed with some caution because of the steady-state assumptions inherent in all ²¹⁰Pb_excess methods and the potential for fluctuating accumulation conditions and open system behavior in salt marshes.     

Thermobarometric data from a fossil zircon partial annealing zone in high pressure–low temperature rocks of eastern and central Crete, Greece

A fossil partial annealing zone of fission tracks in zircon is described from high pressure–low temperature (HP–LT) rocks of the Phyllite–Quartzite Unit (PQ) on the island of Crete, Greece. Correlation of regional trends in fission track age populations with independent thermobarometric and microstructural data, and with new experimental annealing results, allows a calibration of this low temperature thermochronological method to a degree hitherto not available from other field examples.The zircon fission track (FT) ages of samples from the PQ across Crete range from original detrital signature through reduced to completely reset. The annealing is the result of a single heating period related to the HP–LT metamorphism with near-peak temperatures lasting for only a few million years some time between 24±1 and 20±1 Ma. In eastern Crete, where rocks have experienced temperatures of 300±50 °C and pressures of 0.8±0.3 GPa, zircon FT ages range from 414±24 to 145±10 Ma. Ages above 300 Ma occur mostly near the east coast of the island in rocks which have not been heated to above ca. 280 °C and probably represent a pre-Variscan source. Track lengths are already indicative of a substantial annealing at this temperature. Most of the zircon FT ages from eastern Crete scatter within error around the stratigraphic age. Samples with apparent zircon FT ages significantly younger than the depositional age are only observed in areas where temperatures exceeded ca. 320 °C. Towards the west, a sudden decrease to very young ages ranging from 17±2 to 18±1 Ma reflects a complete resetting at ca. 350 °C. Short tracks, however, are still observed. Throughout the central and western part of the island, ages are consistently below 22 Ma. Thermobarometric data for this area indicate maximum temperatures of 400±50 °C and pressures of 1±0.3 GPa. Only samples from western Crete, which have been exposed to 400±50 °C, show exclusively long tracks. Consequently, the high temperature limit of the zircon partial annealing zone (ZPAZ) appears to be between 350 and 400 °C.A significant influence of elevated confining pressure on the stability of fission tracks in zircon is ruled out by the results of annealing experiments at 0.5 GPa and at different temperatures, which fit the curves previously obtained by other authors at ambient pressure.     

Oxygen isotopic composition and U-Pb discordance in zircon

We have investigated U-Pb discordance and oxygen isotopic composition of zircon using high-spatial resolution δ¹⁸O measurement by ion microprobe. δ¹⁸O in both concordant and discordant zircon grains provides an indication of the relationship between fluid interaction and discordance. Our results suggest that three characteristics of zircon are interrelated: (1) U-Pb systematics and concomitant age discordance, (2) δ¹⁸O and the water-rock interactions implied therein, and (3) zircon texture, as revealed by cathodoluminescence and BSE imaging. A key observation is that U-Pb-disturbed zircons are often also variably depleted in ¹⁸O, but the relationship between discordance and δ¹⁸O is not systematic. δ¹⁸O values of discordant zircons are generally lighter but irregular in their distribution. Textural differences between zircon grains can be correlated with both U-Pb discordance and δ¹⁸O. Discordant grains exhibit either a recrystallized, fractured, or strongly zoned CL texture, and are characteristic of ¹⁸O depletion. We interpret this to be a result of metamictization, leading to destruction of the zircon lattice and an increased susceptibility to lead loss. Conversely, grains that are concordant have less-expressed zoning and a smoother CL texture and are enriched in ¹⁸O. From this it is apparent that various stages of water-rock interaction, as evidenced by systematic variations in δ¹⁸O, leave their imprint on both the texture and U-Pb systematics of zircon.