Climate change in the Atlantic Ocean since the 1940's and the effects on the global climate

En aquest projecte s’ha estudiat la relació entre els canvis en les temperatures superficials de l’Oceà Atlàntic i els canvis en la circulació atmosfèrica en el segle XX. Concretament s’han analitzat dos períodes de estudi: el primer des del 1940 al 1960 i el segon des del 1980 fins al 2000. S’ha posat especial interès en les anomalies en les temperatures superficials del mar en la regió tropical de l’Oceà Atlàntic i la possible interconnexió amb els canvis climàtics observats i predits. Per a la realització de l’estudi s’han dut a terme una sèrie d’experiments utilitzant el model climàtic elaborat a la universitat d’UCLA (UCLA‐AGCM model). Els resultats obtinguts han estat analitzats en forma de mapes i figures per a cada variable d’estudi. També s’ha fet una comparació entre els resultats obtinguts i altres trobats en altres treballs publicats sobre el mateix tema de recerca. Els resultats obtinguts són molt amplis i poden tenir diverses interpretacions. Tot i així algunes de les conclusions a les quals s’ha arribat són: les diferències més significatives per a les variables estudiades i trobades a partir dels resultats obtinguts del model per als dos períodes d’estudi són en els mesos d’hivern i a la zona dels tròpics; concretament a parts del nord de sud Amèrica i a parts del nord d’Àfrica. S’han trobat també canvis significatius en els patrons de precipitació sobre aquestes mateixes zones. També s’ha observant un moviment cap al nord de la zona d’interconvergència tropical i pot ser degut a l’anòmal gradient trobat a la zona equatorial en les temperatures superficial de l’Oceà. Tot i així per a una definitiva discussió i conclusions sobre els resultats dels experiments, seria necessari un estudi més ampli i profund.

Early Cretaceous life, climate and anoxia

Early Cretaceous life and the environment were strongly influenced by the accelerated break up of Pangaea, which was associated with the formation of a multitude of rift basins, intensified spreading, and important volcanic activity on land and in the sea. These processes likely interacted with greenhouse conditions, and Early Cretaceous climate oscillated between "normal" greenhouse, predominantly arid conditions, and intensified greenhouse, predominantly humid conditions. Arid conditions were important during the latest Jurassic and early Berriasian, the late Barremian, and partly also during the late Aptian. Humid conditions were particularly intense and widespread during shorter episodes of environmental change (EECs): the Valanginian Weissert, the latest Hauterivian Faraoni, the latest Barremian earliest Aptian Taxy, the early Aptian Selli, the early late Aptian Fallot and the late Aptian-early Albian Paquier episodes. Arid conditions were associated with evaporation, low biogeochemical weathering rates, low nutrient fluxes, and partly stratified oceans, leading to oxygen depletion and enhanced preservation of laminated, organic-rich mud (LOM). Humid conditions enabled elevated biogeochemical weathering rates and nutrient fluxes, important runoff and the buildup of freshwater lids in proximal basins, intensified oceanic and atmospheric circulation, widespread upwelling and phosphogenesis, important primary productivity and enhanced preservation of LOM in expanded oxygen-minimum zones. The transition of arid to humid climates may have been associated with the net transfer of water to the continent owing to the infill of dried-out groundwater reservoirs in internally drained inland basins. This resulted in shorter-term sea-level fall, which was followed by sea-level rise. These sea-level changes and the influx of freshwater into the ocean may have influenced oxygen-isotope signatures. Climate change preceding and during the Early Cretaceous EECs may have been rapid, but in general, the EECs had a "pre"-history, during which the stage was set for environmental change. Negative feedback on the climate through increased marine LOM preservation was unlikely, because of the low overall organic-carbon accumulation rates during these episodes. Life and climate co-evolved during the Early Cretaceous. Arid conditions may have affected continental life, such as across the Tithonian/Berriasian boundary. Humid conditions and the corresponding tendency to develop dys- to anaerobic conditions in deeper ocean waters led to phases of accelerated extinction in oceans, but may have led to more luxuriant vegetation cover on continents, such as during the Valanginian, to the benefit of herbivores. During Early Cretaceous EECs, reef systems and carbonate platforms in general were particularly vulnerable. They were the first to disappear and the last to recover, often only after several million years. (C) 2011 Elsevier Ltd. All rights reserved.

