Three time-slice simulations of the Eemian (eem), the mid-Holocene (hol), and a pre-industrial (pri) control run with the coupled atmosphere-ocean-model ECHAM5/MPI-OM

Orbital forcing does not only exert direct insolation effects, but also alters climate indirectly through feedback mechanisms that modify atmosphere and ocean dynamics and meridional heat and moisture transfers. We investigate the regional effects of these changes by detailed analysis of atmosphere and ocean circulation and heat transports in a coupled atmosphere-ocean-sea ice-biosphere general circulation model (ECHAM5/JSBACH/MPI-OM). We perform long term quasi equilibrium simulations under pre-industrial, mid-Holocene (6000 years before present - yBP), and Eemian (125 000 yBP) orbital boundary conditions. Compared to pre-industrial climate, Eemian and Holocene temperatures show generally warmer conditions at higher and cooler conditions at lower latitudes. Changes in sea-ice cover, ocean heat transports, and atmospheric circulation patterns lead to pronounced regional heterogeneity. Over Europe, the warming is most pronounced over the north-eastern part in accordance with recent reconstructions for the Holocene. We attribute this warming to enhanced ocean circulation in the Nordic Seas and enhanced ocean-atmosphere heat flux over the Barents Shelf in conduction with retreat of sea ice and intensified winter storm tracks over northern Europe.

Solar sensitivity experiments for the pre-industrial time period using the comprehensive global climate model CCSM3 (Community Climate System Model 3)

The Southern Hemisphere Westerly Winds (SWW) constitute an important zonal circulation that influences large-scale precipitation patterns and ocean circulation. Variations in their intensity and latitudinal position have been suggested to exert a strong influence on the CO2 budget in the Southern Ocean, thus making them a potential factor affecting the global climate. The possible influence of solar forcing on SWW variability during the Holocene is addressed. Solar sensitivity experiments with a comprehensive global climate model (CCSM3) are carried out to study the response of SWW to solar variability. In addition, It is shown that a high-resolution iron record from the Chilean continental slope (41° S), which is interpreted to reflect changes in the position of the SWW, is significantly correlated with reconstructed solar activity during the past 3000 years. Taken together, the proxy and model results suggest that centennial-scale periods of lower (higher) solar activity caused equatorward (southward) shifts of the annual mean SWW.

Tropical climate and vegetation simulations during the Heinrich event 1 using an Earth System Model of Intermediate Complexity (EMIC) - the University of Victoria Earth System-Climate Model (UVic ESCM)

Abrupt climate changes from 18 to 15 thousand years before present (kyr BP) associated with Heinrich Event 1 (HE1) had a strong impact on vegetation patterns not only at high latitudes of the Northern Hemisphere, but also in the tropical regions around the Atlantic Ocean. To gain a better understanding of the linkage between high and low latitudes, we used the University of Victoria (UVic) Earth System-Climate Model (ESCM) with dynamical vegetation and land surface components to simulate four scenarios of climate-vegetation interaction: the pre-industrial era, the Last Glacial Maximum (LGM), and a Heinrich-like event with two different climate backgrounds (interglacial and glacial). We calculated mega-biomes from the plant-functional types (PFTs) generated by the model to allow for a direct comparison between model results and palynological vegetation reconstructions. Our calculated mega-biomes for the pre-industrial period and the LGM corresponded well with biome reconstructions of the modern and LGM time slices, respectively, except that our pre-industrial simulation predicted the dominance of grassland in southern Europe and our LGM simulation resulted in more forest cover in tropical and sub-tropical South America. The HE1-like simulation with a glacial climate background produced sea-surface temperature patterns and enhanced inter-hemispheric thermal gradients in accordance with the "bipolar seesaw" hypothesis. We found that the cooling of the Northern Hemisphere caused a southward shift of those PFTs that are indicative of an increased desertification and a retreat of broadleaf forests in West Africa and northern South America. The mega-biomes from our HE1 simulation agreed well with paleovegetation data from tropical Africa and northern South America. Thus, according to our model-data comparison, the reconstructed vegetation changes for the tropical regions around the Atlantic Ocean were physically consistent with the remote effects of a Heinrich event under a glacial climate background.

