Are sea-level-rise trends along the coasts of the north Indian Ocean consistent with global estimates?

Mean-sea-level data from coastal tide gauges in the north Indian Ocean were used to show that low-frequency variability is consistent among the stations in the basin. Statistically significant trends obtained from records longer than 40 years yielded sea-level-rise estimates between 1.06–1.75 mm yr^− 1, with a regional average of 1.29 mm yr^− 1, when corrected for global isostatic adjustment (GIA) using model data. These estimates are consistent with the 1–2 mm yr^− 1 global sea-level-rise estimates reported by the Intergovernmental Panel on Climate Change.     

Reconstruction of Mediterranean sea level fields for the period 1945–2000

The distribution of sea level in the Mediterranean Sea is recovered for the period 1945–2000 by using a reduced space optimal interpolation analysis. The method involves estimating empirical orthogonal functions from satellite altimeter data spanning the period 1993–2005 that are then combined with tide gauge data to recover sea level fields over the period 1945–2000. The reconstruction technique is discussed and its robustness is checked through different tests. For the altimetric period (1993–2000) the prediction skill is quantified over the whole domain by comparing the reconstructed fields with satellite altimeter observations. For past times the skill can only be tested locally, by validating the reconstruction against independent tide gauge records. The reconstructed distribution of sea level trends for the period 1945–2000 shows a positive peak in the Ionian Sea (up to 1.5 mm yr^− 1) and a negative peak of − 0.5 mm yr^− 1 in a small area to the south-east of Crete. Positive trends are found nearly everywhere, being larger in the western Mediterranean (between 0.5 and 1 mm yr^− 1) than in the eastern Mediterranean (between 0 and 0.5 mm yr^− 1). The estimated rate of mean sea level rise for the period 1945–2000 is 0.7 ± 0.2 mm yr^− 1, i.e. about a half of the rate estimated for global mean sea level. These overall results do not appear to be very sensitive to the distribution of tide gauges. The poorest results are obtained in open-sea regions with intense mesoscale variability not correlated with any tide gauge station, such as the Algerian Basin.     

Cryospheric change in China

This paper provides an overview of the current status of the cryosphere in China and its changes. Up-to-date statistics of the cryosphere in China are summarized based on the latest available data. There are 46,377 glaciers in China, covering an area of 59,425 km². The glacier ice reserve is estimated to be about 5600 km³ and the annual glacier runoff is about 61.6 × 10⁹ m³. The continuous snow cover extent (> 60 days) in China is about 3.4 × 10⁶ km² and the maximum water equivalent is 95.9 × 10⁹ m³ yr^− 1. The permafrost area in China is about 1.72 × 10⁶ km². The total ground ice reserve on the Qinghai–Tibetan Plateau is estimated to be about 10,923 km³. Recent investigations indicated that glacier areas in China have shrunk about 2–10% over the past 45 yr. Total glacier area has receded by about 5.5%. Snow mass has increased slightly. Permafrost is clearly degrading, as indicated by shrinking areas of permafrost, increasing depth of the active layer, rising of lower limit of permafrost, and thinning of the seasonal frost depth. Some models predict that glacier area shrinkage could be as high as 26.7% in 2050, with glacier runoff increasing until its maximum in about 2030. Although snow mass shows an increasing trend in western China, in eastern China the trend is toward decreasing snow mass, with increasing interannual fluctuations. Permafrost degradation is likely to continue, with one-third to one-half of the permafrost on the Qinghai–Tibetan Plateau anticipated to degrade by 2100. Most of the high-temperature permafrost will disappear by then. The permafrost in northeastern China will retreat further northward.     

Comparison of Mediterranean sea level fields for the period 1961–2000 as given by a data reconstruction and a 3D model

