Diagnosing GCM errors over West Africa using relaxation experiments. Part I: summer monsoon climatology and interannual variability

The CNRM atmospheric general circulation model Arpege-Climat is relaxed towards atmospheric reanalyses outside the 10°S?C32°N 30°W?C50°E domain in order to disentangle the regional versus large-scale sources of climatological biases and interannual variability of the West African monsoon (WAM). On the one hand, the main climatological features of the monsoon, including the spatial distribution of summer precipitation, are only weakly improved by the nudging, thereby suggesting the regional origin of the Arpege-Climat biases. On the other hand, the nudging technique is relatively efficient to control the interannual variability of the WAM dynamics, though the impact on rainfall variability is less clear. Additional sensitivity experiments focusing on the strong 1994 summer monsoon suggest that the weak sensitivity of the model biases is not an artifact of the nudging design, but the evidence that regional physical processes are the main limiting factors for a realistic simulation of monsoon circulation and precipitation in the Arpege-Climat model. Sensitivity experiments to soil moisture boundary conditions are also conducted and highlight the relevance of land?Catmosphere coupling for the amplification of precipitation biases. Nevertheless, the land surface hydrology is not the main explanation for the model errors that are rather due to deficiencies in the atmospheric physics. The intraseasonal timescale and the model internal variability are discussed in a companion paper.     

Interannual variability and predictability of African easterly waves in a GCM

The interannual variability of African Easterly Waves (AEWs) is assessed with the help of spatio-temporal spectral analysis (STSA) and complex empirical orthogonal functions methods applied to the results of ten-member multiyear ensemble simulations. Two sets of experiments were conducted with the Météo-France ARPEGE-Climat GCM, one with interactive soil moisture (control), and the other with soil moisture relaxed towards climatological monthly means calculated from the control. Composites of Soudano–Sahelian AEWs were constructed and associated physical processes and dynamics were studied in the frame of the waves. It is shown that the model is able to simulate realistically some interannual variability in the AEWs, and that this dynamical aspect of the West African climate is potentially predictable (i.e. signal can be extracted from boundary conditions relatively to internal error of the GCM), especially along the moist Guinean coast. Compared with ECMWF 15-year reanalysis (ERA15), the maximum activity of AEWs is located too far to the South and is somewhat too zonal, but the main characteristics of the waves are well represented. The major impact of soil moisture relaxation in the GCM experiments is to reduce the seasonal potential predictability of AEWs over land by enhancing their internal variability.     

Solitary rossby waves over variable relief and their stability. Part II: numerical experiments

The barotropic, quasi-geostrophic vorticity equation describing large scale, rotating flows over zonal relief supports nonlinear permanent form solutions, namely nonlinear topographic Rossby waves. Through an analytical theory, these solutions have been shown to be neutrally stable to infinitesimal perturbations.Numerical algorithms, which necessarily truncate the infinite number of degrees of freedom of any continuum model to a finite number, are capable of reproducing the numerical equivalent of these form-preserving solutions. Moreover, these numerical solutions are shown to preserve their shape throughout the numerical experiment not only in the limit of small amplitude, but also for high amplitude (Rossby number → O(1)).Through numerical simulation, the stability analysis is carried far beyond the analytical limit of infinitesimal perturbations. The solutions maintain their stability in agreement with the analytical theory, up to perturbations having intensities almost of the same order as the solutions themselves. For higher-amplitude perturbations, the solutions break up and typical turbulent behavior ensues. The passage from wave-like to turbulent behavior, upon surpassing a critical perturbation value, can be observed in the sudden loss of phase locking of the permanent solution Fourier modes.     

Interactions between synoptic, intraseasonal and interannual convective variability over Southern Africa

After removing the annual cycle, a principal component analysis is applied to the daily outgoing longwave radiation anomaly field, used here as a proxy for atmospheric convection. The analysis is carried out over the southern African region (7.5°E–70°E, 10°S–40°S) for austral summer (November through February) for the period 1979–1980 to 2006–2007. The first five principal components (PC) are retained. The first two PCs describe spatial patterns oriented north-west to south-east from tropical southern Africa (SA) to the mid-latitudes. They are interpreted to be different possible locations for synoptic-scale tropical–temperate troughs (TTT), one dominant rainfall-producing synoptic system in the region. The phase relationship between these two PCs describes a tendency for these TTT to propagate eastwards from SA to the Mozambique Channel and southern Madagascar. The next three PCs describe convective fluctuations, respectively, located over the north-west, the south and the centre of SA. Their time series are significantly associated with Madden–Julian oscillation (MJO) activity in the tropics. However, we find that TTT systems are statistically independent of the MJO, i.e. they are equally liable to occur during any phase of the MJO. Three PCs out of five also show a significant association with El Niño southern oscillation, confirming that El Niño years mostly coincide with suppressed convection at the intraseasonal time-scales, a result consistent with its impact on seasonal averages diagnosed in previous studies.     

