Oscillatory dynamics do not mask linear trends in the timing of ice breakup for Northern Hemisphere lakes from 1855 to 2004

Our analyses partition the relative influence of progressive climate change and large-scale climate drivers that can be associated with the Quasi-Biennial Oscillation (QBO), El Niño Southern Oscillation (ENSO), North Atlantic Oscillation (NAO), solar sunspot cycle, and multi-decadal oscillations on lake ice breakup dates for thirteen Northern Hemisphere lakes. Oscillatory dynamics explain 26 % of the total variance in the time series compared with 15 % for linear trends, leaving 60 % unexplained and likely attributable, in part, to local weather. Significant oscillatory dynamics include frequencies in 2–3 year periods (9.4 % of the total variance), 3–6 year periods (8.2 %), 10–12 year periods (1.6 %) and various multidecadal periods (0.4–1.3 %). All 13 study lakes, although widely scattered in the Northern Hemisphere, had similar oscillatory dynamics and linear trends, emphasizing that global processes influence lake ice breakup locally. We illustrate that while quasi-periodic dynamics associated with large-scale climate drivers are important, they do not mask the clear evidence for progressive climate change.     

A brief introduction to the issue of climate and marine fisheries

Climatic variability has profound effects on the distribution, abundance and catch of oceanic fish species around the world. The major modes of this climate variability include the El Niño-Southern Oscillation (ENSO) events, the Pacific Decadal Oscillation (PDO) also referred to as the Interdecadal Pacific Oscillation (IPO), the Indian Ocean Dipole (IOD), the Southern Annular Mode (SAM) and the North Atlantic Oscillation (NAO). Other modes of climate variability include the North Pacific Gyre Oscillation (NPGO), the Atlantic Multidecadal Oscillation (AMO) and the Arctic Oscillation (AO). ENSO events are the principle source of interannual global climate variability, centred in the ocean–atmosphere circulations of the tropical Pacific Ocean and operating on seasonal to interannual time scales. ENSO and the strength of its climate teleconnections are modulated on decadal timescales by the IPO. The time scale of the IOD is seasonal to interannual. The SAM in the mid to high latitudes of the Southern Hemisphere operates in the range of 50–60 days. A prominent teleconnection pattern throughout the year in the Northern Hemisphere is the North Atlantic Oscillation (NAO) which modulates the strength of the westerlies across the North Atlantic in winter, has an impact on the catches of marine fisheries. ENSO events affect the distribution of tuna species in the equatorial Pacific, especially skipjack tuna as well as the abundance and distribution of fish along the western coasts of the Americas. The IOD modulates the distribution of tuna populations and catches in the Indian Ocean, whilst the NAO affects cod stocks heavily exploited in the Atlantic Ocean. The SAM, and its effects on sea surface temperatures influence krill biomass and fisheries catches in the Southern Ocean. The response of oceanic fish stocks to these sources of climatic variability can be used as a guide to the likely effects of climate change on these valuable resources.     

全球海洋蒸发量年代际变化归因：动力因子分析

本文基于动力调整方法,利用客观分析海气通量（OAFlux）资料研究了1958～2016年全球海洋蒸发量变化及其动力作用和辐射强迫分量的变化,发现海洋蒸发量及其动力作用分量具有一致性年代际变化特征,特别是在20世纪70年代及90年代末期存在明显的年代际转折。进一步分析发现:主要动力因子有太平洋—北美遥相关型（PNA）、北极涛动（AO）、北大西洋涛动（NAO）、厄尔尼诺—南方涛动（ENSO）和阿留申低压（AL）,并受到太平洋年代际振荡（PDO）的影响,其中,1970年代末期的转折与PNA、PDO、ENSO和AL密切相关,而1990年代末期的转折还与NAO变化有关。动力作用分量的前六个模态解释方差达到67.5%,其中,低纬北太平洋和印度洋蒸发异常主要与海表温度（SST）及其引起的环流异常有关,南太平洋、中纬北太平洋和北大西洋蒸发异常与环流异常直接相关。ENSO与PDO在全球海洋蒸发量上的影响要大于NAO。单因子相关分析发现南方涛动指数（SOI）、NAO和PDO与海洋蒸发年代际变化密切相关。总体来说,动力作用分量在海洋蒸发的年代际变化中起主导作用,其中,以ENSO、NAO和PDO的影响最大。     