Contribution pollinique à l'étude du pléistocène de la région lémanique (Suisse) : paléoenvironnement des derniers 800'000 ans

Abstact : Pollinic contribution to a Pleistocene study of the Lake Geneva region (Switzerland) Palaeoenvironment of the last 800,000 years. Although the "Plateau Romand" (tableland of French-speaking Switzerland) has been marked by the repetitive action of different Quaternary glacial thrusts, numerous preserved outcrops here and there, thanks to favourable topographical conditions, provide sufficient pollinic content to reconstitute the significant botanical fluctuations and allow bio- and chrono-stratigraphical correlations with the long European sequences. The present study examines around ten Pleistocene deposits rich in pollens, whose presumed ages cover a range of 800,000 years. - The Early and Early-Middle Pleistocene are examined at Ecoteaux (VD), undeniably the oldest site in this study. The "Lower Formation of Ecoteaux", whose low pollinic content reveals an environment of the cold desert type, shows an inverse palaeomagnetic remanence. An age earlier than 780,000 years BP (Matuyama period or earlier) is proposed. The "Upper Formation of Ecoteaux" contains over eight major fluctuations in the vegetal covering, of which at least four temporal phases and one climatic optimum. The presence of botanical genera, relics of the Tertiary in particular Pterocarya and Carya, associated with the vegetal dynamic made up of short oscillations, contrives to link this formation to the Cromerian Complex. The thermomeres recorded could correspond to one or several interglacial periods from Middle Pleistocene. Possible correlations with marine isotopic stages 15, 17 and 19 are discussed. - The Middle and Late Pleistocene is approached through the sediments of Onnens (VD), Creux d'Enfer (FR), Port-Valais (VS), La Dénériaz (VD) and Cortaillod (NE). Each site provides a temperate interglacial flora, whose pollinic content is evaluated according to the different conceivable bio-stratigraphic correlations. - Finally the Würm is studied thanks to the outcrop of Marly (FR), together with the excavations at Villars-sous-Yens (VD), Versoix (GE) and Ollon (VD). Together, these observations tend to show that the west of the Swiss Plateau is a zone of transition between the meridian vegetation of the lower Rhone valley and that of beyond the Jura. At the outlet of the Rhone glacier, and on the borders of the prealpine and Jurassian glaciers, the region is situated at the crossroads of post-glacial botanical migrations, between the major axis of the Rhone valley and the Eastern refuge zones; this situation confers particular importance on it. Although believed to be deeply scraped by ice, one notices today that the "Plateau Romand" harbours overtwenty interglacial and interstadial sites. The abrasive action of the Rhone glacier has thus been more modest than previously estimated. The effective extension of glacial tongues, together with atmospheric circulation, are doubtless varying from one interglacial to another. Finally, the speeds of glacial movements, and the unbelievable vegetal vitality, are important parameters, capable of leading to deposits, very early enriched in pollens. RÉSUMÉS : Contribution pollinique à l'étude du Pléistocène de la région lémanique (Suisse) - Paléoenvironnement des derniers 800'000 ans Bien que le Plateau romand soit marqué par l'action répétitive des différentes poussées glaciaires quaternaires, de nombreux affleurements préservés ça et là, à la faveur de conditions topographiques favorables, livrent un contenu pollinique suffisant pour reconstituer des fluctuations botaniques significatives et autoriser des corrélations bioet chronostratigraphiques avec les longues séquences européennes. Le présent travail étudie une dizaine de dépôts pléistocènes riches en pollens, dont les âges présumés balaient un éventail de plus de 800'000 ans. - Le Pléistocène Ancien et Moyen inférieur, sont examinés à Ecoteaux (VD); indéniablement le site le plus ancien de cette étude. La "Formation inférieure d'Ecoteaux", dont le maigre contenu pollinique traduit un environnement de type désertique froid, présente une rémanence magnétique inverse. Un âge antérieur à 780'000 ans BP (période de Matuyama ou antérieure) est proposé. La "Formation supérieure d'Ecoteaux" comprend plus de huit fluctuations majeures du couvert végétal, parmi lesquelles on compte au moins quatre épisodes tempérés et un optimum climatique. La présence de genres botaniques reliques du Tertiaire, en particulier Pterocarya et Carya, associée à la dynamique végétale faite de courtes oscillations, concourent à lier cette formation au Complexe Cromérien. Les thermomères enregistrées pourraient correspondre à un- ou plusieurs interglaciaires du Pléistocène moyen. Les corrélations possibles avec les stades isotopiques marins, MIS 15, 17 et 19 sont discutées. - Le Pléistocène Moyen et Récent est abordé à travers les sédiments d'Onnens (VD), du Creux d'Enfer (FR), de Port-Valais (VS), de La Dénériaz (VD) et de Cortaillod (NE). Chaque site livre une flore tempérée interglaciaire, dont le contenu pollinique est évalué en fonction des différentes corrélations bio-stratigraphiques envisageables. - Enfin le Würm est étudié grâce à l'affleurement de Marly (FR), accompagné des forages de Villars-sous-Yens (VD), Versoix (GE) et Ollon (VD). L'ensemble de ces observations tend à montrer que l'ouest du Plateau suisse est une zone charnière entre la végétation méridionale de la basse vallée du Rhône et celle d'au-delà du Jura. Au débouché du glacier du Rhône et aux confins des glaciers préalpins et jurassiens, cette région est située au carrefour des migrations botaniques post-glaciaires, entre l'axe majeur de la vallée du Rhône et les zones refuges de l'Est. Cette situation lui confère une importance toute particulière. Alors qu'on pensait le Plateau romand profondément raboté par les glaces, on constate aujourd'hui qu'il recèle plus d'une vingtaine de sites interglaciaires et interstadiaires. L'action abrasive du glacier rhodanien a donc été plus modeste qu'envisagée jusqu'ici. L'extension effective des langues glaciaires, tout comme la circulation atmosphérique, diffèrent sans doute d'un interglaciaire à l'autre. Enfin la vitesse des mouvements glaciaires et l'incroyable vitalité végétale sont des paramètres importants, susceptibles d'avoir entraîné des dépôts, très précocement enrichis en pollens.

