Present and future near-surface wind climate of Greenland from high resolution regional climate modelling

The present and twenty-first century near-surface wind climate of Greenland is presented using output from the regional atmospheric climate model RACMO2. The modelled wind variability and wind distribution compare favourably to observations from three automatic weather stations in the ablation zone of southwest Greenland. The Weibull shape parameter is used to classify the wind climate. High values (κ > 4) are found in northern Greenland, indicative of uniform winds and a dominant katabatic forcing, while lower values (κ < 3) are found over the ocean and southern Greenland, where the synoptic forcing dominates. Very high values of the shape parameter are found over concave topography where confluence strengthens the katabatic circulation, while very low values are found in a narrow band along the coast due to barrier winds. To simulate the future (2081–2098) wind climate RACMO2 was forced with the HadGEM2-ES general circulation model using a scenario of mid-range radiative forcing of +4.5 W m^?2 by 2100. For the future simulated climate, the near-surface potential temperature deficit reduces in all seasons in regions where the surface temperature is below the freezing point, indicating a reduction in strength of the near-surface temperature inversion layer. This leads to a wind speed reduction over the central ice sheet where katabatic forcing dominates, and a wind speed increase over steep coastal topography due to counteracting effects of thermal and katabatic forcing. Thermally forced winds over the seasonally sea ice covered region of the Greenland Sea are reduced by up to 2.5 m s^?1.     

A 40-year accumulation dataset for Adelie Land, Antarctica and its application for model validation

The GLACIOCLIM-SAMBA (GS) Antarctic accumulation monitoring network, which extends from the coast of Adelie Land to the Antarctic plateau, has been surveyed annually since 2004. The network includes a 156-km stake-line from the coast inland, along which accumulation shows high spatial and interannual variability with a mean value of 362?mm water equivalent a^?1. In this paper, this accumulation is compared with older accumulation reports from between 1971 and 1991. The mean and annual standard deviation and the km-scale spatial pattern of accumulation were seen to be very similar in the older and more recent data. The data did not reveal any significant accumulation trend over the last 40?years. The ECMWF analysis-based forecasts (ERA-40 and ERA-Interim), a stretched-grid global general circulation model (LMDZ4) and three regional circulation models (PMM5, MAR and RACMO2), all with high resolution over Antarctica (27–125?km), were tested against the GS reports. They qualitatively reproduced the meso-scale spatial pattern of the annual-mean accumulation except MAR. MAR significantly underestimated mean accumulation, while LMDZ4 and RACMO2 overestimated it. ERA-40 and the regional models that use ERA-40 as lateral boundary condition qualitatively reproduced the chronology of interannual variability but underestimated the magnitude of interannual variations. Two widely used climatologies for Antarctic accumulation agreed well with the mean GS data. The model-based climatology was also able to reproduce the observed spatial pattern. These data thus provide new stringent constraints on models and other large-scale evaluations of the Antarctic accumulation.     

Turbulence characteristics in the atmospheric surface layer during summer monsoon of 1997 over a semi-arid location in India

A land surface processes experiment (LASPEX) was conducted in the semi-arid region of Northwest India during January 1997–February 1998. Analysis of turbulent components of wind and air temperature collected in the surface layer (SL) at Anand (22°35′N, 72°55′E) during the Indian summer monsoon season from June to September 1997 is presented. Turbulent fluctuation of wind components and air temperature observed at Anand varied as a function of terrain features and stability of the surface layer. Under neutral conditions, the standard deviation of vertical velocity (σ _w ) and temperature (σ _T ) were normalized using respective surface layer scaling parameter u _* and T _* which fitted the expressions σ _w /u _* = 1.25 and σ _T /T _* ≈ 4. Micrometeorological spectrum of wind and temperature at 5 m above ground level (AGL) at Anand showed peaks at time scale of 1–3 min at the low-frequency end. The inertial sub-range characteristics (?2/3 slope) of the spectrum are exhibited mostly. However, in some occasions, slope of ?1 denoting brown noise was depicted by the wind and temperature spectrum, which indicated anisotropy in turbulence.     