Two Mediterranean sea level distributions spanning the last decades are examined. The first one is a reconstruction of sea level obtained by a reduced-space optimal interpolation applied to tide gauge and altimetry data. The second distribution is obtained from a 3D (baroclinic) regional circulation model. None of the two representations includes the mechanical atmospheric forcing. Results are presented for two different periods: 1993–2000 (for which altimetry data are available) and 1961–2000 (the longest period common to both distributions).The first period is examined as a test period for the model, since the reconstruction is very similar to altimetry observations. The modelled sea level is in fair agreement with the reconstruction in the Western Mediterranean and in the Aegean Sea (except in the early nineties), but in the Ionian Sea the model departs from observations. For the whole period 1961–2000 the main feature is a marked positive trend in the Ionian Sea (up to 1.8 mm yr^− 1), observed both in the reconstruction and in the model. Also the distribution of positive trends in the Western Mediterranean (mean value of 1.1 mm yr^− 1) and the smaller trends in the Aegean Sea (0.5 mm yr^− 1) are similar in the reconstruction and in the model, despite the first implicitly accounts for sea level variations due to remote sources such as ice melting and the second does not. The interannual sea level variability associated with key regional events such as the Eastern Mediterranean Transient is apparently captured by the reconstruction but not by the model (at least in its present configuration). Hence, the reconstruction can be envisaged as a useful tool to validate further long-term numerical simulations in the region.     

A multitracer study of peat profiles from Tunguska, Siberia

Two peat columns from Tunguska (Siberia) were analysed for pollen, spores, charcoal, trace elements and γ-emitters in order to identify the fingerprints of the impact of a still unidentified cosmic body (TCB), which occurred in the summer of 1908, and the level of environmental pollution in a background area of central Siberia. Peat layers were subject to non-destructive γ-ray spectrometry to derive radiochronology by the excess ²¹⁰Pb method. The age-to-depth relationship was crosschecked by using both 1963 horizon of ¹³⁷Cs associated to maximum global fallout deposition and palynological data profiles. Vertical distributions of trace elements in the peat columns were obtained by PIXE multielemental analysis allowing determination of the levels of environmental contamination in a background region of the Siberian taiga.The association of heavy metals such as Ni, Co and Cu in the profiles suggests the connection of the area with mining and metal smelting activity in the north of the region through atmospheric circulation. As concerns global scale contamination, the inventory of the artificial radionuclide ¹³⁷Cs (4.6 kBq m^− 2) shows a value typical of remote slightly contaminated areas resulting from global scale redistribution of radioactive fallout from Cold War nuclear weapon testing. The atmospheric inventory of the natural radionuclide ²¹⁰Pb, for which a mean annual flux of 200 Bq m^− 2 yr^− 1 has been calculated, is typical of continental regions.The influence of Tunguska Cosmic Body in the peat is recognizable by a large discontinuity in the palynological profile of the peat monolith at a depth coinciding with the 1908 layer as determined by the ²¹⁰Pb technique, showing a large peak of total pollen counting attributed to the impact of the shockwave on the area in which huge tree stands were destroyed. Following the event, tree pollen concentration decreases abruptly showing the temporary inception of a mire environment with an increase of Sphagnum spore concentrations. Results of elemental analysis so far available do not show anomalies in the concentration profiles at depths coinciding with the Tunguska event layer indicating the need for pre-concentration technique enabling the detection of element associations typical of extraterrestrial materials.     

Climate forced atmospheric CO2 variability in the early Holocene: A stomatal frequency reconstruction

The dynamic climate in the Northern Hemisphere during the early Holocene could be expected to have impacted on the global carbon cycle. Ice core studies however, show little variability in atmospheric CO₂. Resolving any possible centennial to decadal CO₂ changes is limited by gas diffusion through the firn layer during bubble enclosure. Here we apply the inverse relationship between stomatal index (measured on sub-fossil leaves) and atmospheric CO₂ to complement ice core records between 11,230 and 10,330 cal. yr BP. High-resolution sampling and radiocarbon dating of lake sediments from the Faroe Islands reconstruct a distinct CO₂ decrease centred on ca. 11,050 cal. yr BP, a consistent and steady decline between ca. 10,900 and 10,600 cal. yr BP and an increased instability after ca. 10,550 cal. yr BP. The earliest decline lasting ca. 150 yr is probably associated with the Preboreal Oscillation, an abrupt climatic cooling affecting much of the Northern Hemisphere a few hundred years after the end of the Younger Dryas. In the absence of known global climatic instability, the decline to ca. 10,600 cal. yr BP is possibly due to expanding vegetation in the Northern Hemisphere. The increasing instability in CO₂ after 10,600 cal. yr BP occurs during a period of increasing cooling of surface waters in the North Atlantic and some increased variability in proxy climate indicators in the region.The reconstructed CO₂ changes also show a distinct similarity to indicators of changing solar activity. This may suggest that at least the Northern Hemisphere was particularly sensitive to changes in solar activity during this time and that atmospheric CO₂ concentrations fluctuated via rapid responses in climate.     

