Were last glacial climate events simultaneous between Greenland and western Europe?

During the last glacial period, several large abrupt climate fluctuations took place on the Greenland ice cap and elsewhere. Often these Dansgaard/Oeschger events are assumed to have been synchronous, and then used as tie-points to link chronologies between the proxy archives. However, if temporally separate events are lumped into one illusionary event, climatic interpretations of the tuned events will obviously be flawed. Here, we compare Dansgaard/Oeschger-type events in a well-dated record from south-eastern France with those in Greenland ice cores. Instead of assuming simultaneous climate events between both archives, we keep their age models independent. Even these well-dated archives possess large chronological uncertainties, that prevent us from inferring synchronous climate events at decadal to multi-centennial time scales. If possible, tuning of proxy archives should be avoided.

Denudation rates of small transient catchments controlled by former glaciation: The Hörnli nunatak in the northeastern Swiss Alpine Foreland

Extensive glaciers repeatedly occupied the northern Alpine Foreland during the Pleistocene and left a strongly glacially overprinted low slope landscape. Only few islands appeared as nunataks standing above the surface of the large piedmont glacier lobes. These nunatak areas kept their original shape, manifested in steep catchments with mean slopes up to 33 . Even though not glaciated, these catchments where significantly affected by base-level changes occurring as a consequence of phases of glacier advances and retreats. Both domains, the glacially eroded and non-eroded, are therefore prone to different mechanisms and time-scales of fluvial and colluvial re-adjustment. In this study we investigate these effects by exploring the spatial distribution and magnitude of denudation in the Hörnli region of the eastern Swiss Alpine Foreland in the present Interglacial. The area represents both domains in a relatively small area with largely uniform tectonic, lithologic and climatic conditions. The differences in Holocene andscape evolution are investigated using topographic analyses and catchment-averaged denudation rates derived from 10Be concentrations in fluvial quartz sand. We find that in formerly non-glaciated, fluvially dominated catchments close hillslope-channel coupling prevails and that these catchments yield high average denudation rates of 350 mm/ka. Glacially overprinted catchments yielded catchment-wide denudation rates an order of magnitude lower. These low denudation rates are hypothesized to be the consequence of both (i) a dominance of slow hillslope processes and (ii) admixture of high concentration, pre-LGM glacial sediment. This suggests that a) a careful field investigation must accompany the denudation rate studies and b) that the concept of area-weighted cosmogenic nuclide denudation rates must be considered in light of the predominant catchment processes.

Observations of the northern seasonal polar cap on Mars II: HiRISE photometric analysis of evolution of northern polar dunes in spring

We present an overview of our analyses of HiRISE observations of spring evolution of selected dune areas of the north polar erg. The north polar erg is covered annually by seasonal volatile ice layer, a mixture of CO2 and H2O with mineral dust contamination. In spring, this layer sublimes creating visually enigmatic phenomena, e.g. dark and bright fan-shaped deposits, dark–bright–dark bandings, dark down-slope streaks, and seasonal polygonal cracks. Similar phenomena in southern polar areas are believed to be related to the specific process of solid-state greenhouse effect. In the north, it is currently unclear if the solid-state greenhouse effect is able to explain all the observed phenomena especially because the increased influence of H2O on the time scales of this process has not yet been quantified. HiRISE observations of our selected locations show that the ground exhibits a temporal behaviour similar to the one observed in the southern polar areas: a brightening phase starting close to the spring equinox with a subsequent darkening towards summer solstice. The resolution of HiRISE enabled us to study dunes and substrate individually and even distinguish between different developments on windward and slip face sides of single dunes. Differences in the seasonal evolution between steep slip faces and flatter substrate and windward sides of dunes have been identified and compared to CRISM data of CO2 and H2O distributions on dunes. We also observe small scale dark blotches that appear in early observations and tend to sustain a low reflectivity throughout the spring. These blotches can be regarded as the analogue of dark fan deposits in southern polar areas, leading us to the conclusion that both martian polar areas follow similar spring evolutions.