The seasonal footprinting mechanism in CFSv2: simulation and impact on ENSO prediction

The seasonal footprinting mechanism (SFM) is thought to be a pre-cursor to the El Nino Southern Oscillation (ENSO). Fluctuations in the North Pacific Oscillation (NPO) impact the ocean via surface heat fluxes during winter, leaving a sea-surface temperature (SST) “footprint” in the subtropics. This footprint persists through the spring, impacting the tropical Pacific atmosphere–ocean circulation throughout the following year. The simulation of the SFM in the National Centers for Environmental Prediction (NCEP)/Climate Forecast System, version 2 (CFSv2) is likely to have an impact on operational predictions of ENSO and potentially seasonal predictions in the United States associated with ENSO teleconnection patterns. The ability of the CFSv2 to simulate the SFM and the relationship between the SFM and ENSO prediction skill in the NCEP/CFSv2 are investigated. Results indicate that the CFSv2 is able to simulate the basic characteristics of the SFM and its relationship with ENSO, including extratropical sea level pressure anomalies associated with the NPO in the winter, corresponding wind and SST anomalies that impact the tropics, and the development of ENSO-related SST anomalies the following winter. Although the model is able to predict the correct sign of ENSO associated with the SFM in a composite sense, probabilistic predictions of ENSO following a positive or negative NPO event are generally less reliable than when the NPO is not active.     

The retrospective prediction of ENSO from 1881 to 2000 by a hybrid coupled model: (II) Interdecadal and decadal variations in predictability

In this study, the retrospective predictions of ENSO (El Niño and Southern Oscillation) were performed for the period from 1881 to 2000 using a hybrid coupled model, which is an ocean general circulation model coupled to a linear statistical atmospheric model, and using a newly developed initialization scheme of SST assimilation by Ensemble Kalman Filter. With the retrospective predictions of the past 120 years, some important issues of ENSO predictability (measured by correlation and RMSE skills of NINO3 sea surface temperature anomaly index) were studied including decadal/interdecadal variations in ENSO predictability and the mechanisms responsible for these variations. Emphasis was placed on investigating the relationship between ENSO predictability and various characteristics of ENSO system such as the signal strength, the irregularity of periodicity, the noise and the nonlinearity. It is found that there are significant decadal/interdecadal variations in the prediction skills of ENSO during the past 120 years. The ENSO events were more predictable during the late nineteenth and the late twentieth centuries. The decadal/interdecadal variations of prediction skills are strongly related to the strength of sea-surface temperature anomaly (SSTA) signals, especially to the strength of SSTA signals at the frequencies of 2–4 year periods. The SSTA persistence, dominated by SSTA signals at frequencies over 4-year periods, also has a positive relationship to prediction skills. The high-frequency noise, on the other hand, has a strong inverse relationship to prediction skills, suggesting that it also probably plays an important role in ENSO predictability.     

Impacts of the central and eastern Pacific types of ENSOon sea surface temperature in the South Pacific

The present study examines the relationship between two types of El Niño–Southern Oscillation (ENSO), the central Pacific (CP) ENSO and the eastern Pacific (EP) ENSO, and the sea surface temperature (SST) variability over the South Pacific (SP) (20° S–60° S, 145° E–70° W) using NOAA OI SST for the period 1982–2006. The SP SST variability associated with the two types of ENSO varies with season. These two types of ENSO can excite different atmospheric patterns associated with the Pacific–South American mode, through which they influence the SP SST variability. Both the surface turbulent air–sea heat fluxes and the heat advection by Ekman currents (i.e., Ekman heat fluxes) have an important impact on the SST variability. An analysis of the surface mixed layer heat budget indicates that the heat fluxes (the sum of turbulent heat fluxes and Ekman heat fluxes) can effectively explain much of the SST variability related to the two types of ENSO.     

Quantifications of the Two “Flavours” of El Niño using Upper-Ocean Heat Content

