{"@context":{"@language":"en","Affiliation":"http:\/\/vivoweb.org\/ontology\/core#departmentOrSchool","AggregatedSourceRepository":"http:\/\/www.europeana.eu\/schemas\/edm\/dataProvider","Campus":"https:\/\/open.library.ubc.ca\/terms#degreeCampus","Creator":"http:\/\/purl.org\/dc\/terms\/creator","DateAvailable":"http:\/\/purl.org\/dc\/terms\/issued","DateIssued":"http:\/\/purl.org\/dc\/terms\/issued","Degree":"http:\/\/vivoweb.org\/ontology\/core#relatedDegree","DegreeGrantor":"https:\/\/open.library.ubc.ca\/terms#degreeGrantor","Description":"http:\/\/purl.org\/dc\/terms\/description","DigitalResourceOriginalRecord":"http:\/\/www.europeana.eu\/schemas\/edm\/aggregatedCHO","FullText":"http:\/\/www.w3.org\/2009\/08\/skos-reference\/skos.html#note","Genre":"http:\/\/www.europeana.eu\/schemas\/edm\/hasType","IsShownAt":"http:\/\/www.europeana.eu\/schemas\/edm\/isShownAt","Language":"http:\/\/purl.org\/dc\/terms\/language","Program":"https:\/\/open.library.ubc.ca\/terms#degreeDiscipline","Provider":"http:\/\/www.europeana.eu\/schemas\/edm\/provider","Publisher":"http:\/\/purl.org\/dc\/terms\/publisher","Rights":"http:\/\/purl.org\/dc\/terms\/rights","ScholarlyLevel":"https:\/\/open.library.ubc.ca\/terms#scholarLevel","Title":"http:\/\/purl.org\/dc\/terms\/title","Type":"http:\/\/purl.org\/dc\/terms\/type","URI":"https:\/\/open.library.ubc.ca\/terms#identifierURI","SortDate":"http:\/\/purl.org\/dc\/terms\/date"},"Affiliation":[{"@value":"Science, Faculty of","@language":"en"},{"@value":"Earth, Ocean and Atmospheric Sciences, Department of","@language":"en"}],"AggregatedSourceRepository":[{"@value":"DSpace","@language":"en"}],"Campus":[{"@value":"UBCV","@language":"en"}],"Creator":[{"@value":"Ye, Zhengqing","@language":"en"}],"DateAvailable":[{"@value":"2011-02-10T20:50:01Z","@language":"en"}],"DateIssued":[{"@value":"2007","@language":"en"}],"Degree":[{"@value":"Doctor of Philosophy - PhD","@language":"en"}],"DegreeGrantor":[{"@value":"University of British Columbia","@language":"en"}],"Description":[{"@value":"The observed El Nin\u0303o\/Southern Oscillation (ENSO) has shown eastward displacements of the zonal wind stress anomalies and surface heat flux anomalies in the equatorial Pacific during El Nin\u0303o episodes in the 1981-1995 regime relative to the 1961-1975 regime. Numerical experiments with a modified Zebiak-Cane intermediate complexity coupled model generally reproduced such displacements when the model climatological fields were replaced by the observed and simulated climatologies for the 1981-1995 regime. The later regime climatology resulted in eastward shifts in the ENSO system during El Nin\u0303o but not La Nin\u0303a, through the eastward shift of the atmosphere convergence heating rate in the coupled model. The ENSO period and ENSO predictability were also enhanced in the coupled model under the later regime climatology. A similar behavior was found in the Lorenz (1963) 3-component chaos system, i.e. in both the ENSO and Lorenz systems, stronger nonlinearity appears to lengthen the period thereby enhancing the predictability. Using nonlinear principal component analysis, we demonstrated that the leading ENSO mode had changed on a physical basis since the late 1970s. The ENSO modes resembled the destabilized 'ocean basin mode' and 'the recharge-mode' in the pre- and post-1980s regimes, respectively. The surface zonal current acted as an intensifier of ENSO in the earlier regime and played a role in the transition of the ENSO cycle in the later regime. With data from an ensemble of coupled general circulation models in the IPCC-AR4 project, climate under year 2000 greenhouse gas (GHG) level was compared with climate under pre-industrial conditions. In the tropical Pacific, the warming in the mean sea surface temperatures (SST) was found to have an El Nin\u0303o-like pattern, in agreement with the observed SST data (1900-1999). The models showed that both the equatorial zonal overturning circulation and the shallow meridional overturning circulation weakened under increased GHG forcing. For ENSO, the asymmetry in the SST anomalies between El Nin\u0303o and La Nin\u0303a was found to be enhanced under increased GHG, for both the ensemble model data and the observed data. The enhanced asymmetry was mainly caused by the intensified vertical nonlinear dynamic heating. Furthermore, the enhanced GHG simulations showed that the asymmetry between El Nin\u0303o and La Nin\u0303a increased in the zonal wind stress anomalies, the equatorial undercurrent anomalies and the shallow meridional overturning circulation anomalies.","@language":"en"}],"DigitalResourceOriginalRecord":[{"@value":"https:\/\/circle.library.ubc.ca\/rest\/handle\/2429\/31187?expand=metadata","@language":"en"}],"FullText":[{"@value":"Changes in the El Nino-Southern Oscillation under climate regime shift and increased greenhouse gases by ZHENGQING Y E B.Sc, The Nanjing Institute of Meteorology, 1989 M.Sc, The Nanjing Institute of Meteorology, 1992 A THESIS SUBMITTED IN PARTIAL F U L F I L M E N T OF T H E REQUIREMENTS FOR T H E D E G R E E OF DOCTOR OF PHILOSOPHY in T H E F A C U L T Y OF G R A D U A T E STUDIES (Oceanography) The University of British Columbia April, 2007 \u00a9 ZHENGQING Y E , 2007 Abstract The observed E l Nino\/Southern Oscillation (ENSO) has shown eastward displacements of the zonal wind stress anomalies and surface heat flux anomalies in the equatorial Pacific during E l Nino episodes in the 1981-1995 regime relative to the 1961-1975 regime. Numerical experiments with a modified Zebiak-Cane intermediate complexity coupled model generally reproduced such displacements when the model climatological fields were replaced by the observed and simulated climatologies for the 1981-1995 regime. The later regime climatology resulted in eastward shifts in the ENSO system during E l Nino but not La Nina, through the eastward shift of the atmosphere convergence heating rate in the coupled model. The ENSO period and ENSO predictability were also enhanced in the coupled model under the later regime climatology. A similar behavior was found in the Lorenz (1963) 3-component chaos system, i.e. in both the ENSO and Lorenz systems, stronger nonlinearity appears to lengthen the period thereby enhancing the predictability. Using nonlinear principal component analysis, we demonstrated that the leading ENSO mode had changed on a physical basis since the late 1970s. The ENSO modes resembled the destabilized 'ocean basin mode' and 'the recharge-mode' in the pre- and post-1980s regimes, respectively. The surface zonal current acted as an intensifier of ENSO in the earlier regime and played a role in the transition of the ENSO cycle in the later regime. With data from an ensemble of coupled general circulation models in the IPCC-AR4 project, climate under year 2000 greenhouse gas (GHG) level was compared with climate ii Abstract under pre-industrial conditions. In the tropical Pacific, the warming in the mean sea sur-face temperatures (SST) was found to have an E l Nino-like pattern, in agreement with the observed SST data (1900-1999). The models showed that both the equatorial zonal overturning circulation and the shallow meridional overturning circulation weakened under increased GHG forcing. For ENSO, the asymmetry in the SST anomalies between E l Nino and La Nina was found to be enhanced under increased GHG, for both the ensemble model data and the observed data. The enhanced asymmetry was mainly caused by the intensified vertical nonlinear dynamic heating. Furthermore, the enhanced G H G simulations showed that the asymmetry between El Nino and La Nina increased in the zonal wind stress anomalies, the equatorial undercurrent anomalies and the shallow meridional overturning circulation anomalies. iii \\ Table of Contents Abstract ii Table of Contents iv List of Tables vii List of Figures viii Lis t of Acronyms xx Acknowledgements xxi Chapter 1 Introduction 1 1.1 Scientific background: the physics of ENSO 1 1.2 ENSO decadal variability 6 1.3 Interdecadal ENSO nonlinear ity and predictability 7 1.4 Response of ENSO to global warming 8 1.5 Objectives of the thesis 9 1.6 References 11 Chapter 2 The Influence of Climate Regime Shift on E N S O 15 2.1 Introduction 15 iv Table of Contents 2.2 Data and methods 17 2.2.1 Data 17 2.2.2 Nonlinear principal component analysis 18 2.3 Coupled model 20 2.4 Results 28 2.4.1 Observational results 28 2.4.2 Model simulation results 32 2.5 Discussion 51 2.6 References 55 Chapter 3 Changes in the Nonlineari ty and Predictabi l i ty 59 3.1 Introduction 59 3.2 Intermediate coupled results 60 3.3 Lorenz attractor data 66 3.4 Summary and conclusion 75 3.5 References 75 Chapter 4 Changes in the leading E N S O modes: Role of surface zonal current 78 4.1 Introduction 78 4.2 Combined N L P C A of the sea level and current anomalies 79 4.3 Summary 84 4.4 References 86 Chapter 5 Changes in E N S O and associated overturning circulations from enhanced greenhouse gases 89 v Table of Contents 5.1 Introduction 89 5.2 Observed SST 93 5.3 Model data 94 5.4 Climate change simulated 100 5.5 Changes in the ENSO SST and zonal WS 108 5.6 Changes in the ENSO ocean circulation 112 5.7 Diagnostic analysis of the surface temperature equation 116 5.8 Conclusion 121 5.9 References 123 Chapter 6 Summary and Conclusions 128 6.1 Summary and Conclusions 128 6.2 References 132 Appendix A The Zebiak-Cane coupled model 133 A . l Coupled Model Equations 133 A.2 Computational Procedures 138 A. 3 References 144 Appendix B Nonlinear principal component analysis ( N L P C A ) 145 B. l Open curves 145 B.2 Closed curves 152 B.3 References 153 v i List of Tables 2.1 The eastward shift (in degrees longitude) of the SST and zonal WS anomaly centers in the post-shift regime (1981-1995 for A, B, C, D; 1998-2003 for E) relative to the pre-shift regime (1961-1975 for A, B, C, D; 1992-1998 for E) during strong E l Nino and strong La Nina, as determined from N L P C A . In set A, atmospheric and oceanic monthly climatology were prescribed for both regimes, in B only the atmospheric, and in C only the oceanic. See text for set D and E. The quantized nature of the shift values arose from the model atmospheric zonal grid being 5.625\u00b0 39 3.1 Ensemble mean results for Lorenz systems with different parameter pairs. . 70 5.1 The models used in present Chapter, including their configurations near the equator 99 vii List of Figures List of Figures 2.1 A schematic diagram illustrating the neural network(NN) model for perform-ing the N L P C A . The model is a standard feedforward N N (i.e. multi-layer perceptron) (Bishop, 1995), with 3 'hidden' layers of variables or 'neurons' (denoted by circles) sandwiched between the input layer x on the left and the output layer x' on the right. Next to the input layer is the encoding layer, followed by the 'bottleneck' layer with a single neuron u, then the decoding layer, and finally the output layer, i.e. a total of 4 layers of transfer func-tions are needed to map from the inputs to the outputs. A neuron V{ at the ith layer receives its value from the neurons Vj_i in the preceding layer, i.e. Vi = \/i(wj \u2022 Vj_i + b), where Wj is a vector of weight parameters and b a bias parameter, and the transfer functions fa and fa are the hyperbolic tangent functions, while fa and fa are simply the identity functions. Effectively, a nonlinear function u = F(x.) maps from the higher dimension input space to the lower dimension bottleneck space, followed by an inverse transform x' = G(u) 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 = (||x-x'|| 2) (i.e. the mean square error, MSE) is minimized. Through the optimization, the values of the weight and bias parameters are solved. See Appendix B for details 19 ix List of Figures 2.2 Relations between P C I of the T S U D anomalies and (a) P C I and (b) PC2 of the h anomalies, as found by nonlinear regression via neural networks (over-lapping circles), and linear regression (dashed lines), with the data shown as dots. The corresponding relations for PC2 of the T s u b anomalies are shown in (c) and (d) 23 2.3 Observed SST (\u00b0C) anomalies along equator (averaged in 5\u00b0S-5\u00b0N). The con-tour interval is 1 \u00b0C with positive values shaded 25 2.4 Zonal WS (dyn c m - 2 ) and SST (\u00b0C) anomalies along the equator (averaged over 5\u00b0S-5\u00b0N) in the coupled model from the original model (top panels) and the new model with the neural network parameterization (bottom panels). Positive values are shaded 26 2.5 Thermocline depth (10 m) and the T s u b (\u00b0C) anomalies in the coupled model along the equator (averaged over 5\u00b0S-5\u00b0N) from the original model (top pan-els) and the new model (bottom panels). Positive values are shaded 27 2.6 The decadal differences in the (a) SST (\u00b0C) and (b) WS (dyn cm\"2) over the tropical Pacific. Shaded regions indicate 95% significance by the corrected t test (Zwiers and Von Storch, 1995) 29 x List of Figures 2.7 The observed zonal wind stress WSx (dyn c m - 2 ) anomaly patterns from the leading N L P C A mode for the 1961-1975 regime ((a) and (b)) and for the 1981-1995 regime ((c) and (d)), with panels in the left column showing the anomalies when the N L P C u assumes its maximum value (i.e. during strong El Nino), and the right column, the minimum value (strong La Nina). Regions with values greater than 0.1 dyn c m - 2 or less than -0.1 dyn c m - 2 are shaded. \"L\" and \"H\" mark the location of the lowest and highest values, respectively 2.8 The net shortwave radiation flux (W m - 2 ) anomaly pattern ((a) - (d)) and latent heat flux (W m - 2 ) anomaly pattern ((e) - (h)) at the sea surface from the leading N L P C A mode in the 1961-1975 regime (top) and in the 1979-1993 regime (bottom) when the N L P C u assumes its maximum value (strong El Nino) and minimum value (strong La Nina). Regions with values greater than 10 W m~ 2 or less than -10 W m - 2 are shaded. \"L\" and \"H\" mark the location of the lowest and highest values, respectively 2.9 Set A experiments using the climatology of the oceanic and the atmospheric fields for the 1961-1975 regime (top panels (a), (b), (e) and (f)) and for the 1981-1995 regime (bottom panels ((c), (d), (g) and (h)): The SST anomaly pattern ((a)-(d)) and zonal WS anomaly pattern ((e)-(h)) from the leading N L P C A mode when the N L P C u takes its maximum value (strong El Nino) and minimum value (strong La Nina). Regions with anomaly magnitude exceeding 1.0 \u00b0C or 0.2 dyn c m - 2 are shaded. \"L\" and \"H\" mark the location of the lowest and highest values, respectively 2.10 As in Fig. 2.9, but for set B experiment List of Figures 2.11 As in Fig. 2.9, but for set C experiment 42 2.12 As in Fig. 2.9, but for set E experiment 44 2.13 The model convergence heating rate field Qc (m 2s~ 3) from the leading N L P C A mode when the N L P C takes its maximum value (i.e. strong E l Nino, left panels) and its minimum value (i.e. strong La Nina, right panels) when the 1961-1975 (top panels) and 1981-1995 (bottom panels) climatologies were used. Regions with anomaly magnitude exceeding 0.01 m 2 s - 3 are shaded. . 47 2.14 As in Fig. 2.13, but for Qs 48 2.15 As in Fig. 2.9, but without Qc 49 2.16 Predictability of the Nino 3.4 SST anomaly index, as given by the cross-validated correlation between the predicted and actual index in the coupled model using climatologies from the pre-shift and post-shift regimes. Cross-validation was performed by dividing the 150 years into five segments, where for each segment chosen to test the forecast correlation skills, the other four were used to build the forecast models 52 2.17 Spectrum of the Nino 3.4 SST anomaly index in the coupled model using climatologies from the pre-shift regime (solid curve) and post-shift regime (dashed curve) 53 3.1 The SST anomaly (\u00b0C ) pattern from the leading N L P C A mode when the N L P C takes its maximum value (strong E l Nino) and minimum value (strong La Nina) for the last 100-year model output in the experiments using the climatology of the oceanic and the atmospheric fields for the 1961-1975 regime ((a)-(b)) and for the 1981-1995 regime ((c)-(d)). \"L\" and \"H\" mark the location of the lowest and highest values, respectively 62 List of Figures 3.2 Histogram of the first PCs for all 21 members (a) for pre-shift climate state run and (b) for post-shift climate state run 64 3.3 As in Fig. 3.1, but for the ensemble mean 65 3.4 Ensemble mean predictability (top) and MSE (bottom) of the Nino 3.4 SST anomaly index, as given by the cross-validated correlation between the pre-dicted and actual index in the coupled model using climatologies from the 1961-75 regime (solid lines) and the 1981-95 regime (dashed lines). The thick lines are from nonlinear regression, the thin lines, linear regression. Cross-validation was performed by dividing each 100-year data record into five seg-ments, where for each segment chosen to test the forecast correlation skills, the other four were used to build the forecast model. Error bars indicate \u00b11 standard error of the ensemble mean 67 3.5 Spectrum of the Nino 3.4 SST anomaly index in the coupled model using climatologies from the pre-shift regime (solid curve) and post-shift regime (dashed curve) 68 3.6 Normalized time series of x (top panel) and z (bottom panel) in the Lorenz model for (a) Case 1, (b) Case 2, (c) Case 3 and (d) Case 4, where the time series for different cases are vertically shifted by multiples of 5 for legibility. 72 3.7 Ensemble mean power spectrum of the normalized (a) x and (b) z compo-nents. Dash, solid, dash-dot and dotted lines represent Case 1, 2, 3 and 4, respectively 73 xiii List of Figures 3.8 Ensemble mean predictability in terms of correlation skill (top) and MSE (bottom) of z in the Lorenz system, as given by the cross-validated correlation and MSE between the predicted and actual time series. Dash, solid, dash-dot and dotted lines represent Case 1 (a), 2 (b), 3 (c) and 4 (d), respectively, 4.1 The first and second modes of the SST (top), sea level height (middle) and surface zonal current (bottom) anomalies obtained from a Combined Prin-cipal Component Analysis (CPCA). The percentage of overall variance ex-plained is 36% by the first C P C A mode and 14% by the second mode. The percentage of variance explained for each variable is given in the title of each panel 81 4.2 Trajectory plot of the principal components of the first and second PCs for the pre-1980s and the post-1980s, \"x\" and \"*\" mark the mean location in the PC space of the pre-1980s and the post-1980s data, respectively 82 4.3 Phase-longitude section of the SST (top panels), sea level height (center panels), and the surface zonal current (bottom panels) anomalies over the equatorial band (5\u00b0S to 5\u00b0N) associated with the first N L P C A mode, (a) to 5.1 (a) The observed mean SST (in \u00b0C) during 1950-1999, and (b) the mean SST of 1950-1999 minus that of 1900-1949, with the 5% and 10% significance levels based on the i-test (von Storch and Zwiers, 1999, p. 115) shown in with error bars showing \u00b11 standard error 74 (c) for the pre-1980s and (d) to (f) for the post-1980s 85 dark and light gray respectively. 95 xiv List of Figures 5.2 The trend in monthly mean SST anomalies for 1900 to 1999 in degrees Celsius per century (from NOAA ERSST) 96 5.3 Composite maps of SST anomalies (in \u00b0C) during E l Nino (left column) and during La Nina (right column), with (a) and (b) for the 1900-1949 regime, and (c) and (d) for 1950-1999. \"H\" and \"L\" mark the location of the highest and lowest values, respectively, and the 5% and 10% significance levels based on the t-test are shown in gray 97 5.4 Schematic view of upper ocean (0-500 m depth) circulations in the Pacific from 0\u00b0 to 30\u00b0N 102 5.5 The multi-model ensemble mean for the Commit runs, its difference from the mean for the PIcntrl runs, and the standard deviation of the difference: (a) mean SST (\u00b0C), (b) the difference in the mean SST, (c) the standard deviation of the SST difference, (d) the mean zonal WS (0.01 Pa), (e) the difference in the mean zonal WS, and (f) the standard deviation of the WS difference. The 5% and 10% significance levels from the t-test are shown in gray in (b) and (e) 103 5.6 SST mean (\u00b0C) for the Commit runs (left) and its difference from the mean in the PIcntrl runs (right) for 12 individul models 104 5.7 The zonal mean velocity (averaged over 2\u00b0S and 2\u00b0N in cm s _ 1 ) for the Commit runs (left) and its difference from that in the PIcntrl runs (right) for 12 individual models 109 xv List of Figures 5.8 The multi-model ensemble mean of velocity (in cm s _ 1 ) for the Commit runs and its difference from that in the PIcntrl runs: (a) the mean zonal velocity (averaged over 2\u00b0S and 2\u00b0N), (b) the difference in the mean zonal velocity, (c) the mean meridional velocity averaged over 100 - 300 m depth, and (d) the difference in the mean meridional velocity. The 5% and 10% significance levels from the i-test are shown in gray in (b) and (d) I l l 5.9 Multi-model ensemble mean of SST (\u00b0C) and zonal WS (0.01 Pa) composites for E l Nino (left panels) and La Nina (right panels), with (a), (b), (e), (f) for the PIcntrl runs and (c), (d), (g), (h) for the Commit runs. \"H\" and \"L\" mark the location of the highest and lowest values, respectively, and the 5% and 10% significance levels from the i-test are shown in gray 113 5.10 Vertical section along the equator showing the multi-model ensemble mean of zonal velocity anomalies (cm s _ 1 ) (averaged between 2\u00b0S and 2\u00b0N) for the composite E l Nino (left panels) and La Nina (right panels), with (a) and (b) for the PIcntrl runs, and (c) and (d) for the Commit runs. The 5% and 10% significance levels from the i-test are shown in gray 115 5.11 Multi-model ensemble mean of subsurface meridional velocity anomalies (cm s _ 1 ) (averaged over 100-300 m depth) for the composite E l Nino (left panels) and La Nina (right panels), with (a) and (b) for the PIcntrl runs, and (c) and (d) for the Commit runs. The 5% and 10% significance levels from the i-test are shown in gray 117 xvi List of Figures 5.12 Multi-model ensemble mean of vertical temperature gradient (\u00b0C\/50m, a and b), vertical velocity (10~6 m s - 1 , c and d), vertical nonlinear dynamic heating anomalies (\u00b0C month - 1 , e and f) and total nonlinear dynamic heating anomalies (\u00b0C month - 1 , g and h) for the Commit runs (left panels) and the change from the PIcntrl runs (right panels), in the ocean surface-layer. The 5% and 10% significance levels from the t-test are shown in gray 119 5.13 Multi-model ensemble mean of surface-layer nonlinear dynamic heating (\u00b0C month - 1) for the composite E l Nino (left panels) and La Nina (right panels), with (a), (b) for the PIcntrl runs, and (c), (d) for the Commit runs. The 5% and 10% significance levels from the t-test are shown in gray 120 5.14 Multi-model ensemble mean of the (a-d) surface vertical nonlinear dynamic heating (\u00b0C month - 1 ), (e-h) vertical temperature gradient anomalies (\u00b0C\/50m), and (i-1) vertical velocity anomalies (10 - 6 m s - 1 ) during E l Nino and La Nina for the Commit runs (left panels) and the difference between the Commit and PIcntrl runs (right panels). The 5% and 10% significance levels from the t-test are shown in gray 122 A . l Ocean model levels 135 A.2 The SST anomalies at the equator at different model times (0, 2, 4, 18 months) with (a) as =0 day - 1 ; (b) as =(125 days)\" 1 143 xvii List of Figures B . l (a) A schematic diagram of the NN model for calculating the N L P C A . There are 3 layers of hidden neurons sandwiched between the input layer x on the left and the output layer x' on the right. Next to the input layer is the encoding layer, followed by the 'bottleneck' layer (with a single neuron u), which is then followed by the decoding layer. A nonlinear function maps from the higher dimension input space to the 1-dimension bottleneck space, followed by an inverse transform mapping from the bottleneck space back to the original space represented by the outputs, which are to be as close to the inputs as possible by minimizing the cost function J \u2014 (||x \u2014 x'||2). Data compression is achieved by the bottleneck, with the bottleneck neuron giving u, the nonlinear principal component (NLPC). (b) A schematic diagram of the NN model for calculating the N L P C A with a circular node at the bot-tleneck (NLPCA(cir)). Instead of having one bottleneck neuron u, there are now two neurons p and q constrained to lie on a unit circle in the p-q plane, so there is only one free angular variable A, the N L P C . This network is suited for extracting a closed curve solution 146 B.2 First N L P C A mode for the combined SST and WS anomalies in the tropical Pacific (125\u00b0E-70\u00b0W, 30\u00b0S-30\u00b0N) during 1961-1995. The first N L P C A mode is indicated by the (overlapping) small circles, with the input data shown as dots. The input data were the first six PCs from combined PC A of the SST and WS anomalies (with each variable first normalized by its standard deviation). The N L P C A solution is shown projected onto (a) the PC1-PC2 plane, (b) the PC1-PC3 plane, (c) the PC1-PC4 plane, and (d) the PC1-PC5 plane 151 xviii List of Figures B.3 First NLPCA(cir) mode from the combined SST, zonal surface current and sea level height anomalies in the tropical Pacific (120\u00b0E-70\u00b0W, 20\u00b0S-20\u00b0N) during 1958-2001. The first NLPCA(cir) mode is indicated by the (over-lapping) small circles, with the input data shown as dots. The input data were the first six PCs from a combined P C A of the SST, zonal surface cur-rent and sea level height anomalies (with each variable first normalized by its standard deviation). The NLPCA(cir) solution is a closed curve in the six-dimensional PC space. The NLPCA(cir) solution is shown projected onto (a) the PC1-PC2 plane, (b) the PC1-PC3 plane, (c) the PC1-PC4 plane, and (d) the PC1-PC5 plane 154 xix List of Acronyms ENSO E l Nino\/Southern Oscillation EOF Empirical Orthogonal Function ERSST Extended Reconstructed Sea Surface Temperature , GHG greenhouse gas ICOADS International Comprehensive Ocean-Atmosphere Data Set IPCC-AR4 the 4th Assessment Report of Intergovernmental Panel on Climate Change MSE Mean Squared Error NLCCA nonlinear canonical correlation analysis NLPCA Nonlinear Principal Component Analysis NN Neural Network OCEE Oceanic Convergence zone - Eastern Edge PCA Principal Component Analysis PC Principal Component PDF Probobility Density Function SLP Sea Level Pressure SST Sea Surface Temperature WS Wind Stress Acknowledgements I would like to thank my supervisor Prof. William Hsieh for his help and guidance during the course of this work. Thanks are also due to Dr. Soon-Il An, who collaborated much of the work presented in Chapter 4. I am grateful to Dr. Stephen Zebiak for providing the Zebiak-Cane coupled model. As well, I am sincerely thankful for the constructive suggestion and comments from the members of my committee - Dr. Phil Austin and Dr. Rich Pawlowicz. I also acknowledge the modeling groups for providing their data for analysis, the Program for Climate Model Diagnosis and Intercomparison (PCMDI) for collecting and archiving the model output, and the JSC\/CLIVAR Working Group on Coupled Modelling (WGCM) for organizing the model data analysis activity. The multi-model data archive is supported by the Office of Science, U.S. Department of Energy. I thank Dr. Aiming Wu and Dr. Youmin Tang for many fruitful and stimulating discussions which greatly benefit the work. Most of all, I would like to express my love and gratitude to my family. Without their understanding and support, it would be impossible for me to finish the thesis. xxi Chapter 1 Introduction 1.1 Scientific background: the physics of E N S O The E l Nino-Southern Oscillation (ENSO) phenomenon is the most important inter-annual component in the tropical atmosphere-ocean coupled system. It has important global effects on climate (Wallace and Gutzler, 1981). A better understanding of the ENSO physics could therefore have useful societal payoffs in addition to being a scientifically interesting problem. Much progress has been made in ENSO theory and predictability (Neelin et al. 1998; Battisti et al, 1995; Latif et al. 1998) during the past decades. Several hypotheses about the ENSO mechanism have been proposed: Bjerknes hypothesis: The earliest hypothesis about an ENSO mechanism dates back to a theory of Bjerknes (Bjerknes, 1969). The essence of Bejerknes postulate is that the trade winds and zonal sea surface temperature (SST) affect each other in the tropical Pacific. A change towards a steeper (flatter) zonal pressure gradient strengthens (slackens) the Walker Circulation which results in an increase (decrease) of equatorial easterly winds. An increase (decrease) of equatorial easterly winds produces an increase (decrease) in upwelling at the eastern equatorial Pacific and hence provides for an increase (decrease) in east-west SST contrast which then strengthens (decrease) the zonal pressure gradient. This hypothesis is able to 1 Chapter 1. Introduction describe the self-sustainment or growth in the warm or cold phase of ENSO. It cannot explain the turnabout between phases. The canonical E N S O : Another ENSO mechanism hypothesis, called the \"canonical ENSO\" by Cane (Cane, 1983), describes the different ENSO phases by summarizing a composite picture (Rasmusson and Carpenter, 1982). A canonical ENSO episode develops as follows: The prelude stage: This stage includes a long period (more than 18 months) of \"build-up\" in the equatorial Pacific. During this period, stronger than normal easterly winds move the water from the eastern tropical Pacific to the western tropical Pacific. There also is a stronger than normal upwelling in eastern equatorial Pacific and hence a stronger than normal zonal SST gradient. The sea level increases in the western Pacific and decreases in the eastern Pacific. Thermocline at the equator is deepened in the western Pacific and depressed in the eastern Pacific. The onset stage: The onset occurs during the boreal fall before a warm event. In this phase the rapid collapse of the easterly wind stress excites an equatorial Kelvin wave in September or October. The easterly wind stress is replaced by westerly wind stress in October or November, leading to a relaxation of the sea level slope. The event stage: Warming off the coast of South America begins in December or January and continues on until June of the following year. Accompanying the positive SST anomalies in the central equatorial Pacific are large positive SST anomalies in the equatorial eastern Pacific area which reach their maxima in April, May or June. Positive SST anomalies extend westward or northwestward from the area off the coast of South America into the central equatorial Pacific in subsequent months. During this time sea level in the equatorial eastern Pacific is anomalously high and the westerly trade winds are intensified in the central 2 Chapter 1. Introduction equatorial Pacific. The mature stage: The positive SST anomalies concentrate on the central equatorial Pacific during the winter and reach their peaks. The mature phase is followed, one season later, by a rapid decline to cool SSTs , with the central Pacific cooling. Delayed oscillator: Delayed oscillator theory (Suarez and Schopf, 1988; Battisti and Hirst, 1989) explains the observed ENSO phenomena with Kelvin and Rossby Wave propagation across the Pacific. It assumes the strongest coupling between the wind stress anomaly and the equatorial SST pattern takes place in the central equatorial Pacific. An initial westerly wind stress anomaly induced by a positive SST perturbation in the central Pacific leads to the generation of an eastward propagating downwelling Kelvin wave with maximum signal on the equator and an upwelling Rossby wave with maximum signal symmetrically located on either side of the equator. The Kelvin wave causes the thermocline to deepen in the eastern equatorial Pacific, thereby increases the warming. In the meantime, the Rossby wave traveling westward is reflected off the western boundary as an upwelling Kelvin wave. The reflected upwelling Kelvin wave returns the SST and thermocline back to normal, terminating the warm event. Western Pacific oscillator: The western Pacific oscillator hypothesis (Weisberg and Wang, 1997) obviates the need for a reflecting weastern boundary by emphasizing the atmosphere-ocean interaction in the western Pacific. An initial westerly wind stress anomaly induced by a positive SST anomaly leads to off-equator cyclones in the western central equatorial Pacific. The westerly anomaly increases the thermocline depth and SST, resulting in positive feedback in the eastern and central equatorial Pacific. In the meantime, the enhanced westerly anomaly produces a high sea level pressure (SLP) over the western Pacific through Ekman pumping, hence 3 Chapter 1. Introduction off-equator anti-cyclones appear in the western tropical Pacific. Then the off-equator anti-cyclones cause thermocline and SST to decrease in the western Pacific. As a result of the off-equator anti-cyclones, an equatorial easterly wind is generated in the western Pacific. The upwelling region initiated by the easterly wind evolves eastward along the equator, providing a negative feedback in the eastern equatorial Pacific area. Recharge-discharge oscillator: The recharge-discharge oscillator model does not emphasize the explicit role of wave propagation (Jin, 1997). The tropical Pacific heat content is of essential importance in the ENSO dynamical adjustment. The tropical subsurface ocean heat content is recharged after a mature ENSO cold phase and is discharged after a mature ENSO warm phase. An initial westerly anomaly wind stress induced by a SST anomaly in the western central equatorial Pacific leads to a positive SST anomaly and a positive thermocline depth anomaly in the central and eastern Pacific. Such a process involves a positive feedback between a westerly wind stress anomaly in the central equatorial Pacific and a positive SST anomaly in the central and eastern equatorial Pacific. The positive feedback continues until the ENSO oscillation develops into its mature warm phase. At the same time, the wind stress also gradually reduces the thermocline depth in the western Pacific, leading to a negative zonal mean thermocline depth anomaly across the Pacific. The tropical subsurface heat content is discharged during this period, and is only replenished when the ENSO oscillation develops into its mature cold phase. Advective-reflective oscillator: In this model (Picaut, 1997), an initial westerly wind anomaly drives a local zonal current that advects the oceanic convergence zone-eastern edge (OCEE) of the warm pool against the weaker mean zonal current. The eastward warm pool displacement decreases the 4 Chapter 1. Introduction east-west SST gradient and strengthens the equatorial westerly wind anomaly. The positive feedback increases the SST in the eastern Pacific. ENSO is in its warm phase. At the same time, the westerly wind anomaly also excites Rossby waves and Kelvin waves. Downwelling Kelvin waves propagate eastward and upwelling Rossby waves propagate westward. There are three sets of currents in the equatorial Pacific. One set is generated by the Kelvin waves coming from the western boundary. One set is generated by the Rossby waves coming from the eastern boundary. Another set is the mean zonal current. The combination of currents in the equatorial Pacific stops the movement of the OCEE. Eventually, the OCEE moves westward and finally moves back toward the middle of the Pacific. Once the midpoint is crossed, the wind shifts from westerly to easterly and E l Nino turns into La Nina; then La Nina turns into E l Nino and the ENSO phases repeat indefinitely. Unified oscillator: The delayed oscillator, the western Pacific oscillator, the recharge-discharge oscillator, and the advective-reflective oscillator have previously been proposed to interpret the oscilla-tory nature of ENSO. Motivated by the existence of these different ENSO oscillators, Wang (Wang, 2001) developed a unified ENSO theory that includes the physics of the previous ENSO oscillators. In this model, if the western Pacific wind-forced response is neglected, thermocline and zonal wind stress anomalies in the western Pacific are decoupled from the coupled system, and the unified oscillator reduces to the delayed oscillator. If wave reflection at the western boundary is neglected, the unified oscillator reduces to the western Pacific oscillator. Most of the physics of the advective-reflective oscillator are implicitly included in the unified oscillator. Linear system with the stochastic forcing: Another theory assumes that the ENSO is generated by a damping linear system with 5 Chapter 1. Introduction stochastic atmospheric forcings (Penland and Sardeshmukh, 1995; Moore and Kleeman, 1999; Thompson and Battisti, 2001). Thompson and Battisti (2001) stochastically forced a linear model that originally did not sustain internally generated variability and found that it reproduces all of the basic features of ENSO, suggesting that the essential characteristics of ENSO appear to be governed by linear processes. 1.2 E N S O decadal variability Besides interannual variability, the Pacific Ocean has decadal\/interdecadal climate variabil-ity, leading to long-term modulation of ENSO (Zhang et al, 1997; Wang and An, 2001; Rodgers et al, 2004). Zhang et al. (1997) argued that besides the well-known abrupt cli-mate regime shifts towards a warmer tropical eastern Pacific and a colder extratropical central North Pacific in the mid 1970s, the global SST also has experienced an analogous shift during 1920s and an opposite polarity shift during the 1940s. The rebound of the climatological state in the Pacific Ocean since 1998 can be viewed as another opposite shift relative to the shift in the mid 1970s (McPhaden and Zhang, 2004). Like the interannual ENSO dynamics mentioned in the previous section, the fundamental dynamics of ENSO decadal variability is still not settled. We do not know what causes the ENSO decadal variability. The current literature includes a number of studies proposing mechanisms for this. The prevailing mechanisms for ENSO decadal variability can be separated into two broad categories: (1) tropical-extratropical interactions and (2) purely tropical processes. Within these two broad categories there are a number of competing hypotheses. For ex-ample, regarding the mechanism for the first category, there have been a number of studies emphasizing the importance of oceanic processes (Gu and Philander, 1997; Kleeman et al, 1999; McPhaden and Zhang, 2002). Atmospheric tropical-extratropical bridge processes 6 Chapter 1. Introduction have also been considered as a potential mechanism for ENSO decadal variability (Barnett et al, 1999). For the second category mechanism, some studies suggested that the decadal variability in the tropical Pacific can be generated by purely tropical nonlinear processes (Knutson and Manabe, 1998; Jin, 2001), while others proposed it comes from atmospheric noise and unrelated to low-frequency changes in the tropical Pacific mean state (Thompson and Battisti, 2001; Fliigel et al, 2004). Deser et al (2004) also examined the linkage between the tropical interdecadal variability and the interdecadal climate fluctuations over the North Pacific during boreal winter based on the analysis of a wide variety of observational datasets in the 20th century. Their study showed that the interdecadal variability over the tropical Pacific can impact the interdecadal variability over the North Pacific. Thus the interdedacal variability over the tropical Pacific and that over the North Pacific may influence each other. 1.3 Interdecadal E N S O nonlinearity and predictability The topic of nonlinearity and predictability in ENSO has advanced significantly in recent years. Predicting the future state with a given initial state is a fundamental problem in climate research. The nonlinearity and predictability of ENSO is quite different in the numerical models at the decadal\/interdecadal time scale (Ji et al, 1996; Chen et al, 2004). Chen et al (2004) made retrospective 148-year SST hindcasts using an intermediate coupled model with a specific lead time in the tropical Pacific. The model successfully predicted all prominent E l Nino episodes within 148 years at a lead time up to 2 years. The predictability of the Nino3.4 index (averaged SST over 5\u00b0S-5\u00b0N, 170\u00b0W-120\u00b0W) varied decade by decade in their coupled model. An (2004) mentioned that the interdecadal change in the predictability 7 Chapter 1. Introduction was related to the interdecadal change in ENSO asymmetry (between the warm El Nino and the cold La Nina states). Eccles and Tziperman (2004) found that the ENSO period was affected by ENSO's nonlinearity in a delayed oscillator model. They calculated the dependence of the ENSO period on its amplitude in both the strongly nonlinear regime and the weakly nonlinear regime. The ENSO period was found to decrease with increasing ENSO amplitude in the strongly nonlinear regime. In contrast, in the weakly nonlinear regime, the period tended to increase with the amplitude. They concluded that the weakly nonlinear case reflected \"the dependence of the period on the ocean-atmosphere coupling strength rather than the dependence on the amplitude.\" 1 . 4 R e s p o n s e o f E N S O t o g l o b a l w a r m i n g How ENSO responds to global warming is still an open question. Some climate models suggest that global warming will foster stronger and more frequent E l Nifios (Timmermann et al., 1999; Collins 2000). Collins (2000) argued that increased greenhouse gases may result in ENSO events becoming larger in amplitude and more frequent than present episodes in the climate model simulations based on version 2 of the Hadley Centre coupled model (HadCM2). However, in a subsequent analysis based on version 3 of the Hadley model (HadCM3), Collins et al (2005) found that he could not detect a change in magnitude or frequency of ENSO as greenhouse gases increased, thus contradicting the results of his earlier study. These differences highlight the level of uncertainty associated with ENSO under global warming. On the other hand, the circulations associated with ENSO over the tropical Pacific has changed with global warming (Vecchi et al, 2006). McPhaden and Zhang (2002) also pointed out that the shallow meridional overturning circulation has slowed down in the 8 Chapter 1. Introduction recent decades. 1.5 Objectives o f the thesis The first thesis objective is to investigate ENSO decadal\/interdecadal variability and its possible cause in the tropical Pacific. As the tropical Pacific climatology is calculated from individual ENSO episodes, it is affected by changes in the properties of ENSO. On the other hand, does the climatology have an influence on the properties of individual ENSO episodes? Some studies have shown that ENSO properties such as the period and amplitude have changed with different climatologies. Fedorov and Philander (2001) found that the changes in the mean state played critical roles in determining the stability of the tropical ocean-atmosphere coupling modes. Their theoretical studies displayed the changes in ENSO properties caused by background changes in the wind stress intensity, thermocline depth and temperature difference across the thermocline. They concluded that the observed ENSO property changes since the 1960s were consistent with their theory. In this thesis, how the climatological state impacts the ENSO characteristic properties is illustrated based on observational data and outputs from an intermediate coupled model. The nonlinearity and predictability of ENSO are investigated as well. The role of heat flux is examined in both the observational data and modeling results. Next, the role of ocean surface current in ENSO is also studied under the 2 different climatological states (before and after the climate regime shift of the mid 1970s). As aforementioned, ENSO changes its characteristic properties (such as the period, in-tensity, nonlinearity and predictability) on the decadal\/interdecadal time scale. Our studies show that changes in the climatology can, at least in part, cause changes in the ENSO prop-erties. The next step is to investigate the mean state change in the tropical Pacific due to 9 Chapter 1. Introduction global warming. Thus the second objective of this thesis is to study how global warming affects the tropical mean state and ENSO properties. The thesis is structured as follows: Chapter 2 illustrates the impact of regime shift in the mean state on ENSO properties (including predictability). Chapter 3 compared the nonlinearity and predictability in the ENSO system with that of the Lorenz (1963) nonlinear system. The role of surface currents in ENSO is discussed over the different climatological states based on the observed data in Chapter 4. The effects of increased greenhouse gas on ENSO and associated ocean circulations are shown in Chapter 5. Summary and discussion are given in Chaper 6. 1 0 Chapter 1. Introduction 1.6 References Barnett T. P., D. W. Pierce, M . Latif, D. Dommenget and R. Saravanan, 1999: Inter-decadal interactions between the tropics and midlatitudes in the Pacific basin. Geo-phys. Res. Lett, 26, 615-618. DDI: 10.1029\/1999GL900042. Bat t i s t i D. S., and A. C. Hirst, 1989: Interannual variability in the tropical atmosphere-ocean system: Influence of the basic state, ocean geometry, and nonlinearity, J. Atmos. Sci., 46, 1687-1712. Ba t t i s t i D. S., and E. S. Sarachik, 1995: Understanding and predicting ENSO. Rev. Geo-phys., (Suppl.), 1367-1376. Bjerknes J., 1969: Atmospheric teleconnections from the equatorial Pacific, Mon. Wea. Rev., 97, 163-172. Cane M . A., 1983: Oceanographic events during E l Nino, Science, 222, 1189-1195. Chen D, M . A. Cane, A. Kaplan, S. E. Zebiak, and D. Huang, 2004: Predictability of E l Nino in the past 148 years. Nature 428, 733-736. Coll ins M . , 2000: The E l Nino Southern Oscillation in the second Hadley center coupled model and its response to greenhouse warming. J. Climate, 13, 1299-1312. Coll ins M . and the CMIP Modelling Groups, 2005: E l Nino- or La Nina-like climate change?. Climate Dyn., 24, 89-104. Deser C, A. S. Phillips and J. W. Hurrell, 2004: Pacific interdecadal climate variability: linkages between the Tropics and the North Pacific during boreal winter since 1900. J. Climate, 17, 3109-3124. Chapter 1. Introduction Eccles F. and E. Tziperman, 2004: Nonlinear effects on ENSO's period. J. Atmos. Sci., 61, 474-482. Fedorov A. V . and S. G. H. Philander, 2001: A stability analysis of tropical ocean-atmosphere interactions: Bridging measurements and theory for E l Nino). J. Climate, 14, 3086-3101. Fl i ige l M . , P. Chang and C. Penland, 2004: The role of stochastic forcing in modulating ENSO predictability. J. Climate, 17, 3125-3140. G u D. and S. G. H. Philander, 1997: Interdecadal climate fluctuations that depend on exchanges between the tropics and extratropics. Science, 275, 805-807. J i M . , A. Leetmaa and V. E. Kousky, 1996: Coupled model predictions of ENSO during the 1980s and the 1990s at the National Centers for Environmental Prediction, J. Climate, 9, 3105-3120. J i n F.-F., 1997: An Equatorial Ocean Recharge Paradigm for ENSO Part I: Conceptual Model, J. Atmos. Sci., 54, 811-829. J i n F.-F., 2001: Low-frequency modes of tropical ocean dynamics, J. Clim., 14, 3874-3881. Kleeman R., J. P. McCreary and B. A. Klinger, 1999: A mechanism for generating ENSO decadal variability. Geophys. Res. Lett., 26, 1743-1746. Knutson T. R. and S. Manabe, 1998: Model assessment of decadal variability and trends in the tropical Pacific Ocean. J. Climate, 11, 2273-2296. La t i f M . , and T. P. Barnett, 1996: Decadal climate variability over the North Pacific and North America: Dynamics and predictability. J. Climate, 9, 2407-2423. 12 Chapter 1. Introduction Lat i f M . , D. Anderson, T. P. Barnett, M . A. Cane, R. Kleeman, A . Leetmaa, J. J. O'Brien, A. Rostati, and E. K. Schneider, 1998: Predictability and Prediction. J. Geophys. Res., 103, 14375-14394. M c P h a d e n M . J., and D. Zhang, 2002: Slowdown of the meridional verturning circulation in the upper Pacific Ocean. Nature, 415, 606-608. M c P h a d e n M . J. and D. Zhang, 2004: Pacific Ocean circulation rebounds, Geophys. Res. Lett., 31, L18301. D01:10.1029\/2004GL020727. Moore A. M . and R. Kleeman, 1999: Stochastic forcing of ENSO by the intraseasonal oscillation. J. Climate 12, 1199-1220. Neel in J. D., D. S. Battisti, A. C. Hirst, F.-F. Jin, Y . Wakata, T. Yamagata, and S. E. Zebiak, 1998: ENSO theory. J. Geophys. Res., 103, 14261-14292. Penland C. and P. D. Sardeshmukh, 1995: The optimal-growth of tropical sea-surface temperature anomalies, J. Climate, 8, 1999-2024. Picaut J., F. Masia, Y . Dupenhoat, 1997: An advective-reflective conceptual-model for the oscillatory nature of the ENSO, Science, 277, 663-666. D01:10.1126\/science. 277.5326.663. Rasmusson E. M . and T. H. Carpenter, 1982: Variations in tropical sea surface tempera-ture and surface wind fields associated with the Southern Oscillation\/El Nino, Mon. Wea. Rev., 110, 345-384. Rodgers K. B., P. Friederichs and M . Latif, 2004: Tropical Pacific decadal variability and its relation to decadal modulations of ENSO. J. Climate, 17, 3761-3774. 13 Chapter 1. Introduction Suarez M . J. and P. S. Schopf, 1988: A delayed action oscillator for ENSO, J. Atmos. Sci., 45, 3283-3287. Thompson C. J. and D. S. Battisti, 2001: A linear stochastic dynamical model of ENSO. Part II: Analysis. J. Climate 14, 445-466. Timmermann A, Oberhuber, J., Bacher, A., Esch, M . , Latif, M . and Roeckner, E. 1999: Increased E l Nino frequency in a climate model forced by future greenhouse warming. Nature, 398, 694-696. Vecchi G. A., B. J. Soden, A. T. Wittenberg, I. M . Held, A. Leetmaa and M . J. Harrison, 2006: Weakening of tropical Pacific atmospheric circulation due to anthropogenic forcing. Nature, 44, 73-76. Wallace J. M . and D. S. Gutzler, 1981: Teleconnections in the geopotential height field during the Northern Hemisphere winter. Mon. Wea. Rev., 109, 784-812. Wang B. and S.-I. An, 2001: Why the properties of E l Nino changed during the late 1970s. Geophys. Res. Lett, 28, 3709-3712. Wang C., 2001: A unified oscillator model for the E l Nino-Southern Oscillation, J. Climate, 14, 98-115. Weisberg R. H., C. Wang, 1997: A western Pacific oscillator paradigm for the El Nino-Southern Oscillation. Geophys. Res. Lett., 24, 779-782. Zhang Y , J. M . Wallace and D. S. Battisti, 1997: ENSO-like interdecadal variability: 1900-93. J. Climate, 10, 1004-1020. 14 Chapter 2 The Influence of Climate Regime Shift on ENSO 1 2.1 Introduction The characteristics of the El Nino\/Southern Oscillation (ENSO) have experienced an abrupt change or 'regime shift' in the mid 1970s, which is usually attributed to decadal\/interdecadal climate variability in the Pacific Ocean (Zhang et al, 1997, Wang and An, 2001). Recent studies (McPhaden and Zhang, 2004) found that the Pacific Ocean has rebounded since 1998, where the rebound can be viewed as an opposite shift relative to the shift in the mid 1970s. There is no consensus on the fundamental dynamics of ENSO decadal variability. Some modeling results suggested that the decadal variability in the tropical Pacific origi-nates from mid-latitude variability (Latif and Barnett, 1996), while others proposed it comes from tropical internal nonlinear instability (Knutson and Manabe, 1998) or uncoupled at-mospheric noise (Thompson and Battisti, 2001; Fliigel et al, 2004). Two hypotheses have been proposed to explain the linkage between the extratropical decadal variability and the tropical decadal variability, one by an oceanic teleconnection pathway (Gu and Philander, 1997; Kleeman et al, 1999), the other by an atmospheric bridge (Barnett et al, 1999). *A version of this chapter has been published. Ye, Z. and Hsieh, W . W . (2006) The influence of climate regime shift on E N S O . Climate Dyn., 26: 823-833. D01:10.1007\/s00382-005-0105-5. 15 Chapter 2. The Influence of Climate Regime Shift on ENSO Since a decadal\/interdecadal change in the tropical Pacific amounts to a change in the background mean state for the ENSO phenomenon, how the background state affects ENSO properties has become an area of active investigation. Studies have shown that changes in the background state would alter ENSO in period, intensity, direction of propagation and spatial structure. For instance, Fedorov and Philander (2001) studied theoretically the changes in ENSO properties caused by background changes in the wind stress intensity, thermocline depth and temperature difference across the thermocline, and concluded that the observed ENSO property changes since the 1960s are consistent with their theory. Wang and An (2002) proposed that the changes in the background winds and the corresponding ocean dynamic fields played the dominant role in the abrupt change in ENSO properties since the late 1970s. A caveat is that the periods used to calculate the mean states are relatively short, and the presence of an extra E l Nino or a very strong E l Nino during a particular decadal period could introduce a difference in the calculated mean states, i.e. the mean decadal states are themselves influenced by ENSO episodes. Recently, the nonlinearity and asymmetry of decadal variability in the tropical Pacific were recognized (Rodgers, 2004). Jin et al (2003) suggested that nonlinear dynamic heating in the tropical Pacific ocean heat budget is essential in generating the nonlinearity and asymmetry of the decadal-scale ENSO. Changes in the asymmetry and nonlinearity of ENSO over different climate regimes have also been studied. Wu and Hsieh (2003) examined that the leading nonlinear canonical correlation analysis (NLCCA) mode between the tropical Pacific sea surface temperature (SST) and wind stress (WS), which showed that after the climate regime shift of the mid 1970s, there was an eastward shift in the westerly anomalies during warm ENSO episodes, but no shift in the easterly anomalies during the cold ENSO episodes, thereby enhancing the asymmetry between E l Nino and La Nina in the post-shift 16 Chapter 2. The Influence of Climate Regime Shift on ENSO regime compared to the pre-shift regime. ENSO predictability has been found to be higher during the decades when the ENSO asymmetry is large (An, 2004). In this chapter, by analyzing a coupled model of intermediate complexity and observed data, we investigate how changes in the background climate state affect properties of ENSO, especially the asymmetry between E l Nino and La Nina, zonal shifts in the observed net shortwave radiation flux and latent heat flux, shifts in the model convergent heating rate, and ENSO predictability. This chapter is organized as follows: In Section 2.2, we introduce the datasets and the nonlinear data analysis method. Section 2.3 describes the revised coupled model. In Section 2.4, observed data are first analyzed for the 'regime shift' in the asymmetry of ENSO, then the intermediate coupled model is run to see how different background states (of surface wind, heat fluxes, and sea surface temperatures) change the ENSO properties. A discussion follows in Section 2.5. 2.2 Data and methods 2.2.1 Data The National Oceanic and Atmospheric Administration (NOAA) extended reconstructed SST dataset was used in this study. The extended reconstructed sea surface temperature (ERSST) was constructed using the most recently available international comprehensive ocean-atmosphere data set (ICOADS) SST data and improved statistical methods that allow stable reconstruction using sparse data. This monthly analysis begins in January 1854, but because of sparse data the analyzed signal was heavily damped before 1880 (Smith and Reynolds, 2004). The monthly WS was objectively analyzed at Florida State University (Bourassa et al, 2001), while the sea surface heat flux data is from da Silva et al 17 Chapter 2. The Influence of Climate Regime Shift on ENSO (1994) (which terminated after 1993). The climatological monthly mean was removed from the data to obtain the monthly anomaly. 2.2.2 N o n l i n e a r p r i n c i p a l c o m p o n e n t a n a l y s i s For the time series data x(i) = [xi,...,xi\\, where each variable Xi is a time series containing n observations, the principal components analysis (PCA) method looks for u, a linear combination of the Xi, and an associated vector a, with u(t) = a \u2022 x(t) so that J = ( | | x ( i ) - a u ( t ) | | 2 ) is minimized. Here (...) denotes a sample or time mean, u, the first principal component (PC), is a time series, while a, the first empirical orthogonal function (EOF), describes a spatial pattern. The fundamental difference between nonlinear principal components analysis (NLPCA) and P C A is that N L P C A allows a nonlinear mapping from x to u and a nonlinear inverse mapping from u to x', whereas P C A only allows linear mapping (Kramer, 1991; Hsieh, 2004). For N L P C A , u = .F(x), x' = G{u), where F and G denote the nonlinear mapping and inverse mapping functions respectively (Fig. 2.1). For more details, see Appendix B. The N L P C A method is used to analyze the nonlinear modes of SST and WS in the tropical Pacific. Prior to N L P C A , combined P C A , also called combined EOF analysis, 18 Chapter 2. The Influence of Climate Regime Shift on ENSO Figure 2.1: A schematic diagram illustrating the neural network(NN) model for performing the N L P C A . The model is a standard feedforward NN (i.e. multi-layer perceptron) (Bishop, 1995), with 3 'hidden' layers of variables or 'neurons' (denoted by circles) sandwiched be-tween the input layer x on the left and the output layer x' on the right. Next to the input layer is the encoding layer, followed by the 'bottleneck' layer with a single neuron u, then the decoding layer, and finally the output layer, i.e. a total of 4 layers of transfer functions are needed to map from the inputs to the outputs. A neuron at the zth layer receives its value from the neurons V j _ i in the preceding layer, i.e. vi = \/i(wj \u2022 V j _ i + b), where W j is a vector of weight parameters and b a bias parameter, and the transfer functions fa and fa are the hyperbolic tangent functions, while fa and fa are simply the identity functions. Effectively, a nonlinear function u = F(x) maps from the higher dimension input space to the lower dimension bottleneck space, followed by an inverse transform x' \u2014 G(u) 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 \u2014 (||x \u2014 x'| | 2) (i.e. the mean square error, MSE) is minimized. Through the optimization, the values of the weight and bias parameters are solved. See Appendix B for details. 19 Chapter 2. The Influence of Climate Regime Shift on ENSO is performed on the SST and WS anomalies in order to reduce the data dimension - we performed P C A analysis on the combined field of SST anomalies and WS anomalies, with each variable first normalized by its standard deviation, then, the leading 6 PCs are used as the inputs to the N L P C A model, which finds a curve in the 6-dimensional PC space. 2.3 Coupled model The coupled model is based on the Zebiak-Cane model (Zebiak and Cane, 1987), where the atmospheric component consists of steady-state, linear shallow-water equations on an equatorial beta plane with a nonlinear convergence feedback, while the oceanic component is a linear reduced-gravity model (see Appendix A for details). With a time step of (about) 10 days, the ocean model domain covers the tropical Pacific from 29\u00b0S to 29\u00b0N, and from 124\u00b0E to 80\u00b0W, where the grid for ocean dynamics is 2\u00b0 longitude by 0.5\u00b0 latitude, and the grid for SST physics and the atmospheric model is 5.625\u00b0 longitude by 2\u00b0 latitude. The temperature equation for the surface layer follows Zebiak and Cane (1987): dT \u2014 dT \u2014 = - u i - v ( r + T ) - u i - \\7T-{M(ws + ws)-M(ws)} x \u2014 -M(ws + w s ) ^ ^ - a s T - K t ^ - ^ + A h A h T , (2.1) Hi Hi where we have added the last two terms to simulate vertical mixing and horizontal diffusion in the temperature equation (Boulanger and Menkes, 2001). Here ui(a:,y, t) and ws(x,y, t) are the prescribed climatological monthly mean horizontal current and upwelling in the surface layer respectively, T(x,y,t) is the prescribed mean SST, dT(x)\/dz the prescribed mean vertical temperature gradient, the mean surface layer depth Hi \u2014 50 m, the diffusion coefficient as = (125 day) - 1 , Kt - 2.5 x 1 0 - 5 m s - 1 , Ah = 2000m 2s - 1 , the function M is 20 Chapter 2. The Influence of Climate Regime Shift on ENSO denned by 0. x < 0 M(x) = { (2.2) x x>0, and the entrainment velocity is wa = H1( dui dx + dvi dy (2.3) The entrainment temperature anomaly, Te, is given by T e - 7 T s u b + ( l - 7 ) T , (2.4) The subsurface temperature anomaly T s ub is assumed to arise from depressing or raising the thermocline according to the anomalous inflow or outflow of mass in the upper layer. This corresponds roughly to considering vertical displacements of the mean subsurface tem-perature profile in accordance with the anomaly in the upper layer depth. To quantify this, the mean vertical temperature profile is approximated by the analytic form of a hyperbolic tangent function. Thus, T s ub originally has the empirical parameterization form (Zebiak and Cane, 1987) : with h(x) the prescribed mean upper layer depth, 7 = 0.75, T\\ = 28\u00b0C, T2 \u2014 \u2014 40\u00b0C, b\\ = (80 m)\" 1 and b2 = (33 m)\" 1 . Here an empirical parameterization scheme was developed to estimate T s ub in terms of other ocean dynamical variables in the temperature equation (2.1). With the ocean model forced by the observed WS during 1964-1998, The equation (2.1) was inverted for T s ub using the observed SST for T and the model anomalous currents (Zhang et al, 2005). P C A was separately applied to the simulated upper layer depth anomalies h, and to T s ub to extract h>0 h<0 (2.5) 21 Chapter 2. The Influence of Climate Regime Shift on ENSO the leading PCs. The nonlinear regression relationships between the first three PCs of h (the predictors) and the first three PCs of T s u b (the predictands) were obtained by a neural network method (Tang and Hsieh, 2002; L i et al, 2005). Then the nonlinear regression relationships were used to replace equation (2.5) in the coupled model. Nonlinear relations between the leading PCs of T s u b and those of h can be seen in Fig. 2.2. L i et al (2006) found that the nonlinear regression was better than the linear regression, although Fig. 2.2 does not exhibit strong nonlinear relationship between h and T s u b- The nonlinear regression improves the SST simulation in the coupled model. The observed SST averaged in equatorial area (5\u00b0S-5\u00b0N) showed that during almost each El Nino between 1960-2000 the warm maximum SST anomalies occupied the whole eastern Pacific. The E l Nino of 1997\/1998 was the only one where the maximum SST anomalies were located in the eastern boundary area (Fig. 2.3). During La Nina, the observed minimum SST anomalies were located further west than the maximum SST anomalies during El Nino. Fig. 2.4 gives an example of the simulated SST and WS anomalies during model years 280-300 in the coupled model to compare the original parameterization with the new one. The simulated SST anomalies in the equatorial area (5\u00b0S-5\u00b0N) were always manifested mainly in the eastern area in the original coupled model during both E l Nino and La Nina (Fig. 2.4b). However, the simulated maximum SST anomalies for the new parameterization extended much further west into the central equatorial Pacific (Fig. 2.4d), close to the observations (Fig. 2.3). The maximum westerly WS anomalies (Fig. 2.4a) for the original parameterization were located around 160\u00b0 W-140\u00b0 W during E l Nino, exhibiting an eastward shift relative to the observed WS anomalies located around 180\u00b0W-160\u00b0W (not shown). Under the new parameterization, the model maximum WS has shifted to the correct location of 180\u00b0W-160\u00b0W during E l Nino (Fig. 2.4c). The negative (i.e. easterly) WS anomalies 22 Chapter 2. The Influence of Climate Regime Shift on ENSO ( a ) 30 T - 10 o CL \u2022\u00a7-10 -30 -50 40 20 CNJ o 0 Q -XI 5-20 -40 -60 -50 -30 -10 10 30 50 h PC1 (C) -60 -40 -20 0 20 40 60 h PC1 -30 -20 - 10 0 10 20 30 40 h PC2 (d) 40 20 CM o 0 CL XI 5-20 r --40 -60 -30 -20 -10 0 10 20 30 40 h PC2 Figure 2.2: Relations between PCI of the T s ub anomalies and (a) P C I and (b) PC2 of the h anomalies, as found by nonlinear regression via neural networks (overlapping circles), and linear regression (dashed lines), with the data shown as dots. The corresponding relations for PC2 of the T s ub anomalies are shown in (c) and (d). 23 Chapter 2. The Influence of Climate Regime Shift on ENSO in the model with the new parameterization have also strengthened and shifted westward compared to the original parameterization (2.4a and c). Thus the new parameterization has corrected the original coupled model bias where the variability (for both SST and WS) was located too far east. Note that there are unrealistic significant easterly wind stress anomalies in the equatorial area near 130\u00b0W in the model with the new parameterization. Such easterly wind stress anomalies are absent in the observations (Deser and Wallace, 1990). These spurious anomalies of zonal wind are a known problem of the Gil l atmospheric model that ties atmospheric heating too closely to SST anomalies. To analyze the cause of shifts in SST and WS anomalies in the new coupled model, both the simulated subsurface temperature and thermocline anomalies with the original and the new parameterization in the coupled model are shown in Fig. 2.5. In the new coupled model, the thermocline depth anomalies extend further west during E l Nino than in the original model (2.5a and c). The maximum thermocline anomalies were located between longitude 180\u00b0 and 120\u00b0W under the new parameterization. The corresponding positive subsurface temperature anomalies extended to the central equatorial Pacific (Fig. 2.5d), whereas in the original model, the maximum subsurface temperature anomalies appeared further east during El Nino with double the amplitude. The change in subsurface temperature simulation in the new model led to the better SST simulation in the coupled model, since the subsurface temperature directly affected the SST via the upwelling in the eastern equatorial Pacific. 24 Chapter 2. The Influence of Climate Regime Shift on ENSO Figure 2.3: Observed SST (\u00b0C) anomalies along equator (averaged in 5\u00b0S-5\u00b0N). The contour interval is 1 \u00b0C with positive values shaded. 25 Chapter 2. The Influence of Climate Regime Shift on ENSO Figure 2.4: Zonal WS (dyn c m - 2 ) and SST (\u00b0C) anomalies along the equator (averaged over 5\u00b0S-5\u00b0N) in the coupled model from the original model (top panels) and the new model with the neural network parameterization (bottom panels). Positive values arc shaded. 26 Chapter 2. The Influence of Climate Regime Shift on ENSO (a) Thermocline depth (b) Subsurface temperature (c) Thermocline depth (d) Subsurface temperature 0298- 0298-q ^ ' z : . . - , * noQfi 0296-0294-0294-J = = l , . : . o = ? g ^ v ^ - - ^ \u00abS 0292 \u2022 V >> - H 0290 -4) \u2022d \u00b0 0288\u2022 0286-\u00ab0292-0) l>\u00bb - 0290 -V \"O \u00b00288-S 0286-><: \u00ab===\u2022! ,_\u201e o\u2014 0284-0282 0284-0282-r l r~ ' 1 ' where (...) denotes the temporal mean, @y measures the size of the nonlinear term relative to the size of the linear terms in (3.3), and similarly, Bz for (3.4). We next perform four numerical experiments by varying the parameters (a,b,c): Case 1 uses the same parameters as in Lorenz (1963) (Table 3.1); Case 2 multiplies the parameters 66 Chapter 3. Changes in the Nonlinearity and Predictability (a) Lead months (b) 0.3 r \u2022 0 5 10 15 Lead months Figure 3.4: Ensemble mean predictability (top) and MSE (bottom) of the Nino 3.4 SST anomaly index, as given by the cross-validated correlation between the predicted and actual index in the coupled model using climatologies from the 1961-75 regime (solid lines) and the 1981-95 regime (dashed lines). The thick lines are from nonlinear regression, the thin lines, linear regression. Cross-validation was performed by dividing each 100-year data record into five segments, where for each segment chosen to test the forecast correlation skills, the other four were used to build the forecast model. Error bars indicate \u00b1 1 standard error of the ensemble mean. 67 Chapter 3. Changes in the Nonlinearity and Predictability x 107 I 1 1 1 1 1 1 1 r 1\/52 I I t :': 1\/49 i! I E -*\u2014\u00bb o CD O. 0.015 0.016 0.017 0.018 0.019 0.02 0.021 0.022 0.023 0.024 0.025 F r e q u e n c y ( 1 \/ m o n t h s ) Figure 3.5: Spectrum of the Nino 3.4 SST anomaly index in the coupled model using climatologies from the pre-shift regime (solid curve) and post-shift regime (dashed curve). 68 Chapter 3. Changes in the Nonlinearity and Predictability of Case 1 by the uniform factor of 1.3; Case 3 uses the parameters from Basu and Foufoula-Georgiou (2002); and Case 4, from Eisner and Tsonis (1992). The nature of the oscillations varies among the four cases (Fig. 3.6). For each case, 21 ensemble runs were made by adding one percent stochastic noise onto the initial data (i.e., the initial data x0,y0: zQ were changed to ( l + 0.01ai)xo, ( l + 0.01a2)yo, (l+0.01a3).zo, where a i , a2 and a 3 were Gaussian random numbers with zero mean and unit standard deviation). Table 3.1 lists j3v and 3Z calculated from (3.5) for the four cases. The main spectral periods Tx and Tz determined from the power spectra of x and z, respectively (Fig. 3.7) are also listed in Table 3.1. The spectral behavior of y (not shown) is basically the same as as that for x, but z shows sharper spectral peaks than x (Fig. 3.7). Proceeding inversely from Case 4 to Case 1, one finds a progressive increase in By and 8Z, and in Tx and Tz, indicating that an increase in the nonlinearity of the system coincides with an increase in the spectral period of the oscillations. 69 Case 1 (Lorenz, 1963) 2 3 (Basu, 2002) 4 (Eisner & Tsonis, 1992) param(a, b, c) (10, 28, 8\/3) (13, 36.4, 3.47) (16, 45.92, 4.0) (16, 120.1, 4.0) nonlinearity (3y 0.033 0.026 0.020 0.008 nonlinearity \/?z 0.400 0.307 0.271 0.260 x period (step) 15.2 12.4 11.2 10.4 z period (step) 15.2 11.3 9.2 5.4 Table 3.1: Ensemble mean results for Lorenz systems with different parameter pairs. Chapter 3. Changes in the Nonlinearity and Predictability Again (x, y, z) at time t were used as predictors in nonlinear regression models to separately predict these three variables at time t+lead time (where the lead time ranged from 0 to 20 time steps). Fig. 3.8 shows enhanced predictability (with higher correlation coefficient and lower MSE) for z as one proceeds inversely from Case 4 to Case 1, where the nonlinearity of the system increased. This behavior in predictability was also found for x (not shown), though for a given level of correlation skill, x tends to have longer lead times than z. Since relative to x, the predictability lead time for z tends to be shorter, and its main spectral period Tz is also < Tx and more sensitive to changes in the Lorenz system pa-rameters (Table 3.1), we now focus on z. Its predictability lead time (for 0.80 correlation skill) divided by the main spectral period Tz is 0.69, 0.69, 0.70 and 0.65 for Cases 1 to 4, respectively, i.e. in all 4 cases, the nonlinear regression model when predicting ahead by 0.7 times the period in z can attain a correlation skill of 0.8 for all four cases, thereby demonstrating that the enhanced predictability was basically due to the lengthened period of the oscillations when the nonlinearity of the Lorenz system was enhanced from Case 4 to Case 1. For the correlation skill of 0.95, its predictability lead time divided by the main spectral period Tz is 0.51, 0.50, 0.50 and 0.46 for Cases 1 to 4, respectively. So a similar conclusion can be drawn regardless of the correlation skill level chosen. A caveat on the Lorenz system: The system is chaotic only over a range of parameter values. Our numerical runs were done over a range where the Lorenz system was chaotic. If the parameters were increased or decreased by the right amount, the Lorenz system could simply converge from the initial state to a single point in the (x, y, z) space. 71 Chapter 3. Changes in the Nonlinearity and Predictability 0 10 20 30 40 50 60 T ime Figure 3.6: Normalized time series of x (top panel) and z (bottom panel) in the Lorenz model for (a) Case 1, (b) Case 2, (c) Case 3 and (d) Case 4, where the time series for different cases are vertically shifted by multiples of 5 for legibility. 72 Chapter 3. Changes in the Nonlinearity and Predictability Figure 3.7: Ensemble mean power spectrum of the normalized (a) x and (b) z components. Dash, solid, dash-dot and dotted lines represent Case 1, 2, 3 and 4, respectively. 73 Chapter 3. Changes in the Nonlinearity and Predictability i . i i i 1 0 4 8 12 16 20 Lead time (x 0.05) Figure 3.8: Ensemble mean predictability in terms of correlation skill (top) and MSE (bot-tom) of z in the Lorenz system, as given by the cross-validated correlation and MSE between the predicted and actual time series. Dash, solid, dash-dot and dotted lines represent Case 1 (a), 2 (b), 3 (c) and 4 (d), respectively, with error bars showing \u00b1 1 standard error. 74 Chapter 3. Changes in the Nonlinearity and Predictability 3 . 4 S u m m a r y a n d c o n c l u s i o n Numerical coupled model results show that the predictability of ENSO is closely related to its nonlinearity and period. Under the post-1980s climatology, ENSOs have stronger nonlinearity and longer period. The longer period enhance the system persistence, leading to better predictability. This behavior was also found in the Lorenz chaotic system. It is important to emphasize that it is not the stronger nonlinearity that causes the better predictability, nor vice versa. The dynamics have to adjust both the nonlinearity and the predictability to be consistent with each other according to the nonlinearity-predictability relation found here. The ENSO system under different climate states and the Lorenz system with different parameters studied here were all in the chaotic domain. We did not discuss the 'extreme' nonlinear state like the limit cycle, where, the system becomes periodic. 3 . 5 R e f e r e n c e s An S.-L, 2004: Interdecadal changes in the E l Nino La-Niha asymmetry. Geophys. Res. Lett, 31, L23210. DOI:10.1029\/2004GL021699. An S.-L, and F.-F. Jin, 2004: Nonlinearity and asymmetry of ENSO. J. Climate, 17, 2399-2412. An S.-I. and B. Wang, 2000: Interdecadal change of the structure of the ENSO mode and its impact on the ENSO frequency. J. Climate 13, 2044-2055. Basu S. and E. Foufoula-Georgiou, 2002: Detection of nonlinearity and chaoticity in time series using the transportation distance function, Phys. Lett. A, 301, 413-423. 75 Chapter 3. Changes in the Nonlinearity and Predictability Bishop C. M . , 1995: Neural Networks for Pattern Recognition, Clarendon Press, 482 pp. Chen D, M . A. Cane, A. Kaplan, S. E. Zebiak, and D. Huang, 2004: Predictability of E l Nino in the past 148 years. Nature 428, 733-736. Eccles F. and E. Tziperman, 2004: Nonlinear effects on ENSO's period. J. Atmos. Sci., 61, 474-482. Eisner J. B. and A. A. Tsonis, 1992: Nonlinear prediction, chaos and noise. Bull Amer Meteor Soc, 73, 49-60. Hsieh W. W., 2004: Nonlinear multivariate and time series analysis by neural network methods. Rev. Geophys., 42, RG1003. DOI: 1 0 . 1 0 2 9 \/ 2 0 0 2 R G 0 0 0 1 1 2 . J i M . , A. Leetmaa and V. E. Kousky, 1996: Coupled model predictions of ENSO during the 1980s and the 1990s at the National Centers for Environmental Prediction, J. Climate, 9, 3105-3120. J i n F.-F., D. Neelin and M . Ghil, 1994: ENSO on the devil's staircase, Science, 264, 70-72. K i r t m a n B. P. and P. S. Schopf, 1998: Decadal variability in ENSO predictability and prediction, J. Climate, 11, 2804-2822. Lorenz E. N. , 1963: Deterministic non-periodic flow. J. Atmos. Sci., 20, 130-141. M i i n n i c h M . , M . A. Cane and S. E. Zebiak, 1991: A study of self-excited oscillations of the tropical ocean-atmosphere system. Part II: Nonlinear cases, J. Atmos. Sci., 48, 1238-1248. Penland C. and P. D. Sardeshmukh, 1995: The optimal-growth of tropical sea-surface temperature anomalies, J. Climate, 8, 1999-2024. 76 Chapter 3. Changes in the Nonlinearity and Predictability Ye Z. and W. W. Hsieh, 2006: The influence of climate regime shift on ENSO. Climate Dyn., 26, 823-833. DOI: 10.1007\/s00382-005-0105-5. Zebiak S. E. and M . A. Cane, 1987: A model E l Nino-Southern Oscillation. Mon. Wea. Rev., 115, 2262-2278. 77 Chapter 4 C h a n g e s i n t h e l e a d i n g E N S O m o d e s : R o l e o f s u r f a c e z o n a l o c u r r e n t \u00b0 4 . 1 I n t r o d u c t i o n In the late 1970s, an abrupt change of the climate state in the midlatitude as well as in the tropical Pacific was observed (Nitta and Yamada, 1989; Trenberth and Hurrell, 1994). Together with this climate shift, the main characteristics of the interannual variability over the tropical Pacific have changed (Wang, 1995; An and Wang, 2000; An et a l , 2005). On one hand, the multi-decadal variation in the climate state causes non-stationary characteristics of the E l Nino-Southern Oscillation (ENSO) (e.g., Fedorov and Philander, 2000; An and Jin, 2000). On the other hand, the multi-decadal variation can also be attributed to the internal nonlinear process of ENSO in a certain parameter regime (Timmermann et al., 2003; An and Jin, 2004). Although it is still uncertain whether the ENSO system is linear or nonlinear (e.g., Philander and Fedorov, 2003; Eisenman et al., 2005), the late 1970s 3 A version of this chapter has been published. An, S.-L, Ye, Z. and Hsieh, W . W . (2006) Changes in the leading E N S O modes associated with the late 1970s climate shift: Role of surface zonal current. Geophys. Res. Lett. 33, L14609, doi:10.1029\/2006GL026604. 78 Chapter 4. Changes in the leading ENSO modes: Role of surface zonal current climate shift appears strong enough to change the leading ENSO mode. An et al. (2005) obtained the nonlinear ENSO cycle in the evolution of tropical Pa-cific thermocline anomalies using a variant of the nonlinear principal component analysis (NLPCA) method, which performs a closed curve fit (see Appendix B.2) of the data in the multi-dimensional PC (principal component) space. Interestingly, the shape of the closed curve, especially its asymmetry, significantly changed since the late 1970s, from which we infer a change of the leading ENSO mode. An and Jin (2000) showed that the interdecadal change could lead to a quantitative change in the leading eigenmode of a stripped-down version of the intermediate ENSO model. In this chapter, observational evidence of the interdecadal change in the leading ENSO mode, especially in the role of the surface zonal currents, is revealed by N L P C A , while An et al (2006) also performed the eigenanalysis of the intermediate ENSO model linearized with respect to the two climatological mean background states. 4.2 Combined N L P C A of the sea level and current anomalies The data used are the University of Maryland Simple Ocean Data Assimilation (SODA) monthly mean sea level height, SST and upper-50m surface zonal current for January 1958-December 2001 (Carton et al., 2000). The anomalies were obtained by removing the cli-matological monthly mean, and then a 7-month running mean was applied. A combined principal component analysis (CPCA) on the monthly mean sea level height, sea surface temperature (SST) and upper-50m surface zonal current anomalies was then performed. In Fig. 4.1, the main features of sea level height, SST and surface zonal current variability are 79 Chapter 4. Changes in the leading ENSO modes: Role of surface zonal current summarized by the two leading modes: the first and second modes of sea level height are respectively a zonal contrast mode (Fig. 4.1c) and an approximate zonal-symmetric mode (Fig. 4. Id) with respect to the equator; those of SST are a horseshoe pattern like the ma-ture phase of ENSO (Fig. 4.1a) and an equatorially-confined pattern like an initial phase of ENSO (Fig. 4.1b), respectively; those of the surface current are a localized zonal\/meridional contrast mode (Fig. 4.1e) and an equatorial basin-wide mode (Fig. 4.If), respectively. Ob-viously, the first (leading) mode is highly correlated with the ENSO index, and the second mode shows a lead\/lag relationship with the ENSO index. The trajectory plot of the two leading PCs associated with C P C A is shown in Fig. 4.2. The mean location of the tracks for the pre-1980s (1958-1979) is located at the first quadrant in the PC1-PC2 plane, while that for the post-1980s (1980-2001) is located at the third quadrant, suggesting a possible regime shift. From the first six PCs of the C P C A modes, we extracted the first combined N L P C A mode, giving a closed curve in the 6-dimensioned PC space. In Fig. 4.3, the time evolution of SST, sea level height and surface zonal current associated with the first N L P C A mode along the equatorial band over one cycle are shown. These modes explain 67.5% and 73.3% of total variance during the pre- and post-1980s, respectively, which are greater than those of the leading C P C A mode (36%). For the pre-shift regime, the standing-type oscillation in the sea level height with its node around the dateline and its center at 120\u00b0W is observed (Fig. 4.3b); the surface zonal current also shows a standing-type oscillation with a basin-wide pattern with its minimum at the boundaries; and SST anomalies (Fig. 4.3a) show a westward propagating tendency. The surface zonal current pattern is very similar to the so-called 'ocean basin mode' (Neelin and Jin, 1993). The ocean basin mode is a natural mode produced by the equatorial trapped waves and their reflection at boundaries: the equatorial 80 Chapter 4. Changes in the leading ENSO modes: Role of surface zonal current (a) SST (mode 1, 58%) (b) SST (mode 2, 5%) Figure 4.1: The first and second modes of the SST (top), sea level height (middle) and surface zonal current (bottom) anomalies obtained from a Combined Principal Component Analysis (CPCA). The percentage of overall variance explained is 36% by the first C P C A mode and 14% by the second mode. The percentage of variance explained for each variable is given in the title of each panel. 81 Chapter 4. Changes in the leading ENSO modes: Role of surface zonal current Figure 4.2: Trajectory plot of the principal components of the first and second PCs for the pre-1980s and the post-1980s. \"x\" and \"*\" mark the mean location in the PC space of the pre-1980s and the post-1980s data, respectively. 82 Chapter 4. Changes in the leading ENSO modes: Role of surface zonal current Kelvin wave travels eastward until it hits the eastern boundary where it is mainly reflected into an equatorially-trapped Rossby wave, which moves westward, reflecting at the western boundary into the Kelvin wave \u2014 thus completing a semi-closed loop by traveling-waves with 8-12 month period in the tropical Pacific basin (with some energy loss at the eastern boundary due to the poleward traveling coastal Kelvin wave; Jin, (1997)). Neelin and Jin (1993) mentioned that this ocean basin mode can be destabilized by the zonal advective feedback (i.e. zonal advection of the mean zonal temperature gradient by the anomalous zonal current). Once the ocean basin mode becomes unstable, it could be changed into the interannual mode, and have larger amplitude than the uncoupled ocean basin mode. An and Jin (2000) showed that the zonal temperature gradient of the mean SST during the 1961-75 is relatively larger than that during 1981-95, and it caused the destabilization of the ocean basin mode during the pre-1980s. In many respects, the leading N L P C A mode of the observations is consistent with such previous studies. The evolution features are quite different for the post-1980s from those for the pre-1980s (Fig. 4.3). The standing structure in the sea level height has turned into eastward propagating patterns (Fig. 4.3e); the basin-wide pattern of the surface current disappears; and the westward moving tendency of SST changes to the eastward moving tendency (e.g., An and Wang 2000). The maximum surface current in the eastern Pacific (Fig. 4.3f) appears when the zonal gradient of the sea level height is at maximum (Fig. 4.3e), indicating that the surface current is mostly due to the local pressure gradient from the zonal mass difference than to the surface wind stress \u2014 the positive SST anomaly and the resulting convective atmospheric heat generate the positive wind stress anomaly in the western part of the center of the SST anomaly, which may enforce the westerly zonal surface current, but the zonal current anomaly is negative as shown here. A n and Jin (2001), under the 83 Chapter 4. Changes in the leading ENSO modes: Role of surface zonal current two-box approximation of the intermediate ENSO model, theoretically proved this point, and mentioned that the two-box approximation could represent the ocean adjustment but not the ocean wave dynamics, thus only the recharge mode is possible. In this regard, the ENSO mode inferred from the N L P C A pattern for the post-shift regime is presumably the recharge mode due to the destabilization of ocean adjustment mode. Note that the lifetime and amplitude for the positive zonal current anomaly and those for the negative zonal current anomaly are not the same (Fig. 4.3f). This asymmetric oscillation was dominant during the post-1980s (An et al., 2005). The phase relationships among the SST, sea level height, and surface zonal current appeared in the N L P C A mode give an idea on the role of the surface current for each decadal period. For the pre-1980s, as shown in Fig. 4.3, the sea level height anomalies over the equatorial eastern Pacific are almost in-phase with the surface zonal current anomalies. Since a positive (negative) surface zonal current anomaly generates a positive (negative) SST anomaly through the warm (cold) surface water advection by the eastward (westward) zonal current, the surface current plays a role in the intensification of this system. For the post-1980s, on the other hands, the surface zonal current anomalies in the equatorial Eastern Pacific lead the sea level height anomalies by about a quarter of a cycle. This suggests that the surface zonal current plays a role in the transition mechanism of this oscillation. Thus, in each climate regime, the role of the surface zonal current in ENSO evolution is different. 4.3 Summary Using the nonlinear principal analysis (NLPCA) technique, we demonstrated that the lead-ing ENSO mode in a physical basis changed since the late 1970s. The ENSO mode during the 84 Chapter 4. Changes in the leading ENSO modes: Role of surface zonal current (a) SST (mode 1, 81%) (d) SST (mode 1, 83%) 150E 180 1S0I 120W 90W 150E 1B0 1S0W 120W 90W (b) SSH (mode 1, 60%) -180 (e) SSH (mode 1, 73%) 130E 180 1S0W 1201 90W -180 150E 180 150W 120W 90V ( c ) Current (mode 1, 47%) 180 -180 (f) Current (mode 1, 52%) 1B0 150E 180 150V 120V 90V -180 150V 120V Figure 4.3: Phase-longitude section of the SST (top panels), sea level height (center panels), and the surface zonal current (bottom panels) anomalies over the equatorial band (5\u00b0S to 5\u00b0N) associated with the first N L P C A mode, (a) to (c) for the pre-1980s and (d) to (f) for the post-1980s. 85 Chapter 4. Changes in the leading ENSO modes: Role of surface zonal current pre-1980s resembles the destabilized 'ocean basin mode', while that during the post-1980s, 'the recharge-mode'. In particular, for the pre-1980s, the surface zonal current mainly acts as an intensifier of the ENSO, while that during post-1980s it plays a role in the transition of the ENSO cycle. 4 . 4 R e f e r e n c e s An S.-L, W. W. Hsieh, and F.-F. Jin, 2005: A nonlinear analysis of the ENSO cycle and its Interdecadal changes. J. Climate, 18, 3229-3239. An S.-L, and F.-F. Jin, 2001: Collective role of thermocline and zonal advective feedbacks in the ENSO mode. J. Climate, 14, 3421-3432. An S.-L, and F.-F. Jin, 2000: An eigen analysis of the Interdecadal changes in the structure and frequency of ENSO mode. Geophy. Res. Lett, 27, 1573-1576. An S.-L, and F.-F. Jin, 2004: Nonlinearity and asymmetry of ENSO. J. Climate, 17, 2399-2412. An S.-L, A. Timmermann, L. Bejarano, F.-F. Jin, F. Justino, Z. Liu, and A. W. Tudhope, 2004: Modeling evidence for enhanced E l Nino-Southern Oscillation amplitude during the Last Glacial Maximum. Paleoceanography, 19, PA4009, DOI: 10.1029\/2004PA001020. An S.-I. and B. Wang, 2000: Interdecadal change of the structure of the ENSO mode and its impact on the ENSO frequency. J. Climate 13, 2044-2055. Carton J., G. Chepurin, X . Cao, and B. Giese, 2000: A simple ocean data assimilation analysis of the global upper ocean 1950-95. Part I: Methodology. J. Phys. Oceanogr., 30, 294-309. 86 Chapter 4. Changes in the leading ENSO modes: Role of surface zonal current Eisenman I., L. Yu, and E. Tziperman, 2005: Westerly wind bursts: ENSO's tail rather than the dog? J. Climate, 18, 5224-5238. Fedorov A. V. and S. G. H. Philander, 2000: Is E l Nino changing?. Science, 228, 1997-2000. J i n F.-F., 1997: An Equatorial Ocean Recharge Paradigm for ENSO Part I: Conceptual Model, J. Atmos. Sci., 54, 811-829. J i n F.-F., and J. D. Neelin, 1993: Modes of interannual tropical ocean-atmosphere inter-action - a unified view. Part I: Numerical results. J. Atmos. Sci., 50, 3477-3503. Neel in J. D., and F.-F. Jin, 1993: Modes of interannual tropical ocean-atmosphere in-teraction - a unified view. Part II: Analytical results in the weak-coupling limit. J. Atmos. Sci., 50, 3504-3522. N i t t a T., and S. Yamada, 1989: Recent warming of tropical sea surface temperature and its relationship to the Northern Hemisphere circulation. J. Meteor. Soc. Japan, 67, 375-383. Philander G. H, 1990: E l Nino, La Nina, and the Southern Oscillation, Academic Press, pp 293. Philander G. H., and A. Fedorov, 2003: Is E l Nino sporadic or cyclic? Annual Review of Earth and Planetary Sciences, 31, 579-594. Timmermann A., F.-F. Jin, and J. Abshagen, 2003: A nonlinear theory for E l Nino bursting. J. Atmos. Sci., 60, 152-165. Trenberth K. E. and J. W. Hurrel, 1994: Decadal atmosphere-ocean variations in the Pacific. Climate Dyn., 9, 303-319. Chapter 4. Changes in the leading ENSO modes: Role of surface zonal current Wang B. and S.-I. An, 2001: Why the properties of E l Nino changed during the late 1970s. Geophys. Res. Lett, 28, 3709-3712. 88 Chapter 5 C h a n g e s i n E N S O a n d a s s o c i a t e d o v e r t u r n i n g c i r c u l a t i o n s f r o m e n h a n c e d g r e e n h o u s e g a s e s 4 5.1 Introduction The global atmospheric and tropical oceanic surface temperatures increased in the past century by > 0.5\u00b0C (Jones and Moberg, 2003; Anthes et al, 2006). This warming is, at least in part, a result of emissions of greenhouse gases (GHGs) from human activities (IPCC, 2001). The climate change caused by the emissions of GHGs from human activities will also change the tropical Pacific climate state. ENSO, the strongest interannual signal in the tropical Pacific, changes its characteristics with the change in the background climatology (Fedorov and Philander, 2001; An and Wang, 2000; Ye and Hsieh, 2006). Our primary interest here is to study how changes in GHGs from the pre-industrial level to the year 2000 level would impact ENSO. There is debate in the climate literature as to whether the low-frequency change of 4 A version of this chapter has been submitted to Journal of Climate for publication. Ye, Z. and Hsieh, W . W . (2006) Changes in E N S O and associated overturning circulations from enhanced greenhouse gases. 89 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases ENSO in the tropical Pacific is generated within the tropics by tropical internal instability (Knutson and Manabe, 1998), or uncoupled atmospheric noise (Thompson and Battisti, 2001; Flugel et al, 2004), or whether it involves the interaction between the extratropics and the tropics (Gu and Philander, 1997; Kleeman et al, 1999; Barnett et al, 1999). The zonal overturning circulation in the equatorial Pacific directly transports water mass and heat energy between the western warm pool and the eastern cool upwelling area in the tropical Pacific, while the shallow meridional overturning circulation conveys water mass and heat energy between tropical Pacific and subtropical Pacific. McCreary and Lu (1994) and Liu (1994) investigated the fundamental dynamics of the shallow meridional overturning circulation using simple layered thermocline models. They postulated that the subtropical water reaches the equator from the northeast subtropics via low-latitude western boundary currents or via an interior path directly linking the northeast subtropics to the central tropics. There is also debate on whether the climate is changing into E l Niho-like warming or La Nina-like cooling in the tropical Pacific. Cane et al (1997) found that the SST had a cooling trend in the eastern tropical Pacific from the 1900-1991 observations. However, the Intergovernmental Panel on Climate Change (IPCC) (2001) found that the trend for SST was E l Nino-like in the tropical Pacific in many models. Meehl and Washington (1996) also noted an E l Niho-like climate change in their coupled general circulation climate model (CGCM) with C 0 2 doubling. This chapter will use the observed SST from 1900-1999 which was constructed using the most recently available International Comprehensive Ocean-Atmosphere Data Set (ICOADS) SST data and improved statistical methods that allow stable reconstruction using sparse data, including a modified historical bias correction for the 1939-1941 period. 90 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases The shallow meridional overturning circulation has been slowing down in the last decades during the period 1950-1999 (McPhaden and Zhang, 2002), though it rebounded during the short period 1998-2003 with respect to the period 1992-1998 (Zhang and McPhaden, 2004). According to Zhang and McPhaden (2006), there appeared to be a strong decadal vari-ability superimposed on a linear weakening trend in the period 1953-2001. They found the shallow meridional overturning circulation in the tropical and subtropical Pacific to be highly anti-correlated with SST anomalies in the central and eastern tropical Pacific in the observational datasets over the period 1953-2001, and explored this relationship in the simulations of 20th century climate from 14 state-of-the-art coupled climate models, though the model results were unclear. The models had difficulty simulating this observed relation because they did not adequately simulate the ocean's response to remote forcing from the subtropical Pacific due to an underestimation of the extratropical winds (Solomon and Zhang, 2006). Using a coupled model, Merryfield and Boer (2005) found that anthro-pogenic forcing may have contributed to the observed slowdown of the shallow meridional overturning circulation. The zonal Walker circulation driven by convection in the western equatorial Pacific and subsidence in the east shows a weakening trend since the mid-19th century due to anthropogenic forcing (Vecchi et al, 2006). Here we will investigate the slowdown in the oceanic zonal and shallow meridional circulations using the multi-model ensemble based on 12 CGCMs in the 4th Assessment Report of the Intergovernmental Panel on Climate Change (IPCC-AR4) . There is much uncertainty in how ENSO will change its characteristics (e.g. the ampli-tude and frequency) under increased GHG in the CGCMs. Guilyardi (2006) assessed the ENSO-mean state-seasonal cycle interactions in 23 coupled ocean-atmosphere models by comparing the pre-industrial control and the stabilized 2xC02 and 4xC02 scenario runs. 91 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases The ENSO amplitude was found to be an inverse function of both the mean trade winds and the relative strength of the seasonal cycle, but the relation was less clear for the ENSO frequency. Van Oldenborgh et al (2005) found very little influence of global warming on ENSO from the results of 17 models. Merryfield (2006) compared the differences in ENSO amplitude, period, and pattern under pre-industrial conditions and under CO2 doubling: The amplitude changes were not strongly related to the magnitude or pattern of surface warming. A narrow (wide) wind stress response was associated with ENSO amplitude decrease (increase). The models exhibited a mean fractional decrease in ENSO period by about 5%, and an increase in the amplitude of SST variations in the central tropical Pacific. On the other hand, Meehl et al (2006) found that ENSO events decreased in magnitude in the future warmer climate (2xC02, 4 x C 0 2 and other scenarios). Collin et al (2005) showed that the most likely scenario is for no large amplitude change toward mean E l Nino or La Niha-like condition in the 20 models submitted to the coupled model inter-comparison project (CMIP) by comparing the 80-year control simulation with fixed CO2 level and the 80-year simulation in which C 0 2 was increased from the control value at a rate of 1% per year compounded. In this chapter, the changes in ENSO under the committed C 0 2 level (i.e. the year 2000 CO2 level) are examined. Finally, we point to the differences between this work and previous studies: While I have also used the datasets from IPCC-AR4, most pre-vious studies trying to detect changes in ENSO used the data from the experiments under C 0 2 doubling (2xC02) or Special Report on Emissions Scenarios A2 (SRES-A2), instead of the PIcntrl and Commit experiments used here. Previous works tended to focus on the magnitude and frequency of ENSO oscillations, while here the focus is on the changes in the nonlinearity and in the tropical zonal and shallow meridional overturning circulations of ENSO. 92 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases ENSO displays considerable asymmetry between its warm phase (El Nino) and its cold phase (La Nina) (An, 2004; Rodgers et al, 2004). The asymmetry is due to the nonlinearity of the ENSO system, as the nonlinear dynamic heating in the tropical Pacific ocean heat budget is essential in producing decadal changes in ENSO nonlinearity and asymmetry (Jin et al, 2003; An and Jin, 2004). Changes in the asymmetry and nonlinearity of ENSO and associated overturning circulations are also addressed in our study. The observed SST and the ensemble model data are described in Sec. 5.2 and Sec. 5.3, respectively. By comparing the ensemble model simulations under the pre-industrial GHG conditions and under the present conditions, we identified climate change in the tropical Pacific (Sec. 5.4), changes in the ENSO SST and zonal wind stress anomalies (Sec. 5.5), and changes in the ENSO ocean circulation (Sec. 5.6). A diagnostic analysis of the ocean surface temperature equation gave some explanation for the enhanced ENSO asymmetry (Sec. 5.7). 5.2 Observed SST The NOAA extended reconstructed SST data (Smith and Reynolds, 2004) were used, with the data divided into the 1900-1949 and 1950-1999 periods. The climatological difference between the two periods suggests that the SST in the eastern-central tropical Pacific has increased by at least 0.3\u00b0C (Fig. 5.1). The climatological difference pattern resembles the E l Nino pattern, with the maximum positive value in the eastern equatorial Pacific. The SST trend during 1900-1999 is also E l Nino-like (Fig. 5.2). Prior to calculating the composite maps for E l Nino and La Nina episodes, a three-month running mean was applied after the climatological seasonal cycle and linear trend were removed in each period. We examined the composite maps of SST anomalies for E l 93 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases Nino and La Nina episodes during each period (Fig. 5.3). E l Nino and La Nina episodes were defined based on the Niho-3.4 index (i.e. the averaged SST anomalies over the region 5\u00b0S-5\u00b0N, 120\u00b0W-170\u00b0W exceeding \u00b1 1 standard deviation). Compared to the 1900-1949 period, a prominent eastward shift of positive SST anomalies appears over the eastern equatorial Pacific during E l Nino in the 1950-1999 period. The maximum anomaly \"H\" has shifted from 120\u00b0W (5.3a) to east of 90\u00b0W (Figs. 5.3c), although the real center of mass of the anomaly pattern does not shift so far. During La Nina, a westward shift of negative SST anomalies is evident in the latter period, i.e. the shift during La Nina is in opposite direction compared to that during E l Nino (5.3b and d). These shifts obviously enhanced the E l Niho-La Nina asymmetry and nonlinearity in the latter period. To derive an index for the asymmetry, the SST composite (Fig. 5.3) for E l Nino and that for La Nina were added together at each grid point, then the root-mean-square (RMS) deviation (averaged over the tropical Pacific) is a measure of the asymmetry, with a zero value indicating the La Nina pattern to be completely symmetrical to the E l Nino pattern. The RMS deviation was 0.11\u00b0C for the period 1900-1949 and 0.19\u00b0C for 1950-1999, confirming the enhanced asymmetry in the latter period. This change in the asymmetry has a longer time scale (50 years) compared to the decadal change in ENSO asymmetry found after the late 1970s (Ye and Hsieh, 2006; Wu and Hsieh, 2003; An, 2004). To find a possible cause for this longer time scale change in ENSO, we will examine the effects of GHG forcing in the following sections. 5.3 M o d e l data Model data from the IPCC-AR4 database were downloaded from the archive hosted by the Program for Climate Model Diagnosis and Intercomparison (PCMDI). The 12 CGCMs 94 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse 150E 180 150W 120W 90W Figure 5.1: (a) The observed mean SST (in \u00b0C) during 1950-1999, and (b) the mean SST of 1950-1999 minus that of 1900-1949, with the 5% and 10% significance levels based on the i-test (von Storch and Zwiers, 1999, p. 115) shown in dark and light gray respectively. 9 5 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases SST Trend (1900-1999) Figure 5.2: The trend in monthly mean SST anomalies for 1900 to 1999 in degrees Celsius per century (from N O A A E R S S T ) . 96 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases Figure 5.3: Composite maps of SST anomalies (in \u00b0C) during E l Nino (left column) and during La Nina (right column), with (a) and (b) for the 1900-1949 regime, and (c) and (d) for 1950-1999. \"H\" and \"L\" mark the location of the highest and lowest values, respectively, and the 5% and 10% significance levels based on the i-test are shown in gray. 97 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases used here are the CCCMA-CGCM3-1 , CNRM-CM3, CSIRO-MK3-0, GFDL-CM2-0, GISS-MODEL-E-R, IAP-FGOALS1-0-G, INMCM3-0, IPSL-CM4, MIROC3-2-MEDRES, MIUB-ECHO-G, MRI-CGCM2-3-2A and U K M O - H A D C M 3 (Table 5.1). The documentation for the models is available on the website http: \/\/www-pcmdi. l l n l . gov\/ipcc\/model_documentation\/ ipcc_model_documentation.php. We consider the pre-industrial control runs (PIcntrl) and the Commit runs (Commit). PIcntrl is the pre-industrial climate simulation with GHG-induced forcing fixed at the level of year 1850, whereas Commit simulates committed climate change using the forcing (GHG + aerosols) at year 2000 levels. For each simulation run, 100 years of data were taken from each of the 12 models. 98 Model Institution (country) Atmosphere resolution Ocean resolution CGCM3.1 (T47) C C C M A (Canada) T47L31 1.85\u00b0xl.85\u00b0L29 CNRM-CM3 Meteo-France\/CNRM (France) T63L45 2\u00b0x0.5\u00b0L31 CSIRO-MK3.0 CSIRO (Australia) T63L18 1.875\u00b0 x0.84\u00b0L31 GFDL-CM2.0 GFDL (USA) 2.5\u00b0x2\u00b0L24 l\u00b0x0.33\u00b0L50 GISS-ER NASA\/GISS (USA) 5\u00b0x4\u00b0L20 5\u00b0x4\u00b0L13 IAP-FGOALS-gl.O L A S G \/ I A P (China) T42L26 l\u00b0xl\u00b0L33 INM-CM3.0 I N M C M (Russia) 5\u00b0x4\u00b0L21 2.5\u00b0x2\u00b0L33 IPSL-CM4 IPSL (France) 2.5\u00b0x3.75\u00b0L19 2\u00b0x0.5\u00b0L31 MIROC3.2 (medres) C C S R \/ N I E S \/ F R C G C (Japan) T42L20 1.4\u00b0x0.5\u00b0L43 MIUB-ECHO-g M I U B \/ (Germany) T30 L19 T42L20 MRI-CGCM2.3.2A MRI (Japan) T42L30 2.5\u00b0x0.5\u00b0L23 UKMO-HadCM3 HadleyCentre (UK) 3.75\u00b0x2.5\u00b0L19 1.25\u00b0xl.25\u00b0L20 Table 5.1: The models used in present Chapter, including their configurations near the equator. Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases For all models, the tropical Pacific SST and zonal wind stress (WS) were analyzed. The zonal current in 2\u00b0S-2\u00b0N and the meridional current in 30\u00b0S-30\u00b0N from the sea surface to 450 m depth were also analyzed. The meridional currents averaged over 100 m to 300 m depth were used to detect the shallow meridional overturning circulation in the region 30\u00b0S-30\u00b0N in the Pacific. The temperature and vertical velocity in the equatorial Pacific region (5\u00b0S-5\u00b0N, 0-50m) were also used in the diagnostic analysis. For each model, the monthly mean SST, WS, ocean upper level velocities were interpo-lated onto a 5\u00b0 longitude by 4\u00b0 latitude regular grid (identical to that used in the GISS-ER model) using bi-linear interpolation. Anomalies were computed with respect to the mean seasonal cycle in each specific model simulation. 5 . 4 Climate change simulated In this section, we focus on the human-induced impacts on the SST in the tropical Pacific and the three large-scale circulations (Fig. 5.4) in the upper tropical and subtropic Pacific from 30\u00b0S to 30\u00b0N area. The three circulations are shown in the schematic figure with shade area standing for the northern Pacific water in upper layer (Fig. 5.4). The climate change induced by the anthropogenic forcing were computed by subtracting the multi-model ensemble mean in the PIcntrl experiment from that in the Commit exper-iment. The warm pool in the western tropical Pacific and the cool tongue in the eastern tropical Pacific were both reproduced well in the mean climate in the Commit experiment (Fig. 5.5a), although there is a cool model bias (up to 1\u00b0C) in the east-central equatorial Pacific Ocean and an exaggerated westward extent of the cold tongue (cf. Fig. 5.1a). These biases have been recognized as a common problem in most of the state-of-the-art coupled models which have a double Inter-Tropical Convergence Zone (ITCZ) (Zhang and Wang, 100 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse 2006; Dai, 2006). Xie et al (2006) found that the interaction between the low-cloud and SST over the eastern tropical Pacific played an important role in generating and maintaining the single ITCZ in their numerical model. The difference between the means of the two exper-iments shows nearly 1\u00b0C warming over all of tropical Pacific (Fig. 5.5b), with the warming pattern vaguely resembling the E l Nino pattern which has the maximum centered in the eastern-central equatorial Pacific. While every model revealed an E l Nino-like pattern of SST change (Fig. 5.6), there is considerable variability in the exact location of, e.g., the maximum warming center which ranged from the eastern to the western equatorial Pacific. Therefore the standard deviation of the SST change in Fig. 5.5c is large. The large spatial spread (in tropical western\/eastern Pacific) in the maximum SST anomalies in the simu-lated ENSO behavior was also found by Saji et al (2006), in the simulated twentieth-century data from 17 models in IPCC-AR4 dataset. Their results showed that the models exhibited poor ability in simulating ENSO's seasonality and in the phase locking to the annual cycle. Only one of 17 models appears to resemble to the observed temporal evolution. The difference between the means of the two experiments also showed a positive zonal WS pattern in the western tropical Pacific (Fig. 5.5e), which also resembles the zonal WS anomaly pattern found during E l Nino (shown later). The enhanced GHG forcing decreased the mean easterlies (Fig. 5.5d) along the western equatorial Pacific (Fig. 5.5e). This agrees with the finding that the Walker circulation has been slowing down in the tropical atmosphere since the mid-19th century (Vecchi et al, 2006), as the trade wind blowing across the tropical Pacific from east to west is associated with the bottom branch of the zonal Walker circulation. Fig. 5.5e also showed a strengthening of the easterlies south of 10\u00b0S. The inter-model variability of the zonal WS change (Fig. 5.5f) revealed that the models were most consistent with each other along the equatorial belt. 101 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse Meridional Overturning 30N West East (Equator) Figure 5.4: Schematic view of upper ocean (0-500 m depth) circulations in the Pacific from 0\u00b0 to 30\u00b0N. 102 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse 20N ION 20S 150E 180 1501 120W 90W 20N ION EQ 10S 20S (d) Zonal WS (Commit) ^ ~ _ 5 \u2014 \" -4\u2014 3 \u2014 -1-' --Z-,-1 -3\u2014. 150E 180 150W 120W 90W 20N ION (b) SST (Commit -PIcntr l ) 20N ION 10S 20S (e) Zonal WS (Commit -PIcntr l ) \" : \u2022 v Oi 20N ION 150E 160 150W 1201 901 (c) Std of SST difference 20S+* l 1 ' 1 \u2022 r- 1 \u2014 1 20N 150E 180 (f) Std of zonal WS difference 205 150E 180 150W 120W 20S-I 90W 150E 1B0 150W 120W 90W Figure 5.5: The multi-model ensemble mean for the Commit runs, its difference from the mean for the PIcntrl runs, and the standard deviation of the difference: (a) mean SST (\u00b0C), (b) the difference in the mean SST, (c) the standard deviation of the SST difference, (d) the mean zonal WS (0.01 Pa), (e) the difference in the mean zonal WS, and (f) the standard deviation of the WS difference. The 5% and 10% significance levels from the i-test are shown in gray in (b) and (e). 103 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases 20N-1 CCCMA-CGCM3.1 (Commit) CCCMA-CGCM3.1 (Commit-PIcntrl) 150E 180 190V 120V CNRM-CM3 (Commit) 90V \"1501 160 150V 120V 90V CNRM-CM3 (Commit-PIcntrl) 150E 160 150V 120V CSIRO-MK3.0 (Commit) 90V ~~I50E 180 150V 120V 90V CSIRO-MK3.0 (Commit-PIcntrl) 20S 150E 20N 160 150V 120V GFDL\u2014CM2.0 (Commit) 2031 90V 150E 180 150V 120V 90V GFDL-CM2.0 (Commit-PIcntrl) 160 150V 120V 90V 150E IAP-FGOALS-gl.O (Commit) 120V 20N EQ 7 ^ Q' 6 _js^==0 0.9_\u2014^\"\"^ 180 150V 90V Figure 5.6: SST mean (\u00b0C) for the Commit runs (left) and its difference from the mean in the PIcntrl runs (right) for 12 individul models. 104 er 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases INMCM3.0 (Commit) 1S0E 160 1501 120\u00a5 IPSL-CM4 (Commit) 20N BO. 20S INMCM3.0 (Commit-PIcntrl) \"'11 901 150E 180 1901 120* 90* IPSL-CM4 (Commit-PIcntrl) 20N MIROC3.2(medres) (Commit) 90V 150E 180 160V 120V 90V : o M MIB0C3.2(medres) (Commit-PIcntrl) 20N 150E 180 150V 120V 90V 1S0E 180 1S0V 120V 90V MIUB-ECHO-g (Commit) _ m MIUB-ECHO-g JCommit-PIcntrl) 90V 150E 20N MRI-CGCM2.3.2A (Commit) 20Ni MRI-CGCM2.3.2A (Commit-PIcntrl) 20S 20N 901 150E 180 150V 120V 90V UKMO-HadCM3 (Commit-PIcntrl) 90V Figure 5.6 (Continued) 105 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse The zonal current along the equator (averaged over 2\u00b0S to 2\u00b0N) showed general con-sistency among the 12 individual models (Fig. 5.7). Although the simulated equatorial undercurrent along the thermocline is not as strong as the observed values (maximum is up to 100 cm s _ 1 in the observations) in most of the models in the subsurface layer, all models have simulated the equatorial undercurrent (left panels of Fig. 5.7). The differences between the Commit run and PIcntrl run also showed general consistency among the 12 models (right panels of Fig. 5.7). The zonal current in all models have eastward zonal cur-rent anomalies near the surface and westward zonal current anomalies near the subsurface. In the other words, the zonal equatorial overturning circulations tend to slow down in all 12 models in the Commit runs compared to the PIcntrl runs. In the multi-model ensemble mean of the zonal current, the equatorial undercurrent along the thermocline in the subsurface layer (Fig. 5.8a) is not as strong as observed cur-rent (with a maximum of about 100 cm s _ 1 ) . The overturning equatorial circulation has weakened as shown in the difference between the means from the Commit and PIcntrl ex-periments (Fig. 5.8b), while the undercurrent has shifted upward in the Commit runs, as Vecchi et al. (2006) have also noted an upward trend in the undercurrent position in the GFDL-CM2.1 simulation data during 1861-2000. As the observed Walker circulation in the tropical atmosphere has been weakening in recent decades, Vecchi et al. felt that global warming was a likely causative factor in the weakening of the wind pattern. Here the sim-ulated results confirm that the overturning oceanic circulation in the equatorial Pacific has also been weakening between PIcntrl and Commit. Next consider the meridional currents which are important for interactions between the tropics and the extratropical Pacific, as studies have shown that extratropical effects influence the tropical ENSO (Gu and Philander, 1997; Kleeman et al, 1999). Observations 106 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases indicated that the cool water in the subtropical Pacific, especially in the eastern area, is subducted down and moved to the tropics by the shallow meridional overturning circulation in the upper Pacific (McPhaden and Zhang, 2004). The Pacific overturning circulation can be presented by the meridional oceanic velocity in the tropical and subtropical Pacific. The difference in the mean subsurface current (Fig. 5.8d) between the Commit and PIcntrl runs shows the Pacific shallow meridional overturning circulation (Fig. 5.8c) tending to slow down, especially in the western tropical Pacific. The recent observed evidence suggests that there does exist such a slowdown of the shallow meridional overturning circulation in the upper Pacific Ocean since the 1970s (McPhaden and Zhang, 2002). Although McPhaden and Zhang (2004) found that the shallow meridional overturning circulation rebounded during 1998-2003, their latest result showed that there was a slowdown trend overlying the decadal variability in the shallow meridional overturning circulation during 1953-2000 (Zhang and McPhaden, 2006). Figs. 5.8c and d indicate that this slowdown in the shallow meridional overturning circulation is at least partially due to the increased G H G forcing. Fig. 5.8d also shows that the slowdown of the overturning circulation occurs along the interior path. In Fig. 5.8d, the changes in the 100-300 m mean meridional velocity within 7\u00b0N\/S are likely in response to changes in equatorial zonal wind (Fig. 5.5e), while the changes in the shallow meridional overturning circulation beyond the equatorial belt are associated with wind changes outside the equatorial belt (Nonaka et al, 2002). The difference in the mean between the Commit runs and PIcntrl runs also show that the GHG-induced changes are E l Nino-like patterns in the tropical Pacific SST. In the next two sections, E l Nino-like patterns are also found in the changes of the zonal WS, zonal overturning circulation and shallow meridional overturning circulation. The observed E l Niho-like pattern in SST (Fig. 5.1b) therefore appears at least in part contributed by the 107 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse GHG forcing. 5 . 5 Changes in the E N S O SST and zonal W S The data from GISS-ER and IAP-FGOALS-gl.O are excluded in the analysis from now on due to no ENSO variability in the GISS-ER model (Guilyardi, 2006) and the very small seasonal cycle with the unrealistic regular interannual cycle in the IAP-FGOALS-gl.O model (Guilyardi, 2006; Van Oldenburg et al, 2005). To identify E l Nino and La Nina episodes in the model data, the simulated SST anomalies in the tropical Pacific (150\u00b0E-80\u00b0W, 20\u00b0S-20\u00b0N) were smoothed by a 25-month running mean, [as the period of the simulated ENSO was greater than 2 years in all 10 coupled models (Guilyardi, 2006)], then the leading PC of the smoothed SST anomalies was defined as a proxy ENSO index for each model run, with E l Nino\/La Nina defined based on the proxy ENSO index exceeding \u00b1 1 standard deviation. Composite maps of E l Nino and La Nina for various variables (without the 25-month running mean smoothing) were computed for each individual model. The multi-model ensemble mean for each variable was then calculated from the individual composite maps. ENSO is an interactive thermodynamic system between the atmosphere and ocean, where the atmosphere dynamically forces the ocean by the surface WS, while the ocean thermally drives the atmosphere by heating\/cooling. Multi-model ensemble means of com-posite maps for SST and zonal WS anomalies over the tropical Pacific are shown in Fig. 5.9. A pronounced shift in the zonal location of positive SST anomalies during E l Nino occurred in the Commit runs relative to the PIcntrl runs (Figs. 5.9a and c). The shift is in the same direction as observed (Figs. 5.3a and c), but there is model bias since the model SST anomalies are located further west than the observed during E l Nino, as noted by Capotondi 108 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases CCCMA-CGCM3.1 (Commit) CCCMA-CGCM3.1 (Commit-PIcntrl) 190V 120V 90V 190E 160 190V 120V 90V CNRM-CM3 (Commit-PIcntrl) 160 190V 120V CSIRO-MK3.0 (Commit) 90V 150E 180 190V 120V CSIBO-MK3.0 (Commit-PIcntrl) 90V 190E -100 -200 -300 ^- 6^ ^ f > J , ^ 4\u2014\u2014\u2014 ^ \/ 160 190V 120V 90V GISS-ER (Commit-PIcntrl) 90V 150E 160 190V 120V IAP-FGOALS-gl.0 (CjmmU-PIcntrl) Figure 5.7: The zonal mean velocity (averaged over 2\u00b0S and 2\u00b0N in cm s l ) for the Commit runs (left) and its difference from that in the PIcntrl runs (right) for 12 individual models. 109 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases INMCM3.0 (Commit) INMCM3.0 (Commit-PIcntrl) 160 150V 120V IPSL-CM4 (Commit) gov \/ -1\" s\u2014 180 150V 120V MIROC3.2(medres) (Commit 90V 150E 100 150V 120V 901 MIROC3.2(medres) (Commit-PIcntrl) 160 150V 120V MIUB-ECHO-g (Commit) | -400 90V 150E 160 150V 120V 90V MIUB-ECHO-g (Commit-PIcntrl) 160 150V 120V MRI-CGCM2.3.2A (Commit) 90V 150E 160 150V 120V 90V UKMO-HadCM3 (Commit-PIcntrl) 90V 150E 90V Figure 5.7 (Continued) 110 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases Figure 5.8: The multi-model ensemble mean of velocity (in cm s _ 1 ) for the Commit runs and its difference from that in the PIcntrl runs: (a) the mean zonal velocity (averaged over 2\u00b0S and 2\u00b0N), (b) the difference in the mean zonal velocity, (c) the mean meridional velocity averaged over 100 - 300 m depth, and (d) the difference in the mean meridional velocity. The 5 % and 10% significance levels from the i-test are shown in gray in (b) and (d). I l l Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse et al (2006). During La Nina, the shift of the negative SST anomalies is undistinguished between the Commit and PIcntrl runs (Figs. 5.9b and d), missing the westward shift of anomalies in the observational data (Figs. 5.3b and d). The RMS index of the asymmetry between E l Nino and La Nina was computed to be 0.045 \u00b0C averaged for the PIcntrl runs and 0.057 \u00b0C for the Commit runs. Although these values are considerably weaker than the observed values given in Sec. 5.2, they indicate the increase of ENSO asymmetry with increased GHG. For the zonal WS during E l Nino, there was eastward shift and strengthening of the westerly anomaly center along the equator (as indicated by \"H\") in the Commit runs relative to the PIcntrl runs (Figs. 5.9e and g), while during La Nina, the easterly WS anomalies along the equator showed no obvious shift (5.9f and h). The contrast in the WS shift between E l Nino and La Nina meant that the asymmetry and nonlinearity of ENSO were enhanced in the Commit runs. The RMS index of asymmetry for the zonal WS was 0.078 (in O.OlPa) for the PIcntrl runs and 0.100 (in 0.01 Pa) for the Commit runs. The change in the asymmetry and nonlinearity of ENSO in these experiments is similar to the observed change after the late 1970s on the decadal time scale (Wu and Hsieh, 2003; An, 2004). 5 . 6 C h a n g e s i n t h e E N S O o c e a n c i r c u l a t i o n In the vertical section showing zonal currents along the equator, during E l Nino the west-ward subsurface zonal current anomaly strengthened slightly and \"L\" shifted eastward and upward in the Commit runs relative to the PIcntrl runs (Figs. 5.10a and c); in contrast, during La Nina, \"H\" in the eastward subsurface zonal current anomaly shifted slightly to the west (Figs. 5.10b and d). Thus the asymmetry in the zonal current anomalies between El Nino and La Nina strengthened in the Commit runs relative to the PIcntrl runs, as the 112 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases Figure 5.9: Multi-model ensemble mean of SST (\u00b0C) and zonal WS (0.01 Pa) composites for E l Nino (left panels) and La Nina (right panels), with (a), (b), (e), (f) for the PIcntrl runs and (c), (d), (g), (h) for the Commit runs. \"H\" and \"L\" mark the location of the highest and lowest values, respectively, and the 5% and 10% significance levels from the t-test are shown in gray. 113 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases Figure 5.9 (Continued) 114 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases Figure 5.10: Vertical section along the equator showing the multi-model ensemble mean of zonal velocity anomalies (cm s _ 1 ) (averaged between 2\u00b0S and 2\u00b0N) for the composite E l Nino (left panels) and La Nina (right panels), with (a) and (b) for the PIcntrl runs, and (c) and (d) for the Commit runs. The 5% and 10% significance levels from the t-test are shown in gray. RMS index of asymmetry was 0.77 cm s - 1 for the PIcntrl runs and 1.01 cm s _ 1 for the Commit runs. In both the E l Nino and La Nina composites (Fig. 5.11), the subsurface meridional currents averaged between 100 - 300 m depth clearly showed the interior path (McCreary and Lu, 1994, Liu, 1994) running from the eastern subtropical Pacific Ocean to the central tropical Pacific Ocean. In the equatorial area, the location of \"L\" (the strongest southward current anomaly) during E l Nino was found to have shifted eastward by 10\u00b0 in the Commit 115 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases runs relative to the PIcntrl runs (Figs. 5.11a and c), whereas during La Nina, \"H\" showed no obvious shift between Figs. 5.11b and d. The RMS index of asymmetry (computed between 5\u00b0S-5\u00b0N) was 0.20 cm s _ 1 for the PIcntrl runs and 0.32 cm s _ 1 for the Commit runs. Hence the asymmetry in the subsurface meridional current anomalies between E l Nino and La Nina was enhanced in the Commit runs relative to the PIcntrl runs. 5 . 7 Diagnostic analysis of the surface temperature equation The ocean surface temperature anomaly equation can be written as where T, u, v, w are, respectively, anomalies of the SST, zonal, meridional and vertical ocean velocities with respect to the climatological mean variables (indicated by the overbar), and \"residuals\" include the SST diffusion, subgrid-scale effects, etc. (An and Jin, 2004). The first bracketed term on the right hand side of (5.1) represents the linear dynamic heating (LDH) in the surface ocean, and the second bracketed term, the nonlinear dynamic heating Since the vertical ocean velocities of the models MIUB-ECHO-g and UKMO-HadCM3 were not available, these two models were excluded in the subsequent analysis. The tem-perature equation was applied to the surface ocean layer, where the average of the SST and velocities values over 0-50 m depth was used to represent the surface layer. The vertical temperature gradient was computed from the difference between the temperature at 0 m and that at 50 m depth. Averaged over time, the SST tendency contributed by the L D H terms is very small compared to the contribution from the N D H terms. Among the three (5.1) (NDH). 116 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse Figure 5.11: Multi-model ensemble mean of subsurface meridional velocity anomalies (cm s _ 1 ) (averaged over 100-300 m depth) for the composite E l Nino (left panels) and La Nina (right panels), with (a) and (b) for the PIcntrl runs, and (c) and (d) for the Commit runs. The 5% and 10% significance levels from the t-test are shown in gray. 117 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases NDH terms, the vertical term - w f j is dominant (Figs. 5.12e and g; f and h). The NDH is about 30% stronger in the Commit runs than in the PIcntrl runs (Figs. 5.12f or h). The enhanced NDH in the Commit runs does not correspond to enhanced upwelling, instead the weakened upwelling (Figs. 5.12c and d) are consistent with the finding in Sec. 5.4 that the equatorial overturning circulation weakened in the Commit runs. The multi-model ensemble composites during E l Nino and La Nina for the NDH in the ocean surface-layer are shown in Fig. 5.13. During both E l Nino and La Nina, the NDH anomalies are positive along the equatorial Pacific (Fig. 5.13). This feature of the NDH anomalies is known to cause the asymmetry between E l Nino and La Nina, as the warm NDH anomalies enhance the warm E l Nino events and weaken the cool La Nina events (An and Jin, 2004). There is very substantial strengthening of the positive NDH anomaly in the equatorial belt and an eastward shift by about 25\u00b0 in the anomaly center (marked by \"H\") in the Commit runs during E l Nino (Figs. 5.13a and c). However, during La Nina, the changes in the NDH anomalies are relatively small between the Commit runs and PIcntrl runs (Figs. 5.13b and d). That the large changes in the magnitude and position of the NDH anomalies between Commit and PIcntrl runs only occurred during E l Nino reflects the nonlinear nature of the change in the NDH under enhanced GHG. The NDH anomalies are mainly caused by the vertical NDH anomalies ( \u2014 ) , as can be seen from the similarity between the total NDH anomaly (Figs. 5.13c and d) and the vertical NDH anomaly (Figs. 5.14a and c) during E l Nino and La Nina. The substantial change in the vertical NDH anomalies during E l Nino in the Commit runs relative to the PIcntrl runs (Figs. 5.14b) is caused by the enhanced vertical temperature gradient anomalies (Figs. 5.14f) and the eastward shift of downwelling anomalies in the central equatorial Pacific (Figs. 5.14i and j). However, during La Nina, the change in the vertical NDH anomalies 118 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse (a) Tz (Commit) (b) Tz (Commit-PIcntrl) 150E 180 1501 120V 90V 150E 180 150V 120V 90V (e) Vertical NDH (Commit) (f) Vertical NDH (Commit-PIcntrl) 150E 180 1501 120V 90V 150E 180 150V 120V 90f (g) NDH (Commit) (h) NDH (Commit-PIcntrl) Figure 5.12: Multi-model ensemble mean of vertical temperature gradient (\u00b0C\/50m, a and b), vertical velocity (10 - 6 m s - 1 , c and d), vertical nonlinear dynamic heating anomalies (\u00b0C month - 1 , c and f) and total nonlinear dynamic heating anomalies (\u00b0C month - 1 , g and h) for the Commit runs (left panels) and the change from the PIcntrl runs (right panels), in the ocean surface-layer. The 5% and 10% significance levels from the t-test are shown in gray. 119 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse Figure 5.13: Multi-model ensemble mean of surface-layer nonlinear dynamic heating (\u00b0C month - 1) for the composite E l Nino (left panels) and La Nina (right panels), with (a), (b) for the PIcntrl runs, and (c), (d) for the Commit runs. The 5% and 10% significance levels from the t-test are shown in gray. -> nn Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse is relatively minor (5.14c and d). The reason for the east-west dipolar pattern seen in | ^ during E l Nino in Figs. 5.14e and f is that most of models have a maximum warm SST anomalies in the western or central equatorial Pacific due to model bias. 5.8 C o n c l u s i o n The E l Niho-like warming in the tropical Pacific seen in the observed SST record (1900-1999) was confirmed by the multi-model results comparing the climate from year 2000 GHG(+aerosol) forcing (Commit) with the pre-industrial climate (PIcntrl). Both the equa-torial zonal overturning circulation and the shallow meridional overturning circulation weak-ened with increased GHG. The idea that the GHGs might play a role on the slowdown of the zonal and shallow meridional overturning circulations received strong support from the multi-model ensemble means. This conclusion also matches the finding that the Walker circulation in the atmosphere has been weakening since the mid-19th century due to the anthropogenic forcing (Vecchi et al, 2006). In the observed SST data, the positive anomalies during E l Nino were located further east in the 1950-1999 composite than in the 1900-1949 composite, whereas the negative anomalies during La Nina were shifted in the opposite direction. From the model data, the positive SST anomalies during E l Nino also shifted eastward in the Commit runs relative to the PIcntrl runs, whereas the negative SST anomalies during La Nina were not shifted zonally. Hence both the model and observed SST results were consistent with increasing asymmetry between El Nino and La Nina as GHG increased. The increase in the asym-metry is associated with an increase in the nonlinearity of ENSO (Jin et al, 2003; An and Jin, 2004), with the nonlinear dynamical heating (NDH) terms producing the asymmetry between E l Nino and La Nina. A diagnostic analysis of the ocean surface temperature 121 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse (a) Vertical NDH (El Nino. Commit) (b) Vertical NDH (El Nino, Commit-PIcntrl) 150E 180 1901 1201 90f (c) Vertical NDH (La Nina. Commit) 150E ISO 1901 1201 90V (d) Vertical NDH (La Nina, Commit-PIcntrl) 190E 100 1901 120V 90W (e) Tz (El Nino. Commit) 190E 100 150V 120V 90V f) Tz (El Nino, Commit-PIcntrl) '''''' ''\/,*,'''!' '\"\"t.\"*'' l'1'1 '' ; ; ' ( 1 5 . 0 8 , 2 - > \" \" \" j ' ' \/ \/ \u2022 EQ '.2 ( . . . . .\"\"i ...... \u2022' 150E 180 1S0V 120V 90V (h) Tz (La Nina. Commit-PIcntrl) 150E 180 190V 1201 901 (1) w \u00bb l e 6 (El Nino, Commit) 150E 180 1501 1201 901 (j) w*le6 (El Nino. Commit-PIcntrl) 150E 180 1501 1201 901 (k) w' leg (La Nina. Commit) 150E 5N (1) w \u00bb l e 6 (La Nina, Commit-PIcntrl) o(y wtf < \u00a3 ^ o < 1501 1201 Figure 5.14: Multi-model ensemble mean of the (a-d) surface vertical nonlinear dynamic heating (\u00b0C month - 1 ), (e-h) vertical temperature gradient anomalies (\u00b0C\/50m), and (i-1) vertical velocity anomalies (10~6 m s _ 1 ) during E l Nino and La Nina for the Commit runs (left panels) and the difference between the Commit and PIcntrl runs (right panels). The 5% and 10% significance levels from the t-test are shown in gray. 122 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases tendency equation revealed the vertical NDH term to be most dominant, with the positive vertical NDH anomalies in the equatorial Pacific enhanced substantially in the Commit runs during E l Nino, caused by the increased vertical temperature gradient anomalies and the eastward shift of downwelling anomalies. Under increased GHG, the enhanced positive NDH anomalies during E l Nino when time-averaged over the whole record would change the SST mean state by an E l Nino-like pattern, which could override other mechanisms trying to induce a La Nina-like change in the mean state (Cane et al., 1997). For the simulated equatorial zonal circulation associated with ENSO, the zonal under-current anomalies during E l Nino strengthened slightly and shifted eastward in the Com-mit runs relative to the PIcntrl runs, whereas the undercurrent anomalies during La Nina shifted slightly westward. For the simulated shallow meridional overturning circulation in the Pacific, the interior path from the eastern subtropical Pacific to the central equatorial region was clearly seen during both E l Nino and La Nina. The enhanced G H G simulations showed that the asymmetry between El Nino and La Nina was enhanced in the anomalies of the SST, the zonal wind stress, the equatorial undercurrent and the shallow meridional overturning circulation. 5.9 References A n S.-L, 2004: Interdecadal changes in the E l Nino La-Nina asymmetry. Geophys. Res. Lett, 31, L23210. DOI:10.1029\/2004GL021699. A n S.-L, and F.-F. Jin, 2004: Nonlinearity and asymmetry of ENSO. J. Climate, 17, 2399-2412. A n S.-I. and B. Wang, 2000: Interdecadal change of the structure of the ENSO mode and 123 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases its impact on the ENSO frequency. J. Climate 13, 2044-2055. Anthes R. A., R. W. Corell, G. Holland, J. W. Hurrell, M . C. MacCracken, and K. E. Tren-berth, 2006: Hurricanes and Global Warming-Potential Linkages and Consequences. Bull. Amer. Meteor. Soc, 87, 623-628. Barnett T. P., D. W. Pierce, M . Latif, D. Dommenget and R. Saravanan, 1999: Inter-decadal interactions between the tropics and midlatitudes in the Pacific basin. Geo-phys. Res. Lett., 26, 615-618. DOI: 10.1029\/1999GL900042. Cane M . A. , A. C. Clement, A. Kaplan, Y . Kushnir, D. Pozdnyakov, R. Seager, S. E. Zebiak and R. Murtugudde, 1997: Twentieth-century sea surface temperature trends. Science, 275, 957-960. Capotondi A. , A. T. Wittenberg, and S. Masina, 2006: Spatial and temporal structure of tropical Pacific interannual variability in 20th century coupled simulations. Ocean Modeling, in press. Coll ins M . and the CMIP Modelling Groups, 2005: E l Nino- or La Nina-like climate change?. Climate Dyn., 24, 89-104. D a i A. , 2006: Precipitation characteristics in eighteen coupled climate models. J. Climate, 19, 4605-4630. Fedorov A. V. and S. G. H. Philander, 2001: A stability analysis of tropical ocean-atmosphere interactions: Bridging measurements and theory for E l Nino). J. Climate, 14, 3086-3101. Fl i ige l M . , P. Chang and C. Penland, 2004: The role of stochastic forcing in modulating ENSO predictability. J. Climate, 17, 3125-3140. 124 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases G u D. and S. G. H. Philander, 1997: Interdecadal climate fluctuations that depend on exchanges between the tropics and extratropics. Science, 275, 805-807. Gui lyard i E., 2006: E l Nino-mean state-seasonal cycle interactions in a multi-model en-semble. Climate Dyn., 26, 329-348. DOI: 10.1007\/s00382-005-0084-6. I P C C TAR (Intergovernmental Panel on Climate Change, Third Assessment Report), 2001: Climate Change 2001: The Scientific Basis, J. T. Houghton, Y . Ding, D. J. Griggs, M . Noguer, P. J. van der Linden, X . Dai, K. Maskell, and C. A. Johnson (Eds.), Cambridge Univ. Press. J i n F.-F., S.-L An, A. Timmermann and J. Zhao, 2003: Strong E l Nino vents and nonlinear dynamical heating. Geophys. Res. Lett. 30, 1120. 01:10.1029\/2002GL016356. Jones P. D., and A. Moberg, 2003: Hemispheric and large-scale surface air temperature variations: an extensive revision and an update to 2001. J. Climate, 16, 206-223. Kleeman R., J. P. McCreary and B. A. Klinger, 1999: A mechanism for generating ENSO decadal variability. Geophys. Res. Lett., 26, 1743-1746. Knutson T. R. and S. Manabe, 1998: Model assessment of decadal variability and trends in the tropical Pacific Ocean. J. Climate, 11, 2273-2296. L i u Z., 1994: A Simple Model of the Mass Exchange between the Subtropical and Tropical Ocean. J. Phys. Oceanogr., 24, 1153-1165. M c C r e a r y J. P., and P. Lu, 1994: Interaction between the Subtropical and Equatorial Ocean Circulations: The Subtropical Cell. J. Phys. Oceanogr., 24, 466-497. M c P h a d e n M . J., and D. Zhang, 2002: Slowdown of the meridional verturning circulation in the upper Pacific Ocean. Nature, 415, 606-608. 125 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse gases M c P h a d e n M . J. and D. Zhang, 2004: Pacific Ocean circulation rebounds, Geophys. Res. Lett, 31, L18301. DOI: 10.1029\/2004GL020727. Meehl G. A., H. Teng, and G. W. Branstator, 2006: Future changes of E l Nino in two global coupled climate models. Climate Dyn., DOI: 10.1007\/s00382-005-0098-0. Merryf ie ld W. J., and G. J. Boer, 2005: Variability of upper Pacific Ocean overturning in a coupled climate model. J. Climate, 18, 666-683. Merryf ie ld W. J., 2006: Changes to ENSO under C02 doubling in a multi-model ensem-ble. J. aim., 19, 4009-4027. Nonaka M . , S.-P. Xie, and J.P. McCreary, 2002: Decadal variations in the Subtropical Cells and equatorial Pacific SST. Geophys. Res. Lett, 29, 10.1029\/2001GL013676. Rodgers K. B., P. Friederichs and M . Latif, 2004: Tropical Pacific decadal variability and its relation to decadal modulations of ENSO. J. Climate, 17, 3761-3774. Saji N . H., S.-P. Xie, and T. Yamagata, 2006: Tropical Indian Ocean variability in the IPCC twentieth-century climate simulations. J. Climate, 19, 4397-4417. Smi th T. M . and R. W. Reynolds, 2004: Improved extended reconstruction of SST, 1854-1997). J. Climate 17, 2466-2477. Solomon A. and D. Zhang, 2006: Pacific Subtropical Cell variability in coupled climate model simulations of the late 19th-20th Century. Ocean Modelling, in press. Thompson C. J. and D. S. Battisti, 2001: A linear stochastic dynamical model of ENSO. Part II: Analysis. J. Climate 14, 445-466. 126 Chapter 5. Changes in ENSO and associated overturning circulations from enhanced greenhouse V a n Oldenborgh G.J., S.Y. Philip, and M . Collins, 2005: E l Nino in a changing climate: a multi-model study. Ocean Sci., 1, 81-95. Vecchi G. A., B. J. Soden, A. T. Wittenberg, I. M . Held, A. Leetmaa and M . J. Harrison, 2006: Weakening of tropical Pacific atmospheric circulation due to anthropogenic forcing. Nature, 44, 73-76. von Storch H. , and F. W. Zwiers, 1999: Statistical Analysis in Climate Research. Cam-bridge Univ. Press, 484 pp. W u A. and W. W. Hsieh, 2003: Nonlinear interdecadal changes of the E l Nino-Southern Oscillation. Climate Dyn., 21, 719-730. DOI: 10.1007\/s00382-003-0361-l. X i e S.-P., T. Miyama, Y . Wang, H. Xu, S.P. de Szoeke, R.J . Small, K . J . Richards, T. Mochizuki and T. Awaji, 2006: A regional ocean-atmosphere model for eastern Pacific climate: Towards reducing tropical biases. J. Climate, in press. Ye Z. and W. W. Hsieh, 2006: The influence of climate regime shift on ENSO. Climate Dyn., 26, 823-833. DOI: 10.1007\/s00382-005-0105-5. Zhang D. and M.J . McPhaden, 2006: Decadal variability of the shallow Pacific meridional overturning circulation: relation to tropical sea surface temperatures in observations and climate change models. Ocean Modelling, DOI: 10.1016\/j .ocemod.2005.12.005. in press. Zhang G. J. and H. Wang, 2006: Toward mitigating the double ITCZ problem in N C A R CCSM3. Geophys. Res. Lett, 33, L06709, DOI: 10.1029\/2005GL025229. 127 Chapter 6 Summary and Conclusions 6.1 Summary and Conclusions This thesis began with a study on the influence of climate regime shift on ENSO. Eastward shifts in the location of the zonal wind stress anomalies were found in the observational data during strong El Nino episodes after the mid 1970s climate regime shift. Similar shifts appeared in the surface shortwave and latent heat fluxes in the observations. The simulated results from a modified Zebiak-Cane coupled model reproduced such shifts when a 1981-1995 climatology was used instead of a 1961-1975 climatology in the numerical model. Further sensitive experiments illustrated that the change in the atmospheric climatology, which caused the eastwards shifts in the westerly anomalies and convergence heating rate anomalies, is the major cause, with the change in the upper ocean climatology being of minor importance. The analysis of the model convergence heating rate Qc (which depends on the mean surface wind convergence and links to the latent heating rate) indeed showed that the strong shift found in Qc was the probable cause for the changes in ENSO. These results are consistent with the observed heat flux. While the change in the surface background wind resulted in the eastward shifts in the surface wind anomalies and convergence heating rate anomalies during strong El Nino, during strong La Nina the eastward shifts were not evident in the 1981-1995 regime. Thus 128 Chapter 6. Summary and Conclusions the asymmetry between the E l Nino state and the La Nina state has been enhanced in the 1981-1995 regime, which also implies enhanced nonlinearity in the ENSO system. Also consistent with observations are the enhanced predictability and lengthened ENSO period found when the model was run with the 1981-1995 climatology. A l l these agreements with observations suggest that despite the uncertainties in using climatologies calculated from relatively short records, using the 1981-1995 climatology in the coupled model did induce changes in ENSO properties consistent with the observed changes. The predictability analysis of the 2100-yr simulated SST and wind stress (WS) in the intermediate coupled model reveals that changes in ENSO properties under climate regime shifts are closely linked to changes in ENSO nonlinearity. The change in ENSO predictabil-ity occurs with the change in ENSO nonlinearity as well as with the change in ENSO period. Stronger nonlinearity is associated with longer period and better predictability. In the Lorenz 3-component chaos system, by adjusting the model parameters to change the nonlinearity of the system, a similar behavior was found. Thus in both the ENSO and Lorenz systems, stronger nonlinearity corresponds to longer period and enhanced predictability. The role of the surface zonal current was clearly distinct in the N L P C A modes associated with the two regimes: for the pre-1980s, the surface zonal current tended to be an intensifier of ENSO, while that for the post-1980s played a role in ENSO transition. An analysis of historical SST provided evidence for the E l Nino-like warming trend in the tropical Pacific during 1900-1999. This E l Niho-like warming trend was confirmed by the IPCC AR4 multi-model data by comparing the climate from year 2000 GHG(+aerosol) forcing (Commit) with the pre-industrial climate (PIcntrl). Analysis showed that both the equatorial zonal overturning circulation and the shallow meridional overturning circulation tended to slow down with increased GHG, as expected since the Walker circulation in the 129 Chapter 6. Summary and Conclusions atmosphere has been weakening since the mid-19th century due to anthropogenic forcing (Vecchi et al, 2006). It is also consistent with the conclusion that the shallow meridional overturning circulation had a slowdown trend during 1953-2000 (McPhaden and Zhang, 2002; Zhang and McPhaden, 2006) The historical SST records showed that the positive SST anomalies during E l Nino appeared further east in the 1950-1999 composite than in the 1900-1949 composite, whereas the negative SST anomalies during La Nina were shifted in the opposite direction. From the ensemble mean of the IPCC AR4 multi-model data, the positive SST anomalies during E l Nino also shifted eastward in the Commit runs relative to the PIcntrl runs, whereas the negative SST anomalies during La Nina were not shifted zonally. Hence both the model and observed SST results were consistent with increasing asymmetry between E l Nino and La Nina as GHG increased. The increase in the asymmetry is associated with an increase in the nonlinearity of ENSO (Jin et al, 2003; An and Jin, 2004), with the nonlinear dynamical heating (NDH) terms producing the asymmetry between E l Nino and La Nina. A diagnostic analysis of the ocean surface temperature tendency equation revealed the vertical NDH term to be most dominant, with the positive vertical NDH anomalies in the equatorial Pacific enhanced substantially in the Commit runs during E l Nino, caused by the enhanced vertical temperature gradient anomalies and the eastward shift of downwelling anomalies. The enhanced GHG simulations showed that the asymmetry between El Nino and La Nina was enhanced in the anomalies of the SST, the zonal wind stress, the equatorial undercurrent and the shallow meridional overturning circulation. There is much scope for further research: Consideration of double or even quadruple CO2 forcing in IPCC scenario experiments would further illustrate the potential for future climate change. The nonlinear nature of the coupled ocean-atmosphere system means that 130 Chapter 6. Summary and Conclusions changes of this magnitude on the ENSO system would be much larger than the Commit runs. We have only speculated on the mechanism responsible for the shift in the surface net shortwave radiation flux and in the latent heat flux. Lack of data is the primary reason for our poor understanding of the mechanism during decadal climate variability. Greater availability of long satellite datasets plus other measurements would allow us to study more detailed decadal variability of surface net shortwave radiation flux, longwave radiation flux, latent heat flux, and sensible heat flux over the tropical Pacific, and to understand the underlying mechanism. The main interest of this thesis focuses on the tropical Pacific; further work could extend to the Indian Ocean and the Atlantic Ocean. 131 Chapter 6. Summary and Conclusions 6.2 References A n S.-L, and F.-F. Jin, 2004: Nonlinearity and asymmetry of ENSO. J. Climate, 17, 2399-2412. J i n F.-F., S.-I. An, A. Timmermann and J. Zhao, 2003: Strong El Nino vents and nonlinear dynamical heating. Geophys. Res. Lett. 30, 1120. OL10.1029\/2002GL016356. M c P h a d e n M . J., and D. Zhang, 2002: Slowdown of the meridional verturning circulation in the upper Pacific Ocean. Nature, 415, 606-608. Vecchi G. A. , B. J. Soden, A. T. Wittenberg, I. M . Held, A. Leetmaa and M . J. Harrison, 2006: Weakening of tropical Pacific atmospheric circulation due to anthropogenic forcing. Nature, 44, 73-76. Zhang D. and M.J . McPhaden, 2006: Decadal variability of the shallow Pacific meridional overturning circulation: relation to tropical sea surface temperatures in observations and climate change models. Ocean Modelling, DOI: 10.1016\/j . ocemod. 2005.12.005. in press. 132 Appendix A The Zebiak-Cane coupled model A . l Coupled M o d e l Equations Atmosphere The dynamics follow Gill (1980), i.e., steady-state, linear shallow-water equations on an equatorial beta plane. Linear dissipation in the form of Rayleigh friction and Newtonian cooling is used. The circulation is forced by a heating anomaly that depends partly on local heating associated with SST anomalies and partly on the low-level moisture convergence (parameterized in terms of the surface wind convergence). The convergence feedback is incorporated into the model using an iterative procedure in which the heating at each iteration depends on the convergence field from the previous iteration. The feedback is nonlinear because the moisture-related heating is operative only when the total wind field is convergent, and this depends not only on the calculated convergence anomaly, but also on the specified mean convergence. On an equatorial \/3-plane the equations for the atmosphere (at iteration n) are as follows (Zebiak and Cane, 1987): euna-BoVvna = -ev2 + B0yu2 = -1 < V Po dx 1 dpn Po dy (A.l) (A.2) -\u00b1] = -Qs-Q? (A.3) 133 Appendix A. The Zebiak-Cane coupled model where Qs = (aT)exp[(T-30\u00b0C)\/16 .rC], (A.4) Qi = 8[M(c + cn)-M(c)}, (A.5) dx dy ' M(x) = 0, x < 0 x, x > 0. u a and u a are the zonal and meridional component of surface wind anomaly, respectively, p is the surface pressure anomaly, p0 is the air reference density, c a = 60 m s - 1 , and e = (2 days) - 1 is the damping coefficient (taken to be the same for momentum and thermal damping). T(x,y,t) is the prescribed monthly mean sea surface temperature (SST), and T the SST anomalies. c(x, y, t) is the prescribed monthly mean surface wind convergence, and c\" the anomalous convergence at iteration n. (x, y) is the horizontal spatial position and t is time. BQ is latitudinal derivative of the Coriolis parameter, a and 3 are both empirical coefficients, a = 0.031 m 2 s - 3 \u00b0 C - 1 , 8 = 1.6 x 104 m 2 s - 2 . Ocean The model ocean basin is rectangular and extends from 120\u00b0E to 80\u00b0W and from 29\u00b0N to 29\u00b0S. The dynamics of the model begin with the linear reduced-gravity equations. Since seasonal responses and large-scale features of the wind field are of interest, the low-frequency, long zonal-scale approximation is made: d\/dt - d\/dx = 0(e0), e0 \u00ab 1. Scaling analysis reveals that v = 0(e0), h = u = 0(u) and the governing equations for the ocean are 134 Appendix A. The Zebiak-Cane coupled model Figure A . l : Ocean model levels du ,dh . . T t - B 0 y v = - g - + \u2014 - r u , (A.6) dh B0yu = -g'\u2014 + \u2014 - r v , (A.7) dy pH dh TT,du dv, , . A \u201e . m+H[d-x + dy] = ~rh> (A-8) where u = H-1[H1u1 + H2u2], with the subscripts 1 and 2 refer to the surface layer and underlying layer, respectively. g' = gAp\/p, Ap is the density difference between the two layers, p is the mean density in the upper layer, and Ap\/p = 0.0057 (i.e., an equivalent depth He = 86 cm). Hi is the surface layer depth, H is the averaged upper layer thickness, and h is the sea level heigh (see Fig. A . l ) . u and v are the zonal and meridional components of the velocity vector u, respectively, and are the zonal and meridional components of the surface wind stress anomalies, respectively. Rayleigh friction damping with decay time r - 1 has been included. The equations governing the shear between layers 1 and 2 are 135 Appendix A. The Zebiak-Cane coupled model rsus-B0yvs = \u2014\u2014, pHi T(v) rsvs + 30yus = \u2014\u2014, pHi where the horizontal currents shear is defined by u s = (us,va) = u i - u 2 . The temperature equation for the surface layer is dT \u2014 dT T \u2014 T \u2014 = -u1-v(T + T)-u1-vT-{M{ws + ws)-M(ws)}x\u2014-M(ws + ws) e - a s T , at oz ii\\ where the entrainment velocity is calculated by ws = Hi dui dvi dx dy ui(x,y,t) and ws(x,y,t) are, respectively, the mean horizontal currents and upwelling ve-locity in the surface layer. T(x,y,t) is the prescribed mean SST, is the prescribed mean vertical temperature gradient, and cts =(125 days) - 1 is the diffusion coefficient. The entrainment temperature anomaly, T e , is defined by Te = lTsub + (1 - 7 ) T , where, in the original Zebiak and Cane (1987) model, Tsut, has the form given by (2.5). In this thesis, the relation between Tsub and h is replaced by a neural network model. First Zebiak and Cane ocean model was forced by the observed Florida State University (FSU) WS (Bourassa et al, 2001) from 1964 to 1998, and the model thermocline depth anomalies h and current anomalies were extracted. T s ub was inversely calculated from the temperature equation (2.1) using the observed SST (Smith and Reynolds, 2004) for T and 136 Appendix A. The Zebiak-Cane coupled model the model anomalous currents. P C A was separately applied to the simulated upper layer depth anomalies h(x,y,t), and to Tsub(a;,y, t) to extract the leading PCs. h(x,y,t) = y_]an(t)\/n(a;,y), n TSub(x,y,t) = ^2cn(t) en(x,y), n where n indicates the nth mode. For the thermocline depth anomaly h, the first three P C A modes accounted for 63% (32%, 22%, 9% for the three leading P C A modes, respectively) of the variance, while for the T s ub, the first three modes accounted for 64% (36%, 15%, 13% for the three leading P C A modes, respectively) of the variance. The N N used has three input neurons, namely the first three principal components (PCs) an(t) for h, and the single output neuron is one of the three leading PCs Cn(i) for T s ub, i.e., a different NN model was used to predict each predictand Cn, using three an as predictors, with no time lag between the predictors and the predictand. Data from 1964 to 1998 was used to train the N N models, The optimal N N T s ub models were obtained by randomly dividing the data into training and test data in 30 times. Then the NN models were used to replace equation (2.5) in the coupled model. Coupl ing The ocean component is forced by surface wind stress anomalies. A standard bulk formula is used to generate stress anomalies from the combination of surface wind anomalies produced by the atmosphere model and the background mean winds. The ocean dynamics time step is 10 days. The atmosphere model is steady-state and in the coupled model, the wind field is deter-mined at 10-day increments. Time dependence is allowed only in the moisture convergence component of the heating. The change in moisture convergence heating is computed at each 137 Appendix A. The Zebiak-Cane coupled model time step, and the assumed background convergence is the total convergence at the previous time step, rather than just the mean convergence (i.e., at time step m, c m = Cm-i + c m _ i in the equation (A.5)). Because of the nonlinearity of the heating parameterization, this time-marching procedure allows the development of small-scale anomalies that can persist and become unrelated to subsequent SST anomaly pattern. To eliminate them, the heat-ing is recalculated periodically using the steady model formulation, based on the current SST anomaly field. (This strategy is similar to the periodical restarts often used with the leapfrog scheme to suppress splitting of the solution.) To summarize, the calculation of the atmospheric heating is split into two parts. The portion related directly to SST gives a wind response in equilibrium with the SST field on a time scale of 10 days. The portion of the heating related to internal moisture convergence feedback operates in a time-stepping sense, and so forces a wind field adjustment on a somewhat shorter time scale. Parameter values used for the coupled simulation are as follows: r =(2.5 years) - 1, c = {g'H)1!2 = 2.9 m s - 1 , H =150 m, HY =50 m, rs =(2 days) - 1 . The standard grid for ocean dynamics is 2\u00b0 longitude by 0.5\u00b0 latitude; the standard grid for SST physics and atmosphere model is 5.625\u00b0 longitude by 2\u00b0 latitude. The intermediate coupled model domain covers the tropical Pacific from 29\u00b0S to 29\u00b0N, and from 124\u00b0E to 80\u00b0W. A . 2 C o m p u t a t i o n a l P r o c e d u r e s Atmosphere The equations (A.l) to (A.3 ) are solved by using spectral decomposition in be, and finite differencing in y. Writing 138 Appendix A. The Zebiak-Cane coupled model (Ua,Va,P,Q3,Ql) = (Ua,Va,P,Qs,Ql) exp(ifcr) and substituting into the equations (A.l) to (A.3 ) gives eV? + 0oyUl eUr:-pQyVra n n (A.10) (A.9) a Po dy Qs Q\\ (A.11) For variable \\I> (e.g., U\u00a3, V\u2122, Pn, Qs or Q\u2122), we use spatial differencing in the y direction where j is the grid index in the y direction, and Ay = 2\u00b0 latitude. The boundary conditions for the equations are that Va = 0 at the poles. In practice the solution decays rapidly outside the forcing region (which is tropical) so the solution is insensitive to the exact location of the zonal walls (There are two walls at 80\u00b0S and 80\u00b0N in the atmosphere model). The domain is from 80\u00b0S to 80\u00b0N in y direction. The spectral decomposition is done using a 64-point fast Fourier transform (FFT), and the domain is the full 360\u00b0 longitude, giving resolution (in x) of 5.625\u00b0. The numerical procedure for solving equations (A.6) to (A.8) is given in detail in Cane and Patton (1984). The solution is obtained by treating the Kelvin wave portion of the solution analytically and calculating the remainder of the solution by finite difference meth-ods on a staggered grid, taking advantage of the fact that this part of the solution has no eastward propagation. The procedure allows large time steps and is limited only by dy dy2 Ocean 139 Appendix A. The Zebiak-Cane coupled model accuracy considerations. The scheme introduces no spurious sources or sinks of mass or energy. The nondimensionalized equations can be obtained from equations (A.6) to (A.8): du dh ^ . . 'dt~yv + di = F~ru> (A-12) yu+^- = G-rv, (A.13) dy dh du dv , . m + d-x + cTy = Q~rh- (A-14) The equations have been nondimensionalized in the usual equatorial way: the time scale is (c80)~1\/2 and length scale is (c\/ \/^) 1 \/ 2 , where c \u2014 {g'H)1\/2 = 2.9 m s - 1 is a scale for wave speeds. The external forcings F, G and Q are due to the zonal wind stress, meridional wind stress and buoyancy flux, respectively. Rayleigh friction damping with decay time r _ 1 has been included; note that, with the definition (u*,v*, h*,F*,G\\ Q*) = exp(-rt)(u, v, h, F, G, Q), the equations for the starred variables are formally the same as (A. 12) to (A. 14) with r = 0. The numerical model treats friction in this way: assume r = 0 for the Kelvin wave analytical solution of the equations (A. 12) to (A. 14). The solution to the equations (A. 12) to (A. 14) may be written as {u, v, h) = ak{x, i )[*-i(y), 0, #-i(y)] + [u'(x, y, t), v'(x, y, t), h'(x, y, t)] where the prime quantities are long Rossby and anti-Kelvin modes, and the first term is the Kelvin wave, *_ 1 (y) = C 7 - 1 e x p ( - 2 \/ 7 2 ) with C2 = \/ exp(-y2)dy, fVN 'ys 140 Appendix A. The Zebiak-Cane coupled model where ys and yjv are the southern boundary and northern boundary, respectively. Then the primed part of the solution, consisting of long Rossby and anti-Kelvin modes, propagates only to the west. The equation governing ak is found by projecting the Kelvin mode on to the equations (A.12) to (A.14): dak dak inhere ( F , G , O ) ' ( * - i , 0 , * - i ) IK = (*_ 1 ,0, tf_ 1 )-(1 '_i ,0, tf_ 1 ) 1 fVN 2 rVN \/ (F + Q^-xdy. (A.15) The differential equations governing the remainder of the solution are found by sub-tracting the Kelvin wave from the equations (A. 12) to (A. 14) (with r = 0) ^ - Y V ' + ^ - = F ' = F - \/ K * - I , (A.16) dt ox dti yu' + ^ - = G, (A.17) dy dh' du' dv' . _ \/ \u00bb i o\\ -^T + TT + IT = Q=Q-fK*-i- (A.18) dt dx dy The finite-difference version is defined on a staggered grid. With n,i,j indexing time, longitude and latitude, respectively, the equations of motion take the form _ L . i V 1 +\u201e\u2122 _ \u201e\u00ab- i _ ,,\u2122-M - - i , - ( v n ~ 1 \/ 2 4-vn~1\/2 ) i^_\/tn-l , L T I tn N p\" -1\/2 ! _ \/ n-1\/2 _ n-1\/2 x_ ^n - 1 \/ 2 + ^ 2 \/ ^i+l \/2j + l\/2 Wi+l\/2,i-l\/2i - ^i+l\/2j-141 Appendix A. The Zebiak-Cane coupled model The difference scheme is second order accurate in time and space. In the absence of forcing and damping, the finite-difference equations identically conserve mass and energy. The boundary conditions for the ocean are those of no flow perpendicular to the basin boundaries. Although the model is designed to allow irregular geometry, a rectangular basin is chosen here, with boundaries at 29.75\u00b0N, 29.75\u00b0S, 124\u00b0E, 80\u00b0W. The temperature anomaly equation is evaluated on a spatially uniform grid using forward time differencing. Since the temperature does not influence the dynamical calculation, the accuracy or stability of that calculation is not affected by the temperature prediction scheme. Horizontal advection terms are evaluated using the modified upwind difference scheme. This has the form: dT, | + ui,j)(Ti,j - T i j - i ) \/ 2 A a ; , , uitj > 0 ,ui,j + ui,j+i)(Ti,j+i - Tij)\/2Ax, Uitj < 0 (A.19) 3 T N I dy (A.20) (vi+1,j + vid)(Ti+1j - Tid)\/2Ay, V i J > 0 (vi-ij + Vij)(Tij - Ti-ij)\/2Ay, vitj < 0 where i and j are the grid indices in the y and x directions, respectively. The diffusive character of this scheme can be shown with a simple temperature equation: 9T \u201e - = - u 1 . V T - a s T . The initial condition has Gaussian SST anomalies (T(x,y) = 5.0 exp(\u2014(x \u2014 XQ)2 \u2014 (y\u2014yo)2)) with Tii = 10 cm\/s in the zonal direction. Fig. A.2 shows that the diffusion due to the advection term scheme (A.19) is much weaker than that due to the diffusive term \u2014asT. The difference scheme is second-order accurate, and is transportable and conservative \u2014 it advects only in the direction of the flow, and does not produce artificial sources. These features of the scheme are important because in part of the domain (that is, away from 142 Appendix A. The Zebiak-Cane coupled model (a) I \/ I A 100 150 200 Longitude 250 300 100 150 200 Longitude 250 300 Figure A.2: The SST anomalies at the equator at different model times (0, 2, 4, 18 months) with (a) a s =0 day - 1 ; (b) as =(125 days) - 1 . the upwelling zones) there are no local advective processes, and therefore are sensitive to deficiencies in the numerical scheme. Equation (A.19) and (A.20) are numerically stable only if At ) +b(x)). (B.2) The transfer function fa is generally nonlinear (usually the hyperbolic tangent or the sig-moidal function, though the exact form is not critical), while f2 is usually taken to be the identity function. Next, a transfer function fa maps from u to the final hidden layer (the decoding layer) b>), 4 U ) = M w ( u W b<\">)fc), (B.3) (fc \u2014 1,..., m); followed by fa mapping from bJ\") to x', the output column vector of length 147 Appendix B. Nonlinear principal component analysis (NLPCA) I, with x'i = \/4((W___(U> + b[u))l). (BA) The cost function J = (||x \u2014 x' | | 2) is minimized by finding the optimal values of W^', b(*), w^, b{x\\ w W , b(\"), W(\") and b{u). The MSE between the N N output x' and the original data x is thus minimized. The N L P C A was implemented using the hyperbolic tangent function for fx and \/ 3 , and the identity function for \/ 2 and f\\, so that u = w ( x ) . h ( x ) + ^ ( R 5 ) x'i = (W____h____ + b{u))l. (B.6) Furthermore, the normalization conditions is adopted: (u) = 0 and (u2) = 1. These conditions are approximately satisfied by modifying the cost function to J = ( | | x - X ' | | 2 ) + ( U ) 2 + ( ( \u00ab 2 ) - l ) 2 . (B.7) The total number of (weight and bias) parameters used by the N L P C A is 2lm + Am + 1 + 1, though the number of effectively free parameters is two less due to the constraints on (u) and (u2). The choice of m, the number of hidden neurons in both the encoding and decoding layers, follows a general principle of parsimony. A larger m increases the nonlinear modelling capability of the network, but could also lead to overfitted solutions (i.e. wiggly solutions which fit to the noise in the data). If fi is the identity function, and m \u2014 1, then (B.6) implies that all x \u2022 are linearly related to a single hidden neuron, hence there can only be a linear relation between the x- variables. Thus, for nonlinear solutions, m > 2 is needed. It is also possible to have more than one neuron at the bottleneck layer. For instance, with two bottleneck neurons, the mode extracted will span a 2-D surface instead of a 1-D curve. 148 Appendix B. Nonlinear principal component analysis (NLPCA) The nonlinear optimization was carried out by the MATLAB function 'fminu', a quasi-Newton algorithm. Because of local minima in the cost function, there is no guarantee that the optimization algorithm reaches the global minimum. Hence a number of runs with random initial weights and bias parameters was made. Also, 20% of the data was randomly selected as validation data and withheld from the training of the NNs. Runs where the MSE was larger for the validation dataset than for the training dataset were rejected to avoid overfitted solutions. Then the run with the smallest MSE was selected as the solution. In general, the most serious problem with N L P C A is the presence of local minima in the cost function. As a result, optimizations started from different initial parameters often converge to different minima, rendering the solution unstable or nonunique. Regularization of the cost function by adding weight penalty terms is an answer. The purpose of the weight penalty terms is to limit the nonlinear power of the N L P C A , which came from the nonlinear transfer functions in the network. The transfer function tanh has the property that given x in the interval [\u2014L,L], one can find a small enough weight w, so that tanh(wx) \u00ab wx, i.e. the transfer function is almost linear. Similarly, one can choose a large enough w, so that tanh approaches a step function, thus yielding Z-shaped solutions. If the use of excessive weights is penalized, the degree of nonlinearity can be limited in the N L P C A solution. This is achieved with a modified cost function J = (||x - x ' | | 2 ) + (u)2 + ((u2) - I ) 2 + Pj2(wif)2 , (B.8) ki where P is the weight penalty parameter. A large P increases the concavity of the cost function, and forces the weights W ( x ) to be small in magnitude, thereby yielding smoother and less nonlinear solutions than when P is small or zero. Hence, increasing P also reduces the number of effectively free parameters of the model. 149 Appendix B. Nonlinear principal component analysis (NLPCA) The percentage of the variance explained by the N L P C A mode is simply 100% x 1 - ( ! ix -x ' | |2 ) ( | | X -X| |2) ) (B.9) with x being the mean of x. In effect, the linear relation u(t) \u2014 a-x(i) in P C A is now generalized toti = \/(x), where \/ can be any nonlinear continuous function representable by a feed-forward N N mapping from the input layer to the bottleneck layer; and instead of (||x(t) \u2014 au(i)|| 2), (||x \u2014 g(u)||2) is minimized. The residual, x \u2014 g(w), can be input into the same network to extract the second N L P C A mode, and so on for the higher modes. That the classical P C A is indeed a linear version of this N L P C A can be readily seen by replacing all the transfer functions with the identity function, thereby removing the nonlinear modelling capability of the N L P C A . Then the forward map to u involves only a linear combination of the original variables as in the P C A . With the N L P C A , for a given value of the N L P C u, one can map from u to the PCs. This is done by assigning the value u to the bottleneck neuron and mapping forward using the second half of the network in Fig. B . la . Each of the PCs can be multiplied by its associated P C A (spatial) eigenvector, and the three added together to yield the spatial pattern for that particular value of u. Unlike P C A which gives the same spatial anomaly pattern except for changes in the amplitude as the PC varies, the N L P C A spatial pattern generally varies continuously as the N L P C changes. Fig. B.2 is an example of N L P C A for the combined SST and WS anomalies in the tropical Pacific. The detailed procedure was given in Section 2.4.1. 150 Appendix B. Nonlinear principal component analysis (NLPCA) Figure B.2: First N L P C A mode for the combined SST and WS anomalies in the tropical Pacific (125\u00b0E-70\u00b0W, 30\u00b0S-30\u00b0N) during 1961-1995. The first N L P C A mode is indicated by the (overlapping) small circles, with the input data shown as dots. The input data were the first six PCs from combined P C A of the SST and WS anomalies (with each variable first normalized by its standard deviation). The N L P C A solution is shown projected onto (a) the PC1-PC2 plane, (b) the PC1-PC3 plane, (c) the PC1-PC4 plane, and (d) the PC1-PC5 plane. 151 Appendix B. Nonlinear principal component analysis (NLPCA) B . 2 C l o s e d c u r v e s The N L P C A is capable of finding a continuous open curve solution, but there are many geophysical phenomena involving waves or quasi-periodic fluctuations, which call for a con-tinuous closed curve solution. Kirby and Miranda (1996) introduced a N L P C A with a circular node at the network bottleneck [henceforth referred to as the NLPCA(cir)], so that the nonlinear principal component (NLPC) as represented by the circular node is an an-gular variable A, and the NLPCA(cir) is capable of approximating the data by a closed continuous curve. Fig. B . lb shows the NLPCA(cir) network, which is almost identical to the N L P C A of Fig. B . la , except at the bottleneck, where there are now two neurons p and q constrained to lie on a unit circle in the p-q plane, so there is only one free angular variable A, the N L P C . At the bottleneck in Fig. B. lb , analogous to u in (B.5), the pre-states p0 and qQ is calculated by p0 = \u2022 h ( x ) + b(x), and qQ = w ( x ) \u2022 h w + Ux), (B.10) where w^) , are weight parameter vectors, and b^ and are bias parameters. Let r = (pl + ql)ll\\ (B . l l ) then the circular node is defined with p = Po\/r, and q = q0\/r, (B.12) satisfying the unit circle equation p2 + q2 \u2014 1. Thus, even though there are two variables p and q at the bottleneck, there is only one angular degree of freedom from A (Fig. B.lb), due to the circle constraint. The mapping from the bottleneck to the output proceeds as in Appendix A, with (B.3) replaced by h[u) = tanh((w____p + w W j + b>>)fc). (B.13) 152 Appendix B. Nonlinear principal component analysis (NLPCA) When implementing NLPCA(cir), Hsieh [2001a] found that there are actually two pos-sible configurations: (i) A restricted configuration where the constraints (p) \u2014 0 = (q) are applied, and (ii) a general configuration without the constraints. With (i), the constraints can be satisfied approximately by adding the extra terms (p)2 and (q)2 to the cost function. If a closed curve solution is sought, then (i) is better than (ii) as it has effectively two fewer parameters. However, (ii), being more general than (i), can more readily model open curve solutions like a regular N L P C A . The reason is that if the input data mapped onto the p-q plane covers only a segment of the unit circle instead of the whole circle, then the inverse mapping from the p-q space to the output space will yield a solution resembling an open curve. Hence, given a dataset, (ii) may yield either a closed curve, or an open curve solution. Its generality comes with a price, namely that there may be more local minima to contend with. The number of parameters is 2lm + 6m + 1 + 2, though under configuration (i), the number of effectively free parameters is two less, due to the imposed constraints. Unlike N L P C A which reduces to P C A when only linear transfer functions are used, NLPCA(cir) does not appear to have a linear counterpart. The NLPCA(cir) analysis is shown in Fig. B.3 for the combined SST, zonal surface current and sea level height in the tropical Pacific. The detailed procedure was given in Section 4.2. B.3 References B i s h o p C. M . , 1995: Neural Networks for Pattern Recognition, Clarendon Press, 482 pp. Hsieh W. W., 2004: Nonlinear multivariate and time series analysis by neural network methods. Rev. Geophys., 42, RG1003. DOI: 10.1029\/2002RG000112. K i r b y M . J. and R. Miranda, 1996: Circular nodes in neural networks. Neural Comp., 8, 153 Appendix B. Nonlinear principal component analysis (NLPCA) ( a ) 60 i \u2022 \u2022 40 -80 1 ' ' ' \u2022 1 -150 -100 -50 0 50 100 -30 -40 ' ' ' ' ' -150 -100 -50 0 50 100 P C 1 (b) o Q. 40 30 20 0 . \u2022 9 . \u2022 \u2022\u2022 soooopo 0 ' - V - - . \u2022 \u2022 -a.-.. -10 o \u2022 0 -10 \u2022 -:>\u2022; \/\"-V-.- : -20 -30 -150 -100 -50 0 50 100 -20 -30 _ 4 0 l , , . , 1 -150 -100 -50 0 50 100 P C 1 Figure B.3: First NLPCA(cir) mode from the combined SST, zonal surface current and sea level height anomalies in the tropical Pacific (120\u00b0E-70\u00b0W, 20\u00b0S-20\u00b0N) during 1958-2001. The first NLPCA(cir) mode is indicated by the (overlapping) small circles, with the input data shown as dots. The input data were the first six PCs from a combined P C A of the SST, zonal surface current and sea level height anomalies (with each variable first normalized by its standard deviation). The NLPCA(cir) solution is a closed curve in the six-dimensional PC space. The NLPCA(cir) solution is shown projected onto (a) the PC1-PC2 plane, (b) the PC1-PC3 plane, (c) the PC1-PC4 plane, and (d) the PC1-PC5 plane. 154 Appendix B. Nonlinear principal component analysis (NLPCA) 390-402. Kramer M . A. , 1991: Nonlinear principal component analysis using autoassociative neural networks. AIChE. J., 37, 233-243. 155 ","@language":"en"}],"Genre":[{"@value":"Thesis\/Dissertation","@language":"en"}],"IsShownAt":[{"@value":"10.14288\/1.0053262","@language":"en"}],"Language":[{"@value":"eng","@language":"en"}],"Program":[{"@value":"Oceanography","@language":"en"}],"Provider":[{"@value":"Vancouver : University of British Columbia Library","@language":"en"}],"Publisher":[{"@value":"University of British Columbia","@language":"en"}],"Rights":[{"@value":"For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https:\/\/open.library.ubc.ca\/terms_of_use.","@language":"en"}],"ScholarlyLevel":[{"@value":"Graduate","@language":"en"}],"Title":[{"@value":"Changes in the El Nin\u0303o-Southern Oscillation under climate regime shift and increased greenhouse gases","@language":"en"}],"Type":[{"@value":"Text","@language":"en"}],"URI":[{"@value":"http:\/\/hdl.handle.net\/2429\/31187","@language":"en"}],"SortDate":[{"@value":"2007-12-31 AD","@language":"en"}],"@id":"doi:10.14288\/1.0053262"}__