Weibull-k Revisited: “Tall” Profiles and Height Variation of Wind Statistics

The Weibull distribution is commonly used to describe climatological wind-speed distributions in the atmospheric boundary layer. While vertical profiles of mean wind speed in the atmospheric boundary layer have received significant attention, the variation of the shape of the wind distribution with height is less understood. Previously we derived a probabilistic model based on similarity theory for calculating the effects of stability and planetary boundary-layer depth upon long-term mean wind profiles. However, some applications (e.g. wind energy estimation) require the Weibull shape parameter (k), as well as mean wind speed. Towards the aim of improving predictions of the Weibull- \(k\) profile, we develop expressions for the profile of long-term variance of wind speed, including a method extending our probabilistic wind-profile theory; together these two profiles lead to a profile of Weibull-shape parameter. Further, an alternate model for the vertical profile of Weibull shape parameter is made, improving upon a basis set forth by Wieringa (Boundary-Layer Meteorol, 1989, Vol. 47, 85–110), and connecting with a newly-corrected corollary of the perturbed geostrophic-drag theory of Troen and Petersen (European Wind Atlas, 1989, Risø National Laboratory, Roskilde). Comparing the models for Weibull-k profiles, a new interpretation and explanation is given for the vertical variation of the shape of wind-speed distributions. Results of the modelling are shown for a number of sites, with a discussion of the models’ efficacy and applicability. The latter includes a comparative evaluation of Wieringa-type empirical models and perturbed-geostrophic forms with regard to surface-layer behaviour, as well as for heights where climatological wind-speed variability is not dominated by surface effects.     

An assessment of oceanic variability for 1960–2010 from the GFDL ensemble coupled data assimilation

The Geophysical Fluid Dynamics Laboratory has developed an ensemble coupled data assimilation (ECDA) system based on the fully coupled climate model, CM2.1, in order to provide reanalyzed coupled initial conditions that are balanced with the climate prediction model. Here, we conduct a comprehensive assessment for the oceanic variability from the latest version of the ECDA analyzed for 51 years, 1960–2010. Meridional oceanic heat transport, net ocean surface heat flux, wind stress, sea surface height, top 300 m heat content, tropical temperature, salinity and currents are compared with various in situ observations and reanalyses by employing similar configurations with the assessment of the NCEP’s climate forecast system reanalysis (Xue et al. in Clim Dyn 37(11):2511–2539, 2011). Results show that the ECDA agrees well with observations in both climatology and variability for 51 years. For the simulation of the Tropical Atlantic Ocean and global salinity variability, the ECDA shows a good performance compared to existing reanalyses. The ECDA also shows no significant drift in the deep ocean temperature and salinity. While systematic model biases are mostly corrected with the coupled data assimilation, some biases (e.g., strong trade winds, weak westerly winds and warm SST in the southern oceans, subsurface temperature and salinity biases along the equatorial western Pacific boundary, overestimating the mixed layer depth around the subpolar Atlantic and high-latitude southern oceans in the winter seasons) are not completely eliminated. Mean biases such as strong South Equatorial Current, weak Equatorial Under Current, and weak Atlantic overturning transport are generated during the assimilation procedure, but their variabilities are well simulated. In terms of climate variability, the ECDA provides good simulations of the dominant oceanic signals associated with El Nino and Southern Oscillation, Indian Ocean Dipole, Pacific Decadal Oscillation, and Atlantic Meridional Overturning Circulation during the whole analyzed period, 1960–2010.     

Snow accumulation and its moisture origin over Dome Argus, Antarctica

