Search

Given the high population density and serious air pollution problems, understanding and predicting air pollution in East Asia are of great importance. Here, we show that the day-to-day variability of East Asia air pollution in winter is remotely controlled by the convection over the tropical Indian Ocean and western Pacific, the so-called Madden–Julian oscillation (MJO), through its extratropical teleconnections. In particular, the concentration of particulate matter with aerodynamic diameter less than 10 micron (PM10) becomes significantly high when the tropical convection is suppressed over the Indian Ocean (MJO phases 5–6). In contrast, PM10 concentration becomes significantly low when the convection is enhanced there (MJO phase 1–2). The station-averaged PM10 difference between the two MJO phases reaches up to 47% of the daily PM10 variability, indicating that the MJO is a primary source of wintertime subseasonal variability of East Asia PM10 concentration. We also show that PM10 anomaly typically lags the tropical convection by one to two weeks. This opens a new window of opportunity for subseasonal PM10 prediction in East Asia.

The stratospheric quasi-biennial oscillation (QBO) and the tropospheric Madden–Julian oscillation (MJO) are strongly linked in boreal winter. In this Review, we synthesize observational and modelling evidence for this QBO–MJO connection and discuss its effects on MJO teleconnections and subseasonal-to-seasonal predictions. After 1980, observations indicate that, during winters when lower-stratospheric QBO winds are easterly, the MJO is ~40% stronger and persists roughly 10 days longer compared with when QBO winds are westerly. Global subseasonal forecast models, in turn, show a 1-week improvement (or 25% enhancement) in MJO prediction skill in QBO easterly versus QBO westerly phases. Despite the robustness of the observed QBO–MJO link and its global impacts via atmospheric teleconnections, the mechanisms that drive the connection are uncertain. Theories largely centre on QBO-related temperature stratification effects and subsequent impacts on deep convection, although other hypotheses propose that cloud radiative effects or QBO impacts on wave propagation might be important. Most numerical models, however, are unable to reproduce the observed QBO–MJO relationship, suggesting biases, deficiencies or omission of key physical processes in the models. While future work must strive to better understand all aspects of the QBO–MJO link, focus is needed on establishing a working mechanism and capturing the connection in models. The quasi-biennial oscillation and Madden–Julian oscillation are strongly connected during boreal winter, typified by a more active and persistent Madden–Julian oscillation during easterly quasi-biennial oscillation phases. This Review outlines the characteristics and potential mechanisms of this coupling, as well as the implications for seasonal-to-seasonal prediction.

Producing accurate weather prediction beyond two weeks is an urgent challenge due to its ever-increasing socioeconomic value. The Madden-Julian Oscillation (MJO), a planetary-scale tropical convective system, serves as a primary source of global subseasonal (i.e., targeting three to four weeks) predictability. During the past decades, operational forecasting systems have improved substantially, while the MJO prediction skill has not yet reached its potential predictability, partly due to the systematic errors caused by imperfect numerical models. Here, to improve the MJO prediction skill, we blend the state-of-the-art dynamical forecasts and observations with a Deep Learning bias correction method. With Deep Learning bias correction, multi-model forecast errors in MJO amplitude and phase averaged over four weeks are significantly reduced by about 90% and 77%, respectively. Most models show the greatest improvement for MJO events starting from the Indian Ocean and crossing the Maritime Continent., The Madden-Julian Oscillation (MJO) is a crucial component of the tropical weather system, but forecasting it has been challenging. Here, the authors present a deep learning bias correction method that significantly improves multi-model forecasts of the MJO amplitude and phase for up to four weeks.

The prediction of wintertime extratropical cyclone activity (ECA) on subseasonal time scales by models participating in the Subseasonal Experiment (SubX) and the Seasonal to Subseasonal Prediction (S2S) is assessed. Consistent with a previous study that investigated the S2S models, the SubX models have skillful predictions of ECA over regions from central North Pacific across North America to western North Atlantic, as well as East Asia and northern and southern part of eastern North Atlantic at 3–4 weeks lead time. SubX provides daily mean data, while S2S provides instantaneous data at 0000 UTC each day. This leads to different variance of ECA. Different S2S and SubX models have different reforecast initialization times and reforecast time periods. These factors can all lead to differences in prediction skill. To fairly compare the prediction skill between different models, we develop a novel way to evaluate the prediction of individual model across the two ensembles by comparing every model to the Climate Forecast System, version 2 (CFSv2), as CFSv2 has 6-hourly output and forecasts initialized every day. Among the S2S and SubX models, the European Centre for Medium-Range Weather Forecasts model exhibits the best prediction skill, followed by CFSv2. Our results also suggest that while the prediction skill is sensitive to forecast lead time, including forecasts up to 4 days old into the ensemble may still be useful for weeks 3–4 predictions of ECA.

There is a growing demand for understanding sources of predictability on subseasonal to seasonal (S2S) time scales. Predictability at subseasonal time scales is believed to come from processes varying slower than the atmosphere such as soil moisture, snowpack, sea ice, and ocean heat content. The stratosphere as well as tropospheric modes of variability can also provide predictability at subseasonal time scales. However, the contributions of the above sources to S2S predictability are not well quantified. Here we evaluate the subseasonal prediction skill of the Community Earth System Model, version 1 (CESM1), in the default version of the model as well as a version with the improved representation of stratospheric variability to assess the role of an improved stratosphere on prediction skill. We demonstrate that the subseasonal skill of CESM1 for surface temperature and precipitation is comparable to that of operational models. We find that a better-resolved stratosphere improves stratospheric but not surface prediction skill for weeks 3–4.

In an assessment of 29 global climate models (GCMs), Part I of this study identified biases in boreal winter MJO teleconnections in anomalous 500-hPa geopotential height over the Pacific–North America (PNA) region that are common to many models: an eastward shift, a longer persistence, and a larger amplitude. In Part II, we explore the relationships of the teleconnection metrics developed in Part I with several existing and newly developed MJO and basic state (the mean subtropical westerly jet) metrics. The MJO and basic state diagnostics indicate that the MJO is generally weaker and less coherent and propagates faster in models compared to observations. The mean subtropical jet also exhibits notable biases such as too strong amplitude, excessive eastward extension, or southward shift. The following relationships are found to be robust among the models: 1) models with a faster MJO propagation tend to produce weaker teleconnections; 2) models with a less coherent eastward MJO propagation tend to simulate more persistent MJO teleconnections; 3) models with a stronger westerly jet produce stronger and eastward shifted MJO teleconnections; 4) models with an eastward extended jet produce an eastward shift in MJO teleconnections; and 5) models with a southward shifted jet produce stronger MJO teleconnections. The results are supported by linear baroclinic model experiments. Our results suggest that the larger amplitude and eastward shift biases in GCM MJO teleconnections can be attributed to the biases in the westerly jet, and that the longer persistence bias is likely due to the lack of coherent eastward MJO propagation.