Climate changes and recent glacier behaviour in the Chilean Lake District

Atmospheric temperatures measured at the Chilean Lake District (38°–42°S) showed contrasting trends during the second half of the 20th century. The surface cooling detected at several meteorological stations ranged from − 0.014 to − 0.021 °C a^− 1, whilst upper troposphere (850–300 gpm) records at radiosonde of Puerto Montt (41°26′S/73°07′W) revealed warming between 0.019 and 0.031 °C a^− 1. Regional rainfall data collected from 1961 to 2000 showed the overall decrease with a maximum rate of − 15 mm a^− 2 at Valdivia st. (39°38′S/73°05′W). These ongoing climatic changes, especially the precipitation reduction, seem to be related to El Niño–Southern Oscillation (ENSO) phenomena which has been more frequent after 1976. Glaciers within the Chilean Lake District have significantly retreated during recent decades, in an apparent out-of-phase response to the regional surface cooling. Moreover, very little is known about upper troposphere changes and how they can enhance the glacier responses. In order to analyse their behaviour in the context of the observed climate changes, Casa Pangue glacier (41°08′S/71°52′W) has been selected and studied by comparing Digital Elevation Models (DEMs) computed at three different dates throughout the last four decades. This approach allowed the determination of ice elevation changes between 1961 and 1998, yielding a mean thinning rate of − 2.3 ± 0.6 m a^− 1. Strikingly, when ice thinning is computed for the period between 1981 and 1998, the resulting rate is 50% higher (− 3.6 ± 0.6 m a^− 1). This enhanced trend and the related area loss and frontal retreat suggests that Casa Pangue might currently be suffering negative mass balances in response to the upper troposphere warming and decreased precipitation of the last 25–30 yr, as well as debris cover would not prevent the glacier from a fast reaction to climate forcing. Most of recent glaciological studies regarding Andean glaciers have concentrated on low altitude changes, namely frontal variations, however, in order to better understand the regional glacier changes, new data are necessary, especially from the accumulation areas.     

Surface processes and tectonics: Forcing of continental subduction and deep processes

It is now well accepted that surface processes provide a critical feedback on the surface tectonic deformation, whatever it is, orogenic building or basin evolution. However, the idea that the influence of these processes may go below the crustal levels, is less common. In this preliminary study, we use coupled thermo-mechanical numerical models to investigate the possible influence of surface processes on the styles of continental collision, in particular, continental subduction. For that, we further exploit the recent successful model of continental subduction of the early stages of India–Asia collision by Toussaint et al. [Toussaint G., Burov, E., and J.-P. Avouac, Tectonic evolution of a continental collision zone: a thermo-mechanical numerical model, Tectonics, 23, TC6003, doi:10.1029/2003TC001604, 2004b.]. On the example of India–Asia-like settings, we show that not only the surface topography but also the total amount of subduction may largely vary as function of denudation rate (controlled by the coefficient of erosion, k). Erosion provides a dynamic discharge of the hanging wall of the major thrust zone, whereas the sedimentation increases loading on the footwall and this helps down-thrusting of the lower plate. Both processes reduce the resistance of the major thrust and subduction channel to subduction. However, very strong or very slow erosion/sedimentation enhance the possibility of plate coupling and promote whole-scale thickening or buckling. The maximal amount of subduction is thus achieved for some intermediate values of erosion rates when the tectonic uplift rate is fine-balanced by the denudation rate. In our case the optimal balance is reached for the values of k on the order of 3000 m²/yr. We then extended our model beyond the conditions of India–Asia collision, in terms of the tested range of k and convergence rates. The experiments suggest that for provided settings, both extra slow (k < 50–100 m²/yr) and extra rapid erosion (k > 6000–8000 m²/yr) limit, by up to 50%, the total amount of subduction, if not totally prevent it. The model demonstrates the large capability of surface processes to adopt to different deformation styles: the orogenic building and subduction successfully develop (subduction number, S > 0.5) in the range of k between 500 m²/yr and 6000 m²/yr at convergence rates ranging from 1 cm/yr to 6 cm/yr. Within this range, some peculiar features of orogenic style such as the geometry of the accretion prism, amount of upper crustal subduction, horizontal progression of the mountain range/thrust fault and the amount of exhumation of metamorphic facies are sometimes quite different. We conclude that surface processes may control deep, mantle level tectonic evolution.     

