Gilbert P. Compo

A practical guide to wavelet analysis

A practical step-by-step guide to wavelet analysis is given, with examples taken from time series of the El Nino–Southern Oscillation (ENSO). The guide includes a comparison to the windowed Fourier transform, the choice of an appropriate wavelet basis function, edge effects due to finite-length time series, and the relationship between wavelet scale and Fourier frequency. New statistical significance tests for wavelet power spectra are developed by deriving theoretical wavelet spectra for white and red noise processes and using these to establish significance levels and confidence intervals. It is shown that smoothing in time or scale can be used to increase the confidence of the wavelet spectrum. Empirical formulas are given for the effect of smoothing on significance levels and confidence intervals. Extensions to wavelet analysis such as filtering, the power Hovmoller, cross-wavelet spectra, and coherence are described. The statistical significance tests are used to give a quantitative measure of change...

Feasibility of a 100-Year Reanalysis Using Only Surface Pressure Data

Climate variability and global change studies are increasingly focused on understanding and predicting regional changes of daily weather statistics. Assessing the evidence for such variations over the last 100 yr requires a daily tropospheric circulation dataset. The only dataset available for the early twentieth century consists of error-ridden hand-drawn analyses of the mean sea level pressure field over the Northern Hemisphere. Modern data assimilation systems have the potential to improve upon these maps, but prior to 1948, few digitized upper-air sounding observations are available for such a “reanalysis.” We investigate the possibility that the additional number of newly recovered surface pressure observations is sufficient to generate useful weather maps of the lower-tropospheric extratropical circulation back to 1890 over the Northern Hemisphere, and back to 1930 over the Southern Hemisphere. Surprisingly, we find that by using an advanced data assimilation system based on an ensemble Kalman filte...

Reanalysis without Radiosondes Using Ensemble Data Assimilation

Abstract Studies using idealized ensemble data assimilation systems have shown that flow-dependent background-error covariances are most beneficial when the observing network is sparse. The computational cost of recently proposed ensemble data assimilation algorithms is directly proportional to the number of observations being assimilated. Therefore, ensemble-based data assimilation should both be more computationally feasible and provide the greatest benefit over current operational schemes in situations when observations are sparse. Reanalysis before the radiosonde era (pre-1931) is just such a situation. The feasibility of reanalysis before radiosondes using an ensemble square root filter (EnSRF) is examined. Real surface pressure observations for 2001 are used, subsampled to resemble the density of observations we estimate to be available for 1915. Analysis errors are defined relative to a three-dimensional variational data assimilation (3DVAR) analysis using several orders of magnitude more observati...

Changes of Probability Associated with El Niño

Abstract Away from the tropical Pacific Ocean, an ENSO event is associated with relatively minor changes of the probability distributions of atmospheric variables. It is nonetheless important to estimate the changes accurately for each ENSO event, because even small changes of means and variances can imply large changes of the likelihood of extreme values. The mean signals are not strictly symmetric with respect to El Nino and La Nina. They also depend upon the unique aspects of the SST anomaly patterns for each event. As for changes of variance and higher moments, little is known at present. This is a concern especially for precipitation, whose distribution is strongly skewed in areas of mean tropospheric descent. These issues are examined here in observations and GCM simulations of the northern winter (January–March, JFM). For the observational analysis, the 42-yr (1958–99) reanalysis data generated at NCEP are stratified into neutral, El Nino, and La Nina winters. The GCM analysis is based on NCEP atmo...

Removing ENSO-Related Variations from the Climate Record

Abstract An important question in assessing twentieth-century climate change is to what extent have ENSO-related variations contributed to the observed trends. Isolating such contributions is challenging for several reasons, including ambiguities arising from how ENSO itself is defined. In particular, defining ENSO in terms of a single index and ENSO-related variations in terms of regressions on that index, as done in many previous studies, can lead to wrong conclusions. This paper argues that ENSO is best viewed not as a number but as an evolving dynamical process for this purpose. Specifically, ENSO is identified with the four dynamical eigenvectors of tropical SST evolution that are most important in the observed evolution of ENSO events. This definition is used to isolate the ENSO-related component of global SST variations on a month-by-month basis in the 136-yr (1871–2006) Hadley Centre Sea Ice and Sea Surface Temperature dataset (HadISST). The analysis shows that previously identified multidecadal v...

The Asian Monsoon, the Tropospheric Biennial Oscillation, and the Indian Ocean Zonal Mode in the NCAR CSM*

Abstract The interaction of the Indian Ocean dynamics and the tropospheric biennial oscillation (TBO) is analyzed in the 300-yr control run of the National Center for Atmospheric Research (NCAR) Climate System Model (CSM). Sea surface temperature (SST) anomalies and equatorial ocean dynamics in the Indian Ocean are associated with the TBO and interannual variability of Asian–Australian monsoons in observations. The air–sea interactions involved in these processes in the coupled ocean–atmosphere model are analyzed, so as to diagnose the causes of the SST anomalies and their role in the development of a biennial cycle in the Indian–Pacific Ocean region. By using singular value decomposition (SVD) analysis, it is found that the model reproduces the dominant mechanisms that are involved in the development of the TBOs influence on the south Asian monsoon: large-scale forcing from the tropical Pacific and regional forcing associated with both the meridional temperature gradient between the Asian continent and ...

Trends and low frequency variability of extra-tropical cyclone activity in the ensemble of twentieth century reanalysis