Lateral sediment sources and knickzones as controls on spatio-temporal variations of sediment transport in an Alpine river

Modern mixed alluvial-bedrock channels in mountainous areas provide natural laboratories for understanding the time scales at which coarse-grained material has been entrained and transported from their sources to the adjacent sedimentary sink, where these deposits are preserved as conglomerates. This article assesses the shear stress conditions needed for the entrainment of the coarse-bed particles in the Glogn River that drains the 400 km2 Val Lumnezia basin, eastern Swiss Alps. In addition, quantitative data are presented on sediment transport patterns in this stream. The longitudinal stream profile of this river is characterized by three ca 500 m long knickzones where channel gradients range from 0·02 to 0·2 m m−1, and where the valley bottom confined into a <10 m wide gorge. Downstream of these knickzones, the stream is flat with gradients <0·01 m m−1 and widths ≥30 m. Measurements of the grain-size distribution along the trunk stream yield a mean D84 value of ca 270 mm, whereas the mean D50 is ca 100 mm. The consequences of the channel morphology and the grain-size distribution for the time scales of sediment transport were explored by using a one-dimensional step-backwater hydraulic model (Hydrologic Engineering Centre – River Analysis System). The results reveal that, along the entire trunk stream, a two to 10 year return period flood event is capable of mobilizing both the D50 and D84 fractions where the Shields stress exceeds the critical Shields stress for the initiation of particle motion. These return periods, however, varied substantially depending on the channel geometry and the pebble/boulder size distribution of the supplied material. Accordingly, the stream exhibits a highly dynamic boulder cover behaviour. It is likely that these time scales might also have been at work when coarse-grained conglomerates were constructed in the geological past.

Two Stable Equilibria of the Atlantic Subpolar Gyre

The cyclonic circulation of the Atlantic subpolar gyre is a key mechanism for North Atlantic climate variability on a wide range of time scales. It is generally accepted that it is driven by both cyclonic winds and buoyancy forcing, yet the individual importance and dynamical interactions of the two contributions remain unclear. The authors propose a simplified four-box model representing the convective basin of the Labrador Sea and its shallow and deep boundary current system, the western subpolar gyre. Convective heat loss drives a baroclinic flow of relatively light water around the dense center. Eddy salt flux from the boundary current to the center increases with a stronger circulation, favors the formation of dense waters, and thereby sustains a strong baroclinic flow, approximately 10%–25% of the total. In contrast, when the baroclinic flow is not active, surface waters may be too fresh to convect, and a buoyancy-driven circulation cannot develop. This situation corresponds to a second stable circulation mode. A hysteresis is found for variations in surface freshwater flux and the salinity of the near-surface boundary current. An analytical solution is presented and analyzed.

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Anisotropic diffusion at the field scale in a 4-year multi-tracer diffusion and retention experiment – I: Insights from the experimental data