Reconstructing the Younger Dryas ice dammed lake in the Baltic Basin: Bathymetry, area and volume

A digital 3D-reconstruction of the Baltic Ice Lake's (BIL) configuration during the termination of the Younger Dryas cold phase (ca. 11 700 cal. yr BP) was compiled using a combined bathymetric–topographic Digital Terrain Model (DTM), Scandinavian ice sheet limits, Baltic Sea Holocene bottom sediment thickness information, and a paleoshoreline database maintained at the Lund University. The bathymetric–topographic DTM, assembled from publicly available data sets, has a resolution of 500 × 500 m on Lambert Azimuthal Equal Area projection allowing area and volume calculations of the BIL to be made with an unprecedented accuracy. When the damming Scandinavian ice sheet margin eventually retreated north of Mount Billingen, the high point in terrain of Southern central Sweden bordering to lower terrain further to the north, the BIL was catastrophically drained resulting in a 25 m drop of the lake level. With our digital reconstruction, we estimate that approximately 7800 km³ of water drained during this event and that the ice dammed lake area was reduced by ca. 18%. Building on previous results suggesting drainage over 1 to 2 years, our lake volume calculations imply that the freshwater flux to the contemporaneous sea in the west was between about 0.12 and 0.25 Sv. The BIL reconstruction provides new detailed information on the paleogeography in the area of southern Scandinavia, both before and after the drainage event, with implications for interpretations of geological records concerning the post-glacial environmental development.     

Recent evolution and mass balance of Cordón Martial glaciers, Cordillera Fueguina Oriental

Past and present glacier changes have been studied at Cordón Martial, Cordillera Fueguina Oriental, Tierra del Fuego, providing novel data for the Holocene deglaciation history of southern South America and extrapolating as well its future behavior based on predicted climatic changes. Regional geomorphologic and stratigraphic correlations indicate that the last glacier advance deposited the ice-proximal (“internal”) moraines of Cordón Martial, around 330 ¹⁴C yr BP, during the Late Little Ice Age (LLIA). Since then glaciers have receded slowly, until 60 years ago, when major glacier retreat started. There is a good correspondence for the past 100 years between the surface area variation of four small cirque glaciers at Cordón Martial and the annual temperature and precipitation data of Ushuaia. Between 1984 and 1998, Martial Este Glacier lost 0.64 ± 0.02 × 10⁶ m³ of ice mass (0.59 ± 0.02 × 10⁶ m³ w.e.), corresponding to an average ice thinning of 7.0 ± 0.2 m (6.4 ± 0.2 m w.e), according to repeated topographic mapping. More detailed climatic data have been obtained since 1998 at the Martial Este Glacier, including air temperature, humidity and solar radiation. These records, together with the monthly mass balance measured since March 2000, document the annual response of the Martial Este Glacier to the climate variation. Mass balances during hydrological years were positive in 2000, negative in 2001 and near equilibrium in 2002. Finally, using these data and the regional temperature trend projections, modeled for different future scenarios by the Atmosphere-Ocean Model (GISS-NASA/GSFC), potential climatic-change effects on this mountain glacier were extrapolated. The analysis shows that only the Martial Este Glacier may survive this century.     

The effect of climate change on carbon in Canadian peatlands

Peatlands, which are dominant features of the Canadian landscape, cover approximately 1.136 million km², or 12% of the land area. Most of the peatlands (97%) occur in the Boreal Wetland Region (64%) and Subarctic Wetland Region (33%). Because of the large area they cover and their high organic carbon content, these peatlands contain approximately 147 Gt soil carbon, which is about 56% of the organic carbon stored in all Canadian soils.A model for estimating peatland sensitivity to climate warming was used to determine both the sensitivity ratings of various peatland areas and the associated organic carbon masses. Calculations show that approximately 60% of the total area of Canadian peatlands and 51% of the organic carbon mass in all Canadian peatlands is expected to be severely to extremely severely affected by climate change.The increase in average annual air temperature of 3–5 °C over land and 5–7 °C over the oceans predicted for northern Canada by the end of this century would result in the degradation of frozen peatlands in the Subarctic and northern Boreal wetland regions and severe drying in the southern Boreal Wetland Region. In addition, flooding of coastal peatlands is expected because of the predicted rise in sea levels. As a result of these changes, a large part of the carbon in the peatlands expected to be severely and extremely severely affected by climate change could be released into the atmosphere as carbon dioxide (CO₂) and methane (CH₄), which will further increase climate warming.     