The spatial and temporal variability of snow accumulation near Dome Argus, Antarctica, is assessed using new snow pit and stake measurement data together with existing snow pit, ice core and automatic weather station records. Snow accumulation rate shows large inter-annual variations, but stable multi-decadal levels over the last seven centuries. Spatial variations in snow accumulation within the space of 50 km of Dome Argus are relatively small, probably thanks to the smooth topography. A comparison of theses accumulation observations with ECMWF reanalyses (ERA-40 and ERA-Interim) suggests ECMWF reanalysis captures the seasonal variations, but underestimates the overall snow accumulation at Dome Argus by ~50 %. The moisture sources for precipitation over Dome Argus are examined by means of a Lagrangian moisture source diagnostic, based on the tracing of specific humidity changes along air parcel trajectories, for the period 2000–2004 using operational ECMWF analysis data. Dome Argus mainly receives moisture from the mid-latitude (46 ± 4°S) South Indian Ocean, with a seasonal latitudinal shift of about 6°. Compared to other central East Antarctic deep ice core sites such as Dome F, Dome C, Vostok, and EPICA Dronning Maud Land, Dome Argus has a more southerly moisture origin, probably due to topographic influences on the moisture transport paths. These results have important implications for the interpretation of future ice cores at Dome Argus.     

Sources of uncertainty in projections of twenty-first century westerly wind changes over the Amundsen Sea,West Antarctica,in CMIP5 climate models

The influence of changes in winds over the Amundsen Sea has been shown to be a potentially key mechanism in explaining rapid loss of ice from major glaciers in West Antarctica, which is having a significant impact on global sea level. Here, Coupled Model Intercomparison Project Phase 5 (CMIP5) climate model data are used to assess twenty-first century projections in westerly winds over the Amundsen Sea (U _AS ). The importance of model uncertainty and internal climate variability in RCP4.5 and RCP8.5 scenario projections are quantified and potential sources of model uncertainty are considered. For the decade 2090–2099 the CMIP5 models show an ensemble mean twenty-first century response in annual mean U _AS of 0.3 and 0.7 m s^?1 following the RCP4.5 and RCP8.5 scenarios respectively. However, as a consequence of large internal climate variability over the Amundsen Sea, it takes until around 2030 (2065) for the RCP8.5 response to exceed one (two) standard deviation(s) of decadal internal variability. In all scenarios and seasons the model uncertainty is large. However the present-day climatological zonal wind bias over the whole South Pacific, which is important for tropical teleconnections, is strongly related to inter-model differences in projected change in U _AS (more skilful models show larger U _AS increases). This relationship is significant in winter (r = ?0.56) and spring (r = ?0.65), when the influence of the tropics on the Amundsen Sea region is known to be important. Horizontal grid spacing and present day sea ice extent are not significant sources of inter-model spread.     

GPCC's new land surface precipitation climatology based on quality-controlled in situ data and its role in quantifying the global water cycle

In 1989, the need for reliable gridded land surface precipitation data sets, in view of the large uncertainties in the assessment of the global energy and water cycle, has led to the establishment of the Global Precipitation Climatology Centre (GPCC) at Deutscher Wetterdienst on invitation of the WMO. The GPCC has calculated a precipitation climatology for the global land areas for the target period 1951–2000 by objective analysis of climatological normals of about 67,200 rain gauge stations from its data base. GPCC's new precipitation climatology is compared to several other station-based precipitation climatologies as well as to precipitation climatologies derived from the GPCP V2.2 data set and from ECMWF's model reanalyses ERA-40 and ERA-Interim. Finally, how GPCC's best estimate for terrestrial mean precipitation derived from the precipitation climatology of 786 mm per year (equivalent to a water transport of 117,000 km³) is fitting into the global water cycle context is discussed.     

Downscaling over Vietnam using the stretched-grid CCAM: verification of the mean and interannual variability of rainfall