The sea surface temperature (SST) or sea level pressure (SLP) has usually been used to measure the strength of El Niño–Southern Oscillation (ENSO) events. In this study, two new indices, based on the upper-ocean heat content (HC), are proposed to quantify the two “flavours” of El Niño (i.e., the Cold Tongue El Niño (CTE) and Warm Pool El Niño (WPE)). Compared with traditional SST or SLP indices, the new HC-based indices can distinguish CTE and WPE events much better and also represent the two leading modes of the interannual variability of the atmosphere–ocean coupled system in the tropical Indo-Pacific region. The two leading modes are obtained by performing multivariate Empirical Orthogonal Function analysis on two oceanic variables (SST and HC) over the tropical Pacific (30°S–30°N, 120°E–80°W) and six atmospheric variables (outgoing longwave radiation, SLP, streamfunction, and velocity potential at 850?hPa and 200?hPa) over the tropical Indo-Pacific region (30°S–30°N, 60°E–80°W) for the period 1980–2010. Because the two new HC-based indices are capable of better depicting coherent variations between the ocean and atmosphere, they can provide a supplementary tool for ENSO monitoring of and climate research into the two flavours of El Niño.     

Role of stochastic forcing in ENSO in observations and a coupled GCM

A procedure is presented to estimate the role of atmospheric stochastic forcing (SF) in El Ni?o–Southern Oscillation (ENSO) simulated by a coupled ocean–atmosphere general circulation model (CGCM), in direct comparison to observations represented by a global reanalysis product. SF is extracted from the CGCM and reanalysis as surface wind anomalies linearly independent of the sea-surface temperature anomalies. Madden–Julian Oscillation (MJO) is isolated from SF to quantify its role in ENSO. A coupled ocean–atmosphere model of intermediate complexity is forced with SF, as well as its MJO and non-MJO components, from the reanalysis and CGCM. The role of SF is estimated by comparing the original ENSO in observations and the CGCM with that reproduced by the intermediate model. ENSO statistics in both reanalysis and CGCM are better reproduced when the intermediate model is tuned to be weakly stable than unstable. The intermediate model driven by SF from the reanalysis reproduces most characteristics of observed ENSO, such as its spectrum, seasonal phase-locking, fast decorrelation of ENSO SST during boreal spring, and its lag-correlation with SF. In contrast, not all characteristics of ENSO in the CGCM are reproduced by the intermediate model when SF from the CGCM is used. The seasonal phase-locking of ENSO in the CGCM is not reproduced at all. ENSO, therefore, appears to be driven by SF to a lesser degree in the CGCM than in observations. Characteristics of observed ENSO reproduced by the intermediate model (driven by SF) can be largely attributed to the MJO; which, for instance, is responsible for the fast decorrelation of ENSO SST during boreal spring in both reanalysis and CGCM. The non-MJO component seems to be more responsible than the MJO for erroneous features of ENSO in the CGCM.     

Uncertainty of ENSO-amplitude projections in CMIP5 and CMIP6 models

There is a long-standing debate on how the El Niño/Southern Oscillation (ENSO) amplitude may change during the twenty-first century in response to global warming. Here we identify the sources of uncertainty in the ENSO amplitude projections in models participating in the Coupled Model Intercomparison Phase 5 (CMIP5) and Phase 6 (CMIP6), and quantify scenario uncertainty, model uncertainty and uncertainty due to internal variability. The model projections exhibit a large spread, ranging from increasing standard deviation of up to 0.6 °C to diminishing standard deviation of up to − 0.4 °C by the end of the twenty-first century. The ensemble-mean ENSO amplitude change is close to zero. Internal variability is the main contributor to the uncertainty during the first three decades; model uncertainty dominates thereafter, while scenario uncertainty is relatively small throughout the twenty-first century. The total uncertainty increases from CMIP5 to CMIP6: while model uncertainty is reduced, scenario uncertainty is considerably increased. The models with “realistic” ENSO dynamics have been analyzed separately and categorized into models with too small, moderate and too large ENSO amplitude in comparison to instrumental observations. The smallest uncertainties are observed in the sub-ensemble exhibiting realistic ENSO dynamics and moderate ENSO amplitude. However, the global warming signal in ENSO-amplitude change is undetectable in all sub-ensembles. The zonal wind-SST feedback is identified as an important factor determining ENSO amplitude change: global warming signal in ENSO amplitude and zonal wind-SST feedback strength are highly correlated across the CMIP5 and CMIP6 models.

     