We propose a set of MJO teleconnection diagnostics that enables an objective evaluation of model simulations, a fair model-to-model comparison, and a consistent tracking of model improvement. Various skill metrics are derived from teleconnection diagnostics including five performance-based metrics that characterize the pattern, amplitude, east–west position, persistence, and consistency of MJO teleconnections and additional two process-oriented metrics that are designed to characterize the location and intensity of the anomalous Rossby wave source (RWS). The proposed teleconnection skill metrics are used to compare the characteristics of boreal winter MJO teleconnections (500-hPa geopotential height anomaly) over the Pacific–North America (PNA) region in 29 global climate models (GCMs). The results show that current GCMs generally produce MJO teleconnections that are stronger, more persistent, and extend too far to the east when compared to those observed in reanalysis. In general, models simulate more realistic teleconnection patterns when the MJO is in phases 2–3 or phases 7–8, which are characterized by a dipole convection pattern over the Indian Ocean and western to central Pacific. The higher model skill for phases 2, 7, and 8 may be due to these phases producing more consistent teleconnection patterns between individual MJO events than other phases, although the consistency is lower in most models than observed. Models that simulate realistic RWS patterns better reproduce MJO teleconnection patterns.

The impact of the stratospheric quasi‐biennial oscillation (QBO) on the prediction of the tropospheric Madden‐Julian oscillation (MJO) is evaluated in reforecasts from nine models participating in subseasonal prediction projects, including the Subseasonal Experiment (SubX) and Subseasonal to Seasonal (S2S) projects. When MJO prediction skill is analyzed for December to February, MJO prediction skill is higher in the easterly phase of the QBO than the westerly phase, consistent with previous studies. However, the relationship between QBO phase and MJO prediction skill is not statistically significant for most models. This insignificant QBO‐MJO skill relationship is further confirmed by comparing two subseasonal reforecast experiments with the Community Earth System Model v1 using both a high‐top (46‐level) and low‐top (30‐level) version of the Community Atmosphere Model v5. While there are clear differences in the forecasted QBO between the two model top configurations, a negligible change is shown in the MJO prediction, indicating that the QBO in this model may not directly control the MJO prediction and supporting the insignificant QBO‐MJO skill relationship found in SubX and S2S models., Key Points The QBO‐MJO prediction skill relationship is mostly insignificant in nine subseasonal reforecastsThe QBO‐MJO prediction skill relationship is sensitive to the choice of QBO eventsThe insignificant QBO‐MJO skill relationship is supported by comparing reforecasts by a high‐top and low‐top version of CESM1

There has been an accelerating interest in forecasting the weather and climate within the subseasonal time range. The Madden–Julian oscillation (MJO), an organized envelope of tropical convection, is recognized as one of the leading sources of subseasonal predictability. This review synthesizes the latest progress regarding the MJO predictability and prediction. During the past decade, the MJO prediction skill in dynamical prediction systems has exceeded the skill of empirical predictions. Such improvement has been mainly attributed to more observations and computer resources, advances in theoretical understanding, and improved numerical models aided in part by multinational efforts through field campaigns and multimodel experiments. The state-of-the-art dynamical forecasts have shown MJO prediction skill up to 5 weeks. Prediction skill can be extended by improving the ensemble generation approach tailored for MJO prediction and by averaging multiensembles or multimodels. MJO prediction skill can be influenced by the tropical mean state and low-frequency climate mode variations, as well as by the extratropical circulation. MJO prediction skill is proven to be sensitive to model physics, ocean–atmosphere coupling, and quality of initial conditions, while the impact of the model resolution seems to be marginal. Remaining challenges and recommendations on new research avenues to fully realize the predictability of the MJO are discussed.

We investigate how the Madden‐Julian Oscillation (MJO), the dominant mode of tropical subseasonal variability, modulates the lifecycle of cool‐season North Pacific atmospheric rivers (ARs). When the enhanced (suppressed) convection center is located over the Indian Ocean (western Pacific), more AR events originate over eastern Asia and with fewer over the subtropical northern Pacific. When the enhanced (suppressed) convection is over the western Pacific (Indian Ocean), the opposite changes occur, with more AR events originate over the subtropical northern Pacific and fewer over eastern Asia. Dynamical processes involving anomalous MJO wind and seasonal mean moisture are found to be the dominant factors impacting these variations in AR origins. The MJO‐related anomalous geopotential height patterns are also shown to modulate the propagation of the AR events. These MJO–AR lifecycle relationships are further supported by model simulations.