Rainfall over Vietnam is highly variable from north to south, due to the interaction of the monsoonal winds with the terrain. There is high rainfall from April to September, and little rainfall from October to March (except along the central Vietnam coast). In order to study the ability of the Commonwealth Scientific and Industrial Research Organisation stretched-grid Conformal Cubic Atmospheric Model (CCAM) to capture the climatic and interannual variability of rainfall, downscaled simulations at approximately 20 km horizontal resolution over the region were produced for the period 1979–2001. A scale-selective digital filter was used to force the winds, temperature and sea-level pressure from the ERA-Interim reanalysis for length scales greater than about 700 km. For wind and temperature, the forcing is applied for pressure-sigma levels above about 0.9. ERA-Interim sea surface temperatures were used over the oceans. The simulations were primarily validated against the gridded Asian Precipitation Highly Resolved Observational Data Integration Toward Evaluation of the Water Resources rainfall dataset and station observations using standard statistical methods. It was found that CCAM reproduces well the amount and spatial variability of rainfall, with an area-averaged bias for the entire study domain of less than 1 mm day^?1; CCAM is also able to capture the rainfall pattern under different El Niño Southern Oscillation phases reasonably well for the dry season. For interannual variability, the simulation generally performed better for North and Central Vietnam than for South Vietnam, where rainfall variability was overestimated.     

Cold air outbreaks along a non-frozen sea channel: effects of wind on snow bands

Wintertime cold air outbreaks along a non-frozen sea channel or a long lake can become destructive if the related bands of heavy snowfall hit onto land. The forcing for such bands is studied with a 2D numerical model set across an east–west sea channel at 60^oN (‘Gulf of Finland’), varying the basic geostrophic wind V _g. Without any V _g opposite coastal land breezes emerge with convergence. This results in a quasi-steady rising motion w _max ~ 7.5 cm/s at 600 m in the middle of the gulf, which can force a snow band. During weak V _g, the rising motion is reduced but least so for winds from 60^o to 80^o (~ENE), when modest alongshore bands could exist near the downstream (Estonian) coast. During V _g of 4–6 m/s from any direction, the land breezes and rising motions are reduced more effectively, so snow bands are not expected during moderate basic flow. In contrast, during a strong V _g of 20–25 m/s from 110^o to 120^o (~ESE) the land breeze perturbations are intense with w _max up to 15–18 cm/s. The induced alongshore bands of heavy snowfall are located in these cases at the sea but quite close to the downstream (Finnish) coast. They can suddenly make a landfall if the basic wind turns clockwise.     

Decadal variation of surface solar radiation in the Tibetan Plateau from observations, reanalysis and model simulations

In this study, the annual and seasonal variations of all-sky and clear-sky surface solar radiation (SSR) in the eastern and central Tibetan Plateau (TP) during the period 1960–2009 are investigated, based on surface observational data, reanalyses and ensemble simulations with the global climate model ECHAM5-HAM. The mean annual all-sky SSR series shows a decreasing trend with a rate of ?1.00 Wm^?2 decade^?1, which is mainly seen in autumn and secondly in summer and winter. A stronger decrease of ?2.80 Wm^?2 decade^?1 is found in the mean annual clear-sky SSR series, especially during winter and autumn. Overall, these results confirm a tendency towards a decrease of SSR in the TP during the last five decades. The comparisons with reanalysis show that both NCEP/NCAR and ERA-40 reanalyses do not capture the decadal variations of the all-sky and clear-sky SSR. This is probably due to a missing consideration of aerosols in the reanalysis assimilation model. The SSR simulated with the ECHAM5-HAM global climate model under both all-sky and clear-sky conditions reproduce the decrease seen in the surface observations, especially after 1980. The steadily increasing aerosol optical depth (AOD) at 550 nm over the TP in the ECHAM5-HAM results suggests transient aerosol emissions as a plausible cause.     

Influence of non-flat terrain and wind direction shear on canopy turbulence