Development of a Technique for Estimating Rain Parameters from Rain Noise Measurements

The measurement of global precipitation is of great importance in climate modeling since the release of latent heat associated with tropical convection is one of the pricipal driving mechanisms of atmospheric circulation.Knowledge of the larger-scale precipitation field also has important potential applications in the generation of initial conditions for numerical weather prediction models Knowledge of the relationship between rainfall intensity and kinetic energy, and its variations in time and space is important for erosion prediction. Vegetation on earth also greatly depends on the total amount of rainfall as well as the drop size distribution (DSD) in rainfall.While methods using visible,infrared, and microwave radiometer data have been shown to yield useful estimates of precipitation, validation of these products for the open ocean has been hampered by the limited amount of surface rainfall measurements available for accurate assessement, especially for the tropical oceans.Surface rain fall measurements(often called the ground truth)are carried out by rain gauges working on various principles like weighing type,tipping bucket,capacitive type and so on.The acoustic technique is yet another promising method of rain parameter measurement that has many advantages. The basic principle of acoustic method is that the droplets falling in water produce underwater sound with distinct features, using which the rainfall parameters can be computed. The acoustic technique can also be used for developing a low cost and accurate device for automatic measurement of rainfall rate and kinetic energy of rain.especially suitable for telemetry applications. This technique can also be utilized to develop a low cost Disdrometer that finds application in rainfall analysis as well as in calibration of nozzles and sprinklers. This thesis is divided into the following 7 chapters, which describes the methodology adopted, the results obtained and the conclusions arrived at.

A Study on Aerosol Transport Dynamics over the Indian Landmass and the Adjoining Oceanic Regions Using Satellite Data

The present study brings out the influence of transport dynamics on the aerosol distribution over the Indian region at a few selected geographically distinct locations. Over the Bay of Bengal the dominant pathway of aerosol transport during the pre-monsoon period is through higher altitudes (~ 3 km); directed from the Indian main land. In contrast, the aerosol pathways over the Arabian Sea during the same period are quite complex. They are directed from geographically different environments around the ocean through different altitudes. However in general, the day-to-day variability of AOD at both these regions is significantly influenced by the features of atmospheric circulation especially, the wind convergence at higher altitudes (around 3 km). Over the Ganga Basin during the winter period, the wind convergence at lower altitudes (< I km) govems the shon term variations in AOD, while the mean AOD distribution at this location is mainly governed by the local anthropogenic sources.

How may tropical cyclones change in a warmer climate?