Initialized Earth System predictions are made by starting a numerical prediction model in a state as consistent as possible to observations and running it forward in time for up to 10 years. Skilful predictions at time slices from subseasonal to seasonal (S2S), seasonal to interannual (S2I) and seasonal to decadal (S2D) offer information useful for various stakeholders, ranging from agriculture to water resource management to human and infrastructure safety. In this Review, we examine the processes influencing predictability, and discuss estimates of skill across S2S, S2I and S2D timescales. There are encouraging signs that skilful predictions can be made: on S2S timescales, there has been some skill in predicting the Madden–Julian Oscillation and North Atlantic Oscillation; on S2I, in predicting the El Niño–Southern Oscillation; and on S2D, in predicting ocean and atmosphere variability in the North Atlantic region. However, challenges remain, and future work must prioritize reducing model error, more effectively communicating forecasts to users, and increasing process and mechanistic understanding that could enhance predictive skill and, in turn, confidence. As numerical models progress towards Earth System models, initialized predictions are expanding to include prediction of sea ice, air pollution, and terrestrial and ocean biochemistry that can bring clear benefit to society and various stakeholders. The foundations of this Review emerged from a workshop held by the National Academies of Sciences, Engineering and Medicine in 2015 at Woods Hole, MA, USA, and the authors gratefully acknowledge support from A. Purcell and N. Huddleston. G.A.M., J.H.R. and N.R. were supported in part by the Regional and Global Model Analysis (RGMA) component of the Earth and Environmental System Modeling Program of the US Department of Energy’s Office of Biological & Environmental Research (BER) via National Science Foundation IA 1844590, and by the National Center for Atmospheric Research (NCAR), which is a major facility sponsored by the National Science Foundation (NSF) under Cooperative Agreement No. 1852977. M.H.E. acknowledges support from the Australian Research Council (Grant CE170100023). H.T. was partly supported by DOE/BER RGMA HiLAT-RASM. M.E.M. was supported by a grant from the NSF Paleoclimate Program #1748097. F.D.-R. and M.G.D. were supported by the H2020 EUCP project under Grant agreement no. 776613, and M.G.D also by the Ramón y Cajal 2017 grant reference RYC-2017-22964. A.C. acknowledges support from the National Oceanic and Atmospheric Administration (NOAA) Climate Program Office’s Modeling Analysis, Prediction and Projections (MAPP) Program and from the NOAA Climate Program Office’s Climate Variability and Predictability (CVP) Program. A.C.S. acknowledges support from the NOAA Climate Variability and Predictability Program (Award NA18OAR4310405) and NOAA-MAPP (NA17OAR4310106). N.S.L. is grateful for support from the NSF (OCE-1752724). D.T. acknowledges support from the NCAR Advanced Study Program and NSF (OCE-1931242). S.C.S was supported by the Joint Institute for the Study of the Atmosphere and Ocean (JISAO) Postdoctoral Fellowship. A.A.S. and D.S were supported by the Met Office Hadley Centre Climate Programme funded by the Department for Business, Energy & Industrial Strategy (BEIS) and Department for Environment, Food and Rural Affairs (Defra) and by the European Commission Horizon 2020 EUCP project (GA 776613). Peer Reviewed "Article signat per 26 autors/es: Gerald A. Meehl, Jadwiga H. Richter, Haiyan Teng, Antonietta Capotondi, Kim Cobb, Francisco Doblas-Reyes, Markus G. Donat, Matthew H. England, John C. Fyfe, Weiqing Han, Hyemi Kim, Ben P. Kirtman, Yochanan Kushnir, Nicole S. Lovenduski, Michael E. Mann, William J. Merryfield, Veronica Nieves, Kathy Pegion, Nan Rosenbloom, Sara C. Sanchez, Adam A. Scaife, Doug Smith, Aneesh C. Subramanian, Lantao Sun, Diane Thompson, Caroline C. Ummenhofer & Shang-Ping Xie"

Since its discovery in the early 1970s, the crucial role of the Madden-Julian Oscillation (MJO) in the global hydrological cycle and its tremendous influence of high-impact climate and weather extremes have been well recognised. The MJO also serves as a primary source of predictability for global Earth system variability on subseasonal time scales. The MJO remains poorly represeted in our state-of-the-art climate and weather forecasting models, however. Moreover, despite the advances made in recent decades, theories for the MJO still disagree at a fundamental level. The problems of understanding and modeling the MJO have attracted significant interest from the research community. As part of the AGU's Centennial collection, this article provides a review of recent progress, particularly over the last decade, in observational, modeling and theoretical study of the MJO. A brief outlook for near-future MJO research directions is also provided.

In this study, the intraseasonal variations in storm-track activity, surface air temperature, and precipitation over North America associated with the Madden–Julian oscillation (MJO) in boreal winter (November–April) are investigated. A lag composite strategy that considers different MJO phases and different lag days is developed. The results highlight regions over which the MJO has significant impacts on surface weather on intraseasonal time scales. A north–south shift of storm-track activity associated with the MJO is found over North America. The shift is consistent with the MJO-related surface air temperature anomaly over the eastern United States. In many regions over the western, central, and southeastern United States, the MJO-related precipitation signal is also consistent with nearby storm-track activity. An MJO-related north–south shift of precipitation is also found near the west coast of North America, with the precipitation over California being consistent with the MJO-related storm-track activity over the eastern Pacific. MJO-related temperature and storm-track anomalies are also found near Alaska. Further analyses of streamfunction anomalies and wave activity flux show clear signatures of Rossby wave trains excited by convection anomalies related to MJO phases 3 and 8. These wave trains propagate across the Pacific and North America, bringing an anticyclonic (cyclonic) anomaly to the eastern part of North America, shifting the westerly jet to the north (south), thereby modulating the surface air temperature and storm-track activity over the continent. Rossby waves associated with phases 2 and 6 are also found to impact the U.S. West Coast.

This study uses a global land–atmosphere coupled model, the land–atmosphere component of the Global Seasonal Forecast System version 5, to quantify the degree to which soil moisture initialization could potentially enhance boreal summer surface air temperature forecast skill. Two sets of hindcast experiments are performed by prescribing the observed sea surface temperature as the boundary condition for a 15-year period (1996–2010). In one set of the hindcast experiments (noINIT), the initial soil moisture conditions are randomly taken from a long-term simulation. In the other set (INIT), the initial soil moisture conditions are taken from an observation-driven offline Land Surface Model (LSM) simulation. The soil moisture conditions from the offline LSM simulation are calibrated using the forecast model statistics to minimize the inconsistency between the LSM and the land–atmosphere coupled model in their mean and variability. Results show a higher boreal summer surface air temperature prediction skill in INIT than in noINIT, demonstrating the potential benefit from an accurate soil moisture initialization. The forecast skill enhancement appears especially in the areas in which the evaporative fraction—the ratio of surface latent heat flux to net surface incoming radiation—is sensitive to soil moisture amount. These areas lie in the transitional regime between humid and arid climates. Examination of the extreme 2003 European and 2010 Russian heat wave events reveal that the regionally anomalous soil moisture conditions during the events played an important role in maintaining the stationary circulation anomalies, especially those near the surface.

Predictions of changes of the land monsoon rainfall (LMR) in the coming decades are of vital importance for successful sustainable economic development. Current dynamic models, though, have shown little skill in the decadal prediction of the Northern Hemisphere (NH) LMR (NHLMR). The physical basis and predictability for such predictions remain largely unexplored. Decadal change of the NHLMR reflects changes in the total NH continental precipitation, tropical general circulation, and regional land monsoon rainfall over northern Africa, India, East Asia, and North America. Using observations from 1901 to 2014 and numerical experiments, it is shown that the decadal variability of the NHLMR is rooted primarily in (i) the north–south hemispheric thermal contrast in the Atlantic–Indian Ocean sector measured by the North Atlantic–south Indian Ocean dipole (NAID) sea surface temperature (SST) index and (ii) an east–west thermal contrast in the Pacific measured by an extended El Niño–Southern Oscillation (XEN) index. Results from a 500-yr preindustrial control experiment demonstrate that the leading mode of decadal NHLMR and the associated NAID and XEN SST anomalies may be largely an internal mode of Earth’s climate system, although possibly modified by natural and anthropogenic external forcing. A 51-yr, independent forward-rolling decadal hindcast was made with a hybrid dynamic conceptual model and using the NAID index predicted by a multiclimate model ensemble. The results demonstrate that the decadal changes in the NHLMR can be predicted approximately a decade in advance with significant skills, opening a promising way forward for decadal predictions of regional land monsoon rainfall worldwide.