An objective cyclone tracking algorithm is applied to twentieth century reanalysis (20CR) 6-hourly mean sea level pressure fields for the period 1871–2010 to infer historical trends and variability in extra-tropical cyclone activity. The tracking algorithm is applied both to the ensemble-mean analyses and to each of the 56 ensemble members individually. The ensemble-mean analyses are found to be unsuitable for accurately determining cyclone statistics. However, pooled cyclone statistics obtained by averaging statistics from individual members generally agree well with statistics from the NCEP-NCAR reanalyses for 1951–2010, although 20CR shows somewhat weaker cyclone activity over land and stronger activity over oceans. Both reanalyses show similar cyclone trend patterns in the northern hemisphere (NH) over 1951–2010. Homogenized pooled cyclone statistics are analyzed for trends and variability. Conclusions account for identified inhomogeneities, which occurred before 1949 in the NH and between 1951 and 1985 in the southern hemisphere (SH). Cyclone activity is estimated to have increased slightly over the period 1871–2010 in the NH. More substantial increases are seen in the SH. Notable regional and seasonal variations in trends are evident, as is profound decadal or longer scale variability. For example, the NH increases occur mainly in the mid-latitude Pacific and high-latitude Atlantic regions. For the North Atlantic-European region and southeast Australia, the 20CR cyclone trends are in agreement with trends in geostrophic wind extremes derived from in-situ surface pressure observations. European trends are also consistent with trends in the mean duration of wet spells derived from rain gauge data in Europe.

The International Atmospheric Circulation Reconstructions over the Earth (ACRE) Initiative

In 2006, climate applications scientists in Queensland, Australia, asked the lead author if a longer and more complete historical weather record could be created and fed directly into various crop, pasture, and production models. Existing dynamical reanalyses were steps toward such a product, but they spanned only the last six decades and had well-known shortcomings. To meet the needs of application scientists, new reanalyses would have to extend much further back in time while maintaining accuracy with limited observations. They would also need to be disseminated in a way that is easy to use directly and to downscale to small regions.

Continental heat anomalies and the extreme melting of the Greenland ice surface in 2012 and 1889

Recent decades have seen increased melting of the Greenland ice sheet. On 11 July 2012, nearly the entire surface of the ice sheet melted; such rare events last occurred in 1889 and, prior to that, during the Medieval Climate Anomaly. Studies of the 2012 event associated the presence of a thin, warm elevated liquid cloud layer with surface temperatures rising above the melting point at Summit Station, some 3212 m above sea level. Here we explore other potential factors in July 2012 associated with this unusual melting. These include (1) warm air originating from a record North American heat wave, (2) transitions in the Arctic Oscillation, (3) transport of water vapor via an Atmospheric River over the Atlantic to Greenland, and (4) the presence of warm ocean waters south of Greenland. For the 1889 episode, the Twentieth Century Reanalysis and historical records showed similar factors at work. However, markers of biomass burning were evident in ice cores from 1889 which may reflect another possible factor in these rare events. We suggest that extreme Greenland summer melt episodes, such as those recorded recently and in the late Holocene, could have involved a similar combination of slow climate processes, including prolonged North American droughts/heat waves and North Atlantic warm oceanic temperature anomalies, together with fast processes, such as excursions of the Arctic Oscillation, and transport of warm, humid air in Atmospheric Rivers to Greenland. It is the fast processes that underlie the rarity of such events and influence their predictability.

Modulation of equatorial subseasonal convective episodes by tropical‐extratropical interaction in the Indian and Pacific Ocean regions

Composite relationships among outgoing longwave radiation, sea level pressure, surface winds, and upper tropospheric circulation are examined for northern winter during subseasonal episodes of eastward progression of convection from the Indian Ocean to the western Pacific. This evolution often culminates with westerly wind burst events and strong air-sea interaction associated with regional-scale convective blowups in the western equatorial Pacific. We first document some of these interactions in the composites for two timescales, the submonthly (6–30 days) and that of the Madden-Julian Oscillation (MJO) timescale (30–70 days). We then analyze the December 1992 period during the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE) to illustrate how these composite relationships are manifested in a case study. Convection in the Indian Ocean for the composites is shown to be associated with a northern hemisphere wave train at 200 mbar, arcing through the midlatitudes, that can contribute to convective blowups farther east on the submonthly (6–30 days) timescale in the Intertropical Convergence Zone (ITCZ) in the eastern Pacific. The eastern Asian trough that is part of this wave train is associated with pressure surges from the northern hemisphere and subsequent convection over Southeast Asia. As the MJO convective envelope moves east to Australasia, midlatitude wave trains in either hemisphere include upper level troughs east of Asia and Australia and pressure surges from either hemisphere that contribute to pressure rises over the Indonesian region and a subsequent shift of the convective envelope to the western Pacific. The vertical wind structure for the December 1992 case study is consistent with the composite surface and upper level winds and also shows strong vertical wind shear in the boundary layer, a sharply defined westerly maximum near 700 mbar and an intensification of the upper level easterlies near 100 mbar. Very deep westerlies (to 200 mbar) are confined to shorter timescales. The case study illustrates the various time and space scale interactions noted in the composites. Reciprocal interactions between the tropics and the midlatitudes on the submonthly and MJO timescales in both the composites and the case study involve pressure surges and wave interaction that influence subsequent convection as the convective envelope migrates eastward from the tropical Indian to Pacific Ocean region.