Influences of local weather, large-scale climatic drivers, and the ca. 11 year solar cycle on lake ice breakup dates; 1905–2004

We investigate the temporal patterns in inter-annual variability in ice breakup dates for Lakes Mendota and Monona, Wisconsin, between 1905 and 2004. We analyze the contributions of long-term trends attributed to climate change, local weather, indices of sunspots, and large-scale climatic drivers, such as the North Atlantic Oscillation (NAO) and El Niňo Southern Ocean Index (ENSO) on time series of lake-ice breakup. The relative importance of the aforementioned explanatory variables was assessed using linear regression and variation partitioning models accounting for cyclic temporal dynamics as represented by Moran Eigenvector Maps (MEM). Model results explain an average of 58 % of the variation in ice breakup dates. A combination of the long-term linear trends, rain and snowfall in the month prior to breakup, air temperature in the winter prior to breakup, cyclic dynamics associated with sunspot numbers, ENSO, and for Lake Mendota, NAO, all significantly influence the timing of ice breakup. Significant cycle lengths were 3.5, 9, 11, and 50 years. Despite their proximity, Lakes Mendota and Monona exhibit differences in how and which explanatory variables were incorporated into the models. Our results indicate that lake ice dynamics are complex in both lakes and multiple interacting processes explain the residuals around the linear warming trends that characterize lake ice records.     

Irregularity and decadal variation in ENSO: a simplified model based on Principal Oscillation Patterns

A new method of estimating the decay time, mean period and forcing statistics of El Niño-Southern Oscillation (ENSO) has been found. It uses a two-dimensional stochastically forced damped linear oscillator model with the model parameters estimated from a Principal Oscillation Pattern (POP) analysis and associated observed power spectra. It makes use of extended observational time series of 150 years of sea surface temperature (SST) and sea level pressure (SLP) as well as climate model output. This approach is motivated by clear physical relationships that SST and SLP POP patterns have to the ENSO cycle, as well as to each other, indicating that they represent actual physical modes of the climate system. Moreover, the leading POP mode accounts for 20–50 % of the variance on interannual time scales. The POP real part is highly correlated with several Niño indices near zero lag while the imaginary part exhibits a 6–9 month lead time and thus is a precursor. The observed POP power spectra show markedly different behavior for the peak and precursor, the former having more power at ENSO frequencies and the latter dominating at low frequencies. The results realistically suggest a period of oscillation of 4–6 years and a decay time of 8 months, which corresponds to the practical ENSO prediction limit. A fundamental finding of this approach is that the difference between the observed peak and precursor spectra at low frequencies can be related to the forcing statistics using the simple model, as well as to the difference between patterns of decadal and interannual variability in the Pacific.     

南极海冰涛动与ENSO的关系

对近30年南极海冰密集度资料的EOF和SVD分析，发现南极地区在罗斯海外围和别林斯高晋海的海冰密集度场存在着“翘翘板”的变化特征，并与ENSO有密切联系。由此定义两个海冰关键区的差值为南极海冰涛动指数（ASOI），ASOI超前SOI和Nino3指数2个月时，其正、负相关系数达到最大，并通过α=0.001的信度检验。ASOI高、低指数阶段对应的南半球海平面气温、气压场和风场的合成分析表明，海冰关键区的异常变化可能引起温度、气压、风场的响应而影响南太平洋的洋流，进而对ENSO的发生、发展产生影响。     

Atmospheric bridge on orbital time scales

Large-scale atmospheric patterns are examined on orbital timescales using a climate model which explicitly resolves the atmosphere–ocean–sea ice dynamics. It is shown that, in contrast to boreal summer where the climate mainly follows the local radiative forcing, the boreal winter climate is strongly determined by modulation of circulation modes linked to the Arctic Oscillation/North Atlantic Oscillation (AO/NAO) and the Northern/Southern Annular Modes. We find that during a positive phase of the AO/NAO the convection in the tropical Pacific is below normal. The related atmospheric circulation provides an atmospheric bridge for the precessional forcing inducing a non-uniform temperature anomalies with large amplitudes over the continents. We argue that this is important for mechanisms responsible for multi-millennial climate variability and glacial inception.     