Prediction skills of the wintertime atmospheric rivers (ARs) and moisture flux over the Northeast Pacific in response to El Nino Southern Oscillation (ENSO) forcing are evaluated from the North American Multi-Model Ensemble (NMME) hindcasts (CFSv2, CCSM4, CanCM3, CanCM4, GEOS5, and GFDL CM2.1). The skill is estimated for the active AR season, December–February (DJF) with initial conditions around early November. Models underestimate the climatological moisture flux to different extents corresponding with various climatological biases in predictions of sea surface temperature (SST) and large-scale atmospheric circulation fields. The anomalous moisture flux and AR frequency over the Northeast Pacific are predicted in the models but in weaker amplitude than the reanalysis. Significant regional biases are shown in the anomalous landfalling AR frequency corresponding with ENSO, underlining the challenge in regional precipitation forecasts.

The relative contributions of atmospheric energy transport (via heat and moisture advection) and sea ice decline to recent Arctic warming were investigated using high-resolution reanalysis data up to 2017. During the Arctic winter, a variation of downward longwave radiation (DLR) is fundamental in modulating Arctic surface temperature. In the warm Arctic winter, DLR and precipitable water (PW) are increasing over the entire Arctic; however, the major drivers for such increases differ regionally. In areas such as the northern Greenland Sea, increasing DLR and PW are caused mainly by convergence of atmospheric energy transport from lower latitudes. In regions of maximum sea ice retreat (e.g., northern Barents–Kara Seas), continued sea ice melting from previous seasons drive the DLR and PW increases, consistent with the positive ice–insulation feedback. Distinct local feedbacks between open water and ice-retreat regions were further compared. In open water regions, a reduced ocean–atmosphere temperature gradient caused by atmospheric warming suppresses surface turbulent heat flux (THF) release from the ocean to the atmosphere; thus, surface warming cannot accelerate. Conversely, in ice-retreat regions, sea ice reduction allows the relatively warm ocean to interact with the colder atmosphere via surface THF release. This increases temperature and humidity in the lower troposphere consistent with the positive ice–insulation feedback. The implication of this study is that Arctic warming will slow as the open water fraction increases. Therefore, given sustained greenhouse warming, the roles of atmospheric heat and moisture transport from lower latitudes are likely to become increasingly critical in the future Arctic climate.

The impact of the seasonal mean background state on the physical processes responsible for better Madden-Julian Oscillation (MJO) propagation and prediction was investigated using 20-year ECMWF reforecast datasets. While the reforecast captures the MJO propagating signal in the circulation field, the OLR and column-integrated moist static energy (MSE) anomalies show a weaker signal than the reanalysis when the MJO detours and propagates in the southern Maritime Continent (MC). By comparing individual MSE budget terms, it is found that the predicted positive MSE tendency to the east of the MJO convection center is weaker than the reanalysis due to the weaker horizontal MSE advection. The dry bias in the seasonal mean lower tropospheric moisture is a key factor that weakens the horizontal MSE advection, and thus deteriorates the MJO propagation in the reforecasts. Therefore, improvements in the mean state could further extend the MJO prediction and associated subseasonal weather phenomena over the globe.

An interdecadal weakening in the North Atlantic storm track (NAST) and a poleward shift of the North Pacific storm track (NPST) are found during October–March for the period 1979–2015. A significant warming of surface air temperature (Ts) over northeastern North America and a La Niña–like change in the North Pacific under the background of Arctic amplification are found to be the contributors to the observed changes in the NAST and the NPST, respectively, via modulation of local baroclinicity. The interdecadal change in baroclinic energy conversion is consistent with changes in storm tracks with an energy loss from eddies to mean flow over the North Atlantic and an energy gain over the North Pacific. The analysis of simulations from the Community Earth System Model Large Ensemble project, although with some biases in storm-track and Ts simulations, supports the observed relationship between the NAST and Ts over northeastern North America, as well as the link between the NPST and El Niño–Southern Oscillation. The near-future projections of Ts and storm tracks are characterized by a warmer planet under the influence of increasing greenhouse gases and a significant weakening of both the NAST and the NPST. The potential role of the NAST in redistributing changes in Ts over the surrounding regions is also examined. The anomalous equatorward moisture flux associated with the weakening trend of the NAST would enhance the warming over its upstream region and hinder the warming over its downstream region via modulation of the downward infrared radiation.

The year-to-year changes in atmospheric rivers (ARs) and moisture transport over the northeast Pacific and western North America are investigated during December to February (DJF) from 1979/80 to 2015/16. Changes in AR frequency, intensity, and landfall characteristics are compared between three ENSO phases: central Pacific El Nino (CPEN), eastern Pacific El Nino (EPEN), and La Nina (NINA). During EPEN events, the subtropical jet extends to the south and east with an anomalous cyclonic flow around a deeper Aleutian Low. More moisture is transported towards North America and AR frequency is increased over western North America. In CPEN events, the Aleutian low shifts further southward relative to its position in EPEN, resulting in an increase in the frequency and intensity of landfalling ARs over the southwestern US. In NINA events, the landfalling AR frequency is reduced associated with anomalous anticyclonic circulation over the eastern North Pacific. We diagnose the contribution of multiple factors to the seasonal mean moisture transport using moisture budgets. During the three ENSO phases, the change in low-frequency circulation (dynamical process) is the leading contributor to the seasonal mean moisture flux divergence, while the contributions of the synoptic anomalies and the change in moisture anomaly (thermodynamic process) are not significant along the west coast of North America.