We have analyzed eddy covariance data collected within open canopy to investigate the influence of non-flat terrain and wind direction shear on the canopy turbulence. The study site is located on non-flat terrain with slopes in both south-north and east-west directions. The surface elevation change is smaller than the height of roughness element such as building and tree at this site. A variety of turbulent statistics were examined as a function of wind direction in near-neutral conditions. Heterogeneous surface characteristics results in significant differences in measured turbulent statistics. Upwind trees on the flat and up-sloping terrains yield typical features of canopy turbulence while upwind elevated surface with trees yields significant wind direction shear, reduced u and w skewness, and negligible correlation between u and w. The directional dependence of turbulence statistics is due that strong wind blows more horizontally rather than following terrain, and hence combination of slope related momentum flux and canopy eddy motion decreases the magnitude of Sk _w and r _uw for the downslope flow while it enhances them for the upslope flow. Significant v skewness to the west indicates intermittent downdraft of northerly wind, possibly due to lateral shear of wind in the presence of significant wind direction shear. The effects of wind direction shear on turbulent statistics were also examined. The results showed that correlation coefficient between lateral velocities and vertical velocity show significant dependence on wind direction shear through change of lateral wind shear. Quadrant analysis shows increased outward interaction and reduced role of sweep motion for longitudinal momentum flux for the downslope flow. Multi-resolution analysis indicates that uw correlation shows peak at larger averaging time for the upslope flow than for the downslope flow, indicating that large eddy plays an active role in momentum transfer for the upslope flow. On the other hand, downslope flow shows larger velocity variances than other flows despite similar wind speed. These results suggest that non-flatness of terrain significantly influences on canopy-atmosphere exchange.     

Shapes of annual frequency distributions of wind speed observed on high meteorological masts

Annual wind distributions from six masts are used to investigate the annual variability of hourly-average wind. A variability extremum occurs at the height where the average diurnal cycle of wind reverses its phase from a nocturnal minimum to a nocturnal maximum. A simple profile model shows that in non-complex terrain, this reversal height varies approximately between 50 m at coasts and 90 m inland. The Weibull distribution shape factor k has a maximum at the reversal height, and with decreasing height, k decreases approximately as a linear function of height. Therefore estimation of wind distribution shape from surface data is possible in the surface layer. In the upper PBL, however, such estimation is not very feasible, because no simple methods are available to estimate the reversal height from routine surface observations.Also: Dept. of Applied Physics, Technical University, Delft, Netherlands.     

Future surface mass balance of the Antarctic ice sheet and its influence on sea level change, simulated by a regional atmospheric climate model

A regional atmospheric climate model with multi-layer snow module (RACMO2) is forced at the lateral boundaries by global climate model (GCM) data to assess the future climate and surface mass balance (SMB) of the Antarctic ice sheet (AIS). Two different GCMs (ECHAM5 until 2100 and HadCM3 until 2200) and two different emission scenarios (A1B and E1) are used as forcing to capture a realistic range in future climate states. Simulated ice sheet averaged 2 m air temperature (T_2m) increases (1.8–3.0 K in 2100 and 2.4–5.3 K in 2200), simultaneously and with the same magnitude as GCM simulated T_2m. The SMB and its components increase in magnitude, as they are directly influenced by the temperature increase. Changes in atmospheric circulation around Antarctica play a minor role in future SMB changes. During the next two centuries, the projected increase in liquid water flux from rainfall and snowmelt, together 60–200 Gt year^?1, will mostly refreeze in the snow pack, so runoff remains small (10–40 Gt year^?1). Sublimation increases by 25–50 %, but remains an order of magnitude smaller than snowfall. The increase in snowfall mainly determines future changes in SMB on the AIS: 6–16 % in 2100 and 8–25 % in 2200. Without any ice dynamical response, this would result in an eustatic sea level drop of 20–43 mm in 2100 and 73–163 mm in 2200, compared to the twentieth century. Averaged over the AIS, a strong relation between $\Updelta$ SMB and $\Updelta\hbox{T}_{2{\rm m}}$ of 98 ± 5 Gt w.e. year^?1 K^?1 is found.     

Using the Twentieth Century Reanalysis to assess climate variability for the European wind industry

We characterise the long-term variability of European near-surface wind speeds using 142 years of data from the Twentieth Century Reanalysis (20CR), and consider the potential of such long-baseline climate data sets for wind energy applications. The low resolution of the 20CR would severely restrict its use on its own for wind farm site-screening. We therefore perform a simple statistical calibration to link it to the higher-resolution ERA-Interim data set (ERAI), such that the adjusted 20CR data has the same wind speed distribution at each location as ERAI during their common period. Using this corrected 20CR data set, wind speeds and variability are characterised in terms of the long-term mean, standard deviation and corresponding trends. Many regions of interest show extremely weak trends on century timescales, but contain large multidecadal variability. Since reanalyses such as ERAI are often used to provide the background climatology for wind farm site assessments, but contain only a few decades of data, our results can be used as a way of incorporating decadal-scale wind climate variability into such studies, allowing investment risks for wind farms to be reduced.     