Large-scale variability modes of freshwater flux and precipitation over the Atlantic

 Precipitation (P) and freshwater (E-P) fluxes at the air-sea interface are investigated in the Atlantic Ocean sector using the reanalyses of the European Centre for Medium Range Weather Forecasts (ERA) and of the National Centers for Environmental Prediction (NCEP). A canonical correlation analysis method between these fields and sea level pressure (SLP) is used to identify patterns. We also test whether precipitation and freshwater fluxes can be reconstructed from SLP data. In the winter months, patterns associated with both the North Atlantic Oscillation (NAO) and the East Atlantic (EA) mode are identified. The signals are strong enough to be reconstructed from the reanalysis fields, and they correspond to a significant part of the variability. The NAO signal is more robust than the EA one. The NAO-related variability mode is also present when the monthly precipitation rate is averaged for the winter season and even for annual averages. However, in the later case, other variability of natural origin (for instance, ENSO variability) or noise from the model and assimilation system prevents the reconstruction of E-P associated with NAO from SLP variability. Difficulties are identified in the tropical Atlantic with a different behaviour of NCEP and ERA precipitation variability, especially near the Inter Tropical Convergence Zone (ITCZ). The ERA patterns suggest a NAO signature in the tropical Atlantic which has clear monthly patterns and indicates a link between the phase of NAO and changes in the position and intensity of ITCZ. However, the analysis of winter rainfall based on satellite and in situ data does not support the monthly tropical pattern of ERA precipitation although it suggests a relation between convection near 15°S and NAO during northern winter. Received: 10 February 2000 / Accepted: 7 May 2001     

A Non-Linear Dynamical–Statistical Model for Reconstruction of the Air–Sea Element Fields in the Tropical Pacific Ocean

Aiming at tackling the difficulty in establishing a sea surface temperature (SST) dynamical model, this study develops a non-linear dynamical–statistical model of SST fields and their correlative factors based on Genetic Algorithms (GA) and the dynamical system reconstruction idea, which greatly improves the El Niño–Southern Oscillation (ENSO) forecast model. Using Hadley SST data, sea surface wind (SSW) and sea level pressure (SLP) data from the National Centers for Environmental Prediction-National Center for Environmental Research (NCEP-NCAR), with empirical orthogonal function (EOF) time-space for reconstruction, we carry out numerical integral forecasting experiments for SST, SSW, and SLP fields. By statistical analysis of the forecasting experiments, we find that forecasts for less than 25 months perform better than longer term forecasts. Based on the model, we forecast SST, SSW, and SLP fields in September, October, and November 2014 and predict a weak La Niña event. This study explores a novel method for the complex atmosphere–ocean system.     

El Niño teleconnections in CMIP5 models

Teleconnections associated with warm El Niño/southern oscillation (ENSO) events in 20 climate model intercomparison project 5 (CMIP5) models have been compared with reanalysis observations. Focus has been placed on compact time and space indices, which can be assigned a specific statistical confidence. Nearly all of the models have surface temperature, precipitation and 250 hPa geopotential height departures in the Tropics that are in good agreement with the observations. Most of the models also have realistic anomalies of Northern Hemisphere near-surface temperature, precipitation and 500 hPa geopotential height. Model skill for these variables is significantly related to the ability of a model to accurately simulate Tropical 250 hPa height departures. Additionally, most models have realistic temperature and precipitation anomalies over North America, which are linked to a model’s ability to simulate Tropical 250 hPa and Northern Hemisphere 500 hPa height departures. The skills of temperature and precipitation departures over the Northern Hemisphere and North America are associated with the ability to realistically simulate realistic ENSO frequency and length. Neither horizontal nor vertical resolution differences for either the model atmosphere or ocean are significantly related at the 95 % level to variations in El Niño simulation quality. Overall, recent versions of earlier models have improved in their ability to simulate El Niño teleconnections. For instance, the average model skills of temperature and precipitation for the Tropics, Northern Hemisphere and North America for 11 CMIP5 models are all larger than those for prior versions.     