Tropical Cyclones (TC) under different climate conditions in the Northern Hemisphere have been investigated with the Max Planck Institute (MPI) coupled (ECHAM5/MPIOM) and atmosphere (ECHAM5) climate models. The intensity and size of the TC depend crucially on resolution with higher wind speed and smaller scales at the higher resolutions. The typical size of the TC is reduced by a factor of 2.3 from T63 to T319 using the distance of the maximum wind speed from the centre of the storm as a measure. The full three dimensional structure of the storms becomes increasingly more realistic as the resolution is increased. For the T63 resolution, three ensemble runs are explored for the period 1860 until 2100 using the IPCC SRES scenario A1B and evaluated for three 30 year periods at the end of the 19th, 20th and 21st century, respectively. While there is no significant change between the 19th and the 20th century, there is a considerable reduction in the number of the TC by some 20% in the 21st century, but no change in the number of the more intense storms. Reduction in the number of storms occurs in all regions. A single additional experiment at T213 resolution was run for the two latter 30-year periods. The T213 is an atmospheric only experiment using the transient Sea Surface Temperatures (SST) of the T63 resolution experiment. Also in this case, there is a reduction by some 10% in the number of simulated TC in the 21st century compared to the 20th century but a marked increase in the number of intense storms. The number of storms with maximum wind speeds greater than 50ms-1 increases by a third. Most of the intensification takes place in 2 the Eastern Pacific and in the Atlantic where also the number of storms more or less stays the same. We identify two competing processes effecting TC in a warmer climate. First, the increase in the static stability and the reduced vertical circulation is suggested to contribute to the reduction in the number of storms. Second, the increase in temperature and water vapor provide more energy for the storms so that when favorable conditions occur, the higher SST and higher specific humidity will contribute to more intense storms. As the maximum intensity depends crucially on resolution, this will require higher resolution to have its full effect. The distribution of storms between different regions does not, at first approximation, depend on the temperature itself but on the distribution of the SST anomalies and their influence on the atmospheric circulation. Two additional transient experiments at T319 resolution where run for 20 years at the end of the 20th and 21st century, respectively using the same conditions as in the T213 experiments. The results are consistent with the T213 study. The total number of tropical cyclones were similar to the T213 experiment but were generally more intense. The change from the 20th to the 21st century was also similar with fewer TC in total but with more intense cyclones.

Trends in aerosol optical depth in the Russian Arctic and their links with synoptic climatology

Temporal and spatial variability of aerosol optical depth (AOD) are examined using observations of direct solar radiation in the Eurasian Arctic for 1940-1990. AOD is estimated using empirical methods for 14 stations located between 66.2 degrees N and 80.6 degrees N, from the Kara Sea to the Chukchi Sea. While AOD exhibits a well-known springtime maximum and summertime minimum at all stations, atmospheric turbidity is higher in spring in the western (Kara-Laptev) part of the Eurasian Arctic. Between June and August, the eastern (East Siberian-Chukchi) sector experiences higher transparency than the western part. A statistically significant positive trend in AOD was observed in the Kara-Laptev sector between the late 1950s and the early 1930s predominantly in spring when pollution-derived aerosol dominates the Arctic atmosphere but not in the eastern sector. Although all stations are remote, those with positive trends are located closer to the anthropogenic sources of air pollution. By contrast, a widespread decline in AOD was observed between 1982 and 1990 in the eastern Arctic in spring but was limited to two sites in the western Arctic. These results suggest that the post-1982 decline in anthropogenic emissions in Europe and the former Soviet Union has had a limited effect on aerosol load in the Arctic. The post-1982 negative trends in AOD in summer, when marine aerosol is present in the atmosphere, were more common in the west. The relationships between AOD and atmospheric circulation are examined using a synoptic climatology approach. In spring, AOD depends primarily on the strength and direction of air flow. Thus strong westerly and northerly flows result in low AOD values in the East Siberian-Chukchi sector. By contrast, strong southerly flow associated with the passage of depressions results in high A OD in the Kara-Laptev sector and trajectory analysis points to the contribution of industrial regions of the sub-Arctic. In summer, low pressure gradient or anticyclonic conditions result in high atmospheric turbidity. The frequency of this weather type has declined significantly since the early 1980s in the Kara-Laptev sector, which partly explains the decline in summer AOD values. (c) 2004 Elsevier B.V. All rights reserved.