APEX-SZ first light and instrument status

The apex-sz instrument is designed for the discovery and study of galaxy clusters at mm-wavelengths using the Sunyaev Zel’dovich effect. The receiver consists of 320 superconducting transition edge sensor (TES) bolometers cooled to 250 mK with the combination of a three stage He sorption fridge and mechanical pulse tube cooler. The detectors are instrumented with a frequency domain multiplexing readout system. The receiver is mounted on the 12 m apex telescope located at 5100 m on the Atacama plateau in Chile. For the first light engineering deployment of December 2005, the receiver was configured with a 55 element wedge of the bolometers and operating in the 150 GHz atmospheric window. During the engineering run we achieved significant milestones in our instrumentation development efforts, including celestial observations with a monolithically fabricated TES bolometer array cooled with a mechanical cooler and successful implementation of a SQUID-based MHz AC-biased readout. These technology demonstrations point the way toward future large TES bolometer array instruments. Here we describe the results of this deployment and future plans for the apex-sz instrument.     

Estimation of net accumulation rate at a Patagonian glacier by ice core analyses using snow algae

Snow algae in a 45.97-m-long ice core from the Tyndall Glacier (50°59′05″S, 73°31′12″W, 1756 m a.s.l.) in the Southern Patagonian Icefield were examined for potential use in ice core dating and estimation of the net accumulation rate. The core was subjected to visual stratigraphic observation and bulk density measurements in the field, and later to analyses of snow algal biomass, water isotopes (¹⁸O, D), and major dissolved ions. The ice core contained many algal cells that belonged to two species of snow algae growing in the snow near the surface: Chloromonas sp. and an unknown green algal species. Algal biomass and major dissolved ions (Na⁺, K⁺, Mg²⁺, Ca²⁺, Cl⁻, SO₄²⁻) exhibited rapid decreases in the upper 3 m, probably owing to melt water elution and/or decomposition of algal cells. However, seasonal cycles were still found for the snow algal biomass, ¹⁸O, D-excess, and major ions, although the amplitudes of the cycles decreased with depth. Supposing that the layers with almost no snow algae were the winter layers without the melt water essential to algal growth, we estimated that the net accumulation rate at this location was 12.9 m a^− 1 from winter 1998 to winter 1999, and 5.1 m from the beginning of winter to December 1999. These estimates are similar to the values estimated from the peaks of ¹⁸O (17.8 m a^− 1 from summer 1998 to summer 1999 and 11.0 m from summer to December 1999) and those of D-excess (14.7 m a^− 1 from fall 1998 to fall 1999 and 8.6 m a^− 1 from fall to December 1999). These values are much higher than those obtained by past ice core studies in Patagonia, but are of the same order of magnitude as those predicted from various observations at ablation areas of Patagonian glaciers.     

Modeling mineral dust emissions from Chinese and Mongolian deserts

The present study investigates the frequency and intensity of mineral dust emissions over the deserts of eastern Asia from 1996 to 2001. Mineral dust emissions are simulated using a physical dust emission scheme over a region extending from 35.5°N to 47°N and from 73°E to 125°E. The input parameters required by the dust emission model are (1) surface features data including aerodynamic roughness length, soil dry size distribution and texture; and (2) meteorological surface data, mainly wind speed, soil moisture and snow cover. The way by which these surface features and meteorological data can be assessed is described and discussed. The influence of soil moisture and snow cover is taken into account and their effects on simulated dust emission are quantified.The simulations reproduce on a daily basis the location and intensity of the severe events of April 1998 and spring 2001 as recorded by the meteorological stations and/or described in various studies. Based on 6 yr of simulations, the main dust source regions are identified and their relative contributions to the total dust emissions are quantified.The seasonal cycle of the dust storms frequency is well reproduced with a maximum in spring. The simulations suggest that it is mainly controlled by the emissions occurring in the Taklimakan desert in latter spring and in summer, and by those occurring in the northern deserts of China in winter. The Taklimakan desert appears to be the most frequent and steady source of dust emissions during the studied period. On the other hand, in the Gobi desert, only a few dust emission events are simulated, but the dust amount emitted during each event is generally very large. In the northern deserts of China, dust emissions are frequent and their intensity is variable.These results show an important annual and inter-annual variability of the emitted dust (between 100 Mt yr^− 1 and 460 Mt yr^− 1), mainly controlled by the occurrence of severe events in the Gobi desert and in the northern deserts of China.     