Simulation of West African monsoon using the RegCM3. Part I: Model validation and interannual variability

Summary The West African monsoon oscillates each year with remarkable regularity but the interannual variability associated with the monsoon is not fully understood although much progress has been made in recent years. This study examines and evaluates the mean state and the interannual variability of the West African climate as simulated by the International Centre for Theoretical Physics (ICTP) Regional Climate Model version 3 (RegCM3) over the period 1979 through 1990 using the National Center for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR) reanalysis data as lateral boundary conditions. Our analysis shows that the averaged rainfall over the region is well represented by the model and demonstrates considerable skill in reproducing the extreme rainfall regimes. There is however a tendency to overestimate rainfall amounts along the Guinean coast, particularly around mountainous areas, and to underestimate it over the Soudano-Sahel. The increased rainfall along the coast is due to an enhanced low-level convergence of the moist southwesterly winds along the coast leading to a reduction of the moisture content in the atmosphere. The decrease over the Soudano-Sahel could be associated with the weakening of the land–sea temperature gradient and hence the decrease in the low level southerly flows. The spatial and temporal variations in temperature are well captured by the model except for slightly cold bias over the coastal region due to an overestimation of precipitation.     

Oscillations of the Intertropical Convergence Zone and the genesis of easterly waves Part II: numerical verification

A companion paper (Part I: Toma and Webster 2008), argued that the characteristics of the mean Intertropical Convergence Zone (ITCZ) arise from instabilities associated with the strong cross-equatorial pressure gradient (CEPG) that exists in the eastern Pacific Ocean as a result of the latitudinal sea-surface temperature (SST) gradient. Furthermore, it was argued that instabilities of the mean ITCZ resulted in the in situ development of easterly waves. Thus, in Part I, it was hypothesized that the mean and transient state of eastern Pacific convection was due to local processes and less so to the advection of waves from the North Atlantic Ocean. To test this hypothesis and, at the same time, consider others such as a possible role of the equatorial and subtropical orography in generating local instabilities, a series of controlled numerical experiments are designed using the WRF regional model. The domain of the model was configured to include the western Atlantic Ocean, the Isthmus of Panama and the eastern Pacific Ocean to 155°W. Lateral boundaries were set at 40°N and 40°S, thus containing the mountains of Central America, the Andes and the Sierra Madre of Mexico. In a series of experiments, analysis products were used as boundary conditions that were successively updated four times per day, set as 10-day running average fields or as running mean monthly fields. Finally, the model was run with topography essentially eliminated over the land areas. Although there are differences between the details of the resultant fields, the location of mean convection and the form of the transients remain the same. It is concluded, in support of the theoretical and diagnostic studies of Part I that orographic forcing or waves generated in the North Atlantic Ocean are not the major causes of the mean and transient nature of disturbances in the eastern Pacific.     

Influence of aerosol-radiative forcings on the diurnal and seasonal cycles of rainfall over West Africa and Eastern Atlantic Ocean using GCM simulations

Effects of aerosol radiative forcing on the diurnal and seasonal cycles of precipitation over West Africa and eastern Atlantic Ocean are investigated for the boreal summer season: June–July–August. An eight year (2000–2007) average of GCM simulated rainfall data is compared with the corresponding TRMM rainfall data. The comparison shows that the amplitude of the diurnal cycles of rainfall over land and ocean are reasonably well simulated. Over land, the phase of the simulated diurnal cycle of precipitation peaks several hours earlier than that of the TRMM data. Corresponding differences over the ocean(s) are relatively smaller. Some of the key features of the aerosol induced model simulated field anomalies are: (a) aerosol direct radiative forcing which increases the atmospheric stability and reduces the daytime moist convection and convective precipitation; (b) the aerosol induced changes in the diurnal cycle of precipitation are out of phase with those of the TRMM data over land, but are in-phase over the ocean; (c) aerosols reduce the amplitude of the diurnal cycle of precipitation over land and enhance it over ocean. However, the phase of the diurnal cycle is not affected much by the aerosol radiative forcing both over land and ocean. During the boreal summer, aerosol radiative forcing and induced circulation and precipitation cool the Sahel and the southern part of Sahara desert more than the adjacent areas to the north and south, thereby shifting the peak meridional temperature gradient northward. Consequently, an anomalous easterly jet is found north of its climatological location. This anomalous jet is associated with increased cyclonic circulation to the south of its axis, resulting in an anomalous monsoon rain belt in the Sahel.     

North African climate variability. Part 3: Resource prediction