Model - data comparison of western African rainfall evolutions during the Last Interglacial (130-115 ka)

This dataset characterizes the evolution of western African precipitation indicated by marine sediment geochemical records in comparison to transient simulations using CCSM3 global climate model throughout the Last Interglacial (130-115 ka). It contains (1) defined tie-points (age models), newly published stable isotopes of benthic foraminifera and Al/Si log-ratios of eight marine sediment cores from the western African margin and (2) annual and seasonal rainfall anomalies (relative to pre-industrial values) for six characteristic latitudinal bands in western Africa simulated by CCSM3 (two transient simulations: one non-accelerated and one accelerated experiment).

High-resolution carbon dioxide concentration records from the EPICA Dome C ice core (EDC99)

Changes in past atmospheric carbon dioxide concentrations can be determined by measuring the composition of air trapped in ice cores from Antarctica. So far, the Antarctic Vostok and EPICA Dome C ice cores have provided a composite record of atmospheric carbon dioxide levels over the past 650,000 years. Here we present results of the lowest 200 m of the Dome C ice core, extending the record of atmospheric carbon dioxide concentration by two complete glacial cycles to 800,000 yr before present. From previously published data and the present work, we find that atmospheric carbon dioxide is strongly correlated with Antarctic temperature throughout eight glacial cycles but with significantly lower concentrations between 650,000 and 750,000 yr before present. Carbon dioxide levels are below 180 parts per million by volume (p.p.m.v.) for a period of 3,000 yr during Marine Isotope Stage 16, possibly reflecting more pronounced oceanic carbon storage. We report the lowest carbon dioxide concentration measured in an ice core, which extends the pre-industrial range of carbon dioxide concentrations during the late Quaternary by about 10 p.p.m.v. to 172-300 p.p.m.v.

Stable isotope ratios and chemistry on corals from the Caribbean Sea