Surface mean temperature from the observational stations and multiple reanalyses over the Tibetan Plateau

The Tibetan Plateau (TP), also called the “Third pole”, is sensitive to climate change due to extensive areas at high elevation presently dominated by snow and ice. In this study, observed surface temperature trends at 150 stations over the TP during 1979–2018 are analyzed and compared with surface temperatures from multiple reanalyses (NCEP1, NCEP2, ERA-Interim, MERRA, JRA55). Observed warming at the stations has a mean annual rate of 0.46 °C/decade during 1979–2018. Although all reanalyses underestimate observed temperatures (cold bias), most reproduce much of the inter-decadal variations of surface temperature shown in the observations. Absolute errors of mean surface temperature (reanalysis minus observation) are closely correlated with elevation errors, suggesting that parts of the cold bias can be interpreted by elevation errors of reanalysis. After elevation-temperature correction, about half of the cold bias is typically eliminated, more for both ERA-Interim and JRA55. Compared with the observations, corrected NCEP2 surface temperatures still have larger cold biases, and fail to capture the overall warming over the TP. Since the elevation-temperature correction fails to improve trend magnitudes even when a significant proportion of the bias has been removed, this suggests that a more sophisticated modeling of the lapse rate in each reanalysis is required to realistically model warming trends across complex topography.

     

High-resolution sea wind hindcasts over the Mediterranean area

The goal of this study is to develop a high-resolution atmospheric hindcast over the Mediterranean area using the WRF-ARW model, focusing on offshore surface wind fields. In order to choose the most adequate model configuration, the study provides details on the calibration of the experimental saet-up through a sensitivity test considering the October–December 2001 period (the 2001 super-storm event in the West Mediterranean). A daily forecast outperforms the spectral technique of previous products and the boundary data from ERA-Interim reanalysis produces the most accurate estimates in terms of wind variability and hour-to-hour correspondence. According to the sensitivity test, two data sets of wind hindcast are produced: the SeaWind I (30-km horizontal resolution for a period of 60 years) and the SeaWind II (15-km horizontal resolution for 20 years). The validation of the resulting surface winds is undertaken considering two offshore observational datasets. On the one hand, hourly surface buoy stations are used to validate wind time series at specific locations; on the other hand, wind altimeter satellite observations are considered for spatial validation in the whole Mediterranean Sea. The results obtained from this validation process show a very good agreement with observations for the southern Europe region. Finally, SeaWind I and II are used to characterize offshore wind fields in the Mediterranean Sea. The statistical structure of sea surface wind is analyzed and the agreement with Weibull probability distribution is discussed. In addition, wind persistence and extreme wind speed (50 year return period) are characterized and relevant areas of wind power generation are described by estimating wind energy quantities.     

Air–Sea Interaction Features in the Baltic Sea and at a Pacific Trade-Wind Site: An Inter-comparison Study

A systematic comparison of wind profiles and momentum exchange at a trade wind site outside Oahu, Hawaii and corresponding data from the Baltic Sea is presented. The trade wind data are to a very high degree swell dominated, whereas the Baltic Sea data include a more varied assortment of wave conditions, ranging from a pure growing sea to swell. In the trade wind region swell waves travel predominantly in the wind direction, while in the Baltic, significant cross-wind swells are also present. Showing the drag coefficient as a function of the 10-m wind speed demonstrates striking differences for unstable conditions with swell for the wind-speed range 2 m s^?1 < U ₁₀ < 7 m s^?1, where the trade-wind site drag values are significantly larger than the corresponding Baltic Sea values. In striking contrast to this disagreement, other features studied are surprisingly similar between the two sites. Thus, exactly as found previously in Baltic Sea studies during unstable conditions and swell, the wind profile in light winds (3 m s^?1) shows a wind maximum at around 7–8 m above the water, with close to constant wind speed above. Also, for slightly higher wind speeds (4 m s^?1 < U ₁₀ < 7 m s^?1), the similarity between wind profiles is striking, with a strong wind-speed increase below a height of about 7–8 m followed by a layer of virtually constant wind speed above. A consequence of these wind-profile features is that Monin–Obukhov similarity is no longer valid. At the trade-wind site this was observed to be the case even for wind speeds as high as 10 m s^?1. The turbulence kinetic energy budget was evaluated for four cases of 8–16 30- min periods at the trade-wind site, giving results that agree very well with corresponding figures from the Baltic Sea.     