Author(s): Zhou, Y; Kim, H | Abstract: An atmospheric river (AR) event represents strong poleward moisture transport and is defined as a series of spatiotemporally connected instantaneous AR objects. Utilizing an AR tracking algorithm with a depth-first search (a widely used algorithm in computer science), we examine the life cycle characteristics of AR events that make landfall over the U.S. West Coast by their distinct origin locations. Landfalling AR events from the Northwest Pacific (120°E–170°W, WLAR events) temporally last longer (5.3 versus 3.6 days on average) and have stronger intensity of integrated vapor transport (508 versus 388 kg m−1 s−1 on average) than those originating from the Northeast Pacific (125°W–170°W, ELAR events). A persistent tripole geopotential height anomaly pattern over the North Pacific modulates the origin locations and propagation of landfalling AR events. WLAR events are associated with anomalous highs over northeastern Asia and the Northeast Pacific and an anomalous low over the central North Pacific. This pattern provides favorable conditions for WLAR events to start, propagate northeastward, and make landfall in the northwestern West Coast. WLAR events contribute approximately 25% of the total winter precipitation over Washington and British Columbia. ELAR events are associated with a similar tripole pattern to the WLAR eventsnwith an eastward shift. The anomalous low over the Northeast Pacific helps ELAR events to start, propagate northeastward, and make landfall in the southwestern West Coast. Precipitation induced by ELAR events contributes up to 30% of total winter precipitation over California.

The characteristics of the MJO propagation across the Maritime Continent are investigated using a 20-yr reforecast dataset from the ECMWF ensemble prediction system. Analysis of the MJO events initialized over the Indian Ocean (phase 2) shows that the initial MJO amplitude and prediction skill relationship is not linear, particularly when the predictions start in moderate (between strong and weak) MJO amplitude category. To examine the key factors that determine the prediction skill, reforecasts in the moderate category are grouped into high- and low-skill events, and the differences in their ocean–atmospheric conditions as well as the physical processes during reforecast period are examined. The initial distribution of OLR anomalies in high-skill events shows a clear dipole pattern of convection with an enhanced convective anomalies over the Indian Ocean and strongly suppressed convective anomalies in the western Pacific Ocean. This dipole mode may support the MJO propagation across the Maritime Continent via the Rossby wave response and associated meridional moisture advection. Prominent ocean–atmosphere coupled processes are also simulated during the propagation of high-skill events. However, in low-skill events, the convective signal over the western Pacific is almost absent and less organized, and the ocean–atmosphere coupled processes are not simulated correctly. It is found that in both high- and low-skill events, the amplitude of the convective anomaly decreases significantly after about day 15, possibly due to the systematic mean model bias. A strong wet bias in the vicinity of the Maritime Continent, a cold SST bias in the equatorial Pacific, and associated circulation biases make the west Pacific area unfavorable for MJO propagation, thus limiting its prediction skill.

This study explores the seasonal-to-interannual near-surface air temperature (TAS) prediction skills of state-of-the-art climate models that were involved in phase 5 of the Coupled Model Intercomparison Project (CMIP5) decadal hindcast/forecast experiments. The experiments are initialized in either November or January of each year and integrated for up to 10 years, providing a good opportunity for filling the gap between seasonal and decadal climate predictions. The long-lead multimodel ensemble (MME) prediction is evaluated for 1981–2007 in terms of the anomaly correlation coefficient (ACC) and mean-squared skill score (MSSS), which combines ACC and conditional bias, with respect to observations and reanalysis data, paying particular attention to the seasonal dependency of the global-mean and equatorial Pacific TAS predictions. The MME shows statistically significant ACCs and MSSSs for the annual global-mean TAS for up to two years, mainly because of long-term global warming trends. When the long-term trends are removed, the prediction skill is reduced. The prediction skills are generally lower in boreal winters than in other seasons regardless of lead times. This lack of winter prediction skill is attributed to the failure of capturing the long-term trend and interannual variability of TAS over high-latitude continents in the Northern Hemisphere. In contrast to global-mean TAS, regional TAS over the equatorial Pacific is predicted well in winter. This is mainly due to a successful prediction of the El Niño–Southern Oscillation (ENSO). In most models, the wintertime ENSO index is reasonably well predicted for at least one year in advance. The sensitivity of the prediction skill to the initialized month and method is also discussed.

This study investigates why OLR plays a small role in the Real-time Multivariate (Madden–Julian oscillation) MJO (RMM) index and how to improve it. The RMM index consists of the first two leading principal components (PCs) of a covariance matrix, which is constructed by combined daily anomalies of OLR and zonal winds at 850 (U850) and 200 hPa (U200) in the tropics after being normalized with their globally averaged standard deviations of 15.3 W m−2, 1.8 m s−1, and 4.9 m s−1, respectively. This covariance matrix is reasoned mathematically close to a correlation matrix. Both matrices substantially suppress the overall contribution of OLR and make the index more dynamical and nearly transparent to the convective initiation of the MJO. A covariance matrix that does not use normalized anomalies leads to the other extreme where OLR plays a dominant role while U850 and U200 are minor. Numerous tests indicate that a simple scaling of the anomalies (i.e., 2 W m−2, 1 m s−1, and 1 m s−1) can better balance the roles of OLR and winds. The revised PCs substantially enhance OLR over the eastern Indian and western Pacific Oceans and change it less notably in other locations, while they reduce U850 and U200 only slightly. Comparisons with the original RMM in spatial structure, power spectra, and standard deviation demonstrate improvements of the revised RMM index.

The Madden–Julian oscillation (MJO), the dominant mode of tropical intraseasonal variability, influences weather and climate in the extratropics through atmospheric teleconnection. In this study, two simulations using the Community Atmosphere Model version 5 (CAM5)—one with the default shallow and deep convection schemes and the other with the unified convection scheme (UNICON)—are employed to examine the impacts of cumulus parameterizations on the simulation of the boreal wintertime MJO teleconnection in the Northern Hemisphere. It is demonstrated that the UNICON substantially improves the MJO teleconnection. When the UNICON is employed, the simulated circulation anomalies associated with the MJO better resemble the observed counterpart, compared to the simulation with the default convection schemes. Quantitatively, the pattern correlation for the 300-hPa geopotential height anomalies between the simulations and observation increases from 0.07 for the default schemes to 0.54 for the UNICON. These circulation anomalies associated with the MJO further help to enhance the surface air temperature and precipitation anomalies over North America, although room for improvement is still evident. Initial value calculations suggest that the realistic MJO teleconnection with the UNICON is not due to the changes in the background wind, but rather primarily to the improved tropical convective heating associated with the MJO.