Instrumental climate data are limited in length and only available with low spatial coverage before the middle of the 20th century. This is too short to reliably determine and interpret decadal and longer scale climate variability and to understand the underlying mechanisms with sufficient accuracy. A proper knowledge of past variability of the climate system is needed to assess the anthropogenic impact on climate and ecosystems, and also important with regard to long-range climate forecasting. Highly-resolved records of past climate variations that extend beyond pre-industrial times can significantly help to understand long-term climate changes and trends. Indirect information on past environmental and climatic conditions can be deduced from climate-sensitive proxies. Large colonies of massive growing tropical reef corals have been proven to sensitively monitor changes in ambient seawater. Rapid skeletal growth, typically ranging between several millimeters to centimeters per year, allows the development of proxy records at sub-seasonal resolution. Stable oxygen isotopic composition and trace elemental ratios incorporated in the aragonitic coral skeleton can reveal a detailed history of past environmental conditions, e.g., sea surface temperature (SST). In general, coral-based reconstructions from the tropical Atlantic region have lagged behind the extensive work published using coral records from the Indian and Pacific Oceans. Difficulties in the analysis of previously utilized coral archives from the Atlantic, typically corals of the genera Montastrea and Siderastrea, have so far exacerbated the production of long-term high-resolution proxy records. The objective of this study is the evaluation of massive fast-growing corals of the species Diploria strigosa as a new marine archive for climate reconstructions from the tropical Atlantic region. For this purpose, coral records from two study sites in the eastern Caribbean Sea (Guadeloupe, Lesser Antilles; and Archipelago Los Roques, Venezuela) were examined. At Guadeloupe, a century-long monthly resolved multi-proxy coral record was generated. Results present the first d18O (Sr/Ca)-SST calibration equations for the Atlantic braincoral Diploria strigosa, that are robust and consistent with previously published values using other coral species from different regions. Both proxies reflect local variability of SST on a sub-seasonal scale, which is a precondition for studying seasonally phase-locked climate variations, as well as track variability on a larger spatial scale (i.e., in the Caribbean and tropical North Atlantic). Coral Sr/Ca reliably records local annual to interannual temperature variations and is higher correlated to in-situ air temperature than to grid-SST. The warming calculated from coral Sr/Ca is concurrent with the strong surface temperature increase at the study site during the past decades. Proxy data show a close relationship to major climate signals from the tropical Pacific and North Atlantic (the El Niño Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO)) affecting the seasonal cycle of SST in the North Tropical Atlantic (NTA). Coral oxygen isotopes are also influenced by seawater d18O (d18Osw) which is linked to the hydrological cycle, and capture large-scale climate variability in the NTA region better than Sr/Ca. Results from a quantitative comparison between extreme events in the two most prominent modes of external forcing, namely the ENSO and NAO, and respective events recorded in seasonal coral d18O imply that SST variability at the study site is highly linked to Pacific and North Atlantic variability, by this means supporting the assumptions of observational- and model-based studies which suggest a strong impact of ENSO and NAO forcings onto the NTA region through a modulation of trade wind strength in winter. Results from different spectral analysis tools suggest that interannual climate variability recorded by the coral proxies is II largely dictated by Pacific ENSO forcing, whereas at decadal and longer timescales the influence of the NAO is dominan. tThe Archipelago Los Roques is situated in the southeastern Caribbean Sea, north of the Venezuelan coast. Year-to-year variations in monthly resolved coral d18O of a nearcentury- long Diploria strigosa record are significantly correlated with SST and show pronounced multidecadal variations. About half of the variance in coral d18O can be explained by variations in seawater d18O, which can be estimated by calculating the d18Oresidual via subtracting the SST component from measured coral d18O. The d18Oresidual and a regional precipitation index are highly correlated at low frequencies, suggesting that d18Osw variations are primarily atmospheric-driven. Warmer SSTs at Los Roques broadly coincide with higher precipitation in the southeastern Caribbean at multidecadal time scales, effectively strengthening the climate signal in the coral d18O record. The Los Roques coral d18O record displays a strong and statistically significant relationship to different indices of hurricane activity during the peak of the Atlantic hurricane season in boreal summer and is a particularly good indicator of decadal-multidecadal swings in the latter indices. In general, the detection of long-term changes and trends in Atlantic hurricane activity is hampered due to the limited length of the reliable instrumental record and the known inhomogeneity in the observational databases which result from changes in observing practice and technology over the years. The results suggest that coral-derived proxy data from Los Roques can be used to infer changes in past hurricane activity on timescales that extend well beyond the reliable record. In addition, the coral record exhibits a clear negative trend superimposed on the decadal to multidecadal cycles, indicating a significant warming and freshening of surface waters in the genesis region of tropical cyclones during the past decades. The presented coral d18O time series provides the first and, so far, longest continuous coral-based record of hurricane activity. It appears that the combination of both signals (SST and d18Osw) in coral d18O leads to an amplification of large-scale climate signals in the record, and makes coral d18O even a better proxy for hurricane activity than SST alone. Atlantic hurricane activity naturally exhibits strong multidecadal variations that are associated with the Atlantic Multidecadal Oscillation (AMO), the major mode of lowfrequency variability in the North Atlantic Ocean. However, the mechanisms underlying this multidecadal variability remain controversial, primarily because of the limited instrumental record. The Los Roques coral d18O displays strong multidecadal variability with a period of approximately 60 years that is closely related to the AMO, making the Archipelago Los Roques a very sensitive location for studying low-frequency climate variability in the Atlantic Ocean. In summary, the coral records presented in this thesis capture different key climate variables in the north tropical Atlantic region very well, indicating that fast-growing Diploria strigosa corals represent a promising marine archive for further proxy-based reconstructions of past climate variability on a range of time scales.

Glaciological characteristics and carbon dioxide record of ice cores from Antarctica

We report a decadally resolved record of atmospheric CO2 concentration for the last 1000 years, obtained from the West Antarctic Ice Sheet (WAIS) Divide shallow ice core. The most prominent feature of the pre-industrial period is a rapid ~7 ppm decrease of CO2 in a span of ~20-50 years at ~1600 A.D. This observation confirms the timing of an abrupt atmospheric CO2 decrease of ~10 ppm observed for that time period in the Law Dome ice core CO2 records, but the true magnitude of the decrease remains unclear. Atmospheric CO2 variations over the time period 1000-1800 A.D. are statistically correlated with northern hemispheric climate and tropical Indo-Pacific sea surface temperature. However, the exact relationship between CO2 and climate remains elusive due to regional climate variations and/or uneven geographical data density of paleoclimate records. We observe small differences of 0 ~2% (0 ~ 6 ppm) among the high-precision CO2 records from the Law Dome, EPICA Dronning Maud Land and WAIS Divide Antarctic ice cores. However, those records share common trends of CO2 change on centennial to multicentennial time scales, and clearly show that atmospheric CO2 has been increasing above preindustrial levels since ~1850 A.D.

