Influence of Sea Level Pressure on Inter-Annual Rainfall Variability in Northern Senegal in the Context of Climate Change ()
1. Introduction
The Sahel is an area of the semi-arid expanse of meadows, shrubs and small thorny trees found in the Southern Sahara Desert [1]. This region stretches to about 5000 km through Africa from east to west, more precisely, the term “Sahel” applies to a small region between latitudes 14˚N and 18˚N. It consists of countries like Sudan, Chad, Mali, Senegal, Niger, Mauritania, Burkina Faso (the northern part), and Nigeria (the northern part). Most of the local population of this region live in rural areas and practice agriculture. Thus, the variability of climatic parameters is of great importance for the region. The specific study area, i.e. Senegal (Figure 1), has two anticyclones: the Azores and Saint Helena. Taking monthly time-step, with no connection with ocean temperatures, [2] initiated a series of studies by showing the Azores anticyclone’s inhibitory influence on the Sahelian precipitation. The enhanced effect of their coupling, a strong Azores anticyclone and a weak St. Helena anticyclone leads to an enhanced rainfall deficit over the Sahel and vice versa. [3] showed that during the wet regime, there is a substantial variance in the seasonal modes linked to the Inter-Tropical Convergence Zone (ITCZ) and multi-annual modes succeeding
Figure 1. Topography (m) of West Africa, the Atlantic Ocean and some African continents are shown. Box indicates the study area of the northern of Senegal region.
each other over time (5 - 8 years then 2 - 4 years) in association with temperatures of the North Tropical Atlantic (TNA) then with the Anticyclone of the Azores. [4] showed that about 50% of the variance in July-September rainfall over the Sahel is explained by changes in the high longitude position of the Azores and low pressure in South Asia.
The hot phase of Enoa is manifested by a weakening of the gradient of pressure anomalies between the east and west of the Pacific and a decrease in the east-west flux of the trade winds, phenomena associated with persistent warming of the eastern and central equatorial Pacific. The opposite phenomena manifest the cold phase of the Enoa: increase in the zonal pressure gradient and the flux of trade winds associated with persistent cooling in the central equatorial Pacific. This phenomenon has a considerable impact on temperature and precipitation fields [5] and is particularly noticed around the Pacific and the Indian Ocean. However, the peri-Atlantic regions are also concerned since Enoa is also associated with the variability of sea surface temperatures and trade winds over the tropical Atlantic ( [6] [7] [8] ).
Since temperatures and sea surface pressure are linked, it is necessary to investigate the interaction between precipitation and MSLP to assess their modification in the West African Sahel’s climate system, especially over northern Senegal. [9] showed that the existence of cooler (warmer) TSOs in the Atlantic Ocean is overcome by an increase (decrease) in the pressure field. The choice of SLP comes from the analysis of [10] who showed that pressure is the most sensitive to the influence of Sea Surface Temperature (SST).
In this paper, datasets and methods used are briefly described in Section 2. Section 3 investigates the characteristic of precipitation and the annual cycle associated with MSLP. Finally, Section 4 summarizes the main results and gives some prospects for future work.
2. Data and Methods
2.1. Data
Precipitation data used was derived from the Global Precipitation Climatology Project (GPCP) version 2.2. The GPCP analyses are derived from rain gauge stations, satellites, and sounding observations that have been merged. The data is available from 1979 to the present at 2.5˚ resolution [11]. This dataset consists of a combination of all precipitation data from polar satellites (estimates of The Special Sensor Microwave Imager (SSMI/I) emissions) and the precipitation data centre of geostationary satellites of the GPCP (estimated rain gauge GPI and OPI). This data product has been validated regionally and globally by [12] and [13]. The ERA-Interim (ERA-I) reanalysis of European Centre for Medium-Range Weather Forecasts (ECMWF, [14] ), available from 1979 to present with a spatial resolution of 0.75˚ × 0.75˚ were used. The variables are specific humidity, zonal and meridional wind component. We also use the monthly SLP from the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) Reanalysis Project available for the period 1979 to the present at 2.5˚ resolution.
2.2. Methodology
Wavelet analysis is a standard tool for decomposing a time series into a time-frequency space and detecting time-frequency variations. The wavelet transformation makes it possible to compare a signal with a wavelet function called mother wavelet [15]. After having tested different mother wavelets, the results obtained are similar. As the wavelet transform is a band-pass filter with a known response function (the wavelet function), it is also a powerful filtering technique. To identify the dominant modes, the wavelet analysis method is used during the period 1979-2003 in the study area (Figure 1).
The Orthogonal Empirical Function (EOF), also called Principal Component Analysis (PCA) is a multivariate statistical technique that consists in reducing a dataset containing a large number (K × 1) of data vectors to a dataset containing fewer new (M × 1) variable vectors u, which are linear combinations of the original vectors and which contain most of the information from the original collection of x. The elements of these new vectors u are called the principal components (PC). Most often, the main component is calculated using anomalies. The first PC, u1, is the linear combination of x' with the most significant variance. The following main component um, m = 2, 3, 4, ..., are the linear combinations having the most significant possible variance, provided that they are not correlated with the main components having lower indices.
The mth principal component of the elements um of u is obtained by the projection of the data vector x' on the mth eigenvector, em, of the covariance matrix of x, [S], by the relationship (Wilks, 2011):
(1)
where the transposition operation is denoted by the exponent T, and each of the M eigenvectors contains an element relating to each of the K variables, xm.
Other methods involve analyzes of the main characteristics of spatial variability (generally given by the most important empirical orthogonal function models (EOF)), and large-scale mechanisms controlling the regional climatic variability given by the canonical correlation ( [16] [17] [18] [19] ).
3. Results and Discussion
Figure 1 shows the topography (m) of West Africa, the Atlantic Ocean and some African continents. The rectangle indicates the study area of the northern region of Senegal. From Figure 1, the study area is between 15˚N - 17˚N and 17˚W - 13˚W. The study area is bordered to the west by the Atlantic Ocean; this justifies the study of Sea Level Pressure (SLP). To understand the different levels of sea level pressure used in the literature Figure 2 is shown.
Figure 2 shows the average global surface pressure. The study area, which is located in the north part of West Africa, is surrounded by two anticyclones, between 15˚N - 60˚N and 90˚W - 05˚W (1010 hPa) and 40˚S - 10˚S and 30˚W - 0˚E (1020 hPa). To explain the precipitation trend changes, [20] consider the response precipitation to the NAO, the regression between sea level pressure (SLP) over the Euro-Atlantic area onto precipitation. Figure 1 allows us to show that our study area, northern Senegal, is surrounded by two high sea level pressures, namely the Azores High and the St. Helena High. Since our study focuses on the influence of these two high-pressure areas, this figure allows us to locate the coordinates of these high-pressure areas in relation to our zone.
As the study is based on the influence of sea-level pressure on interannual variability, is interesting to see how surface pressure varies in different months.
Figure 3 shows the monthly average of Sea Level Pressure. It is observed that SLP, varying between 1006 - 1016 hPa over 20˚S and 20˚N and 1018 - 1022 hPa
Figure 2. Global average of Sea Level pressure for the period 1979-2003. between 30˚S - 30˚N and 40˚W - 30˚E. The pressures are in millibars.
Figure 3. Monthly average (shaded, hPa) of the Sea Level Pressure for the 1979-2003 period between 30˚S - 30˚N and 180˚W - 60˚E. The pressures are in millibars.
beyond 20˚N and below 20˚S. The southerly high pressure of 1022 hPa or greater approached 20˚S and sea level longitude (located between the continents of Africa and America) from June to October. This period corresponds to the rain belt defined by [21]. While the northern one approaches the 20˚N zone and the sea longitude (located between the two continents Africa and America) from June to July. The high-pressure system in the south becomes lower than 1022 hPa from August to October. This figure of the mean monthly pressure indicates the relationship of the Sea Level Pressure with the seasonal cycle of precipitation in our study area. [22] showed that the pressure adjustment mechanism is the main driver of the convergence of meridional surface winds in the eastern tropical Atlantic.
The deficit of rainfall was associated with an increase in the subsidence of the air and a surplus with a more intense monsoon circulation. It is necessary to look at the wind characteristics at 700 hPa and 200 hPa. The wind at these two pressure levels plays an important role in the interannual variability of rainfall. They correspond to the position of the African Easterly Jet (AEJ) and the Tropical Easterly Jet (TEJ). In Figure 3, it can be seen that the SLP to the north and the south accompanies the rainfall belt. It is necessary to represent the Sea Level Pressure in JJAS.
Figure 4(a) shows the 700 hPa average wind (JJAS) and GPCP precipitation (JJAS) while Figure 4(b) shows the divergence of the wind (JJAS) at 700 hPa (vectors) and the Sea Level Pressure (JJAS). The mean wind (JJAS) at 200 hPa (vectors) and precipitation from GPCP (JJAS), the wind divergence (JJAS) at 200 hPa (vectors) and the Sea Level Pressure (JJAS) are shown in Figure 4(c) and Figure 4(d), respectively. Maximum precipitation accompanies the mean east wind of the northern hemisphere, which corresponds to the African Easterly Jet (JEA) (700 hPa) and the Tropical Easterly Jet (JET) (200 hPa). The high pressure located between 25˚N - 40˚N and 50˚W - 20˚W is accompanied by the wind subsidence at 700 hPa and 200 hPa. [23] have shown that the forecast rainfall deficits towards the west are associated with an increase in air subsidence and the surplus with a more intense monsoon circulation and a reinforcement of the Africa Easterly Jet and anomalies in the zonal circulation between the Indian and Atlantic oceans, favouring the subsidence of the air and the transfer of humidity outside the region. The wind convergence is observed at 200 hPa over (15˚N - 17˚N; 17˚W - 12˚W). [24] established a link between low-level divergence of the upper troposphere above West Africa with a strengthening of the Walker-type atmospheric cell above the Atlantic and the occurrence of droughts affecting the whole of West Africa.
Figure 4(e) and Figure 4(f) show the rotational wind at 700 hPa and at 200 hPa. It can be seen that the sea surface pressure (greater than or equal to 1015 hPa) is located between two zones (10˚S - 40˚S and 40˚W - 0˚E, and 15˚N - 40˚N and 50˚W - 25˚W). It is found that the rotational wind at 700 hPa is from the northeast, i.e., of negative sign. This weakness is synonymous with a minimum
Figure 4. (a) Mean wind (vector, m/s) at 700 hPa and precipitation (shaded, mm/day), (b) divergence of wind at 700 hPa (vector, m/s) and SLP (shaded, hPa), (c) mean wind at 200 hPa(vector, m/s) and GPCP precipitation (shaded, mm/day), (d) Wind divergence at 200 hPa (vector, m/s) and SLP (shaded, hPa), (e) Rotational wind (vector, m/s) at 700 hPa and SLP (hPa) and (f) rotational wind at 200 hPa and SLP (in hPa).
vortex. [25] has shown that a maximum of relative vortex marks the increase in cyclonic circulation linked to the continental thermal depression. The figure illustrates the Inter-Tropical Convergence Zone (ITCZ) defined as the zone into which the converging fluxes from the subtropical anticyclones move: The Saint Helena high (southern hemisphere), and the Azores and Saharo-Lybian (northern hemisphere) high. The ITCZ corresponds to the low-pressure area between the subtropical high-pressure belts [26]. It was found that there is probably a relationship between Sea Level Pressure and precipitation. In this study, we aim to highlight the statistical relationship between these two variables.
Figure 5 shows the correlation between SLP and JJAS precipitation between 5˚S - 20˚N and 20˚W - 0˚W. There is a negative correlation over 15˚N - 17˚N and 17˚W - 13˚W. This figure also indicates that, in the study area, SLP is related to precipitation with the opposite sign (i.e. higher values of sea surface pressure in the north are accompanied by a decrease in precipitation). [24] showed that a deficit season in the Sahel occurs when there is a “dipolar” structure of surface pressure (higher values in the North Atlantic, lower values in the South Atlantic). [27] have shown that climate change in the West African monsoon region in autumn is due to increased pressure with abnormally high-pressure circulation over Europe.
The two levels of sea pressure were identified in Figure 4 (10˚S - 40˚S and 40˚W - 0˚E, and 15˚N - 40˚N and 50˚W - 25˚W). The PCs of these SLP are
Figure 5. Correlation between sea surface pressure and JJAS precipitation 1979-2003.
shown in Figures 6-9.
Figure 6 shows the first three PCs of the SLP data between 15˚N - 40˚N and 50˚W - 25˚W during summer (JAS) from 1979 to 2003. PC1 is in a negative phase for the years 1979-1984, 1990-1994, 1996-1997, and 2002. These components account for 62.2% of the total variance. The PC2 of SLP data was generally in its positive phase for the years 1981, 1983, 1986-1987, 1991-1994, 1997 and
Figure 6. The first three PCs of SLP data between 15˚N and 40˚N and 50˚W and 25˚W during June, July and August for 1979-2003. The percentage of the total variance explaining the pattern appears in the upper right corner.
Figure 7. EOF of SLP data on between 15˚N - 40˚N and 50˚W - 25˚W during June, July and August for 1979-2003. The percentage of the total variance explaining the pattern appears in the upper right corner.
Figure 8. The first three PCs of SLP data on 40˚S - 10˚S and 40˚W - 0˚E during June, July and August for 1979-2003. The percentage of the total variance explaining the pattern appears in the upper right corner.
2002 while the PC3 of SLP data was in its negative phase between the years 1979-1981, 1984, 1987-1988, 1990, 1993, 1995-1998 and 2000 and in its positive phase between 1982-1983, 1985-1986, 1989, 1991-1992, 1994, 1999 and 2001-2003. The years 1982-1983 corresponded to years of drought in our study area. [28] showed that any strengthening or displacement towards the south of the Azores anticyclone causes a southward shift of the ITCZ. PC2 and PC3 respectively represent approximately 24.2% and 6.5% of the total variance. The first three corresponding EOFs are shown in Figure 7.
Figure 6 shows the temporal component while Figure 7 shows the spatial component. PC1 is negative over the entire area representing 62.2%. PC2 is negative between 28˚N - 40˚N and positive between 15˚N - 28˚N showing a north-south dipole. PC3 is positive between 50˚W - 25˚W and 15˚N - 24˚N. Our study area is located in this zone. According to Figure 6, PC3 was in its positive 1982-1983; this period corresponds to a period of drought in our study zone. The SLP zone between 15˚N - 40˚N and 50˚W - 25˚W allows us to see its influence on our study area located from 15˚N onwards. We have seen in Figure 2, that the North High is on 15˚N - 40˚N and 50˚W - 25˚W.
Figure 8 shows the first three PCs of Sea Level Pressure (SLP) data between 40˚S - 10˚S and 40˚W - 0˚E during summer (JAS) from 1979 to 2003. PC1 is in a positive phase for the years 1979, 1988, 1990-1992, 1994-1995, 1995 and 2000-2001. It can be noted that the positive phase 2000-2001 corresponds to the return to normal rainfall over this region. These components represent 69.4% of the total variance. The PC2 of the SLP data was generally in its positive phase for the years 1981, 1983-1986, 1988, 1991-1992, 1995-1996, 1999 and 2001-2003, while the PC3 of the SLP data was in its positive phase between the years 1980-1983, 1991-1995, 1997 and 1999-2003 and in its negative phase between 1979, 1984-1990, 1996 and 1998. PC2 and PC3 account for about 12.5% and 8.9% of the total variance, respectively. The first three corresponding EOFs are shown in Figure 9.
Figure 9 shows the temporal component while Figure 10 shows the spatial component. The PC1 is positive or null over the whole area representing 69.4%. The PC2 is negative between 26˚S - 40˚S. PC3 is positive between 30˚S - 10˚S and 40˚W - 0˚E. We will now see the wavelet of the two SLP and precipitation in our study area. Figure 10 shows the average spectrum of wavelets over the period 1979 to 2003 of the SLP data between 15˚N - 40˚N and 50˚W - 25˚W during July August and September 1979-2003. There are peaks in the 4-year band in 1996-1999; oscillations of more than eight years in 1983-1990 and indicate intervals of high variance over the entire study period.
The average wavelet spectrum over the period 1979 to 2003 of the SLP data between 40˚S - 10˚S and 40˚W - 0˚E during July August and September 1979-2003 is shown in Figure 11. A 4-year variability is observed over the period 1989-1993. The figure shows an interval of high variance between 1979-1996 and low variance between 1997-2003. [30] suggested that the Pacific and Atlantic modes are analogous, governed by physics intrinsic to the ITCZ/cold tongue complex. This explains the mode of variability of more than 4 years.
Figure 12 shows the average wavelet spectrum over the period 1979 to 2003 of the GPCP precipitation data between 15˚N - 17˚N and 17˚W - 13˚W during July August and September for 1979-2003. There are fluctuations of 1 to 2 years from 1994-1995 and peaks of 8 years from 1984-1991. Interval of high variance (1984-2003) and low variance (1980-1984) are also noted in the figure. We can
Figure 9. EOFs of SLP data between 40˚S - 10˚S and 40˚W - 0˚E during June, July and August 1979-2003. The percentage of the total variance explaining the pattern appears in the upper right corner.
Figure 10. Average wavelet spectrum over the period 1979 to 2003 from the SLP data between 15˚N - 40˚N and 50˚W - 20˚W during July-August September (JAS) for 1979-2003.
Figure 11. Average wavelet spectrum over the period 1979 to 2003 from SLP data between 40˚S - 10˚S and 20˚W - 10˚E during the months of July August and September 1979-2003.
Figure 12. Average spectrum of wavelets over the period 1979 to 2003 from the GPCP precipitation data between 15˚N - 17˚N and 17˚W - 12˚W during July, August and September 1979-2003.
conclude that the eight-year oscillations observed in 1984-1991 on the Precipitation data are in phase with the same duration’s oscillations with the SLP data of the northern anticyclones. Fontaine et al. (1998) showed that a warm phase of Enoa type tends to move Walker cell eastward, which causes an abnormal upward transport over the Pacific and an abnormal subsidence transport over the peri-Atlantic space. The variability of more than 8 years is due to the El Nino phenomenon which is in phase with the North and South high ( [29] [30] ). [31] showed that there is a significant remote influence of the Pacific ENSO on the variability of the tropical Atlantic.
4. Summary and Conclusions
The most recent research on the interannual variability of the West African monsoon (time scale 1 - 15 years) has focused on oscillation sea surface temperature or pressure. It is shown that the maximum rainfall is accompanying the mean easterly wind of the northern hemisphere corresponding to the African Easterly Jet (AEJ) (700 hPa) and the Tropical Easterly Jet (TEJ) (200 hPa); the high pressure located at 25˚N - 40˚N and 50˚W - 20˚W is escorted from the wind subsidence at 700 hPa and 200 hPa. We also found that sea level pressure (greater than or equal to 1015 hPa) is located between two zones (10˚S - 40˚S and 40˚W - 0˚E and 15˚N - 40˚N and 50˚W - 25˚W). This zone is accompanied to the north by diverging winds at 700 hPa and 200 hPa. It can be observed that in our study area (15˚N - 17˚N and 17˚W - 13˚W) the wind rotation at 700 hPa is north-easterly, i.e. with a negative sign, which is synonymous with a minimum of eddies, and that precipitation in our area is well related to SLP with the opposite sign.
We have shown the first three PC of Sea Level Pressure (SLP) data over 15˚N - 40˚N and 50˚W - 25˚W during summer (JAS) from 1979 to 2003. PC1 accounted for 62.2% of the total variance and was generally in its negative phase between 1979-1984, 90-94, 1996-1997, and 2002. PC2 of SLP were generally in its positive phase over our study area between 1981, 1983, 1986-1987, 1991-1994, 1997 and 2002; while PC3 of pressure was in its negative phase between 1979-1981, 1984, 1987-1988, 1990, 1993, 1995-1998 and 2000 and in its positive phase between 1982-1983, 1985-1986, 1989, 1991-1992, 1994, 1999 and 2001-2003. PC2 and PC3 account for about 54.2% and more than 6.5% of the total difference, respectively. These positive phases of North Sea Level pressures (strengthening or southward shift of the Azores High) correspond to years of drought in the studied area due to a southward shift of the ITCZ.
We have also shown the first three PC of Sea Level Pressure (SLP) data over 40˚S - 10˚S and 40˚W - 0˚E during summer (JAS) from 1979 to 2003. PC1s account for 69.4% of the total variance. PC1s of Sea Level Pressures were generally in its positive phase between 1979, 1988, 1990-1992, 1994-1995, 1995 and 2000-2001; while PC2s of pressures were in its positive phase 1981, 1983-1986, 1988, 1991-1992, 1995-1996, 1999 and 2001-2003. PC2 and PC3 account for about 12.5% and more than 8.9% of the total difference, respectively. These positive phases of southern Sea Level Pressures (strengthening or northward shift of the St. Helena High) correspond to years of high rainfall in the studied due to a northward shift of the ITCZ. The variability lasting one year is the seasonal variability linked to the ITCZ. Its migration depends on the pressure fields’ organisation, itself controlled by the thermal factor (sun) [32]. The role of the Atlantic appears to be the modulation of the north-south excursion of the Intertropical Convergence Zone (ITCZ). Since this zone is at the level of the Atlantic Ocean, it is, therefore, normal to see seasonal variability. This one-year variability related to the ITCZ is related to the temperatures of the North Tropical Atlantic (NTA) and then to the Azores High according to [3].
Wavelet power spectra and variance time series show interdecadal changes in 8-years. Oscillations of eight (8) years on the precipitation data are in phase with oscillations of the same duration on the SLP data from northern highs. The positive phases of North Sea level pressure (15˚N - 40˚N and 50˚W - 25˚W) (strengthening or southward displacement of the Azores High) is an unfavourable condition for rainfall in the northern zone of Senegal. This corresponds to a southward displacement of the ITCZ. On the other hand, the positive phases of Sea Level Pressures in the South (40˚S - 10˚S and 40˚W - 0˚E) (strengthening or northward shift of the St. Helena High) correspond to years of heavy rainfall in the studied area.