Airborne laser altimetry survey of Glaciar Tyndall, Patagonia

The first airborne laser altimetry measurements of a glacier in South America are presented. Data were collected in November of 2001 over Glaciar Tyndall, Torres del Paine National Park, Chilean Patagonia, onboard a Twin Otter airplane of the Chilean Air Force. A laser scanner with a rotating polygon-mirror system together with an Inertial Navigation System (INS) were fixed to the floor of the aircraft, and used in combination with two dual-frequency GPS receivers. Together, the laser–INS–GPS system had a nominal accuracy of 30 cm after data processing. On November 23rd, a total of 235 km were flown over the ablation area of Glaciar Tyndall, with 5 longitudinal tracks with a mean swath width of 300 m, which results in a point spacing of approximately 2 m both along and across track. A digital elevation model (DEM) generated using the laser altimetry data was compared with a DEM produced from a 1975 map (1:50,000 scale — Instituto Geográfico Militar (IGM), Chile). A mean thinning of − 3.1 ± 1.0 m a^− 1 was calculated for the ablation area of Glaciar Tyndall, with a maximum value of − 7.7 ± 1.0 m a^− 1 at the calving front at 50 m a.s.l. and minimum values of between − 1.0 and − 2.0 ± 1.0 m a^− 1 at altitudes close to the equilibrium line altitude (900 m a.s.l.). The thinning rates derived from the airborne survey were similar to the results obtained by means of ground survey carried out at  600 m of altitude on Glaciar Tyndall between 1975 and 2002, yielding a mean thinning of − 3.2 m a^− 1 [Raymond, C., Neumann, T.A., Rignot, E., Echelmeyer, K.A., Rivera, A., Casassa, G., 2005. Retreat of Tyndall Glacier, Patagonia, over the last half century. Journal of Glaciology 173 (51), 239–247.]. A good agreement was also found between ice elevation changes measured with laser data and previous results obtained with Shuttle Radar Topography Mission (SRTM) data. We conclude that airborne laser altimetry is an effective means for accurately detecting glacier elevation changes in Patagonia, where an ice thinning acceleration trend has been observed during recent years, presumably in response to warming and possibly also drier conditions.     

Observations of the temperature and polarization anisotropies with Boomerang 2003

The Boomerang experiment completed its final long duration balloon (LDB) flight over Antarctica in January 2003. The focal plane was upgraded to accommodate four sets of 145 GHz polarization sensitive bolometers (PSBs), identical to those to be flown on the Planck HFI instrument. Approximately, 195 hours of science observations were obtained during this flight, including 75 hours distributed over 1.84% of the sky and an additional 120 hours concentrated on a region covering 0.22% of the sky. We derive the angular power spectra of the cosmic microwave background (cmb) temperature and polarization anisotropies from these data. The temperature anisotropies are detected with high signal to noise on angular scales ranging from several degrees to 10 arcminutes. The curl-free (EE) component is detected at 4.8σ, and a two-sigma upper limit on the curl (BB) component of 8.6 μK² is obtained on scales corresponding to 0.5°. Both the temperature and polarization anisotropies are found to be consistent with a concordance ΛCDM cosmology that is seeded by adiabatic density perturbations. In addition to the cmb observations, Boomerang03 surveyed a 300 square degree region centered on the Galactic plane. These observations represent the first light for polarization sensitive bolometers, which are currently operational in two South-Pole based polarimeters, as well as Planck HFI, at frequencies ranging from 100 to 350 GHz (3 mm to 850 μm).     