Carbon dioxide concentrations for the Last Millennium, Antarctica

The most direct method of investigating past variations of the atmospheric CO2 concentration before 1958, when continuous direct atmospheric CO2 measurements started, is the analysis of air extracted from suitable ice cores. Here we present a new detailed CO2 record from the Dronning Maud Land (DML) ice core, drilled in the framework of the European Project for Ice Coring in Antarctica (EPICA) and some new measurements on a previously drilled ice core from the South Pole. The DML CO2 record shows an increase from about 278 to 282 parts per million by volume (ppmv) between ad 1000 and ad 1200 and a fairly continuous decrease to a mean value of about 277 ppmv around ad 1700. While the new South Pole measurements agree well with DML at the minimum at ad 1700 they are on average about 2 ppmv lower during the period ad 1000-1500. Published measurements from the coastal high-accumulation site Law Dome are considered as very reliable because of the reproducibility of the measurements, high temporal resolution and an accurate time scale. Other Antarctic ice cores could not, or only partly, reproduce the pre-industrial measurements from Law Dome. A comparison of the trends of DML and Law Dome shows a general agreement. However we should be able to rule out co-variations caused by the same artefact. Two possible effects are discussed, first production of CO2 by chemical reactions and second diffusion of dissolved air through the ice matrix into the bubbles. While the first effect cannot be totally excluded, comparison of the Law Dome and DML record shows that dissolved air diffusing to bubbles cannot be responsible for the pre-industrial variation. Therefore, the new record is not a proof of the Law Dome results but the first very strong support from an ice core of the Antarctic plateau.

Trace metals in sediments of the Var Canyon turbidite system at the DYFAMED station in the Ligurian Sea

The magnitude and the chronology of anthropogenic impregnation by Hg and other trace metals of environmental concern (V, Cr, Ni, Cu, Zn, Ag, Cd and Pb, including its stable isotopes) in the sediments are determined at the DYFAMED station, a site in the Ligurian Sea (Northwestern Mediterranean) chosen for its supposed open-sea characteristics. The DYFAMED site (VD) is located on the right levee of the Var Canyon turbidite system, at the end of the Middle Valley. In order to trace the influence of the gravity current coming from the canyon on trace metal distribution in the sediment, we studied an additional sediment core (VA) from a terrace of the Var Canyon, and material collected in sediment traps at the both sites at 20 m above sea bottom. The patterns of Hg and other trace element distribution profiles are interpreted using stable Pb isotope ratios as proxies for its sources, taking into account the sedimentary context (turbidites, redox conditions, and sedimentation rates). Major element distributions, coupled with the stratigraphic examination of the sediment cores point out the high heterogeneity of the deposits at VA, and major turbiditic events at both sites. At the DYFAMED site, we observed direct anthropogenic influence in the upper sediment layer (<2 cm), while on the Var Canyon site (VA), the anthropization concerns the whole sedimentary column sampled (19 cm). Turbiditic events superimpose their specific signature on trace metal distributions. According to the 210Pbxs-derived sedimentation rate at the DYFAMED site (0.4 mm yr-1), the Hg-enriched layer of the top core corresponds to the sediment accumulation of the last 50 years, which is the period of the highest increase in Hg deposition on a global scale. With the hypothesis of the absence of significant post-depositional redistribution of Hg, the Hg/C-org ratio changes between the surface and below are used to estimate the anthropogenic contribution to the Hg flux accumulated in the sediment. The Hg enrichment, from pre-industrial to the present time is calculated to be around 60%, consistent with estimations of global Hg models. However, based on the chemical composition of the trapped material collected in sediment traps, we calculated that epibenthic mobilization of Hg would reach 73%. Conversely, the Cd/C-org ratio decreases in the upper 5 cm, which may reflect the recent decrease of atmospheric Cd inputs or losses due to diagenetic processes.

Isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) of sediment core SO188_342 from the Northern Bay of Bengal for the last 18 ka

The Indian Summer Monsoon (ISM) is a major global climatic phenomenon. Long-term precipitation proxy records of the ISM, however, are often fragmented and discontinuous, impeding an estimation of the magnitude of precipitation variability from the Last Glacial to the present. To improve our understanding of past ISM variability, we provide a continuous reconstructed record of precipitation and continental vegetation changes from the lower Ganges-Brahmaputra-Meghna catchment and the Indo-Burman ranges over the last 18,000 years (18 ka). The records derive from a marine sediment core from the northern Bay of Bengal (NBoB), and are complemented by numerical model results of spatial moisture transport and precipitation distribution over the Bengal region. The isotopic composition of terrestrial plant waxes (dD and d13C of n-alkanes) are compared to results from an isotope-enabled general atmospheric circulation model (IsoCAM) for selected time slices (pre-industrial, mid-Holocene and Heinrich Stadial 1). Comparison of proxy and model results indicate that past changes in the dD of precipitation and plant waxes were mainly driven by the amount effect, and strongly influenced by ISM rainfall. Maximum precipitation is detected for the Early Holocene Climatic Optimum (EHCO; 10.5-6 ka BP), whereas minimum precipitation occurred during the Heinrich Stadial 1 (HS1; 16.9-15.4 ka BP). The IsoCAM model results support the hypothesis of a constant moisture source (i.e. the NBoB) throughout the study period. Relative to the pre-industrial period the model reconstructions show 20% more rain during the mid-Holocene (6 ka BP) and 20% less rain during the Heinrich Stadial 1 (HS1), respectively. A shift from C4-plant dominated ecosystems during the glacial to subsequent C3/C4-mixed ones during the interglacial took place. Vegetation changes were predominantly driven by precipitation variability, as evidenced by the significant correlation between the dD and d13C alkane records. When compared to other records across the ISM domain, precipitation and vegetation changes inferred from our records and the numerical model results provide evidence for a coherent regional variability of the ISM from the Last Glacial to the present.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Re-OS and Pb isotope ratios of sediments from Saanich Inlet, Western Canada

Re-Os and Pb-Pb isotopic analysis of reduced varved sediments cored in the deeper basin of Saanich Inlet (B.C.) are presented. From core top to 61 cm down-core, spanning approximately the last 100 yrs of sedimentation, 187Os/188Os ratio and Os concentration respectively increase from ~0.8 to ~0.9 and from 55 to 60 ppt, whereas Re concentration decreases from 3600 to 2600 ppt. Re correlates with Corg (R2=0.6) throughout the entire section, whereas Os follows Re and Corg trends deeper down-core, suggesting a decoupling of a Re- and Os-geochemistry during burial and/or very early diagenesis. No systematic compositional differences are observed between seasonal laminae. 204Pb-normalized lead isotope ratios increase from sediment surface down to 7 cm down-core, then decrease steadily to pre-industrial levels at ~50 cm down-core. This pattern illustrates the contamination from leaded gasoline until the recent past. The measured Pb isotopic ratios point primarily toward gasoline related atmospheric lead from the USA. The osmium isotopic values measured are significantly lower than those of modern seawater-Os. In comparison with other anoxic environments, the osmium content of Saanich Inlet sediments is low, and its Os isotopic composition suggests significant inputs from unradiogenic sources (detrital and/or dissolved). Ultramafic lithologies in the watershed of the Fraser River are suspected to contribute to sedimentary inputs as well as to the input of dissolved unradiogenic osmium in the water of Saanich Inlet. The presence of some unradiogenic Os from anthropogenic contamination cannot be discounted near the core top, but since deeper, pre-anthropogenic levels also yielded unradiogenic Os results, one is led to conclude that the overall low 187Os/188Os ratios result from natural geochemical processes. Thus, the bulk sediment of Saanich Inlet does not appear to record 187Os/188Os composition of the marine end-member of the only slightly below normal salinity, fjord water. The low seawater-derived Os content of the sediment, coupled with unradiogenic Os inputs from local sources, explains the overall low isotopic values observed. As a consequence, such near-shore anoxic sediments are unlikely to record changes in the past ocean Os isotopic composition.