The North Atlantic Oscillation (NAO) is the dominant large-scale atmospheric oscillation in the North Atlantic and has profound effects on water temperatures in the North Atlantic. In this paper we diagnosed the effects of the NAO on sea surface temperature (SST) in the Northeast US Continental Shelf (NES). Waters in the Gulf of Maine (GOM) originate primarily from Scotian Shelf Water (SSW) at the surface and Labrador Slope Water (LSW) and Warm Slope Water (WSW) at depth through the Northeast Channel. By using a high-resolution SST dataset, we found that the correlation between the NAO and annual mean SST in the GOM is significant and negative at lag of four years. Further spatial correlation analysis shows that the NAO influences SST in the GOM primarily through advection of SSW or shelf water at the surface from the Labrador Sea. Cross-correlation analysis was also applied between the NAO and SSTs in other subregions of the NES (Georges Bank, Southern New England, and Mid-Atlantic Bight), but no statistically significant relationships were found at any lags. Different from temperature at depth in the GOM that is positively influenced by the NAO with a lag of two years, we concluded that the NAO has a significant negative effect on SST in the GOM four years later, while its effects on SSTs in the other three subregions of the NES are negligible. The four-year lagged relationship we found between the NAO and annual mean SST in the GOM provides a robust empirical method to predict the effect of the NAO on annual mean SST in the GOM four years in advance.

The vertically integrated water vapor transport (IVT) over the Pacific–North American sector during three phases of ENSO in boreal winter (December–February) is investigated using IVT values calculated from the Climate Forecast System Reanalysis (CFSR) during 1979–2010. The shift of the location and sign of sea surface temperature (SST) anomalies in the tropical Pacific Ocean leads to different atmospheric responses and thereby changes the seasonal mean moisture transport into North America. During eastern Pacific El Niño (EPEN) events, large positive IVT anomalies extend northeastward from the subtropical Pacific into the northwestern United States following the anomalous cyclonic flow around a deeper Aleutian low, while a southward shift of the cyclonic circulation during central Pacific El Niño (CPEN) events induces the transport of moisture into the southwestern United States. In addition, moisture from the eastern tropical Pacific is transported from the deep tropical eastern Pacific into Mexico and the southwestern United States during CPEN. During La Niña (NINA), the seasonal mean IVT anomaly is opposite to that of two El Niño phases. Analyses of 6-hourly IVT anomalies indicate that there is strong moisture transport from the North Pacific into the northwestern and southwestern United States during EPEN and CPEN, respectively. The IVT is maximized on the southeastern side of a low located over the eastern North Pacific, where the low is weaker but located farther south and closer to shore during CPEN than during EPEN. Moisture enters the southwestern United States from the eastern tropical Pacific during NINA via anticyclonic circulation associated with a ridge over the southern United States.

This study attempts, for the first time, to predict the annual number of tropical cyclones (TCs) affecting New York State (NYS), as part of the effort of the New York State Resiliency Institute for Storms and Emergencies (RISE). A pure statistical prediction model and a statistical–dynamical hybrid prediction model have been developed based on the understanding of the physical mechanism between NYS TCs and associated large-scale climate variability. During the cold phase of El Niño–Southern Oscillation, significant circulation anomalies in the Atlantic Ocean provide favorable conditions for more recurving TCs into NYS. The pure statistical prediction model uses the sea surface temperature (SST) over the equatorial Pacific Ocean from the previous months. Cross validation shows that the correlation between the observed and predicted numbers of NYS TCs is 0.56 for the June 1979–2013 forecasts. Forecasts of the probability of one or more TCs impacting NYS have a Brier skill score of 0.35 compared to climatology. The statistical–dynamical hybrid prediction model uses Climate Forecast System, version 2, SST predictions, which are statistically downscaled to forecast the number of NYS TCs based on a stepwise regression model. Results indicate that the initial seasonal prediction for NYS TCs can be issued in February using the hybrid model, with an update in June using the pure statistical prediction model. Based on the statistical model, for 2014, the predicted number of TCs passing through NYS is 0.33 and the probability of one or more tropical cyclones crossing NYS is 30%, which are both below average and in agreement with the actual activity (0 NYS TCs).

The authors examine the predictability and prediction skill of the Madden–Julian oscillation (MJO) of two ocean–atmosphere coupled forecast systems of ECMWF [Variable Resolution Ensemble Prediction System (VarEPS)] and NCEP [Climate Forecast System, version 2 (CFSv2)]. The VarEPS hindcasts possess five ensemble members for the period 1993–2009 and the CFSv2 hindcasts possess three ensemble members for the period 2000–09. Predictability and prediction skill are estimated by the bivariate correlation coefficient between the observed and predicted Wheeler–Hendon real-time multivariate MJO index (RMM). MJO predictability is beyond 32 days lead time in both hindcasts, while the prediction skill is about 27 days in VarEPS and 21 days in CFSv2 as measured by the bivariate correlation exceeding 0.5. Both predictability and prediction skill of MJO are enhanced by averaging ensembles. Results show clearly that forecasts initialized with (or targeting) strong MJOs possess greater prediction skill compared to those initialized with (or targeting) weak or nonexistent MJOs. The predictability is insensitive to the initial MJO phase (or forecast target phase), although the prediction skill varies with MJO phases. A few common model issues are identified. In both hindcasts, the MJO propagation speed is slower and the MJO amplitude is weaker than observed. Also, both ensemble forecast systems are underdispersive, meaning that the growth rate of ensemble error is greater than the growth rate of the ensemble spread by lead time.

The prediction skill and errors in surface temperature anomalies in initialized decadal hindcasts from phase 5 of the Coupled Model Intercomparison Project (CMIP5) are assessed using six ocean–atmosphere coupled models initialized every year from 1961 to 2008. The initialized hindcasts show relatively high prediction skill over the regions where external forcing dominates, indicating that a large portion of the prediction skill is due to the long-term trend. After removing the linear trend, high prediction skill is shown near the centers of action of the dominant decadal climate oscillations, such as the Pacific decadal oscillation (PDO) and Atlantic multidecadal oscillation (AMO). Low prediction skill appears over the tropical and eastern North Pacific Ocean where the predicted anomaly patterns associated with the PDO are systematically different in model and observations. By statistically correcting those systematic errors using a stepwise pattern projection method (SPPM) based on the data in an independent training period, the prediction skill of sea surface temperature (SST) is greatly enhanced over the North Pacific Ocean. The SST prediction skill over the North Pacific Ocean after the SPPM error correction is as high as that over the North Atlantic Ocean. In addition, the prediction skill in a single model after correction exceeds the skill of the multimodel ensemble (MME) mean before correction, implying that the MME method is not as effective in addressing systematic errors as the SPPM correction.