Recent organic matter accumulation in relation to some climatic factors in ombrotrophic peat bogs near heavy metal emission sources in Finland

Accumulation of organic matter (OM) was studied in four ombrotrophic peat bogs in Finland: Harjavalta (vicinity of a Cu–Ni smelter), Outokumpu (near a closed Cu–Ni mine), Alkkia (Ni-treated site) and Hietajärvi (a pristine site). At each sampling site, two peat cores (15 × 15 × 100 cm) were taken. Age-dating of peat was determined using ²¹⁰Pb method (CRS model). The local annual temperature sum and precipitation for the past 125 years were modeled. The objective was to compare recent net accumulation rates of heavy metal polluted ombrotrophic peat bogs with those of a pristine bog, and to study the relationship between weather and net accumulation rates. Based on ²¹⁰Pb age-dating, the upper 16-cm peat layer at Harjavalta, 35 cm at Outokumpu and 25 cm at Hietajärvi represents 125 years of peat formation, yielding the following average peat accumulation rates: Harjavalta 1.3 mm year^− 1, Outokumpu 2.8 mm year^− 1 and Hietajärvi 2.0 mm year^− 1. At the Alkkia site, the Ni treatment in 1962 had completely stopped the peat accumulation. Net accumulation rates were related to precipitation at Outokumpu, Harjavalta and Hietajärvi sites. In addition, emissions released from the nearby located Cu–Ni smelter could have affected negatively net OM accumulation rate at Harjavalta site.     

Effects of a 2 × CO2 climate on two large lake systems: Pyramid Lake, Nevada, and Yellowstone Lake, Wyoming

The possible effects of trace-gas induced climatic changes on Pyramid and Yellowstone Lakes are assessed using a model of lake temperature. The model is driven by years of hourly meteorological data obtained directly from the output of double-CO₂ experiments (2 × CO₂) conducted with a regional climate model nested in a general circulation model. The regional atmospheric model is the climate version of the National Center for Atmospheric Research/Pennsylvania State University mesoscale model, MM4.Average annual surface temperature of Pyramid Lake for the 2 × CO₂ climate is 15.5 ± 5.4°C (±1 σ), 2.8°C higher than the control. Annual overturn of the lake ceases as a result of these higher temperatures for the 2 × CO₂ climate. Evaporation increases from 1400 mm yr⁻¹ in the control to 1595 mm yr⁻¹ in the 2 × CO₂ simulation, but net water supplied to the Pyramid Lake basin increases from −6 mm yr⁻¹ in the control to +27 mm yr⁻¹ in the 2 × CO₂ simulation due to increased precipitation.For the open water periods, the average annual surface temperature of Yellowstone Lake is 13.2 ± 5.1°C for the 2 × CO₂ climate, a temperature 1.6°C higher than the control. The annual duration of ice cover on the lake is 152 days in the 2 × CO₂ simulation, a reduction of 44 days relative to the control. Warming of the lake for the 2 × CO₂ climate is mostly confined to the near-surface. Simulated spring overturn for the 2 × CO₂ climate occurs earlier in the year and fall overturn later than in the control. Evaporation increases from 544 mm yr⁻¹ to 600 mm yr⁻¹ in the 2 × CO₂ simulation, but net water supplied to the Yellowstone Lake basin increases from +373 mm yr⁻¹ in the control to +619 mm yr⁻¹ due to increased precipitation. The effects of these climatic changes suggest possible deterioration of water quality and productivity in Pyramid Lake and possible enhancement of productivity in Yellowstone Lake.     

Palaeoecological constraints on late Glacial and Holocene ice retreat in the Southern Andes (53°S)