Summary Tropical north Africa depends on rain-fed agriculture as the main economic driver. The variability of climate-sensitive resources is investigated with a goal to develop statistical long-lead prediction models with reasonable skill. Climate data from NCEP is analysed in conjunction with agricultural and economic production in various sectors, in addition to the traditional climatic indices: temperature and rainfall. Key predictors for statistical models include the lower-level zonal wind over the Atlantic and Pacific Oceans. These exhibit a ‘memory’ that is consistent with sea surface temperatures (SST) through equatorial upwelling dynamics. Kinematic predictors outperform SST in hindcast fit by an average 33% with respect to various tropical north African resource indices. A multi-decadal oscillation induces long-term trends in rainfall that contribute to apparently skilful forecasts based on the interaction of Pacific ENSO and the Atlantic zonal overturning circulation. The skill of statistical forecasts is lower when the drying trend is removed.     

Oscillations of the intertropical convergence zone and the genesis of easterly waves. Part I: diagnostics and theory

We examine the mean and transient state of the intertropical convergence zone (ITCZ) by analyzing data and using simple theory. We concentrate on the tropical eastern Pacific Ocean noting that there exists in this region a well-developed mean ITCZ. Furthermore, it is a region where there has been considerable discussion in the literature of whether easterly waves develop in situ or propagate westwards from the Atlantic Ocean. The region is typical of tropical regions where there is a strong cross-equatorial pressure gradient (CEPG): mean convection well removed from the equator but located equatorward of the maximum sea-surface temperature (SST) and minimum sea level pressure (MSLP). Further to the west, near the dateline where the CEPG is much smaller, convection is weaker and collocated with SST and MSLP extrema. It is argued that in regions of significant CEPG that the near-equatorial tropical system is inertially unstable and that the rectification of the instability for a given CEPG determines the location and intensity of the climatological ITCZ. Using simple theoretical arguments, we develop an expression for the mean latitude of the ITCZ as a function of the CEPG. We note on a day-by-day basis that the ITCZ is highly transient state with variability occurring on 3–8 day time scales. Transients with amplitudes about half of the mean ITCZ, propagate northwards from the near-equatorial southern hemisphere as anomalous meridional oscillations, eventually amplifying convection in the vicinity of the mean ITCZ. It is argued that in these longitudes of strong CEPG the mean ITCZ is continually inertially unstable with advections of anticylonic vorticity across the equator resulting in the creation of an oscillating divergence–convergence doublet. The low-level convergence produces convection and the resultant vortex tube stretching generates cyclonic vorticity which counteracts the northward advection of anticylonic vorticity. During a cycle, the mid-troposphere heating near 10oN oscillates between 6 and 12 K/day at the inertial frequency of the latitude of the mean convection. As a result, there exists an anomalous and shallower, oscillating meridional circulation with a magnitude about 50% of the mean ITCZ associated with the stable state following the generation of anticylonic vorticity. Further, it is argued that the instabilities of the ITCZ are directly associated with in situ development of easterly waves which occur with the inertial period of the latitude of the mean ITCZ. The dynamical sequences and the genesis of easterly waves are absent in the regions further to the east where the CEPG is much weaker or absent altogether. In a companion study (Part II), numerical experiments are conducted to test the hypothesis raised in the present study.     

North African climate variability. Part 2: Tropical circulation systems

Summary Tropical North African climate variability is investigated using a Sahel rainfall index and streamflow of the Nile River in the 20^th century. The mechanisms that govern tropical North Africa climate are diagnosed from NCEP reanalysis data in the period 1958–1998: spatially – using composite and correlation analysis, and temporally – using wavelet co-spectral analysis. The Sahelian climate is characterised by a decadal rhythm, whilst the mountainous eastern and equatorial regions exhibit interannual cycles. ENSO-modulated zonal circulations over the Atlantic/Pacific sector are important for decadal variations, and create a climatic polarity between South America and tropical North Africa as revealed through upper-level velocity potential and convection patterns. A more localised N–S shift in convection between the Sahel and Guinea coast is associated with the African Easterly Jet.     

North African climate variability. Part 1: Tropical thermocline coupling

Summary Tropical ocean thermocline variability is studied using gridded data assimilated by an ocean model in the period 1950–2000. The dominant patterns and variability are identified using EOF analysis applied to E–W depth slices of sea temperatures averaged over the tropics. After removing the annual cycle, an east–west ‘see-saw’ with an interannual to decadal rhythm is the leading mode in each of the tropical basins. In the case of the leading mode in the Pacific, the thermocline oscillation forms a dipole structure, but in the (east) Atlantic and (southwest) Indian Ocean there is a single center of action. The interaction of the ocean thermocline and atmospheric Walker circulations is studied through cross-modulus analysis of wavelet-filtered EOF time scores. Our study demonstrates how tropical ocean thermocline variability contributes to zonal circulation anomalies in the atmosphere. The equatorial Pacific thermocline oscillation explains 62 and 53% of the variability of the Pacific and Atlantic zonal overturning circulations, the latter driving convective polarity between North Africa and South America. The Pacific sea-saw leads the Atlantic zonal circulation by a few months.     