Abstract Claystones are considered worldwide as barrier materials for nuclear waste repositories. In the Mont Terri underground research laboratory (URL), a nearly 4-year diffusion and retention (DR) experiment has been performed in Opalinus Clay. It aimed at (1) obtaining data at larger space and time scales than in laboratory experiments and (2) under relevant in situ conditions with respect to pore water chemistry and mechanical stress, (3) quantifying the anisotropy of in situ diffusion, and (4) exploring possible effects of a borehole-disturbed zone. The experiment included two tracer injection intervals in a borehole perpendicular to bedding, through which traced artificial pore water (APW) was circulated, and a pressure monitoring interval. The APW was spiked with neutral tracers (HTO, HDO, H2O-18), anions (Br, I, SeO4), and cations (Na-22, Ba-133, Sr-85, Cs-137, Co-60, Eu-152, stable Cs, and stable Eu). Most tracers were added at the beginning, some were added at a later stage. The hydraulic pressure in the injection intervals was adjusted according to the measured value in the pressure monitoring interval to ensure transport by diffusion only. Concentration time-series in the APW within the borehole intervals were obtained, as well as 2D concentration distributions in the rock at the end of the experiment after overcoring and subsampling which resulted in �250 samples and �1300 analyses. As expected, HTO diffused the furthest into the rock, followed by the anions (Br, I, SeO4) and by the cationic sorbing tracers (Na-22, Ba-133, Cs, Cs-137, Co-60, Eu-152). The diffusion of SeO4 was slower than that of Br or I, approximately proportional to the ratio of their diffusion coefficients in water. Ba-133 diffused only into �0.1 m during the �4 a. Stable Cs, added at a higher concentration than Cs-137, diffused further into the rock than Cs-137, consistent with a non-linear sorption behavior. The rock properties (e.g., water contents) were rather homogeneous at the centimeter scale, with no evidence of a borehole-disturbed zone. In situ anisotropy ratios for diffusion, derived for the first time directly from field data, are larger for HTO and Na-22 (�5) than for anions (�3�4 for Br and I). The lower ionic strength of the pore water at this location (�0.22 M) as compared to locations of earlier experiments in the Mont Terri URL (�0.39 M) had no notable effect on the anion accessible pore fraction for Cl, Br, and I: the value of 0.55 is within the range of earlier data. Detailed transport simulations involving different codes will be presented in a companion paper.

Atmospheric CO₂ response to volcanic eruptions: The role of ENSO, season, and variability

Tropical explosive volcanism is one of the most important natural factors that significantly impact the climate system and the carbon cycle on annual to multi-decadal time scales. The three largest explosive eruptions in the last 50�years�Agung, El Chichón, and Pinatubo�occurred in spring/summer in conjunction with El Niño events and left distinct negative signals in the observational temperature and CO2 records. However, confounding factors such as seasonal variability and El Niño-Southern Oscillation (ENSO) may obscure the forcing-response relationship. We determine for the first time the extent to which initial conditions, i.e., season and phase of the ENSO, and internal variability influence the coupled climate and carbon cycle response to volcanic forcing and how this affects estimates of the terrestrial and oceanic carbon sinks. Ensemble simulations with the Earth System Model (Climate System Model 1.4-carbon) predict that the atmospheric CO2 response is �60 larger when a volcanic eruption occurs during El Niño and in winter than during La Niña conditions. Our simulations suggest that the Pinatubo eruption contributed 11�±�6 to the 25�Pg terrestrial carbon sink inferred over the decade 1990�1999 and �2�±�1 to the 22�Pg oceanic carbon sink. In contrast to recent claims, trends in the airborne fraction of anthropogenic carbon cannot be detected when accounting for the decadal-scale influence of explosive volcanism and related uncertainties. Our results highlight the importance of considering the role of natural variability in the carbon cycle for interpretation of observations and for data-model intercomparison.

An Ensemble Kalman Filter multi-tracer assimilation: Determining uncertain ocean model parameters for improved climate-carbon cycle projections

An Ensemble Kalman Filter is applied to assimilate observed tracer fields in various combinations in the Bern3D ocean model. Each tracer combination yields a set of optimal transport parameter values that are used in projections with prescribed CO2 stabilization pathways. The assimilation of temperature and salinity fields yields a too vigorous ventilation of the thermocline and the deep ocean, whereas the inclusion of CFC-11 and radiocarbon improves the representation of physical and biogeochemical tracers and of ventilation time scales. Projected peak uptake rates and cumulative uptake of CO2 by the ocean are around 20% lower for the parameters determined with CFC-11 and radiocarbon as additional target compared to those with salinity and temperature only. Higher surface temperature changes are simulated in the Greenland–Norwegian–Iceland Sea and in the Southern Ocean when CFC-11 is included in the Ensemble Kalman model tuning. These findings highlights the importance of ocean transport calibration for the design of near-term and long-term CO2 emission mitigation strategies and for climate projections.

Long-term climate change commitment and reversibility: An EMIC Intercomparison