Long-term change, interannual and intra-seasonal variability in climate and glacier mass balance in the Central Greater Caucasus, Russia

Long-term trends, interannual and intra-seasonal variability in the mass-balance record from Djankuat glacier, central Greater Caucasus, Russia, are related to local climate change, synoptic and large-scale anomalies in atmospheric circulation. A clear warming signal emerged in the central Greater Caucasus in the early 1990s, leading to a strong increase in ablation. In the absence of a compensating change in winter accumulation, the net mass balance of Djankuat has declined. The highest value of seasonal ablation on record was registered in the summer of 2000. At the beginning of the 21st century these trends reversed. Ablation was below average even in the summer of 2003, which was unusually warm in western Europe. Precipitation and winter accumulation were high, allowing for a partial recovery of net mass balance. The interannual variability in the components of mass balance is weakly related to the North Atlantic Oscillation (NAO) and the Scandinavian teleconnection patterns, but there is a clear link with the large-scale circulation anomalies represented by the Rossby pattern. Five synoptic categories have been identified for the ablation season of 2005, revealing a strong separation between components of radiation budget, air temperature and daily melt. Air temperature is the main control over melt. The highest values of daily ablation are related to the strongly positive NAO which forces high net radiation, and to the warm and moist advection from the Black Sea.

Changes in the global sulfate burden due to perturbations in global CO2 concentrations

A large ensemble of general circulation model (GCM) integrations coupled to a fully interactive sulfur cycle scheme were run on the climateprediction.net platform to investigate the uncertainty in the climate response to sulfate aerosol and carbon dioxide (CO2) forcing. The sulfate burden within the model (and the atmosphere) depends on the balance between formation processes and deposition (wet and dry). The wet removal processes for sulfate aerosol are much faster than dry removal and so any changes in atmospheric circulation, cloud cover, and precipitation will feed back on the sulfate burden. When CO2 is doubled in the Hadley Centre Slab Ocean Model (HadSM3), global mean precipitation increased by 5%; however, the global mean sulfate burden increased by 10%. Despite the global mean increase in precipitation, there were large areas of the model showing decreases in precipitation (and cloud cover) in the Northern Hemisphere during June–August, which reduced wet deposition and allowed the sulfate burden to increase. Further experiments were also undertaken with and without doubling CO2 while including a future anthropogenic sulfur emissions scenario. Doubling CO2 further enhanced the increases in sulfate burden associated with increased anthropogenic sulfur emissions as observed in the doubled CO2-only experiment. The implications are that the climate response to doubling CO2 can influence the amount of sulfate within the atmosphere and, despite increases in global mean precipitation, may act to increase it.

ENSO at 6ka and 21ka from ocean–atmosphere coupled model simulations

We analyze how the characteristics of El Niño-Southern Oscillation (ENSO) are changed in coupled ocean–atmosphere simulations of the mid-Holocene (MH) and the Last Glacial Maximum (LGM) performed as part of the Paleoclimate Modeling Intercomparison Project phase 2 (PMIP2). Comparison of the model results with present day observations show that most of the models reproduce the large scale features of the tropical Pacific like the SST gradient, the mean SST and the mean seasonal cycles. All models simulate the ENSO variability, although with different skill. Our analyses show that several relationships between El Niño amplitude and the mean state across the different control simulations are still valid for simulations of the MH and the LGM. Results for the MH show a consistent El Niño amplitude decrease. It can be related to the large scale atmospheric circulation changes. While the Northern Hemisphere receives more insolation during the summer time, the Asian summer monsoon system is strengthened which leads to the enhancement of the Walker circulation. Easterlies prevailing over the central eastern Pacific induce an equatorial upwelling that damps the El Niño development. Results are less conclusive for 21ka. Large scale dynamic competes with changes in local heat fluxes, so that model shows a wide range of responses, as it is the case in future climate projections.

Performance of the OPA/ARPEGE-T21 global ocean-atmosphere coupled model