This study assesses the prediction skill of the boreal winter Arctic Oscillation (AO) in the state-of-the-art dynamical ensemble prediction systems (EPSs): the UKMO GloSea4, the NCEP CFSv2, and the NASA GEOS-5. Long-term reforecasts made with the EPSs are used to evaluate representations of the AO, and to examine skill scores for the deterministic and probabilistic forecast of the AO index. The reforecasts reproduce the observed changes in the large-scale patterns of the Northern Hemispheric surface temperature, upper-level wind, and precipitation according to the AO phase. Results demonstrate that all EPSs have better prediction skill than the persistence prediction for lead times up to 3-month, suggesting a great potential for skillful prediction of the AO and the associated climate anomalies in seasonal time scale. It is also found that the deterministic and probabilistic forecast skill of the AO in the recent period (1997-2010) is higher than that in the earlier period (1983-1996).

The relationship between El Niño–Southern Oscillation (ENSO) and tropical storm (TS) activity over the western North Pacific Ocean is examined for the period from 1981 to 2010. In El Niño years, TS genesis locations are generally shifted to the southeast relative to normal years and the passages of TSs tend to recurve to the northeast. TSs of greater duration and more intensity during an El Niño summer induce an increase of the accumulated tropical cyclone kinetic energy (ACE). Based on the strong relationship between the TS properties and ENSO, a probabilistic prediction for seasonal ACE is investigated using a hybrid dynamical–statistical model. A statistical relationship is developed between the observed ACE and large-scale variables taken from the ECMWF seasonal forecast system 4 hindcasts. The ACE correlates positively with the SST anomaly over the central to eastern Pacific and negatively with the vertical wind shear near the date line. The vertical wind shear anomalies over the central and western Pacific are selected as predictors based on sensitivity tests of ACE predictive skill. The hybrid model performs quite well in forecasting seasonal ACE with a correlation coefficient between the observed and predicted ACE at 0.80 over the 30-yr period. A relative operating characteristic analysis also indicates that the ensembles have significant probabilistic skill for both the above-normal and below-normal categories. By comparing the ACE prediction over the period from 2003 to 2011, the hybrid model appears more skillful than the forecast from the Tropical Storm Risk consortium.

The seasonal prediction skill for the Northern Hemisphere winter is assessed using retrospective predictions (1982–2010) from the ECMWF System 4 (Sys4) and National Center for Environmental Prediction (NCEP) CFS version 2 (CFSv2) coupled atmosphere–ocean seasonal climate prediction systems. Sys4 shows a cold bias in the equatorial Pacific but a warm bias is found in the North Pacific and part of the North Atlantic. The CFSv2 has strong warm bias from the cold tongue region of the eastern Pacific to the equatorial central Pacific and cold bias in broad areas over the North Pacific and the North Atlantic. A cold bias in the Southern Hemisphere is common in both reforecasts. In addition, excessive precipitation is found in the equatorial Pacific, the equatorial Indian Ocean and the western Pacific in Sys4, and in the South Pacific, the southern Indian Ocean and the western Pacific in CFSv2. A dry bias is found for both modeling systems over South America and northern Australia. The mean prediction skill of 2 meter temperature (2mT) and precipitation anomalies are greater over the tropics than the extra-tropics and also greater over ocean than land. The prediction skill of tropical 2mT and precipitation is greater in strong El Nino Southern Oscillation (ENSO) winters than in weak ENSO winters. Both models predict the year-to-year ENSO variation quite accurately, although sea surface temperature trend bias in CFSv2 over the tropical Pacific results in lower prediction skill for the CFSv2 relative to the Sys4. Both models capture the main ENSO teleconnection pattern of strong anomalies over the tropics, the North Pacific and the North America. However, both models have difficulty in forecasting the year-to-year winter temperature variability over the US and northern Europe.

The predictability of intraseasonal variation in the tropics is assessed in the present study by using various statistical and dynamical models with rigorous and fair measurements. For a fair comparison, the real-time multivariate Madden–Julian oscillation (MJO) (RMM) index, proposed by Wheeler and Hendon, is used as a predictand for all models. The statistical models include the models based on a multilinear regression, a wavelet analysis, and a singular spectrum analysis (SSA). The prediction limits (correlation skill of 0.5) of statistical models for RMM1 (RMM2) index are at days 16–17 (14–15) for the multiregression model, whereas they are at days 8–10 (9–12) for the wavelet- and SSA-based models. The poor predictability of the wavelet and SSA models is related to the tapering problem for a half-length of the time window before the initial condition. To assess the dynamical predictability, long-term serial prediction experiments with a prediction interval of every 5 days are carried out with Seoul National University (SNU) AGCM and coupled general circulation model (CGCM) for 26 (1980–2005) boreal winters. The prediction limits of RMM1 and RMM2 occur at around 20 days for both AGCM and CGCM. These results demonstrate that the skills of dynamical models used in this study are better than those of the three statistical predictions. The dynamical and statistical predictions are combined using a multimodel ensemble method. The combination provides a superior skill to any of the statistical and dynamical predictions, with a prediction limit of 22–24 days. The dependencies of prediction skill on the initial phase and amplitude of the MJO are also investigated.

The influence of ocean-atmosphere coupling on the simulation and prediction of the boreal winter Madden-Julian Oscillation (MJO) is examined using the Seoul National University coupled general circulation model (CGCM) and atmospheric—only model (AGCM). The AGCM is forced with daily SSTs interpolated from pentad mean CGCM SSTs. Forecast skill is examined using serial extended simulations spanning 26 different winter seasons with 30-day forecasts commencing every 5 days providing a total of 598 30-day simulations. By comparing both sets of experiments, which share the same atmospheric components, the influence of coupled ocean- atmosphere processes on the simulation and prediction of MJO can be studied. The mean MJO intensity possesses more realistic amplitude in the CGCM than in AGCM. In general, the ocean-atmosphere coupling acts to improve the simulation of the spatio-temporal evolution of the eastward propagating MJO and the phase relationship between convection (OLR) and SST over the equatorial Indian Ocean and the western Pacific. Both the CGCM and observations exhibit a near-quadrature relationship between OLR and SST, with the former lagging by about two pentads. However, the AGCM shows a less realistic phase relationship. As the initial conditions are the same in both models, the additional forcing by SST anomalies in the CGCM extends the prediction skill beyond that of the AGCM. To test the applicability of the CGCM to real-time prediction, we compute the Real-time Multivariate MJO (RMM) index and compared it with the index computed from observations. RMM1 (RMM2) falls away rapidly to 0.5 after 17-18 (15-16) days in the AGCM and 18-19 (16- 17) days in the CGCM. The prediction skill is phase dependent in both the CGCM and AGCM.

