Annals 47, 5, 2004 1617 ANNALS OF GEOPHYSICS, VOL. 47, N. 5, October 2004 Key words Mediterranean – precipitation – princi- pal component – canonical correlation – trajectory 1. Introduction The Mediterranean is characterised by hot, dry summers and warm, wet winters. The win- ter precipitation variability patterns of the Mediterranean region, which extends 3600 km from the Iberian Peninsula in the west to Israel in the east, show marked differences between the north-western and the south-eastern areas (e.g., Conte et al., 1989; Douguedroit, 1998). This variability can be described as a telecon- nection; an anti-correlation between the rainfall in Iberia and that in the Near East. What influences the operation of the Mediter- ranean precipitation teleconnection? Conver- gence of moisture and regional precipitation is controlled by the dominant atmospheric pressure structures. The major sources of water close to the Mediterranean region are the Atlantic Ocean, the Mediterranean Sea itself, and further east the Arabian Sea. Consequently pressure systems over the North Atlantic, Europe and North Africa must all be considered in answering the above question. Hurrell (1995) has shown that changes in the North Atlantic pressure system, indicated by extremes of the North Atlantic Oscillation (NAO), strongly alter the moisture transport across the Atlantic Ocean and onto the European landmass. The NAO is a measure of a pressure teleconnection and is defined as a large-scale al- Mediterranean precipitation and its relationship with sea level pressure patterns Roy Thompson and David N. Green School of GeoSciences, The University of Edinburgh, U.K. Abstract The relationship between Mediterranean precipitation and North Atlantic and European sea level pressure fields has been studied using statistical techniques to investigate the variability within the data. A principal component analysis shows the major winter precipitation variability is described by a see-saw fluctuation between the West- ern and Eastern Mediterranean. The pressure-precipitation relationships indicate that a highly variable, pressure region situated to the south of Britain dominates this major precipitation pattern. The large-scale pressure fields which facilitate the precipitation patterns have been isolated using a canonical correlation analysis. Although the well-known major pressure centres of action in the North Atlantic are important, pressure changes in the east are found to also control the transport of moisture across the Mediterranean to a large degree, as the presence of a large high over Kazakhstan causes meridonial flow and impedes the passage of moisture across the Mediter- ranean. The pressure-precipitation relationships are found to be very consistent over multi-decadal, seasonal, monthly and daily time-scales with trajectory analysis confirming many of the features of the average seasonal pressure charts. This steadiness and regularity indicates that the Mediterranean precipitation teleconnection is a robust phenomenon that is affected by large-scale pressure changes to both the east and west. Mailing address: Prof. Roy Thompson, School of Geo- Sciences, The University of Edinburgh, Grant Institute, The King’s Buildings, West Mains Road, Edinburgh EH9 3JW, U.K.; e-mail: roy@ed.ac.uk 1618 Roy Thompson and David N. Green ternation of atmospheric mass between regions of subtropical high pressure (the Azores high) and subpolar low pressure (the Icelandic low) in the North Atlantic (e.g., Lamb and Peppler, 1987). The NAO can be expected to contribute to explaining the dominant variations in the West- ern Mediterranean (e.g., Zorita et al., 1992; Storch et al., 1993; Rodo et al., 1997; Trigo and Palutikof, 2001), however there may well be oth- er pressure systems influencing the precipitation in the Mediterranean Basin. The spatio-temporal precipitation patterns of local regions within the Mediterranean have been described in some detail. Good sum- maries of precipitation and its relationships with circulation are available for the Central (Maheras et al., 1992; Bordi et al., 2001; Bordi and Sutera, 2002), Northwestern (Zhang et al., 1997; Ulbrich et al., 1999) and Eastern (Kuteil and Paz, 1998; Kuteil et al., 2002) Mediter- ranean. These studies point to the importance of both regional and local topographic influences. Here we include analyses of the less well-stud- ied areas of Northern Africa, so increasing the geographical coverage, and concentrate on Mediterranean wide effects. We are primarily concerned with long-term (multi-decadal) changes and variability and so only analyse sites with published precipitation data that extend as far back as at least 1905. The main objective of this paper is to study the extent and temporal be- haviour of Mediterranean precipitation patterns and their relation to sea level pressure situations over Europe and the North Atlantic. The main techniques employed are canonical correlation and principal component analysis. 2. Methodology An overview of the steps employed in our analyses is given in the flow chart of fig. 1. Monthly values of precipitation for stations around the Mediterranean region spanning the time period 1905-1990 were obtained from the Global Historical Climatology Network Data- base (Vose et al., 1992). The winter months, December through to March, were chosen for study as these are the classic months for Mediterranean precipitation and because more pronounced precipitation and pressure variabil- ity occurs in the winter season (Hurrell, 1995). Sites were selected on the basis of complete- ness of the records and on the need to achieve a wide geographical spread. Forty-four satisfac- tory sites were found. We consider this data set to be the best obtainable, providing a good bal- ance between geographical extent, complete- ness and series length. The collection of long-term precipitation data is fraught with difficulties and many records are discontinuous, due to either a break in recording or a change of recording site. Missing data can cause severe problems in cor- relation and time-series analysis. This is espe- cially true for the statistical techniques used in this project, principal component analysis and canonical correlation analysis, as both require complete data sets. 15 100 mean monthly precipitation values were available to us, however on account of the long time period under analysis there were in- evitably many missing values. To form com- Fig. 1. The major stages of the analysis. GHCN = = Global Historical Climatology Network; CRU = Cli- mate Research Unit, University of East Anglia. 1619 Mediterranean precipitation and its relationship with sea level pressure patterns plete records an imputation was executed. The most conservative approach is to use the mean of the available data. However we have fair cor- relation between the mean monthly precipita- tion at nearby sites (r 2 typically 0.55) and so a less conservative approach was preferred. We considered different geographical regions sepa- rately and then used a median polish technique (Velleman and Hoaglin, 1981) to estimate the missing data month by month. In practice for regions with poor site correlation (e.g., North Africa; r 2∼ 0.3) our interpolation method will give results similar to using the time-series mean, while for regions with good spatial co- herence (e.g., Iberia; r 2∼ 0.8) results will be similar to those obtainable by multiple regres- sion. The solid circles in fig. 2c,d are typical ex- amples of our imputation. Missing monthly means for two recent years were imputed for Lisbon, while missing data for eight war years and one recent year were imputed for Limassol. The winter monthly mean-sea-level pres- sure data were obtained on the 5° latitude by 10° longitude grid of Jones (1987). We selected data extending from 50°W to 60°E and from 20°N to 65°N covering the North Atlantic, Eu- Fig. 2a-d. a) The 1st canonical function (precipitation) index. b) The 1st precipitation principal component in- dex. c) The precipitation record at Lisbon, Portugal. d) The precipitation record (scale inverted) at Limassol, Cyprus. a b c d 1620 Roy Thompson and David N. Green rope, North Africa and the western extent of Asia. The data were initially checked for outly- ing values, which were removed. The pressure data grid was almost complete (< 1% data miss- ing) and any outliers were replaced using the imputation scheme outlined above. To allow comparison between individual station records for both sets of data a normali- sation is required to remove the bias due to some regions having larger mean values for the variables than others. As we are only concerned in the variability of precipitation and pressure, dividing each record through by the mean value of that time-series provides a simple normalisa- tion with the advantage that all values are still positive reducing problems in the subsequent analysis. It should be noted that the normalised pressure and precipitation variables still have different ranges. In particular the coefficient of variation (standard deviation of time-series/ arithmetic mean of time-series) is approximate- ly two orders of magnitude larger for the pre- cipitation data compared to that of the pressure data. This causes no mathematical problems in the subsequent analyses except that certain val- ues, as noted in the text below, are not directly comparable. The two data sets are moderately large and therefore for an efficient analysis they required a reduction in size, whilst maintaining the max- imum information content possible. Principal Component Analysis (PCA) provides the ideal candidate for such data compression (e.g., Jol- liffe, 1990) and was applied individually to both the precipitation and pressure data sets. PCA worked well on our two data sets allowing us to study only the first four principal compo- nents of the pressure data and the first three principal components of the precipitation data. The dimensions of the analysis were consider- ably reduced and yet 83.5% of the pressure variance and 51.2% of the precipitation vari- ance was retained. To gain insight into the relationship between the two sets of data, pressure and precipitation, a Canonical Correlation Analysis (CCA) was performed. CCA quantifies the strength of the relationship between two data sets by maximis- ing the correlation between linear combinations of the variables (e.g.,, Hair et al., 1995). The problem is to find the linear combination of pre- cipitation variables that has maximum correla- tion with a linear combination of pressure vari- ables. To simplify proceedings it is beneficial to use the two data sets in their principal compo- nent co-ordinates. From the CCA a time-series can be constructed which emphasises the years at which the pressure and precipitation patterns are highly or inversely correlated. We also found it helpful to consider only the six most extreme years, for both the high and low canonical func- tion values. For these extreme years mean pres- sure maps were constructed. These maps can be easily interpreted, alongside meteorological records, to show which large scale features of the pressure field give rise to the Mediterranean precipitation teleconnection. Storch and Zwiers (1999, p. 382) provide us with a formal defini- tion of a teleconnection pattern. The mean pres- sure maps can also be compared with maps of individual events (see Section 3.3). It is important to emphasise that a statistical analysis based solely on principal component and canonical correlation analysis basically on- ly serves to show that patterns of the analysed fields covary. By itself it does not demonstrate the patterns of major statistical weight are dy- namically significant (being, for example, linear or non-linear normal modes); not that there is a cause and effect connection between the co- varying variables; or that the dominant patterns of co-variation (as identified, by the canonical correlation analysis) are dynamically signifi- cant. So in addition to analysing mean monthly data, daily and sub-daily records were studied for specific synoptic conditions. The daily records chosen for comparison with the mean monthly charts were selected as objectively as possible. The year and month of an event were selected from the exceptionally high and low values of the first CCA timeseries. The day of the event was chosen as that with the highest rainfall at either Lisbon (Portugal) or Paphos (Cyprus), depending on which polarity of the precipitation teleconnection was being consid- ered. Lisbon and Paphos were chosen to repre- sent typical West and East Mediterranean sites respectively. Surface pressure maps were constructed to give a synoptic overview of the conditions over 1621 Mediterranean precipitation and its relationship with sea level pressure patterns Europe, at the time of the precipitation event, and back and forward trajectories were produced us- ing the British Atmospheric Data Centre trajec- tory service to give an indication of the direction and speed of airflow towards and away from the two sites. For this analysis an initial pressure lev- el of 700 mb was chosen over the two recording stations as a compromise corresponding very ap- proximately to an average height of precipitation formation while being removed from surface boundary layer levels. The height along the tra- jectory path was calculated from the pressure level data using the hypsometric equation. 3. Results 3.1. The Mediterranean precipitation teleconnection A spatial view of the Mediterranean precip- itation teleconnection is provided by the 1st PCA loadings map (fig. 3) giving an indication of the geographical extent of the first-order pre- cipitation variability. The zero loading line trends south-west to north-east, with the load- ing size increasing with distance from the zero contour. The time series of the first principal component (fig. 2b) indicates how the variabil- ity in the precipitation across the Mediter- ranean, shown in fig. 3, has changed over the 85 year timespan since 1905. The PCA index is centred on a value of one, therefore high index years such as 1936 show unusually high precip- itation in the west and low rainfall in the east, whereas low index years such as 1989 show the reverse situation. The index shows a slight ten- dency to rising values prior to the 1960s/1970s, and falling values thereafter. 3.2. CCA and the synoptic pressure situations The canonical correlation analysis took the precipitation PCA discussed above and a PCA of the corresponding pressure field and calcu- Fig. 3. The geographical form of the 1st precipitation principal component. The contours represent equal com- ponent loading values and indicate that the regions of maximum variability occur towards the west and south- east extents of the Mediterranean. The smooth contours were drawn by fitting a lowess smoothing function (Cleveland, 1979) to the loadings of the first principal component. Black dots represent the forty-four precipita- tion recording stations used in the analysis. 1622 Roy Thompson and David N. Green lated the strength of the relationship between the two sets. The resulting time series (fig. 2a) is centred on a value of 0.24, and shows the variability of the first canonical function. This CCA time series is extremely well correlated with the precipitation first principal component (r = 0.86), giving an indication that the precipi- tation teleconnection is highly influenced by the pressure situation. The geographical patterns of the CCA show- ing the positions of maximum response between the variability in the pressure and precipitation fields over the 85 year study are shown in fig. 4a,b. The first precipitation canonical function is remarkably similar to the first precipitation prin- cipal component indicating that the predomi- nant Mediterranean precipitation variability is well correlated with the pressure changes over Fig. 4a,b. a) The geographical form of the first precipitation canonical function. The contours indicate equal loading values. Increasing loading values represent increased correlation between the variations of precipitation at that locality and the pressure field. b) The geographical extent of the first pressure canonical function. The contours indicate equal loading values. Increasing loading values represent increased correlation between the variations of pressure at that locality and the whole precipitation field. Thus regions of high loading are partic- ularly responsive. N.B.: the loading values in a) and b) are not directly comparable due to the choice of normal- isation (see text for further details). a b 1623 Mediterranean precipitation and its relationship with sea level pressure patterns the area chosen. The first pressure canonical function map indicates in which regions the pressure variability contributes most to the pre- cipitation changes. The form of the canonical function is at first sight unexpected. Rather than two centres of action over the East and West Mediterranean or a single maximum lying in be- tween the two foci of precipitation variability, major regions of pressure variation well to the north are found. One elongated focus, centred at an approximate position of 80°N and 10°W, ex- tends from South Greenland across the North Pole to Northern Russia. It is similar in shape to the area of high correlation between Portuguese rainfall and North Atlantic sea level pressure as given by Ulbrich et al. (1999), but with a promontory extending towards Siberia. The negative polarity of loadings indicates that pres- sure variations in a high index year (wet in the west/dry in the east) will tend to be unusually high pressure over Greenland and over Northern Russia. The second focus is situated over the ap- proaches to the English Channel (fig. 4a,b) and is of opposite (i.e. positive) polarity indicating that the pressure variations in this area are anti- correlated with those to the north. In high index years higher pressure than usual over the Eng- lish Channel guides maritime air and precipitat- ing weather systems north westwards away from Iberia. The form the canonical function takes can be clarified by examining the six most extreme years for both polarities of the index. Figure 5a,b plots the mean winter sea-level pressure maps for these two sets of extreme conditions. Figure 5a presents the situation causing high precipitation in the Western Mediterranean and low precipitation in the east. The Azores high and Icelandic low pressure systems are dis- placed south of their mean positions causing a series of isobars to run west-east directly across the Iberian Peninsula, providing a route for moisture from the Atlantic Ocean to impinge upon Spain and Portugal. Towards the east of the Mediterranean, a large high-pressure ridge (a western bulge of the Siberian high-pressure Fig. 5a. The mean winter pressure map associated with the six highest canonical function years (precipitation in the west, dry in the east). All contours labelled in millibars. 1624 Roy Thompson and David N. Green system), referred to here as the Kazakhstan high, provides a block causing meridional air- flow. The col on the ridge between the Azores and Kazakhstan highs lies to the south of the Mediterranean over Libya. Figure 5b indicates the situation leading to increased precipitation in the east and dry con- ditions in the west. The Azores high pressure system is displaced to the north and extends over the Iberian Peninsula preventing airflow across the Atlantic and onto the western Mediterranean where below normal precipita- tion is recorded. Consequently tight isobars trend north eastwards from the Atlantic across Northern Scotland and then onwards across Scandinavia and Finland. Situated in the east- ern Mediterranean is a distinct pressure mini- ma, the Cyprus low. This region has been linked to Mediterranean cyclones (Alpert et al., 1990), the associated cyclonic flow causes high pre- cipitation in Israel and the surrounding area (e.g., Zangvil and Druian, 1990). Further east the Kazakhstan high is more subdued and dis- placed somewhat to the south. The ridge be- tween the Azores and Kazakhstan highs now lies to the north of the Mediterranean with a broad col over Hungary and Romania. 3.3. Individual events In addition to examining the 6-year mean winter circulation of the extremes of our Mediterranean precipitation teleconnection, we have also analysed individual precipitation events to see how the sub-daily dynamic situa- tion compares to the long-term average. A high rainfall event recorded in Lisbon on 27th January 1979 was chosen as an example of behaviour in a year of high CCA index. January 1979 was wet across Iberia and dry, or moder- ately dry, in the Eastern Mediterranean (com- pletely dry at Cairo). The synoptic pressure map, trajectory path and trajectory height for this one event can be seen in fig. 6a. The pres- sure map for 26th January 1979 (the 24-h peri- Fig. 5b. The mean winter pressure map associated with the six lowest canonical function years (precipitation in the east, dry in the west). All contours labelled in millibars. (i) (ii) (iii) 1625 Mediterranean precipitation and its relationship with sea level pressure patterns Fig. 6a. i) The pressure map for 26th January 1979. ii) The forward and back trajectories for the 700 mb level at Lisbon (crosses) and Paphos (open circles). The markers are placed at 6 hourly intervals. iii) The trajectory heights for Lisbon (crosses) and Paphos (open circles). Solid symbols 00:00 on 27/01/1979. 1626 Roy Thompson and David N. Green Fig. 6b. i) The pressure map for the 13th February 1990. ii) The forward and back trajectories for the 700 mb level at Lisbon (crosses) and Paphos (open circles). The markers are placed at 6 hourly intervals. iii) The trajec- tory heights for Lisbon (crosses) and Paphos (open circles). Solid symbols 00:00 on 14/02/1990. (i) (ii) (iii) 1627 Mediterranean precipitation and its relationship with sea level pressure patterns od leading up to the precipitation event) shows a well formed high in the east extending from Kazakhstan restricting the Mediterranean air- flow and generating anticyclonic circulation across the Eastern Mediterranean Sea. The Ice- landic low pressure system and the Azores high are displaced much further to the south than their mean positions so allowing depressions to pass across the Atlantic along more southerly tracks than normal onto the Iberian peninsula. The 26th sees lows passing over the Canaries and the Faeroes and a cold front approaching Lisbon. The 3-day back trajectory for Lisbon shows cyclonic circulation and rising air arriv- ing from the SW (ii and iii in fig. 6a). The tra- jectory is consistent with rainfall in Lisbon and with the air passing low over the Atlantic gath- ering moisture before rising to condense and precipitate over the land. The forward trajecto- ry from Lisbon shows the air continuing to as- cend and being carried north-west away from the Mediterranean illustrating how the moisture from the Atlantic does not reach Israel or Cyprus in such high index years. The corresponding trajectory path for Pa- phos, Cyprus for the same six-day period shows dry air circulating around the northern edge of an anticyclone before being transported off the Sahara and passing towards the north-east. The flow corresponds well with the winter trajecto- ries studied by Dayan (1986). Turning to a low canonical function year. Paphos experienced a high rainfall event on the 14th February 1990 (fig. 6b). The pressure map shows a large area of high pressure centred over the Canaries and extending over the Iberian Peninsula with depressions crossing to the north of Scotland and a low pressure system ex- tending into the eastern Mediterranean. There is a high pressure system to the east extending from the Persian Gulf to Kazakhstan, however this is much smaller than in the example above. A series of lows pass through the East Mediter- ranean around this time. One formed over North Italy on the 11th and during the 13th to 15th another passed to the north of Cyprus, an associated cold front producing the torrential rain in Paphos on the 14th. The trajectory paths show air movement consistent with this pres- sure situation. Air crossing the Atlantic is forced around the high pressure system, passing over the Bay of Biscay into the Western Mediterranean. This air has been at a low height across the ocean gathering up water, as it pass- es into the Mediterranean it rises slowly until it reaches approximately 3000 m over Cyprus, be- coming the source for the precipitation event. The air passing over the Iberian Peninsula at this time also travels across the Atlantic, how- ever it moves slowly and descends as it reaches land due to it being trapped within the Azores high pressure system. This descending air has no potential for precipitation. 4. Discussion 4.1. Pressure influences on the teleconnection The principal component analysis shows that the Mediterranean precipitation telecon- nection is the dominant form of precipitation variability across the whole region in the winter months. The canonical correlation analysis pro- vides two synoptic pressure situations that give a first order explanation of this variability. High precipitation in the Western Mediter- ranean occurs due to the moisture transport off the North Atlantic impinging upon the western seaboard of North Africa and Iberia. This is a direct consequence of the southward shift of the Icelandic low and Azores high pressure systems (fig. 5a) comparable to the motions described in the work of Sahsamanoglou (1990). At the same time, the Kazakhstan high provides a block for any westerly airflow across the Mediterranean, forcing meridional flow north- wards over Eastern Europe and Russia and al- lowing the Eastern Mediterranean to remain dry. The meridional flow causes air tempera- tures to be some 1.5°C warmer than average in the region extending from Rumania to Eastern Turkey. At the same time, in North Europe, a zone extending from Iceland across the North- ern Baltic to Northern Russia experiences air temperatures 1.5°C below average. These tem- perature anomalies further demonstrate the re- gional extent and consistency of the pressure pattern producing the Mediterranean precipita- tion teleconnection. 1628 Roy Thompson and David N. Green Conversely with an extended high pressure system over the Iberian Peninsula, as in fig. 5b, Spain and Portugal remain dry as moisture car- ried off the Atlantic is advected north west- wards onto Northern Europe. Air temperatures some 2.0°C above normal are produced over much of Sweden and North Finland by the flow from the south-west. In contrast air tempera- tures over Eastern Turkey are some 2.0°C cool- er than average. The isobars surrounding the extended Azores high exhibit some curvature across Northern Europe leading to the possibil- ity of moisture entering the Eastern Mediter- ranean either by transportation from the At- lantic or from secondary depressions formed within Mediterranean areas in the lee of the Alps (e.g., Buzzi and Tibaldi, 1978). The tra- jectory of fig. 6bii indicates that the former oc- curred on the 13th of February 1979. The Cyprus low in fig. 5b indicates that depressions and cyclonic flow are present in the East which remove considerable heat and moisture from the water (Shay-El and Alpert, 1991). This is il- lustrated by the rising airmass in the trajectory of fig. 6b as it approaches the Cyprus region. The areas of pressure variability which af- fect the precipitation teleconnection most are highlighted by the large loading regions in the first canonical correlation map (fig. 4b). The central area of high variability correlation, cen- tred over the English Channel approaches, is linked to the translation of the Azores high pressure system. The English Channel centre of action is the position that is affected most by this movement. When the Azores high is promi- nent and moves northward from its mean posi- tion (approximately 30°N and 35°W: Sah- samanoglou, 1990) this region underlies high pressure (e.g., fig. 5b) but when the Azores high and the Icelandic low move southwards the re- gion is overlain by a low-pressure system. As explained above this alters the moisture trans- port across the Iberian Peninsula. As this vari- ability affects the movement of moisture into and across the Western Mediterranean, so the pressure variability in Northern Russia affects the eastern end of the teleconnection. As the pressure increases in Northern Russia the Kaza- khstan high blocks any East Mediterranean storm tracks and airflow across the Mediter- ranean (fig. 5a). This coincides with the south- ward migration of the Azores high causing low pressure to be present in the English Channel. This difference in pressure between the western and eastern centres of action is illustrated by the difference in polarity between the canonical function loadings of the two centres of action in fig. 4b. A year in which there is unusually high pressure in Northern Russia corresponds to a year of unusually low pressure in the English Channel approaches. The opposite situation oc- curs when the Kazakhstan high wanes and the Azores high strengthens and pushes northward over the English Channel region. Calculation of the canonical function pair of fig. 4a,b involves a long series of matrix manip- ulations including scaling, extraction of princi- pal components and estimation of canonical loadings. The whole sequence is then followed in reverse to generate the individual canonical components. A very instructive exercise, which can be used to confirm that the above matrix al- gebra has been performed correctly, is to sub- tract the pressure field of fig. 5a (the six highest index years) from that of fig. 5b (the pressure field of the six lowest index years). Here the re- sult, as expected, is an almost identical pattern to that of fig. 4b. Recently Palutikof et al. (1996) and Qua- drelli et al. (2001) have investigated the rela- tionship between Mediterranean precipitation and 500 hPa geopotential height. Although Quadrelli et al.’s (2001) results only spanned 17 years, and so they were in doubt about the long-term robustness of their results, our 85- year sea level pressure/precipitation relation- ships are in excellent agreement with their link- ages to large-scale upper-air circulation anom- alies. By contrast, Palutikof et al. (1996) em- phasise the 500 hPa surface over the Mediter- ranean Sea itself as an indicator of atmospheric pressure which relates to Mediterranean precip- itation variability. 4.2. The second canonical pair The second canonical pair highlights the major correlation between the pressure and rainfall data sets under the constraint that this 1629 Mediterranean precipitation and its relationship with sea level pressure patterns correlation is orthogonal to the first canonical pair. The second precipitation canonical func- tion has high loadings over the North African sites and loadings close to zero over Portugal to the west and Israel to the east. The area of high loadings is centred over the region in which there is most difference between the first prin- cipal component (fig. 3) and the first precipita- tion canonical function (fig. 4a). This can be explained due to the fact that the second canon- ical function, which perturbs the consistency between the canonical function and the princi- pal component is describing a second order ef- fect. The first canonical function has a r 2 value of 0.79 while the second has a r 2 value of only 0.22. The second pressure canonical function shows that the area of large pressure changes which is predominantly correlated with the sec- ond-order precipitation variation in North Africa is situated well to the north of the Mediter- ranean, with a centre over Northern Scandi- navia. The zero contour of the map of this sec- ond pressure canonical function passes along the whole length of the Mediterranean Sea. So again, we find it is the broad regional or far field, pressure variations that principally orches- trate Mediterranean precipitation variability. 4.3. The link between monthly and daily analyses The pressure charts in fig. 5a,b are the aver- age of six mean winter (DJFM) pressure charts chosen due to their high or low canonical index value. Therefore, they are an extremely smoothed version of the actual pressure situation at any particular time within the winter period. Nevertheless the main pressure patterns are ro- bust, representative features. For example the low-pressure centre that lies between Crete and Cyprus, in fig. 5b, occurs as a closed contour feature in three out of the four individual winter months (December through March). A distinc- tive, closed contour, Eastern Mediterranean low is also present in five of the six individual win- ters that were averaged to make up fig. 5b. Fur- thermore the mean pressure charts of fig. 5a,b can be compared with the pressure maps and tra- jectory analyses of fig. 6a,b that shows plots rep- resenting actual events. These type events show broadly similar features to those of the monthly analysis. The most important dynamical features con- sistent with high precipitation in the west and low precipitation in the east are the high over Kazakhstan linked with low pressure over the Atlantic, both of which are present in the daily chart (fig. 6a). This situation causes a strong east-west pressure gradient and air to flow off the ocean onto the Iberian Peninsula, before it passes north as the block creates meridonal flow. The equivalent air mass passing over Cyprus at the same time has its origins in the dry regions of the Northern Sahara. The back trajectory stays at an approximately constant height allowing no advection of moisture. Dayan (1986) classified this type of flow as of ‘inland North African type’. Consistency between the seasonal mean and daily charts is also apparent for the opposite sit- uation in which the Eastern Mediterranean is subject to higher rainfall. The high pressure over Portugal in fig. 6b causes the Atlantic moisture to travel north and then turn around the pressure maximum and back into the Mediterranean. The low pressure system then extends into the eastern region generating cy- clonic flow and transportation of moisture over Israel, Cyprus and coastal regions of Egypt. Thus as discussed by Zangvil and Druian (1990) rainfall over Israel is largely affected by the trajectory of air prior to the rain event. Mar- itime air trajectories are linked with rain events. However, our canonical pressure patterns would suggest that flow from the Atlantic is al- so of importance in addition to Zangvil and Druian’s (1990) postulated source as moisture evaporation from the Mediterranean. At these times there is also airflow across the Atlantic onto Lisbon, however it is travelling extremely slowly within the high pressure system and is sinking as it passes over the land (fig. 6biii), thereby allowing no precipitation. Once again the fact that the same features can be identified in both the seasonal average and daily synoptic charts suggests that a stable recurrent phenom- enon orchestrates the Mediterranean precipita- tion teleconnection. 1630 Roy Thompson and David N. Green 4.4. Mediterranean precipitation and pressure teleconnections It has been well documented by previous workers that the NAO causes precipitation anomalies in the Western Mediterranean area (e.g., Lamb and Peppler, 1987; Hurrell, 1995). However, the present study indicates that pres- sure variability in the east, far removed from the North Atlantic, has a large influence on the pas- sage of moisture through the entire Mediter- ranean region. The relative importance of the Kazakhstan high and Cyprus low have possibly been previously underestimated in terms of their effect on Mediterranean precipitation. The moisture which forms the precipitation comes from the west as seen in all the trajectories (fig. 6a,b), however it is the presence of a block or weakened anticyclone in the east which deter- mines the extent of any meridonial flow, and therefore the extent of any Eastern Mediter- ranean precipitation. The influence of this east- ern precipitation can be seen clearly on the time series (fig. 2a-d). Comparing the CCA index (fig. 2a) to the two individual precipitation records (fig. 2c and 2d) it can be seen that dif- ferent temporal sections have distinct features from one or either of these rainfall time-series. For example the variation in the PCA index (fig. 2b) around 1945 closely follows that of the Lisbon precipitation record whereas the varia- tion around 1970 has more in common with the Paphos time series, these variations can then be followed into the CCA time-series. Rodgers (1997) suggests that the NAO has little influence on the North Atlantic storm tracks and that the major changes in these tracks are consistent with sea level pressure changes in the Bay of Biscay region. This part of the At- lantic is identifiable in fig. 4b as having a large influence on the passage of moisture from the Atlantic into the Mediterranean, as seen by our ‘English Channel’ centre of action centred slightly northwards of Rodgers’ in the Bay of Biscay. When there is low pressure over the Bay of Biscay (fig. 5a), the Azores high is pushed to the south allowing moisture to pass over the At- lantic and onto the Iberian Peninsula causing high precipitation. High pressure over the Bay of Biscay extends across Spain and Portugal re- ducing the precipitation in this area. The storm tracks from the Atlantic in this situation pass to the north of the Bay of Biscay and then turn south into the Mediterranean, as seen by the tra- jectories in fig. 6b. The influence of pressure changes over the Bay of Biscay is described in terms of a series of North Atlantic teleconnec- tions summarised by Rodgers (1990). Recently Kutiel and Benaroch (2002) have identified and described the upper level atmospheric telecon- nection pattern (NCP) at the 500 hPa level be- tween the North Sea and the Caspian areas as the most significant in the European area. The associated circulation patterns are found to have a major and direct impact on the mean surface air temperature over the Eastern Mediterranean. However the impact of the NCP on rainfall is more complex (Kutiel et al., 2002). The strength of our canonical correlation between Mediterranean precipitation variability and the large scale pressure field over an 85 year period, coupled with a very similar rela- tionship on annual, monthly and even daily time-scales, suggest it is a robust phenomenon that would be worth investigating further. In particular more work needs to be done on the influence of pressure variations in the east, in- cluding the Siberian high, before a complete picture of the link between Mediterranean pre- cipitation and sea-level pressure changes is ob- tainable. Winter precipitation persists further east through Iraq, Iran, Turkestan and on through Tashkent. However, the precipitation teleconnection patterns are rather weak through this area. To the south winter precipitation pat- terns change rapidly and there is less scope for additional studies. We were unable to use records on the desert margins (even sites such as Cairo on the Nile Valley proved unworkable) i.e. our analyses in NE Africa were largely re- stricted to the classic Mediterranean zone which only occupies a very narrow coastal strip where the Libyan and Egyptian deserts en- croach on the Mediterranean Sea (Quezel, 1981). The strong canonical pressure loadings over Greenland and Russia, both areas well to the north of the Mediterranean, point to the need of also extending the geographical cover- age of the sea-level pressure variations and pre- cipitation patterns northwards. Indeed we have 1631 Mediterranean precipitation and its relationship with sea level pressure patterns subsequently found that winter precipitation along the eastern coast of the Black Sea corre- lates inversely with precipitation in Northern Spain and Southern France in a very similar fashion to that of the East and West Mediter- ranean. More speculatively, links with other teleconnections could be explored. Interesting- ly we find a significant correlation between the first canonical function and the Indian Ocean Dipole. This link, with a pattern of sea surface temperatures within the Indian Ocean (Saji et al., 1999), further suggests that an interaction from the east plays an important role in creating the correct conditions for the Mediterranean precipitation teleconnection. The well-docu- mented Southern Oscillation shows no link with the first canonical function but does have a significant relationship with the second. This connection, harnessed with the knowledge that precipitation over North-west Africa is correlat- ed with the Southern Oscillation (Ropelwski and Halpert, 1987) and that spring precipitation in Spain shows the most significant relationship with the El Niño-Southern Oscillation in Eu- rope (Rodo et al., 1997), suggests that less pro- nounced Mediterranean precipitation variations may have links with various global pressure phenomena. 4.5. Statistical identification of pressure influences The findings of this paper are based on the identification of correlations between two large data sets. The multivariate statistical methods used provided an extremely powerful tool, al- lowing clear pictures to be obtained of the vari- ability. The consistency of results throughout are remarkable, with multi decadal (86-year), winter, monthly and daily pressure charts all showing the same broad features and highlight- ing the phenomena giving rise to the precipita- tion teleconnection. The consistency is such that even smaller local features such as the Cyprus low are identifiable even on maps of 6- year (24-month) winter averages. Consistency between independent statistical techniques is also found. The first precipitation principal component indicates the major pattern of internal variability within the rainfall data (fig. 3). This produces an extremely similar variability map to that of the first canonical function (fig. 4a), which describes the correla- tion of the variability between the precipitation and pressure data. These two types of analysis do not necessarily produce the same geograph- ical patterns as the former derives from climat- ic variability within a climatic parameter, while the later derives from the relationship between parameters. 5. Conclusions The major pattern of precipitation variabili- ty in the Mediterranean is shown to be orientat- ed such that it is dipolar with maximum fluctu- ations in the north-west and south-east. The ze- ro variability line in the first principal compo- nent strikes in a south-west/north-east direction separating the Iberian Peninsula and Northern Europe from the Eastern Mediterranean and Northern Africa. Canonical correlation analysis shows that the European surface pressure field and Mediter- ranean precipitation field variations are highly correlated with a north-west/south-east precipita- tion see-saw and an anomalous, highly variable, pressure region situated to the south of Britain. Pressure maps for the high index years of the canonical function show that a high pres- sure system over West Kazakhstan causes meridional flow, blocking precipitation from the Eastern Mediterranean. A simultaneous translation of the Azores high to the south al- lows westerly flow from the Atlantic to pass over the Iberian Peninsula increasing precipita- tion. Pressure maps for the low canonical func- tion show the Azores high translated north- wards over the Iberian Peninsula, thereby re- ducing precipitation in the Western Mediter- ranean. These average pressure maps also indi- cate the presence of a depression, the Cyprus low, situated in the Eastern Mediterranean. Its location is in good agreement with the findings 1632 Roy Thompson and David N. Green of Zangvil and Druian (1990) who linked the phenomenon with high precipitation in Israel. The winter pressure and precipitation inter- action in the Mediterranean is correlated with the pressure variations across the North Atlantic i.e. the North Atlantic Oscillation. A further link with the Indian Ocean as described by the Indian Dipole Mode is also found to exist. Our results show excellent consistency be- tween all four individual winter months and the total season (DJFM) over the 86-year period studied, as well as for single precipitation events. The invariability over such a wide range of temporal scales lends confidence to the ro- bustness of the Mediterranean precipitation teleconnection and to its physical link with re- gional pressure patterns. Acknowledgements The data were obtained from databases kindly provided by P. Jones of the Climate Re- search Unit at the University of East Anglia, and R. Vose and colleagues of the Global His- torical Climatology Network. The British At- mospheric Data Centre is thanked for its re- analysis database and trajectory service. The authors are grateful to Keith Weston who pro- vided much meteorological expertise and ad- vice concerning the synoptic analyses, especial- ly the trajectories. REFERENCES ALPERT, P., B.U. NEEMAN and Y. SHAY-EL (1990): Climato- logical analysis of Mediterranean cyclones using ECMWF data, Tellus, 42A, 65-77 BORDI, I. and A. SUTERA (2002): An analysis of drought in Italy in the last fifty years, Il Nuovo Cimento, 25, 185- 206. BORDI, I., S. FRIGIO, P. PARENTI, A. SPERANZA and A. SUTERA (2001): The analysis of the standardized pre- cipitation index in the Mediterranean area: large scale patterns, Ann. Geofis., 44 (5/6), 965-978. BUZZI, A. and S. TIBALDI (1978): Cyclogenesis in the lee of the Alps: a case study, Q. J. R. Meteorol. Soc., 104, 271-287. CLEVELAND, W.S. (1979): Robust locally weighted regres- sion and smoothing scatterplots, J. Am. Stat. Ass., 74, 829-836. CONTE, M., A. GIUFFRIDA and S. TEDESCO (1989): The Mediterranean Oscillation. Impact on precipitation and hydrology in Italy, in Conference on Climate, Water (Publication of the Academy of Finland, Helsinki), 121-137. DAYAN, U. (1986): Climatology of back trajectories from Is- rael based on synoptic analysis. J. Climate Appl. Mete- orol., 25, 591-595. DOUGUEDROIT, A. (1998): Que peut-on dire d’une Oscillation Méditerranéenne?, in «Climate and Environmental Change», Proceedings of the Meeting of the Commission on Climatology of the IGU, Evora, Portugal, 135-136. HAIR, J.F., R.E. ANDERSON, R.L. TATHAM and W.C. BLACK (1995): Multivariate Data Analysis with Readings (Prentice Hall, New Jersey), 4th edition, pp. 745. HURRELL, J.W. (1995): Decadal trends in the North Atlantic Oscillation: regional temperatures and precipitation, Science, 269, 676-679. JOLLIFFE, I.T. (1990): Principal Component Analysis: a be- ginner’s guide, 1. Introduction and application, Weath- er, 45 (10), 375-383. JONES, P.D. (1987): The early twentieth century Arctic High – fact or fiction?, Climate Dynamics, 1, 63-75. KUTIEL, H. and Y. BENAROCH (2002): North Sea-Caspian Pattern (NCP) – An upper level atmospheric telecon- nection affecting the Eastern Mediterranean: identifi- cation and definition, Theor. Appl. climatol., 71 (1/2), 17-28. KUTIEL, H. and S. PAZ (1998): Sea level pressure departures in the Mediterranean and their relationship with monthly rainfall conditions in Israel, Theor. Appl. Cli- matol., 60 (1/4), 93-109. KUTIEL, H., P. MAHERAS, M. TURKES and S. PAZ (2002): North Sea Caspian Pattern (NCP) – An upper level at- mospheric teleconnection affecting the Eastern Mediterranean: implications on the regional climate, Theor. Appl. Climatol., 72 (3/4), 173-192. LAMB, P.J. and R.A. PEPPLER (1987): North Atlantic Oscil- lation: concept and an application, Bull. Am. Meteorol. Soc., 68 (10), 1218-1225. MAHERAS, P., C. BALAFOUTIS and M. VAFIADIS (1992): Pre- cipitation in the Central Mediterranean during the last century, Theor. Appl. Climatol., 45, 209-216. PALUTIKOF, J.P., M. CONTE, J. CASIMIRO MENDES, C.M. GOODESS and F. ESPIRITO SANTO (1996): Climate and climatic change, in Mediterranean Desertification and Land Use, edited by C.J. BRANDT and J.B. THORNES (Wiley, London), 43-86. QUADRELLI, R., V. PAVAN and F. MOLTENI (2001): Winter- time variability of Mediterranean precipitation and its links with large-scale circulation anomalies, Climate Dyn., 17, 457-466. QUEZEL, P. (1981): Floristic composition and phytosocio- logical structure of sclerophyllous matorral around the Mediterranean, in Ecosystems of the World. Mediter- ranean Type-Shrublands, edited by F. DI CASTRI, D.W. GOODALL and R.L. SPECHT (Elsevier Scientific, New York), vol. 11, 107-121. RODGERS, J.C. (1990): Patterns of low-frequency monthly sea level pressure variability (1899-1986) and associat- ed wave cyclone frequencies, J. Climate, 3, 1364-1379. RODGERS, J.C. (1997): North Atlantic storm track variabili- ty and its association to the North Atlantic Oscillation 1633 Mediterranean precipitation and its relationship with sea level pressure patterns and climate variability of Northern Europe, J. Climate, 10, 1635-1647. RODO, X., E. BAERT and F.A. COMIN (1997): Variations in seasonal rainfall in Southern Europe during the present century: relationships with the North Atlantic Oscilla- tion and the El Niño-Southern Oscillation, Climate Dyn., 13, 275-284. ROPELWSKY, C.F. and M.S. HALPERT (1987): Global and re- gional scale precipitation patterns associated with the El Niño/Southern Oscillation, Mon. Weather Rev., 115, 1606-1626. SAHSAMANOGLOU, H.S. (1990): A contribution to the study of action centres in the North Atlantic, Int. J. Climatol., 10, 247-261. SAJI, N.H., B.N. GOSWAMI, P.N. VINAYACHANDRAN and T. YAMAGATA (1999): A dipole mode in the tropical Indi- an Ocean, Nature, 401, 360-363. SHAY-EL, Y. and P. ALPERT (1991) A diagnostic study of winter diabatic heating in the Mediterranean in relation to cyclones, Q. J. R. Meteorol. Soc., 117, 715-747. STORCH, H. VON and F.W. ZWIERS (1999): Statistical Analy- sis in Climate Research (Cambridge University Press, Cambridge), pp. 494. STORCH, H. VON , E. ZORITA and U. CUBASCH (1993): Down- scaling of global climate change estimates to regional scales, an application to Iberian rainfall in wintertime, J. Climate, 6, 1161-1171. TRIGO, R.M. and J.P. PALUTIKOF (2001): Precipitation sce- narios over Iberia: a comparison between direct GCM output and different downscaling techniques, J. Cli- mate, 14, 4422-4446. ULBRICH, U., M. CHRISTOPH, J.G. PINTO and J. CORTE-REAL (1999): Dependence of winter precipitation over Portu- gal on NAO and baroclinic wave activity, Int. J. Cli- matol., 19, 379-390 VELLEMAN, P.F. and D.C. HOAGLIN (1981): Applications, Basics, and Computing of Exploratory Data Analysis (Duxbury, Boston), pp. 354. VOSE, R.S., R.L. SCHMOYER, P.M. STEURER, T.C. PETERSON, R. HEIM, T.R. KARL and J.K. EISCHEID (1992): The Global historical climatologynetwork: long-term monthly tem- perature, precipitation, sea level pressure, and station pressure data, NDP-041, Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, Oak Ridge, Tennessee (on line: http://cdiac.esd.ornl.gov/ghcn/ ghcn.html). ZANGVIL, A. and P. DRUIAN (1990): Upper air trough axis orientation and the spatial distribution of rainfall over Israel. Int. J. Climatol., 10, 57-62 ZHANG, X., X. WANG and J. CORTE-REAL (1997): On the re- lationships between daily circulation patterns and pre- cipitation in Portugal, J. Geophys. Res., 102, 13495- 13507. ZORITA, E., V. KHARIN and H. VON STORCH (1992): The at- mospheric circulation and sea surface temperature in the North Atlantic area in winter: their interaction and relevance for Iberian precipitation, J. Climatol., 5, 1097-1108. (received December 15, 2003; accepted March 4, 2004)