Mineral dust deposition in interglacial and glacial climate conditions: model results

The global aerosol/climate model ECHAM5-HAM is used in order to investigate the dust cycle for four interglacial and one glacial climate conditions. The 20-year time-slices are the pre-industrial control (CTRL), mid-Holocene (6000 years BP), last glacial inception (115000 years BP), Eemian (126000 years BP) and Last Glacial Maximum (LGM) (21000 years BP) time intervals. The study is focused on the Antarctic region. The model is able to reproduce the magnitude order of dust deposition globally for the pre-industial and LGM climates. Correlation coefficient of the natural logarithm of the observed and modeled values is 0.78 for the CTRL and 0.81 for the LGM. For the pre-industrial simulation the model overestimates observed values in Antarctica by a factor of about 2-3 due to overestimation of the Australian dust source and too high wet deposition in the Antarctica interior. In the LGM, the model underestimates dust deposition in eastern Antarctica by a factor of about 4-5 due to underestimation of the South American dust source. More records are needed to validate dust deposition for the past interglacial time-slices. The modeled results show that dust deposition in Antarctica in the past interglacial time-slices is higher than in the CTRL simulation. The largest increase of dust deposition in Antarctica is simulated for the LGM, showing about 10-fold increase compared to CTRL.

Seawater carbonate chemistry in Ischia, Italy, 2008

The atmospheric partial pressure of carbon dioxide (pCO2) will almost certainly be double that of pre-industrial levels by 2100 and will be considerably higher than at any time during the past few million years1. The oceans are a principal sink for anthropogenic CO2 where it is estimated to have caused a 30% increase in the concentration of H+ in ocean surface waters since the early 1900s and may lead to a drop in seawater pH of up to 0.5 units by 2100. Our understanding of how increased ocean acidity may affect marine ecosystems is at present very limited as almost all studies have been in vitro, short-term, rapid perturbation experiments on isolated elements of the ecosystem4, 5. Here we show the effects of acidification on benthic ecosystems at shallow coastal sites where volcanic CO2 vents lower the pH of the water column. Along gradients of normal pH (8.1-8.2) to lowered pH (mean 7.8-7.9, minimum 7.4-7.5), typical rocky shore communities with abundant calcareous organisms shifted to communities lacking scleractinian corals with significant reductions in sea urchin and coralline algal abundance. To our knowledge, this is the first ecosystem-scale validation of predictions that these important groups of organisms are susceptible to elevated amounts of pCO2. Sea-grass production was highest in an area at mean pH 7.6 (1,827 µatm pCO2) where coralline algal biomass was significantly reduced and gastropod shells were dissolving due to periods of carbonate sub-saturation. The species populating the vent sites comprise a suite of organisms that are resilient to naturally high concentrations of pCO2 and indicate that ocean acidification may benefit highly invasive non-native algal species. Our results provide the first in situ insights into how shallow water marine communities might change when susceptible organisms are removed owing to ocean acidification.

Transient Holocene experiments under orbital forcing using the comprehensive global climate model CCSM3 (Community Climate System Model 3)

The Southern Hemisphere Westerly Winds (SWW) have been suggested to exert a critical influence on global climate through wind-driven upwelling of deep water in the Southern Ocean and the potentially resulting atmospheric CO2 variations. The investigation of the temporal and spatial evolution of the SWW along with forcings and feedbacks remains a significant challenge in climate research. In this study, the evolution of the SWW under orbital forcing from the early Holocene (9 kyr BP) to pre-industrial modern times is examined with transient experiments using the comprehensive coupled global climate model CCSM3. Analyses of the model results suggest that the annual and seasonal mean SWW were subject to an overall strengthening and poleward shifting trend during the course of the early-to-late Holocene under the influence of orbital forcing, except for the austral spring season, where the SWW exhibited an opposite trend of shifting towards the equator.