Late Glacial to Holocene ice retreat was investigated along a 120 km long fjord system, reaching from Gran Campo Nevado (GCN) to Seno Skyring in the southernmost Andes (53°S). The aim was to improve the knowledge on regional and global control on glacier recession with special emphasis on latitudinal shifting of the westerlies. The timing of ice retreat was derived from peat and sediment cores, using mineralogical and chemical characteristics, and pollen as proxies. Stratigraphy was based on ¹⁴C-AMS ages and tephrochronology. The ice retreat of the Seno Skyring Glacier lobe is marked by an ice rafted debris layer which was formed around 18,300 to 17,500 cal. yr B.P. Subsequently, fast glacier retreat occurred until around 15,000 to 14,000 cal. yr B.P. during which around 84% of Skyring Glacier were lost. This fast recession was probably also triggered by an increase of the Equilibrium Line Altitude (ELA) from 200 to 300 m. Subsequently, the ice surface was lowered below the ELA in an area that previously made up more than 50% of the accumulation area. Much slower retreat and glacier fluctuations of limited extent in the fjord channel system northeast of GCN occurred between around 14,000 to 11,000 cal. yr B.P. during both the Antarctic Cold Reversal and the Younger Dryas. This slow down of retreat indicates a decline in the general warming trend and/or increased precipitation, due to a southward migration of the westerlies. After around 11,000 cal. yr B.P. pollen distribution shows evolved Magellanic Rainforest and similar climate as at present, which lasted throughout most of the Holocene. Only Late Neoglacial moraine systems were formed in the period 1220–1460 AD, and subsequently in the 1620s AD, and between 1870 and 1910 AD. The results indicate that the Gran Campo Nevado ice cap has reacted more sensitive and partly distinct to climate change, compared to the Patagonian Ice Field.     

The significance of coral reefs as global carbon sinks— response to Greenhouse

Coral reefs are net sinks for C, principally as CaCO₃ accretion. It is possible to predict quite accurately the rate of production, given adequate information about any particular reef environment. The best data set for an extensive region is that for the Great Barrier Reef (GBR). Careful analysis of this region and the incorporation of previously documented present day system calcification rates suggest net production (G) from G = 1 (kg CaCO₃ m⁻² yr⁻¹) for fringing reefs, to G = 1.9 for planar (infiled platform) reefs, G = 3 for ribbon reefs and lagoonal reefs. The 20,055 km² of reefs in the GBR are thus estimated to average G = 2.4, resulting in a total production of 50 million tonnes yr⁻¹. In a 50–100 year Greenhouse scenario of rising sealevel, we predict that recolonisation of present day reef flats will be extensive and prolific. Production will increase substantially, and this could be by as much as 40% (ranging from 0% for deep shoals to 180% for fringing reefs) to give 70 million tonnes yr⁻¹ if the rate of sealevel rise reaches or exceeds 6–8 mm yr⁻¹We estimate 115,000 km² of oceanic atolls worldwide. Drawing on points equivalence from the detailed analysis of the GBR, we estimate the atolls presently produce 160 million tonnes yr⁻¹. We predict that a similar 40% increase could be possible in the next 100 years or so resulting in a production of 220 million tonnes.Accepting an existing estimate of 617,000 km² for reefs worldwide, drawing from our projections from the GBR and the atolls, and making some assumptions about the remaining reef types (we suggest fringing reefs to dominate) we estimate global reef production at the present time to be 900 million tonnes yr⁻¹. Within the next 100 years or so, we suggest this rate could almost double to 1800 million tonnes. In the long term (several centuries) we predict that the continuing trend of recolonisation, particularly of fringing and planar reefs could result in the production of > 3000 million tonnes yr⁻¹ if rates of sealevel rise approaching or exceeding 6–8 mm yr⁻¹ are achieved. Eventually (> 500 yr), reefs could actually “drown” due to inability to match the rate of sealevel increase if that rate significantly exceeds 6–8 mm yr⁻¹.Thus, coral reefs at present act as a sink for 111 million tonnes C yr⁻¹, the equivalent of 2% of present output of anthropogenic CO₂. In the short term Greenhouse scenario (100 yr) we predict this could increase to the equivalent of 4% of the present CO₂ output. In the much longer term (several centuries), if all trends continue, this could increase to the equivalent of as much as 9% of the present CO₂ output.Unfortunately, we also predict that this considerable sink for C will be most likely of negative value in alleviating Greenhouse because of the immediate effect of CaCO₃ precipitation is to raise the P_CO2 of the surface oceans — ie, ot encourage CO₂ efflux to the atmosphere. We do not attempt to quantify this effect.Other Greenhouse changes such as seawater temperature increase, changes in cloud cover, increased rainfall and runoff, increased storm activity, and changes in dissolved CO₂ concentration and surface ocean circulation may complicate the reef response. However, we suggest that sealevel rise will be the dominant influence, at least during the next 100 years or so.