The climatology of the OPA/ARPEGE-T21 coupled general circulation model (GCM) is presented. The atmosphere GCM has a T21 spectral truncation and the ocean GCM has a 2°×1.5° average resolution. A 50-year climatic simulation is performed using the OASIS coupler, without flux correction techniques. The mean state and seasonal cycle for the last 10 years of the experiment are described and compared to the corresponding uncoupled experiments and to climatology when available. The model reasonably simulates most of the basic features of the observed climate. Energy budgets and transports in the coupled system, of importance for climate studies, are assessed and prove to be within available estimates. After an adjustment phase of a few years, the model stabilizes around a mean state where the tropics are warm and resemble a permanent ENSO, the Southern Ocean warms and almost no sea-ice is left in the Southern Hemisphere. The atmospheric circulation becomes more zonal and symmetric with respect to the equator. Once those systematic errors are established, the model shows little secular drift, the small remaining trends being mainly associated to horizontal physics in the ocean GCM. The stability of the model is shown to be related to qualities already present in the uncoupled GCMs used, namely a balanced radiation budget at the top-of-the-atmosphere and a tight ocean thermocline.

General concepts in meteorology and dynamics

The aim of this chapter is to give a general overview of the atmospheric circulation, highlighting the main concepts that are important for a basic understanding of meteorology and atmospheric dynamics relevant to atmospheric data assimilation.

Enhanced signature of solar variability in Eurasian winter climate

We demonstrate that a new geomagnetic index of solar variability exhibits stronger correlations with atmospheric circulation variations than conventional measures. The circulation anomalies are particularly enhanced over the North Atlantic / Eurasian sector, where there are large changes in the occurrence of blocking and the winter mean surface temperature differs by several degrees between high- and low-solar terciles. The relationship is also simpler, being largely linear between high- and low-solar winters. While the circulation anomalies strongly resemble the North Atlantic Oscillation they also extend deeper into Eurasia, in a distinct signature which may be useful for the detection and attribution of observed changes and also the identification of dynamical mechanisms.

Transient teleconnection event at the onset of a planet-encircling dust storm on Mars

We use proper orthogonal decomposition (POD) to study a transient teleconnection event at the onset of the 2001 planet-encircling dust storm on Mars, in terms of empirical orthogonal functions (EOFs). There are several differences between this and previous studies of atmospheric events using EOFs. First, instead of using a single variable such as surface pressure or geopotential height on a given pressure surface, we use a dataset describing the evolution in time of global and fully three-dimensional atmospheric fields such as horizontal velocity and temperature. These fields are produced by assimilating Thermal Emission Spectrometer observations from NASA's Mars Global Surveyor spacecraft into a Mars general circulation model. We use total atmospheric energy (TE) as a physically meaningful quantity which weights the state variables. Second, instead of adopting the EOFs to define teleconnection patterns as planetary-scale correlations that explain a large portion of long time-scale variability, we use EOFs to understand transient processes due to localised heating perturbations that have implications for the atmospheric circulation over distant regions. The localised perturbation is given by anomalous heating due to the enhanced presence of dust around the northern edge of the Hellas Planitia basin on Mars. We show that the localised disturbance is seemingly restricted to a small number (a few tens) of EOFs. These can be classified as low-order, transitional, or high-order EOFs according to the TE amount they explain throughout the event. Despite the global character of the EOFs, they show the capability of accounting for the localised effects of the perturbation via the presence of specific centres of action. We finally discuss possible applications for the study of terrestrial phenomena with similar characteristics.

Decadal-scale rainfall variability in Ethiopia recorded in an annually laminated, Holocene-age, stalagmite

An annually laminated, uranium-series dated, Holocene stalagmite from southeast Ethiopia has been analysed for growth rate and δ13C and δ18O variations at annual to biennial resolution, in order to provide the first long duration proxy record of decadal-scale rainfall variability in this climatically sensitive region. Our study site (10°N) is climatically influenced by both summer (June—August) and spring (March—May) rainfall caused by the annual movement of the Inter-Tropical Convergence Zone (ITCZ) and modulated by large-scale anomalies in the atmospheric circulation and in ocean temperatures. Here we show that stalagmite growth, episodic throughout the last 7800 years, demonstrates decadal-scale (8—25 yr) variability in both growth rate and δ 18O. A hydrological model was employed and indicates that this decadal variability is due to variations in the relative amounts of rainfall in the two rain seasons. Our record, unique in its combination of length (a total of ~1000 years), annual chronology and high resolution δ18O, shows for the first time that such decadal-scale variability in rainfall in this region has occurred through the Holocene, which implies persistent decadal-scale variability for the large-scale atmospheric and oceanic driving factors.