This paper summarizes the results of an intercomparison project with Earth System Models of Intermediate Complexity (EMICs) undertaken in support of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5). The focus is on long-term climate projections designed to 1) quantify the climate change commitment of different radiative forcing trajectories and 2) explore the extent to which climate change is reversible on human time scales. All commitment simulations follow the four representative concentration pathways (RCPs) and their extensions to year 2300. Most EMICs simulate substantial surface air temperature and thermosteric sea level rise commitment following stabilization of the atmospheric composition at year-2300 levels. The meridional overturning circulation (MOC) is weakened temporarily and recovers to near-preindustrial values in most models for RCPs 2.6-6.0. The MOC weakening is more persistent for RCP8.5. Elimination of anthropogenic CO2 emissions after 2300 results in slowly decreasing atmospheric CO2 concentrations. At year 3000 atmospheric CO2 is still at more than half its year-2300 level in all EMICs for RCPs 4.5-8.5. Surface air temperature remains constant or decreases slightly and thermosteric sea level rise continues for centuries after elimination of CO2 emissions in all EMICs. Restoration of atmospheric CO2 from RCP to preindustrial levels over 100-1000 years requires large artificial removal of CO2 from the atmosphere and does not result in the simultaneous return to preindustrial climate conditions, as surface air temperature and sea level response exhibit a substantial time lag relative to atmospheric CO2.

Climatic drivers of hourly to yearly tree radius variations along a 6°C natural warming gradient

Climate affects the timing, rate and dynamics of tree growth, over time scales ranging from seconds to centuries. Monitoring how a tree's stem radius varies over these time scales can provide insight into intra-annual stem dynamics and improve our understanding of climate impacts on tree physiology and growth processes. Here, we quantify the response of radial conifer stem size to environmental fluctuations via a novel assessment of tree circadian cycles. We analyze four years of sub-hourly data collected from 56 larch and spruce trees growing along a natural temperature gradient of ∼6 °C in the central Swiss Alps. During the growing season, tree stem diameters were greatest at mid-morning and smallest in the late evening, reflecting the daily cycle of water uptake and loss. Along the gradient, amplitudes calculated from the stem radius cycle were ∼50% smaller at the upper site (∼2200 m a.s.l.) relative to the lower site (∼800 m a.s.l.). We show changes in precipitation, temperature and cloud cover have a substantial effect on typical growing season diurnal cycles; amplitudes were nine times smaller on rainy days (>10 mm), and daily amplitudes are approximately 40% larger when the mean daily temperature is 15–20 °C than when it is 5–10 °C. We find that over the growing season in the sub-alpine forests, spruce show greater daily stem water movement than larch. However, under projected future warming, larch could experience up to 50% greater stem water use, which may severely affect future growth on already dry sites. Our data further indicate that because of the confounding influences of radial growth and short-term water dynamics on stem size, conventional methodology probably overstates the effect of water-linked meteorological variables (i.e. precipitation and relative humidity) on intra-annual tree growth. We suggest future studies use intra-seasonal measurements of cell development and consider whether climatic factors produce reversible changes in stem diameter. These study design elements may help researchers more accurately quantify and attribute changes in forest productivity in response to future warming.

Vegetation responses to rapid warming and to minor climatic fluctuations during the Late-Glacial Interstadial (GI-1) at Gerzensee (Switzerland)