Interannual variability of winter precipitation over northwest India and adjoining region: impact of global forcings

The role of El Niño/Southern Oscillation (ENSO) and the mechanism through which ENSO influences the precipitation variability over northwest India and the adjoining (NWIA) region is well documented. In this study, the relative role of North Atlantic Oscillation (NAO)/Arctic Oscillation (AO) and ENSO in modulating the Asian jet stream in the Northern Hemisphere winter and their relative impact on the precipitation variability over the region have been estimated through analysis of observed data. It is seen that interannual variations of NWIA precipitation are largely influenced by ENSO. An empirical orthogonal function (EOF) analysis has been carried out to understand dominant modes of interannual variability of zonal wind at 200 hPa of the Northern Hemisphere. The EOF-1 pattern in the tropical region is similar to that of an ENSO pattern, and the principal component (PC) time series corresponds to the ENSO time series. The EOF-2 spatial pattern resembles that of NAO/AO with correlation of PC time series with AO and NAO being 0.74 and 0.62, respectively. The precipitation anomaly time series over the region of interest has marginally higher correlation with the PC-2 time series as compared to that of PC-1. Regression analysis of precipitation and circulation parameters indicates a larger contribution of the second mode to variability of winds and precipitation over the NWIA. Moisture transport from the Arabian Sea during the active phase of NAO/AO and the presence of a cyclonic anomaly lead to higher precipitation over the NWIA region.     

The nonlinear association between ENSO and the Euro-Atlantic winter sea level pressure

A nonlinear projection of the tropical Pacific sea surface temperature anomalies (SSTA) onto the Northern Hemisphere winter sea level pressure (SLP) anomalies by neural networks (NN) was performed to investigate the nonlinear association between El Niño-Southern Oscillation (ENSO) and the Euro-Atlantic winter climate. While the linear impact of ENSO on the Euro-Atlantic winter SLP is weak, the NN projection reveals statistically significant SLP anomalies over the Euro-Atlantic sector during both extreme cold and warm ENSO episodes, suggesting that the Euro-Atlantic climate mainly responds to ENSO nonlinearly. The nonlinear response, mainly a quadratic response to the SSTA, reveals that regardless of the sign of the SSTA, positive SLP anomalies are found over the North Atlantic, stretching from eastern Canada to Europe (with anomaly center located just northwestward of Portugal), and negative anomalies centered over Scandinavia and Norwegian Sea, consistent with the excitation of the positive North Atlantic Oscillation pattern.     

Understanding the SAM influence on the South Pacific ENSO teleconnection

The relationship between the El Niño Southern Oscillation (ENSO) and the Southern Hemisphere Annular Mode (SAM) is examined, with the goal of understanding how various strong SAM events modulate the ENSO teleconnection to the South Pacific (45°–70°S, 150°–70°W). The focus is on multi-month, multi-event variations during the last 50 years. A significant (p < 0.10) relationship is observed, most marked during the austral summer and in the 1970s and 1990s. In most cases, the significant relationship is brought about by La Niña (El Niño) events occurring with positive (negative) phases of the SAM more often than expected by chance. The South Pacific teleconnection magnitude is found to be strongly dependent on the SAM phase. Only when ENSO events occur with a weak SAM or when a La Niña (El Niño) occurs with a positive (negative) SAM phase are significant South Pacific teleconnections found. This modulation in the South Pacific ENSO teleconnection is directly tied to the interaction of the anomalous ENSO and SAM transient eddy momentum fluxes. During La Niña/SAM+ and El Niño/SAM? combinations, the anomalous transient momentum fluxes in the Pacific act to reinforce the circulation anomalies in the midlatitudes, altering the circulation in such a way to maintain the ENSO teleconnections. In La Niña/SAM? and El Niño/SAM+ cases, the anomalous transient eddies oppose each other in the midlatitudes, overall acting to reduce the magnitude of the high latitude ENSO teleconnection.