JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. D15, 4232, 10.1029/2001JD001250, 2002 Representation of the quasi-biennial oscillation in the tropical stratospheric wind by nonlinear principal component analysis Kevin Hamilton Department of Meteorology and International Pacific Research Center, University of Hawaii at Manoa, Honolulu, Hawaii William W. Hsieh Department of Earth and Ocean Sciences, University of British Columbia, Vancouver, British Columbia, Canada Received 27 August 2001; revised 24 December 2001; accepted 27 December 2001; published 3 August 2002. [1] The zonal winds at several levels between 70 and 10 hPa (roughly 20–30 km) measured at near-equatorial stations during 1956–2000 were analyzed to produce a onedimensional approximation. The neural network-based technique applied was the circular nonlinear principal component analysis (NLPCA.cir) designed to characterize quasiperiodic phenomena. The reconstructed height-time series of wind based on the onedimensional NLPCA.cir captures many of the characteristic features of the observed quasi-biennial oscillation (QBO). The nonlinear results were evaluated relative to comparable linear principal component analysis and found to produce a superior onedimensional representation of the data. The NLPCA.cir analysis produces a single time series of QBO phase based on data at all levels. This phase was then applied to demonstrate a strong correlation of the state of the tropical QBO and boreal winter highINDEX TERMS: 3334 Meteorology and Atmospheric latitude stratospheric temperatures. Dynamics: Middle atmosphere dynamics (0341, 0342); 3319 Meteorology and Atmospheric Dynamics: General circulation; 3374 Meteorology and Atmospheric Dynamics: Tropical meteorology; KEYWORDS: quasi-biennial oscillation, neural network, tropical stratosphere 1. Introduction [2] The prevailing winds in the equatorial stratosphere are observed to undergo an oscillation between strong westerly (i.e., eastward) and easterly (i.e., westward) winds, with the transition from one phase to the other occurring roughly every year. This quasi-biennial oscillation (QBO) of the circulation dominates any annual cycle or other variations in the low-latitude stratosphere (see Hamilton [1998a] and Baldwin et al. [2001] for recent reviews). The period of the oscillation has been observed to vary from cycle to cycle between roughly 22 and 33 months, but the mean period over the last half century of observations is close to 28 months. The variability of the period and the fact that the mean period appears to be incommensurate with the annual cycle or any harmonics or subharmonics of the annual cycle indicates that the QBO is driven primarily by the internal fluid mechanics of the atmosphere. Thus it is to be strongly contrasted with the familiar astronomically driven oscillations (seasonal cycle, tides) which are genuinely periodic. However, the stratospheric QBO is much more regular than other extensively studied oscillatory phenomena in the atmosphere (such as the Southern Oscillation or North Atlantic Oscillation). [3] The nearly periodic, but somewhat irregular, behavior makes the QBO one of the most intriguing phenomena in Copyright 2002 by the American Geophysical Union. 0148-0227/02/2001JD001250$09.00 ACL geophysics. A standard height-time record of monthly mean zonal wind constructed from balloon observations up to 10 hPa ($30 km) at near-equatorial stations has been carefully updated over the years [Wallace, 1973; Coy, 1979; Naujokat, 1986; Marquardt and Naujokat, 1997] and has become a familiar and widely analyzed data record. Among the characteristic properties that have been clearly documented in these data are (1) consistently downward phase propagation of wind reversals, (2) a strong tendency for the wind accelerations at any level to be concentrated into limited periods near the wind reversals (leading to a roughly square-wave appearance to the time series at any particular level), (3) a corresponding tendency for the reversals to be associated with the downward passage of strong vertical shear zones, (4) asymmetry between the wind reversals such that the easterly to westerly transition is more rapid than the westerly to easterly transition (and the peak westerly shear is stronger than the peak easterly shear), and (5) fairly regular downward propagation of the westerly transition, but a tendency found in many cycles for the easterly transition to ‘‘stall’’ for several months between about 30 and 50 hPa (this effect accounts for much of the period variability from cycle to cycle). [4] The somewhat complicated vertical structure of the QBO in zonal wind makes it difficult to characterize the ‘‘phase of the QBO’’ in a completely adequate way. The equatorial QBO is known to be correlated with a number of important aspects of atmospheric behavior, such as high-latitude stratospheric circulation [e.g., Holton and 3-1 ACL 3-2 HAMILTON AND HSIEH: REPRESENTATION OF THE QBO Figure 1. A schematic diagram illustrating the NN model for calculating the NLPCA with a circular node at the bottleneck (NLPCA.cir). The model is a standard feedforward NN (i.e., multilayer perception), with three ‘‘hidden’’ layers of variables or ‘‘neurons’’ (denoted by circles) sandwiched between the input layer x on the left and the output layer x0 on the right. Next to the input layer (with l neurons) is the encoding layer (with m neurons), followed by the ‘‘bottleneck’’ layer, then the decoding layer (with m neurons), and finally the output layer (with l neurons), i.e., a total of four layers of transfer functions are needed to map from the inputs to the outputs. A neuron vi at the ith layer receives its value from the neurons viÀ1 in the preceding layer, i.e., vi = fi(wi Á viÀ1 + b), where wi is a vector of weight parameters and b a bias parameter, and the transfer functions f1 and f3 are the hyperbolic tangent functions, while f2 and f4 are simply the identity functions. In NLPCA.cir, the bottleneck contains two neurons p and q confined to lie on a unit circle, i.e., only 1 degree of freedom as represented by the angle q. Effectively, a nonlinear function q = F(x) maps from the higher dimension input space to the lower dimension bottleneck space, followed by an inverse transform x0 = G(q) mapping from the bottleneck space back to the original space, as represented by the outputs. To make the outputs as close to the inputs as possible, the cost function J ¼ hk x À x0 k2 i, (i.e., the mean-square-error, MSE) is minimized (where angle brackets denote a sample or time mean). Through the optimization, the values of the weight and bias parameters are solved. Data compression is achieved by the bottleneck, yielding the nonlinear principal component (NLPC) q. In our application, there are seven inputs and seven outputs, as the wind data are available at seven vertical levels. Tan, 1980, 1982; Dunkerton and Baldwin, 1991; Baldwin and Dunkerton, 1998] and the distribution of stratospheric ozone and other trace constituents [e.g., Tung and Yang, 1994; Randel et al., 1998], and may be connected with such features of tropospheric climate as the Southern Oscillation [e.g., Xu, 1992; Gray et al., 1992; Geller et al., 1997] and even the frequency of occurrence of tropical cyclones [Gray, 1984]. These connections have generally been studied by using the zonal wind at one arbitrarily chosen level to characterize the phase of the QBO, but for many applications a quantity that could characterize the whole vertical structure at a given time would be more relevant. Simplified statistical representations of the QBO wind variations could also be useful in characterizing the modest annual cycle effects apparent in the observed record, and for comparison of observed QBO with corresponding features seen in comprehensive general circulation model (GCM) simulations [e.g., Takahashi, 1996, 1999; Horinouchi and Yoden, 1998; Hamilton et al., 1999, 2001]. Finally, recent work of Salby and Callaghan [2000] suggests that the period of the stratospheric QBO may be modulated by the 11-year solar activity cycle. As Salby and Callaghan noted, such a modulation would have significant implications for understanding and predicting solar effects on the global climate system. Again, Salby and Callaghan performed their analysis using a QBO period estimated by using winds at a single level. A characterization of the QBO phase evolution that used wind information over a deeper layer might be more meaningful. [5] There have been some more sophisticated efforts to obtain a time-series representation of the phase of the equatorial zonal wind QBO (along with an associated characteristic vertical structure). In particular, Wallace et al. [1993] computed the first two principal components (PCs) of the height-time section and found that the time series of the coefficients were quasi-cyclic and in quadrature, allowing the phase of the QBO to be characterized by the arctangent of the ratio of the coefficients. Fraedrich et al. [1993] and Wang et al. [1995] examined the problem via the singular spectrum analysis (often referred to as extended PC analysis) using a window of 40 months. They also found that the leading two extended PCs were in quadrature and described a steadily descending oscillation. [6] Unfortunately, the leading PCs or extended PCs describe only the crudest aspects of the zonal wind QBO and miss many details characteristic of the actual observations. Notably the leading PCs (or extended PCs) tend to describe a near-sinusoidal oscillation rather than the more square-wave behavior apparent in the real QBO. Also the first modes in the PC and extended-PC approaches capture virtually none of the asymmetry between westerly and easterly shear zones (and westerly and easterly acceleration regimes) that is such a prominent characteristic of the observations (in every QBO cycle observed so far). [7] The limitation of the leading PCs (or extended PCs) as descriptions of the full QBO is not surprising, as the PCs are just linear least squares fits to the data. Recently, nonlinear generalizations of the PC approach using neural network techniques have been developed and applied to the analysis of geophysical phenomena such as the El Nin˜o-Southern Oscillation [Monahan, 2001; Hsieh, 2001; Hsieh and Wu, 2002]. The present paper will apply the nonlinear principal component analysis (NLPCA) to the observed equatorial stratospheric wind record. The rich detail characteristic of each observed QBO cycle, along with the slightly irregular cycleto-cycle nature of the QBO, make this a particularly appealing phenomenon for application of the NLPCA technique. The present analysis will show that the nonlinear approach leads to a significantly more complete characterization of the QBO HAMILTON AND HSIEH: REPRESENTATION OF THE QBO Table 1. Root-Mean-Square Error Between Observations and the Reconstructed Component From NLPCA.cir Mode 1 as a Function of the Number m of Hidden Neuron Layersa m 2 3 4 5 6 7 8 a RMSE, m sÀ1 4.93 4.45 3.96 3.84 3.77 3.73 3.69 Results are averaged for all seven pressure levels. behavior in a single time series, than does the comparable linear analysis. 2. Data [8] The data analyzed here were obtained from B. Naujokat of the Free University of Berlin and are monthly means of the zonal wind component measured by twice-per-day ACL 3-3 balloon ascents at: Canton Island (2.8°N) during January 1956 through August 1967, Gan (0.7°S) from September 1967 through December 1975, and Singapore (1.4°N) from January 1976 through December 2000. This represents a continuation of the data record analyzed by Naujokat [1986], Wallace et al. [1993], Fraedrich et al. [1993], Wang et al. [1995], and many other studies. Values at 70, 50, 40, 30, 20, 15, and 10 hPa (i.e., from about 20 to 30 km altitude) are used. A plot of a slightly smoothed version of the monthly data (with 45-year mean value at each level removed) for 20 years is shown in Figure 7b below. 3. Brief Description of the NLPCA.cir Method [9] In principal component analysis (PCA), also known as empirical orthogonal function analysis, a straight-line approximation to the dataset is sought which accounts for the maximum amount of variance in the data. There are now several types of neural network (NN) models which use, instead of the straight line, a continuous curve to approximate Figure 2. The NLPCA.cir mode 1 solution is shown by the (overlapping) circles, while the data are shown as dots. For comparison, the PCA mode 1 solution is shown as a thin straight line. Only three out of seven dimensions are shown, namely, u at the top, middle, and bottom levels (10, 30 and 70 hPa). (a) – (c) The 2-D views; (d) a 3-D view. ACL 3-4 HAMILTON AND HSIEH: REPRESENTATION OF THE QBO Figure 3. The mode 1 nonlinear principal component q time series for (a) m = 2, (b) m = 4, (c) m = 6 (the chosen solution), and (d) m = 8, where m is the number of hidden neurons in the encoding layer. Note q (plotted in p radians) is a cyclic variable bounded between [À1, 1]. The values of q in Figures 3a, 3b, and 3d have each been adjusted by a constant value, so all q have the same starting value. the data. Using a multi-layer perception NN formulation, the nonlinear PCA (NLPCA) method of Kramer [1991] is capable of finding an open curve to approximate the data. However, for periodic or quasiperiodic phenomena, it would be more appropriate to fit the data with a closed curve. Kirby and Miranda [1996] introduced an NLPCA with a circular node at the bottleneck (henceforth NLPCA.cir), capable of extracting closed curve solutions. Applications to the tropical Pacific climate variability have been made with the NLPCA by Monahan [2001] and Hsieh [2001] and with the NLPCA.cir by Hsieh [2001]. [10] The NLPCA.cir used here is described in detail by Hsieh [2001]. The input data are in the form x(t) = [x1,. . ., xl], where each variable xi, (i = 1,. . ., l ), is a time series containing n observations. The information is mapped forward through a bottleneck to the output x0 (Figure 1). The parameters of the network are solved by minimizing the cost function, which is basically the mean square error (MSE) of x0 relative to x. Because of local minima in the cost function, an ensemble of 30 NNs with random initial weights and bias parameters was run. Also, 20% of the data was randomly selected as test data and withheld from the training of the NNs. Runs where the MSE was larger for the test dataset than for the training dataset were rejected to avoid overfitted solutions. Then the NN with the smallest MSE was selected as the solution. As noted by Hsieh [2001], the NLPCA.cir can be set up in a general or a restricted configuration, where in the general configuration, the solution can be either an open or a closed curve, while the restricted configuration favors a closed curve. In this study, the general configuration will be used. 4. Analysis of the Equatorial Wind Data by NLPCA.cir [11] For this analysis the wind values at each of the seven vertical levels had the 45-year mean removed (but the weak seasonal cycle was retained). These zonal wind values are then the seven inputs into the NLPCA.cir (Figure 1), with m, the number of hidden neurons in the encoding layer (and in the decoding layer) varying from 2 to 8. No weight penalty [Hsieh, 2001] was needed in the cost function as the signal to noise ratio was very good. Table 1 gives the root-meansquare error (RMSE) as a function of m. Beyond m = 6 the decrease in RMSE obtained by increasing the hidden layers is very slow, so the m = 6 solution is chosen as the most appropriate one. [12] The NLPCA.cir mode 1 solution gives a closed curve in a seven-dimensional space. The system goes around the closed curve once, as the NLPC q varies through one cycle. Figure 2 shows the solution in three of the seven dimensions, namely, the wind anomalies at 10, 30, and 70 hPa. The HAMILTON AND HSIEH: REPRESENTATION OF THE QBO Figure 4. Contour plot of the NLPCA.cir mode 1 zonal wind anomalies as a function of pressure and q. Contour interval is 5 m sÀ1, with positive contours (i.e., westerly or eastward winds) indicated by solid lines, negative contours indicated by dashed lines, and zero contours indicated by thick lines. NLPCA.cir mode 1 explains 94.8% of the variance. For comparison, the linear PCA yields seven modes explaining 57.8, 35.4, 3.1, 2.1, 0.8, 0.5, and 0.3% of the variance, respectively. Thus the variance explained by the NLPCA.cir mode 1 (94.8%) exceeds the variance explained by PCA 1 and 2 combined (93.2%). The PCA mode 1 is also shown as a straight line in Figure 2. The MSE of NLPCA.cir mode 1 is only 0.12 that of the PCA mode 1, indicating a much better fit to the data by the nonlinear mode. [13] To provide the most appropriate linear approach for comparison with the NLPCA.cir, the time series for the first two linear PCAs were divided by their standard deviation and a phase angle then defined as the arctangent of the ratio of the two normalized PCA coefficients. This was then used as the basis for a one-parameter linear ‘‘circular’’ approximation. This linear model accounts for 83.0% of the variance in the zonal wind time series, versus the 94.8% from the first NLPCA.cir mode. [14] The time series for the nonlinear principal component (NLPC) q is shown in Figure 3 for m = 2, 4, 6, and 8. The results are rather similar for all these values of m, demonstrating the stability of the NLPC approach in this case. The NLPC completes slightly more than 19 cycles over the 45 years considered, for an average period of about 28 months. The NLPCA.cir representation of a single QBO cycle as a function of q is given in Figure 4. However, Figure 3 shows that q does not advance uniformly with time through the QBO cycle (note, for example, the conspicuous slowing down of q progression around 0.7– 0.9). To quantify this effect, all the monthly q values for 1956 –2000 were binned into 20 equal q intervals. The resulting histogram (not shown) has a variation of between 3 and 58 monthly values within the individual bins, confirming the impression in Figure 3 that q consistently advances through some values quickly and through other values slowly. After applying a three-point running mean to the histogram function, the QBO cycle of Figure 4 was replotted with respect to qweighted, which is q weighted by the smoothed histogram distribution (Figure 5). The purpose of this weighting is to produce an independent variable that equally ACL 3-5 Figure 5. Contour plot of the NLPCA.cir mode 1 zonal wind anomalies as a function of pressure and qweighted, where qweighted is more representative of actual time during a cycle than q. Contour interval is 5 m sÀ1. distributes the time the QBO spends in any range. This gives a more accurate picture of how the QBO typically progresses with time over one cycle, and can be regarded as an objectively determined ‘‘composite’’ QBO cycle using the NLPCA.cir approach. Earlier attempts to construct composite QBO cycles [e.g., Naujokat, 1986; Baldwin et al., 2001] were made with the standard ‘‘superposed epoch’’ approach. These earlier composites tended to smooth out some of the features seen consistently in every QBO cycle, notably the strength of the shear zones, and the asymmetry in phase progression and shear strength between the descending easterly and westerly regimes. The present NLPCA.cir ‘‘composite’’ in Figure 5 does a much better job of capturing these subtleties, including the tendency for the descending easterly phase to ‘‘stall’’ between 30 and 50 hPa. The equivalent of Figure 5 was computed using the one-parameter linear circular approximation and is shown here as Figure 6. This is much less able than the NLPCA.cir treatment to capture the characteristic irregularities in QBO wind regime propagation. Figure 6. As in Figure 5, but for the equivalent analysis using a linear circular approach. 3-6 ACL HAMILTON AND HSIEH: REPRESENTATION OF THE QBO Table 2. Correlation and RMSE Between Observations and the Reconstructed Componenta RMSE, m sÀ1 Correlation Pressure, hPa 10 15 20 30 40 50 70 Averageb a b NLRC1 0.966 0.978 0.983 0.985 0.979 0.958 0.851 0.957 RC1 0.840 0.983 0.924 0.512 0.086 0.469 0.780 0.656 RC1+2 0.953 0.984 0.977 0.972 0.973 0.934 0.822 0.945 NLRC1 4.91 4.12 3.56 3.07 3.26 3.76 3.44 3.73 RC1 10.27 3.62 7.47 15.27 15.76 11.59 4.10 9.73 RC1+2 5.75 3.54 4.20 4.17 3.64 4.71 3.73 4.25 From NLPCA.cir mode 1 (NLRC1), from PCA mode 1 (RC1), and from PCA modes 1 and 2 (RC1+2). Average over the seven vertical levels. [15] The NLPCA.cir reconstructed component (NLRC1) was calculated; this represents the approximation of the original seven time series (at the seven vertical levels) by the single NLPCA.cir. The reconstructed time series from linear PCA mode 1 (RC1) and from PCA modes 1 and 2 (RC1+2) were also calculated. The correlation and rootmean-square error (RMSE) between these reconstructed time series and the observations are given in Table 2. The NLRC1 and RC1+2 are plotted for the period 1981 – 2000 in Figure 7, together with the raw observations. Reconstructions using additional linear PCA modes will, of course, fit the observations better, but will also be fitting to the noise in the data. Here RC1 + 2 is already more noisy in appearance than NLRC1. The single RC1 is less Figure 7. (a) The reconstructed component (RC) from NLPCA.cir mode 1, (b) the observed zonal wind anomalies, and (c) the RC from the PCA modes 1 and 2. Only the period 1981 – 2000 is shown. The observed anomalies in Figure 7b were smoothed by a 3-month running mean for better legibility. Negative contours are dashed, while the zero contour is thickened. Contour intervals are 10 m sÀ1. HAMILTON AND HSIEH: REPRESENTATION OF THE QBO Figure 8. The NLPCA.cir phase, qweighted, in each month of the year for each year from 1956 to 1970. The phases are normalized by p and are plotted modulo 2p in a way to show the phase progression continuously each year. noisy than RC1+2, but models the actual data rather poorly, particularly around 30 – 50 hPa (Table 2). The single parameter NLPCA.cir reconstruction is better correlated with the observations and has a lower MSE than the two-mode RC1+2 reconstruction. Inspection of Figure 7 shows that the NLRC1 reconstruction captures much better the strong shear zones and rapid accelerations apparent in the raw data than does the RC1+2. ACL 3-7 [16] A number of earlier studies have addressed the possible role of the annual cycle in the dynamics of the QBO [e.g., Dunkerton, 1983b, 1990]. Wallace et al. [1993] found that their QBO phase angle defined from the first two linear PCAs, on average, advanced most rapidly in April – May and least rapidly in December – January. Figure 8 shows the present NLPCA.cir qweighted plotted each month for the first 15 consecutive years of the record. The tendency for the phase to advance most rapidly in boreal spring is evident in almost all individual years. Figure 9 shows the phase difference between successive months averaged over all years of data for each pair of adjacent calendar months. On average the phase varies almost 3 times as much between May and June as between September and October. This shows a clear dependence of the QBO behavior on season. The rapid phase variation in boreal spring found here is consistent with the earlier study of Wallace et al. However, Figure 9 has a more semiannual nature than the comparable result from Wallace et al., with a secondary maximum appearing in Figure 9 in November– December. [17] The possibility that the QBO period may be subject to some systematic low-frequency modulation was studied originally by Quiroz [1981] and Dunkerton and Delisi [1985]. As noted in section 1, the recent work of Salby and Callaghan [2000], has revived interest in this possibility. Earlier studies have all used a simple characterization of QBO period in terms of the zero crossings of the observed wind at a particular level (generally near 50 hPa). The present objective determination of the phase of the QBO, using data from all the levels, allows this issue to be revisited. A ‘‘local’’ period of the QBO was defined for Figure 9. The average QBO phase progression between successive months for January-to-February (J– F), February-to-March (F –M), through December-to-January (D – J). ACL 3-8 HAMILTON AND HSIEH: REPRESENTATION OF THE QBO Figure 10. Period of the QBO as a function of time. The period was estimated from 2p hÁq/Áti, where angle brackets denote averaging over 28 months. each month proportional to the reciprocal of the rate of phase change, ÁqweightedÁt. The monthly values defined this way are rather noisy, so Figure 10 shows the reciprocal of the 28-month running mean of ÁqweightedÁt. There does seem to be a systematic modulation with the longest periods occurring in the mid to late 1960s, and the mid to late 1980s, while the shortest periods occurred around 1960, the early 1970s and the late 1990s. The reason for such systematic modulation of the QBO period is not known, although possibilities that have been suggested include: variations in volcanic aerosol loading or concentrations of other radiatively active constituents in the stratosphere [Dunkerton, 1983a], variations in tropospheric climate leading to modulation of the waves thought to force the QBO wind accelerations [Geller et al., 1997], or solar cycle variations [Salby and Callaghan, 2000]. The present results in Figure 10 do not seem to show any clear connection with the 11-year solar cycle. The results in the present Figure 10 differ to some extent from earlier time series of QBO periods deduced from simple inspection of zero crossings of the time series at 30 hPa [Geller et al., 1997] or 50 hPa [Salby and Callaghan, 2000]. Features in common in all three studies include a clear lengthening of the period for a 8 – 10 year era centered around 1966 and an era of very short periods in the early 1970s. The analyzed variations in QBO period during the 1980s and 1990s are somewhat different among the three studies, with the 50 hPa estimates of Salby and Callaghan showing a strong maximum centered around 1986 and a deep minimum in the early 1990s, which do not have clear counterparts in the present results or those of Geller et al. 5. Application of NLPCA.cir-Based QBO Phase to Study Extratropical Stratospheric Effects [18] The phase of the tropical QBO is known to be correlated with the extratropical stratospheric circulation, a fact first demonstrated by Holton and Tan [1980, 1982] and now often referred to as the ‘‘Holton-Tan (HT) effect.’’ In particular, HT and later investigators have shown that the polar stratosphere during NH winter is on average warmer, and the westerly vortex circulation is weaker, when the zonal-mean zonal winds at the equator near 40 hPa are easterly than when they are westerly. The HT effect is strongest in November –January and is virtually absent later in the winter. This effect has a plausible physical explanation in terms of the expected modulation of quasi-stationary planetary wave propagation by the mean flow and has been reproduced in various models of the stratospheric circulation including sophisticated GCMs [Hamilton, 1998b; Takahashi, 1999]. [19] The HT effect has generally been studied using the equatorial wind at a single level near 40 hPa to characterize the QBO phase (the 40 hPa level is chosen as this generally is found to maximize correlations with the indicators of extratropical circulation). Physically there is no reason to expect that the winds at any particular level will control the planetary wave modulation. Thus it is of interest to revisit the HT effect with the present NLPCA.cir characterization of the tropical QBO phase. [20] The January-mean North Pole (NP) 30 hPa temperatures from the NCEP/NCAR reanalyses for the 42-year period 1958– 1999 were used to characterize the extratropical winter stratosphere. Then 42 separate (but overlapping) q ranges were defined so that each contained exactly 21 of the q values from the NLPCA.cir analysis for these 42 Januaries. For each of these q ranges the quantity ÁT was computed as the NP temperatures averaged over the 21 Januaries within the range minus that averaged over the 21 Januaries out of the range. The largest ÁT of 5.4°C was obtained for a q range between À0.86675 and 0.44775. This q range then is the analogue of the ‘‘easterly QBO phase’’ as defined using the raw wind value at one level by HT and others. Inspection of Figure 4 shows that this range of q corresponds almost exactly to the period of anomalously easterly wind at 40 hPa in the composite QBO cycle. [21] Table 3 shows the phases for each January defined by this NLPCA.cir criterion. The second column shows the anomaly in the raw equatorial wind time series at 40 hPa. There is general agreement between the two measures of QBO phase, but there are some years where there are differences. For example, January 1968 emerges as ‘‘easterly’’ by the NLPCA.cir definition, but the 40 hPa winds are anomalously westerly, and similarly in January 1995 the 40 hPa winds are westerly, but the NLPCA.cir definition is ‘‘easterly.’’ If the ÁT is computed using the average NP temperatures in the 21 Januaries with most easterly 40 hPa wind minus the average during the 21 most westerly months (essentially following the original HT approach), then a value of 4.0°C is obtained. This suggests that the NLPCA. cir criterion may be better able to characterize the aspect of the tropical QBO that correlates with extratropical NH winter circulation. Another approach based on the linear PCs was also investigated. An angle was defined as the arctangent of the ratio of the amplitude coefficients of the first two linear PCs. The analysis described for the NLPCA.cir q and NP temperature was repeated using this angle. The largest ÁT obtained is 4.2°C. The phase characterization for each January using this procedure is shown in the final column of Table 3. It turns out that this linear PC approach yields identical results to the NLPCA.cir approach in 40 of the 42 years, the exceptions being 1995 and 1996. The situation in January 1995 is particularly HAMILTON AND HSIEH: REPRESENTATION OF THE QBO Table 3. Indications of the Phase of the Tropical QBO in Each January Based on the NLPCA and on the Linear PCAa Year NLPCA Phase 40 hPa u, m sÀ1 PCA Phase 1956 1957 1958 1959 1960 1961 1962 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 W E W E W E W E W W E W E E W E W E W E W E W W E W E W E E W E W W E W E W E E W E W E W 8.3 À24.3 13.6 À19.4 5.6 À20.6 13.2 À19.7 6.6 13.8 À27.5 17.1 6.6 À26.8 11.8 À26.6 15.4 À17.9 13.1 À18.4 16.4 À10.4 14.6 13.6 À21.9 16.5 À12.4 18.1 À6.3 À4.0 17.8 À6.0 11.7 14.4 À26.7 12.6 À13.4 16.0 À0.4 15.9 5.1 À21.3 16.4 À4.3 15.3 W E W E W E W E W W E W E E W E W E W E W E W W E W E W E E W E W W E W E W E W E E W E W a Also given is the anomaly in the January-mean 40 hPa zonal wind, u. W means westerly (i.e., eastward or positive) zonal wind phase. Results shown for 1956 – 2000, although the analysis to define the E and W phase ranges in the NLPCA and linear PCA were defined on the basis of 1958 – 1999 data only. interesting. The linear PC approach classifies this as westerly and the January-mean 40 hPa winds are in fact westerly. However, the 40 hPa winds in January represent something of a one month spike, with the value dropping over February and March, and the sustained westerly phase at 40 hPa beginning in April (see Figure 7b). The NLPCA criterion classifies January 1995 as the end of an easterly phase while the linear PC approach classifies this month as the beginning of a westerly phase. As noted earlier, the NLPCA.cir phase determination for these two Januaries is ‘‘better’’ in the sense that it leads to an enhanced contrast in the average easterly versus westerly NP temperature (5.2°C versus 4.2°C). The statistical significance of the temperature differences was evaluated by comparing to 1000 values computed based on randomly selected sets of 21 years. The 5.2°C threshold was exceeded by only 9 of 1000 of these ACL 3-9 random cases, while the 4.2°C threshold was exceeded by 37 cases. So both approaches lead to results that are 95% significant. 6. Conclusion [22] Many geophysical data sets take the form of time series observations simultaneously at several spatial locations. In such cases it is often desirable to construct lowdimensional approximate representations of the full data set, hence the immense popularity of linear PC analysis. Ideally, such an approximation should fit the individual data values accurately, and in such a way that the basic phenomena are preserved while what may be regarded as noise is minimized. Linear PC analysis is a computationally and conceptually straightforward approach, but it represents just one of many possible approximations that can be used. The results in this paper have shown that application of the NNbased, circular bottleneck NLPCA.cir approach described by Hsieh [2001] produces a one-dimensional approximation to the observed height-time variability of the equatorial stratospheric zonal winds with many desirable features. A single NLPC fits the data more accurately in terms of MSE than a comparable one-parameter linear fit or even the sum of the first two linear PCs. The reconstructed wind field with the single NLPC also shows less apparent noise than the linear PC reconstruction, and better captures some of the important characteristics of the observed QBO wind variations. The NLPC was used to construct a ‘‘typical’’ QBO cycle in the winds (Figure 5) which captures such features as the very strong shear zones and westerly vs. easterly asymmetries better than comparable ‘‘composite’’ QBO cycles that have been constructed using linear approaches. The present typical QBO cycle determined from observations could serve as a useful comparison for a similar analysis of model simulated winds. The NLPC can also be used to study the time-varying modulation in QBO period in a way that uses the observations from all levels simultaneously, in contrast to earlier single-level approaches. Finally, the NLPCA.cir can also serve as the basis to study how the equatorial QBO phase correlates with other aspects of the atmospheric circulation. In this paper it was shown that a stratification of NP 30 hPa January temperature data based on the NLPCA.cir-determined QBO phase resulted in isolation of a stronger QBO-related effect than that obtained using characterization of the QBO phase by either the single level or linear PCA-based approaches. [23] The stratospheric QBO is a good candidate for application of the NLPCA.cir approach, as it is clearly a quasiperiodic variation with subtle characteristic features in each cycle, but one with considerable irregularity in its phase progression. The NLPCA.cir as described here would be less useful for quasiperiodic variations with large modulations in amplitude, and the method would need to be modified to deal with such a case. The NLPCA.cir can be extended in other ways as well. Hsieh and Wu [2002] have developed the NLPCA analogue of linear singular spectrum analysis (also referred to as extended empirical orthogonal function analysis), and the present authors have now applied this technique to analysis of the stratospheric QBO (W.W. Hsieh and K. Hamilton, Nonlinear singular spectrum anal- ACL 3 - 10 HAMILTON AND HSIEH: REPRESENTATION OF THE QBO ysis of the tropical stratospheric wind, submitted to Quarterly Journal of the Royal Meteorological Society, 2001). [24] Acknowledgments. W. Hsieh was supported by research and strategic grants from the Natural Sciences and Engineering Research Council of Canada. The International Pacific Research Center is supported in part by the Frontier Research System for Global Change. The authors thank Barbara Naujokat for supplying the wind data. The authors also acknowledge the helpful comments of the two official reviewers. References Baldwin, M., and T. J. Dunkerton, Quasi-biennial modulation of the Southern Hemisphere stratospheric polar vortex, Geophys. Res. Lett., 25, 3343 – 3346, 1998. Baldwin, M., et al., The quasi-biennial oscillation, Rev. Geophys., 39, 179 – 229, 2001. Coy, L., An unusually large westerly amplitude of the quasi-biennial oscillation, J. Atmos. Sci., 36, 174 – 176, 1979. Dunkerton, T. J., Modification of stratospheric circulation by trace constituent changes, J. Geophys. Res., 88, 10,831 – 10,836, 1983a. Dunkerton, T. J., Laterally-propagating Rossby waves in the easterly acceleration phase of the quasi-biennial oscillation, Atmos. Ocean, 21, 55 – 68, 1983b. Dunkerton, T. J., Annual variation of deseasonalized mean flow acceleration in the equatorial lower stratosphere, J. Meteorol. Soc. Jpn., 68, 499 – 508, 1990. Dunkerton, T. J., and M. P. Baldwin, Quasi-biennial modulation of planetary-wave fluxes in the Northern Hemisphere winter, J. Atmos. Sci., 48, 1043 – 1061, 1991. Dunkerton, T. J., and D. Delisi, Climatology of the equatorial lower stratosphere, J. Atmos. Sci., 42, 1199 – 1208, 1985. Fraedrich, K., S. Pawson, and R. Wang, An EOF analysis of the verticaltime delay structure of the quasi-biennial oscillation, J. Atmos. Sci., 50, 3357 – 3365, 1993. Geller, M. A., W. Shen, M. Zhang, and W.-W. Tan, Calculations of the stratospheric quasi-biennial oscillation for time-varying wave forcing, J. Atmos. Sci., 54, 883 – 894, 1997. Gray, W. M., Atlantic seasonal hurricane frequency, I, El Nin˜o and 30 mb quasi-biennial oscillation influences, Mon. Weather Rev., 112, 1649 – 1668, 1984. Gray, W. M., J. D. Scheaffer, and J. A. Knaff, Hypothesized mechanism for stratospheric QBO influence on ENSO variability, Geophys. Res. Lett., 19, 107 – 110, 1992. Hamilton, K., Dynamics of the tropical middle atmosphere: A tutorial review, Atmos. Ocean, 36, 319 – 354, 1998a. Hamilton, K., An imposed quasi-biennial oscillation in a comprehensive general circulation model: Response of the tropical and extratropical circulation, J. Atmos. Sci., 55, 2393 – 2418, 1998b. Hamilton, K., R. J. Wilson, and R. Hemler, Middle atmosphere simulated with high vertical and horizontal resolution versions of a GCM: Improvement in the cold pole bias and generation of a QBO-like oscillation in the tropics, J. Atmos. Sci., 56, 3829 – 3846, 1999. Hamilton, K., R. J. Wilson, and R. Hemler, Spontaneous stratospheric QBO-like oscillations simulated by the GFDL SKYHI general circulation model, J. Atmos. Sci., 58, 3271 – 3292, 2001. Holton, J. R., and H.-C. Tan, The influence of the equatorial quasi-biennial oscillation on the global circulation at 50 mb, J. Atmos. Sci., 37, 2200 – 2208, 1980. Holton, J. R., and H.-C. Tan, The quasi-biennial oscillation in the Northern Hemisphere lower stratosphere, J. Meteorol. Soc. Jpn., 60, 140 – 148, 1982. Horinouchi, T., and S. Yoden, Wave-mean flow interaction associated with a QBO-like oscillation simulated in a simplified GCM, J. Atmos. Sci., 55, 502 – 526, 1998. Hsieh, W. W., Nonlinear principal component analysis by neural networks, Tellus, Ser. A, 53, 599 – 615, 2001. Hsieh, W. W., and A. Wu, Nonlinear multichannel singular spectrum analysis of the tropical Pacific climate variability using a neural network approach, J. Geophys. Res., 107(C7), 10.1029/2001JC000957, 2002. Kirby, M. J., and R. Miranda, Circular nodes in neural networks, Neural Comput., 8, 390 – 402, 1996. Kramer, M. A., Nonlinear principal component analysis using autoassociative neural networks, AIChE J., 37, 233 – 243, 1991. Marquardt, C., and B. Naujokat, An update of the equatorial QBO and its variability, Tech. Doc. 814, pp. 87 – 90, World Meteorol. Org., Geneva, 1997. Monahan, A. H., Nonlinear principal component analysis: Tropical IndoPacific sea surface temperature and sea level pressure, J. Clim., 14, 219 – 233, 2001. Naujokat, B., An update of the observed quasi-biennial oscillation of the stratospheric winds over the tropics, J. Atmos. Sci., 43, 1873 – 1877, 1986. Quiroz, R. S., Period modulation of the stratospheric quasi-biennial oscillation, Mon. Weather Rev., 109, 665 – 674, 1981. Randel, W. J., J. F. Wu, J. M. Russell III, A. Roche, and J. Waters, Seasonal cycles and QBO variations in stratospheric methane and water vapor observed in UARS HALOE data, J. Atmos. Sci., 55, 163 – 185, 1998. Salby, M., and P. Callaghan, Connection between the Solar Cycle and the QBO: The missing link, J. Clim., 13, 2652 – 2662, 2000. Takahashi, M., Simulation of the stratospheric quasi-biennial oscillation using a general circulation model, Geophys. Res. Lett., 23, 661 – 664, 1996. Takahashi, M., The first realistic simulation of the stratospheric quasi-biennial oscillation in a general circulation model, Geophys. Res. Lett., 26, 1307 – 1310, 1999. Tung, K.-K., and H. Yang, Global QBO in circulation and ozone, I, Reexamination of observational evidence, J. Atmos. Sci., 51, 2699 – 2707, 1994. Wallace, J. M., The general circulation of the tropical lower stratosphere, Rev. Geophys., 11, 191 – 222, 1973. Wallace, J. M., L. Panetta, and J. Estberg, A phase-space representation of the equatorial stratospheric quasi-biennial oscillation, J. Atmos. Sci., 50, 1751 – 1762, 1993. Wang, R., K. Fraedrich, and S. Pawson, Phase-space characteristics of tropical stratospheric quasi-biennial oscillation, J. Atmos. Sci., 52, 4482 – 5000, 1995. Xu, J.-S., On the relationship between the stratospheric quasi-biennial oscillation and the tropospheric Southern Oscillation, J. Atmos. Sci., 49, 725 – 734, 1992. ÀÀÀÀÀÀÀÀÀÀÀ K. Hamilton, Department of Meteorology and International Pacific Research Center, University of Hawaii at Manoa, Honolulu, HI 96822, USA. (kph@soest.hawaii.edu) W. W. Hsieh, Department of Earth and Ocean Sciences, University of British Columbia, Vancouver, B.C., Canada V6T 1Z4.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Faculty Research and Publications /
- Representation of the quasi-biennial oscillation in...
Open Collections
UBC Faculty Research and Publications
Representation of the quasi-biennial oscillation in the tropical stratospheric wind by nonlinear principal… Hamilton, Kevin; Hsieh, William W. Jul 3, 2002
pdf
Page Metadata
Item Metadata
Title | Representation of the quasi-biennial oscillation in the tropical stratospheric wind by nonlinear principal component analysis |
Creator |
Hamilton, Kevin Hsieh, William W. |
Publisher | American Geophysical Union |
Date Issued | 2002-07-03 |
Description | The zonal winds at several levels between 70 and 10 hPa (roughly 20–30 km) measured at near-equatorial stations during 1956–2000 were analyzed to produce a one-dimensional approximation. The neural network-based technique applied was the circular nonlinear principal component analysis (NLPCA.cir) designed to characterize quasiperiodic phenomena. The reconstructed height-time series of wind based on the one-dimensional NLPCA.cir captures many of the characteristic features of the observed quasi-biennial oscillation (QBO). The nonlinear results were evaluated relative to comparable linear principal component analysis and found to produce a superior one-dimensional representation of the data. The NLPCA.cir analysis produces a single time series of QBO phase based on data at all levels. This phase was then applied to demonstrate a strong correlation of the state of the tropical QBO and boreal winter high-latitude stratospheric temperatures. An edited version of this paper was published by AGU. Copyright 2002 American Geophysical Union. |
Genre |
Article |
Type |
Text |
Language | eng |
Date Available | 2016-11-10 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivatives 4.0 International |
DOI | 10.14288/1.0041789 |
URI | http://hdl.handle.net/2429/32834 |
Affiliation |
Science, Faculty of Earth and Ocean Sciences, Department of |
Citation | Hamilton, Kevin, Hsieh, William W. 2002. Representation of the quasi-biennial oscillation in the tropical stratospheric wind by nonlinear principal component analysis. Journal of Geophysical Research Oceans. 107 D15 4232 |
Publisher DOI | 10.1029/2001JD001250 |
Peer Review Status | Reviewed |
Scholarly Level | Faculty |
Copyright Holder | Hsieh, Wiliam W. |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/4.0/ |
AggregatedSourceRepository | DSpace |
Download
- Media
- 52383-Hsieh_AGU_2002_2001JD001250.pdf [ 485.12kB ]
- Metadata
- JSON: 52383-1.0041789.json
- JSON-LD: 52383-1.0041789-ld.json
- RDF/XML (Pretty): 52383-1.0041789-rdf.xml
- RDF/JSON: 52383-1.0041789-rdf.json
- Turtle: 52383-1.0041789-turtle.txt
- N-Triples: 52383-1.0041789-rdf-ntriples.txt
- Original Record: 52383-1.0041789-source.json
- Full Text
- 52383-1.0041789-fulltext.txt
- Citation
- 52383-1.0041789.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.52383.1-0041789/manifest