The influence of sea surface temperature (SST) on the simulation and predictability of the Madden–Julian oscillation (MJO) is examined using the Seoul National University atmospheric general circulation model (SNU AGCM). Forecast skill was examined using serial climate simulations spanning eight different winter seasons with 30-day forecasts commencing every 5 days, giving a total of 184 thirty-day simulations. The serial runs were repeated using prescribing observed SST with monthly, weekly, and daily temporal resolutions. The mean SST was the same for all cases so that differences between experiments result from the different temporal resolutions of the SST boundary forcing. It is shown that high temporal SST frequency acts to improve 1) the MJO activity of 200-hPa velocity potential field over the entire Asian monsoon region at all lead times; 2) the percentage of filtered variance of the two leading EOF modes that explain the eastward propagation of MJO; 3) the power of the wavenumber 1 eastward propagating mode; and 4) the forecast skill of MJO, maintaining it for longer periods. However, the MJO phase relationship between MJO convection and SST, as is often the case with many atmosphere-only models, although well simulated at the beginning of forecast period becomes distorted rapidly as the forecast lead time increases, even with the daily SST forcing case. Comparison of AGCM simulations with coupled GCM (CGCM) integrations shows that ocean–atmosphere coupling improves considerably the phase relationship between SST and convection. The CGCM results reinforce that the MJO is a coupled phenomenon and suggest strongly the need of the ocean–atmosphere coupled processes to extend predictability.

The impact of ocean–atmosphere coupling on the simulation and prediction of the boreal summer intraseasonal oscillation (ISO) has been investigated by diagnosing 22-year retrospective forecasts using the Seoul National University coupled general circulation model (CGCM) and its atmospheric GCM (AGCM) forced with SSTs derived from the CGCM. Numerous studies have shown that the ocean–atmosphere coupling has a significant effect on the improvement of ISO simulation and prediction. Contrary to previous studies, this study shows similar results between CGCM and AGCM, not only in regard to the ISO simulation characteristics but also the predictability. The similarities between CGCM and AGCM include (1) the ISO intensity over the entire Asian-monsoon region; (2) the spatiotemporal evolution of the northward propagating ISO (NPISO); and (3) the potential and practical predictability. A notable difference between CGCM and AGCM is the phase relationship between precipitation and SST anomalies. The CGCM and observation exhibits a near-quadrature relationship between precipitation and SST, with the former lagging about two pentads. The AGCM shows a less realistic phase relationship. The similar structure and propagation characteristics of ISO between the CGCM and AGCM suggest that the internal atmospheric dynamics could be more essential to the ISO than the ocean–atmosphere interaction over the Indian monsoon region.

The reproducibility of boreal summer intraseasonal variability (ISV) and its interannual variation by dynamical models are assessed through diagnosing 21-year retrospective forecasts from ten state-of-the-art ocean–atmosphere coupled prediction models. To facilitate the assessment, we have defined the strength of ISV activity by the standard deviation of 20–90 days filtered precipitation during the boreal summer of each year. The observed climatological ISV activity exhibits its largest values over the western North Pacific and Indian monsoon regions. The notable interannual variation of ISV activity is found primarily over the western North Pacific in observation while most models have the largest variability over the central tropical Pacific and exhibit a wide range of variability in spatial patterns that are different from observation. Although the models have large systematic biases in spatial pattern of dominant variability, the leading EOF modes of the ISV activity in the models are closely linked to the models’ El Nino-Southern Oscillation (ENSO), which is a feature that resembles the observed ISV and ENSO relationship. The ENSO-induced easterly vertical shear anomalies in the western and central tropical Pacific, where the summer mean vertical wind shear is weak, result in ENSO-related changes of ISV activity in both observation and models. It is found that the principal components of the predicted dominant modes of ISV activity fluctuate in a very similar way with observed ones. The model biases in the dominant modes are systematic and related to the external SST forcing. Thus the statistical correction method of this study based on singular value decomposition is capable of removing a large portion of the systematic errors in the predicted spatial patterns. The 21-year-averaged pattern correlation skill increases from 0.25 to 0.65 over the entire Asian monsoon region after applying the bias correction method to the multi-model ensemble mean prediction.

This study examines the influence of the Madden-Julian Oscillation (MJO) on tropical storm (TS) activity in the western North Pacific, using observations and GEOS-5 simulations at 50-km horizontal resolution. While GEOS-5 produces an MJO of faster propagation and weaker amplitude, it nevertheless reproduces the observed modulation of TS activity by the MJO with the highest TS genesis and increased track density in the active phases of MJO. The study suggests that the simulation of the sub-seasonal variability of TS activity could be improved by improving the simulations of the MJO in climate models.

[1] This study assesses the CMIP5 decadal hindcast/forecast simulations of seven state-of-the-art ocean-atmosphere coupled models. Each decadal prediction consists of simulations over a 10 year period each of which are initialized every five years from climate states of 1960/1961 to 2005/2006. Most of the models overestimate trends, whereby the models predict less warming or even cooling in the earlier decades compared to observations and too much warming in recent decades. All models show high prediction skill for surface temperature over the Indian, North Atlantic and western Pacific Oceans where the externally forced component and low-frequency climate variability is dominant. However, low prediction skill is found over the equatorial and North Pacific Ocean. The Atlantic Multidecadal Oscillation (AMO) index is predicted in most of the models with significant skill, while the Pacific Decadal Oscillation (PDO) index shows relatively low predictive skill. The multi-model ensemble has in general better-forecast quality than the single-model systems for global mean surface temperature, AMO and PDO.

During July 2010, a series of monsoonal deluges over northern Pakistan resulted in catastrophic flooding, loss of life and property and an agricultural crisis that may last for years. Was the rainfall abnormal compared to previous years? Furthermore, could a high probability of flooding have been predicted? To address these questions, regional precipitation is analyzed using three dataset sets covering the 1981–2010 time period. It is concluded that the 2010 average May to August (MJJA) rainfall for year 2010 is somewhat greater in magnitude than previous years. However, the rainfall rate of the July deluges, especially in North Pakistan was exceptionally rare as deduced from limited data. The location of the deluges over the mountainous northern part of the country lead to the devastating floods. The European Centre for Medium Range Weather Forecasts (ECMWF) 15-day Ensemble Prediction System (EPS) is used to assess whether the rainfall over the flood affected region was predictable. A multi-year analysis shows that Pakistan rainfall is highly predictable out to 6–8 days including rainfall in the summer of 2010. We conclude that if these extended quantitative precipitation forecasts had been available in Pakistan, the high risk of flooding could have been foreseen. If these rainfall forecasts had been coupled to a hydrological model then the high risk of extensive and prolonged flooding could have anticipated and actions taken to mitigate their impact.