Trends and variability in column-integrated atmospheric water vapor

An analysis and evaluation has been performed of global datasets on column-integrated water vapor (precipitable water). For years before 1996, the Ross and Elliott radiosonde dataset is used for validation of European Centre for Medium-range Weather Forecasts (ECMWF) reanalyses ERA-40. Only the special sensor microwave imager (SSM/I) dataset from remote sensing systems (RSS) has credible means, variability and trends for the oceans, but it is available only for the post-1988 period. Major problems are found in the means, variability and trends from 1988 to 2001 for both reanalyses from National Centers for Environmental Prediction (NCEP) and the ERA-40 reanalysis over the oceans, and for the NASA water vapor project (NVAP) dataset more generally. NCEP and ERA-40 values are reasonable over land where constrained by radiosondes. Accordingly, users of these data should take great care in accepting results as real. The problems highlight the need for reprocessing of data, as has been done by RSS, and reanalyses that adequately take account of the changing observing system. Precipitable water variability for 1988–2001 is dominated by the evolution of ENSO and especially the structures that occurred during and following the 1997–98 El Niño event. The evidence from SSM/I for the global ocean suggests that recent trends in precipitable water are generally positive and, for 1988 through 2003, average 0.40±0.09 mm per decade or 1.3±0.3% per decade for the ocean as a whole, where the error bars are 95% confidence intervals. Over the oceans, the precipitable water variability relates very strongly to changes in SSTs, both in terms of spatial structure of trends and temporal variability (with a regression coefficient for 30°N–30°S of 7.8% K^?1) and is consistent with the assumption of fairly constant relative humidity. In the tropics, the trends are also influenced by changes in rainfall which, in turn, are closely associated with the mean flow and convergence of moisture by the trade winds. The main region where positive trends are not very evident is over Europe, in spite of large and positive trends over the North Atlantic since 1988. A much longer time series is probably required to obtain stable patterns of trends over the oceans, although the main variability could probably be deduced from past SST and associated precipitation variations.     

The Congo basin zonal overturning circulation

The Gulf of Guinea in the equatorial Atlantic is characterized by the presence of strong subsidence at certain times of the year. This subsidence appears in June and becomes well established from July to September. Since much of theWest African monsoon flow originates over the Gulf, Guinean subsidence is important for determining moisture sources for the monsoon. Using reanalysis products, I contribute to a physical understanding of what causes this seasonal subsidence, and how it relates to precipitation distributions across West Africa.There is a seasonal zonal overturning circulation above the Congo basin and the Gulf of Guinea in the ERA-Interim, ERA-40, NCEP2, and MERRA reanalyses. The up-branch is located in the Congo basin around 20°E. Mid-tropospheric easterly flows constitute the returning-branch and sinking over the Gulf of Guinea forms the down-branch, which diverges at 2°W near the surface, with winds to the east flowing eastward to complete the circulation. This circulation is driven by surface temperature differences between the eastern Gulf and Congo basin. Land temperatures remain almost uniform, around 298 K, throughout a year, but the Guinean temperatures cool rapidly from 294 K in May to about 290 K in August. These temperature changes increase the ocean/land temperature contrast, up to 8 K, and drive the circulation.I hypothesize that when the overturning circulation is anomalously strong, the northward moisture transport and Sahelian precipitation are also strong. This hypothesis is supported by ERA-Interim and PERSIANN-CDR (Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Climate Data Record) data.