Composite patterns of easterly disturbances over West Africa and the tropical Atlantic: a climatology from the 1979–95 NCEP/NCAR reanalyses

 The horizontal and vertical structure of the 3–5-day and 6–9-day easterly waves over West Africa and tropical Atlantic are investigated. NCEP/NCAR reanalyses are used for the period 1979–1995 to produce a 17-year climatology of both 3–5-day and 6–9-day easterly waves. Composite patterns of convection, wind, temperature and vertical velocity are analysed with respect to the following: the modulation by 3–5-day and 6–9-day wave regimes; the contrasts between the ITCZ (5°N–10°N) and the Sahelo-Saharan band (15°N–20°N); the difference between land and ocean, and seasonal variations. Similarities and differences in the characteristics of the two wave regimes are identified. Received: 18 August 1999 / Accepted: 14 March 2001     

Variability of seasonal-mean fields arising from intraseasonal variability. Part 3: Application to SH winter and summer circulations

A study has been made, using the National Centers for Environmental Prediction and National Center for Atmospheric Research re-analysis 500 hPa geopotential height data, to determine how intraseasonal variability influences, or can generate, coherent patterns of interannual variability in the extratropical summer and winter Southern Hemisphere atmospheric circulation. In addition, by separating this intraseasonal component of interannual variability, we also consider how slowly varying external forcings and slowly varying (interannual and longer) internal dynamics might influence the interannual variability of the Southern Hemisphere circulation. This slow component of interannual variation is more likely to be potentially predictable. How sea surface temperatures are related to the slow components is also considered. The four dominant intraseasonal modes of interannual variability have horizontal structures similar to those seen in both well-known intraseasonal dynamical modes and statistical modes of intraseasonal variability. In particular, they reflect intraseasonal variability in the high latitudes associated with the Southern Annular Mode, and wavenumber 4 (summer) and wavenumber 3 (winter) patterns associated with south Pacific regions of persistent anomalies and blocking, and possibly variability related to the Madden-Julian Oscillation (MJO). The four dominant slow components of interannual variability, in both seasons, are related to high latitude variability associated with the Southern Annular Mode, El Nino Southern Oscillation (ENSO) variability, and South Pacific Wave variability associated with Indian Ocean SSTs. In both seasons, there are strong linear trends in the first slow mode of high latitude variability and these are shown to be related to similar trends in the Indian Ocean. Once these are taken into account there is no significant sea surface temperature forcing of these high latitude modes. The second and third ENSO related slow modes, in each season, have high correlations with tropical sea surface temperature variability in the Pacific and Indian Oceans, both contemporaneously and at one season lag. The fourth slow mode has a characteristic South Pacific wave structure of either a wavenumber 4 (summer) or wavenumber 3 (winter) pattern, with strongest loadings in the South Pacific sector, and an association simultaneously with a dipole SST temperature gradient in the subtropical Indian Ocean.     

Surface energy balance complexity in GCM land surface models. Part II: coupled simulations

 Global coupled simulations with the Bureau of Meteorology Research Centre climate model and the CHAmeleon Surface Model (CHASM) are used to examine how four general extensions to the representation of the basic land surface energy balance affect simulated land-atmosphere interface variables: evaporation, precipitation, skin temperature and air temperature. The impacts of including separate surface energy balance calculations for: vegetated and non-vegetated portions of the land surface; an explicit parametrisation of canopy resistance; explicit bare ground evaporation; and explicit canopy interception are isolated and quantified. The hypothesis that these aspects of surface energy balance parametrisation do not contain substantial information at the monthly time scale (and are therefore not important to consider in a land surface model) is shown to be false. Considerable sensitivity to each of the four general surface energy balance extensions is identified in average pointwise monthly changes for important land-atmosphere interface variables. Average pointwise changes in monthly precipitation and land evaporation are equal to about 40 and 31–37% of the global-average precipitation and land evaporation respectively. Average pointwise changes for land surface skin temperature and lowest model layer air temperature are about 2 and 0.9 K respectively. The average pointwise change and average pointwise biases are statistically significant at 95% in all cases. Substantial changes to zonally average variables are also identified. We demonstrate how the globally averaged surface resistance parameter can vary from 150 to 25 s/m depending on which aspects of the surface energy balance are treated implicitly. We also show that if interception is treated implicitly, the effective surface resistance must vary geographically in order to capture the behaviour of a model which treats this process explicitly. The implication of these results for the design of land surface models is discussed. Received: 8 July 1999 / Accepted: 1 September 2000