High-resolution pollen analyses made on the same samples on which the ratios of oxygen isotopes were measured that provided the time scale and a temperature proxy after correlation to NorthGRIP. (1) A primary succession: The vegetation responded to the rapid rise of temperatures around 14,685 yr BP, with a primary succession on a decadal to centennial time scale. The succession between ca 15,600 and 13,000 yr BP included: (1.1.) The replacement of shrub-tundra by woodland of Juniperus and tree birch (around 14,665 yr BP) (1.2.) The response of Juniperus pollen to the shift in oxygen isotopes in less than 20 yr, (1.3.) A sequence of population increases of Hippophaë rhamnoides (ca 14,600 yr BP), Salix spp. (ca 14,600 yr BP), Betula trees (ca.14,480 yr BP), Populus cf. tremula (ca. 14,300 yr BP), and Pinus cf. sylvestris (ca. 13,830 yr BP). (2) Biological processes: Plants responded to the rapid increase of summer temperatures on all organisational levels: (2.1) Individuals may have produced more pollen (e.g. Juniperus); (2.2) Populations increased or decreased (e.g. Juniperus, Betula, later Pinus), and (2.3) Populations changed their biogeographical range and may show migrational lags. (2.4) Plant communities changed in their composition because the species pools changed through immigration and (local) extinction. Some plant communities may have been without modern analogue.These mechanisms require increasing amounts of time. (2.5) Processes on the level of ecosystems, with species interactions, may involve various time scales. Besides competition and facilitation, nitrogen fixation is discussed. (3) The minor fluctuations of temperature during the Late-Glacial Interstadial, which are recorded in δ18O, resulted in only very minor changes in pollen during the Aegelsee Oscillation (Older Dryas biozone, GI-1d) and the Gerzensee Oscillation (GI-1b). (4) Biodiversity: The afforestation at the onset of Bølling coincided with a gradual increase of taxonomic diversity up to the time of the major Pinus expansion.

Local Rates of Ice-Sheet Thickness Change in Greenland

The rate of ice-sheet thickness change is calculated for 10 sites in Greenland by comparing measured values of ice vertical velocity and snow-accumulation rate. Vertical velocities are derived from repeat surveys of markers using precision global positioning system techniques, and accumulation rates are determined from stratigraphic analysis of firn cores. The results apply to time-scales covered by the firn-core records, which in most cases are a few decades. A spectrum of thickness-change rates is obtained, ranging from substantial thinning to slow thickening. The sites where ice-sheet thinning is indicated are located near the ice-sheet margin or in outlet glacier catchments. Interior and high-elevation sites are predominantly in balance or thickening slowly. Uncertainties in the rates of thickness change are dominated by errors in the determination of accumulation rates. The results of this work are broadly comparable with regional estimates of mass balance obtained from the analysis of catchment input vs discharge.

Barry Saltzman and the Theory of Climate

Barry Saltzman was a giant in the fields of meteorology and climate science. A leading figure in the study of weather and climate for over 40 yr, he has frequently been referred to as the "father of modern climate theory." Ahead of his time in many ways, Saltzman made significant contributions to our understanding of the general circulation and spectral energetics budget of the atmosphere, as well as climate change across a wide spectrum of time scales. In his endeavor to develop a unified theory of how the climate system works, lie played a role in the development of energy balance models, statistical dynamical models, and paleoclimate dynamical models. He was a pioneer in developing meteorologically motivated dynamical systems, including the progenitor of Lorenz's famous chaos model. In applying his own dynamical-systems approach to long-term climate change, he recognized the potential for using atmospheric general circulation models in a complimentary way. In 1998, he was awarded the Carl-Gustaf Rossby medal, the highest honor of the American Meteorological Society "for his life-long contributions to the study of the global circulation and the evolution of the earth's climate." In this paper, the authors summarize and place into perspective some of the most significant contributions that Barry Saltzman made during his long and distinguished career. This short review also serves as an introduction to the papers in this special issue of the Journal of Climate dedicated to Barry's memory.

Terahertz macrospin dynamics in insulating ferrimagnets

We investigate numerically the excitation of nonlinear magnetic interactions in a ferrite material by an energetic pump pulse of terahertz (THz) radiation. The calculations are performed by solving the coupled Maxwell and Landau-Lifshitz-Gilbert differential equations. In a time-resolved THz pump/THz probe scheme, it is demonstrated that Faraday rotation of a delayed THz probe pulse can be used to map these interactions. Our study is motivated by the ability of soft x-ray free electron lasers to perform time-resolved imaging of the magnetization process at the submicrometer and subpicosecond length and time scales.