Coastal flooding operational forecasting in the US is limited to short-range temporal scales (3–7 days), which limits the response time for emergency preparation and planning. The sub-seasonal prediction project (SubX), which produces weather forecasts with a lead time of up to four weeks, provides an opportunity to assess the potential for creating probabilistic flood forecasts with longer lead times. Using the ADCIRC hydrodynamic model for coastal storm surge, two major hurricanes, Isabel (2003) and Katrina (2005), were used as case studies to test coastal flood predictions induced by wind and pressure fields generated from five global weather models within SubX. The storm surges simulations are forced by Sea Level Pressure (SLP) and 10 m winds fields from SubX models for a lead-time of up to 30 days before storm landfall. The subseasonal surge forecasts are evaluated temporally and spatially at 1–4 weeks lead-time against the NOAA tide gages observations and a verification dataset derived by forcing the storm surge model with wind and pressure fields from the NCEP-Reanalysis. The results are evaluated in terms of lead-time and forecast skill metrics. The storm surge forecast skill is measured using the mean square error skill score (MSESS) relative to the verification dataset and an approximate of the climatology. A skill score greater than 0.55 is considered here useful for flood forecasting. The multi-model ensemble (MME) mean surge forecasts driven by several members of SubX models demonstrate skill greater than 0.55 up to a 4-day and 10-day lead for Katrina and Isabel, respectively. A sharper decrease in MSESS was noted from week 1 to week 3 lead-times for Katrina, in comparison to Isabel. Some ensemble members forecasted hurricanes and storm surges as early as 3–4 weeks lead-time. However, due to the offsets developed in the timing and magnitude of the peak at these lead-times, and based on a sample size of only two events, it is hard to establish the significance of these longer lead-time results. While a follow up study involving flood reforecasts over the entire SubX reforecast period (1999–2015) is needed to support more robust statistics of the forecast skill, our case studies demonstrate the feasibility of probabilistic flood forecasting at subseasonal timescales using the SubX models.

Tropical cyclone (TC) landfalls over the U.S. mid-Atlantic region are very infrequent. However, when they do occur, the resulting human and material losses can be severe, as was the case with Hurricane Sandy in 2012. Therefore, it is important to predict these land-falling events as far in advance as possible. In this study, we investigate the relationship between mid-Atlantic TC landfalls and the Madden-Julian Oscillation (MJO), which is the dominant source of climate variability in the tropics on intraseasonal time scales. This is largely accomplished by using a high-atmospheric-resolution ensemble prediction system based on the European Centre for Medium-Range Weather Forecasts (ECMWF) operational model (Project Minerva) to compile the statistics of these rare events, and the velocity potential MJO (VPM) index to define the phase and amplitude of the MJO. We find that at longer lead times (between 14 and 7 days prior to landfall), statistically significant peak landfall probabilities are present during MJO phases 1 and 7 and, to some extent, phase 8. This result is largely supported by observational data. At shorter lead times (between 6-day lead and landfall), phase 1 is strongly favored with some contribution from phases 2 and 3. These findings suggest a potential for extended-range predictions of the mid-Atlantic TC landfall risk based on the phase of the MJO.

Tropical cyclone (TC) landfalls over the U.S. mid-Atlantic region, which include the so-called Sandy-like, or westward-curving, tracks, are among the most infrequent landfalls along the U.S. East Coast. However, when these events do occur, the resulting economic and societal consequences can be devastating. A recent example is Hurricane Sandy in 2012. Multimodel ensemble seasonal hindcasts conducted with a high-atmospheric-resolution coupled prediction system based on the ECMWF operational model (Project Minerva) are used here to compile the statistics of these rare events. Minerva hindcasts are found to exhibit skill in reproducing climatological characteristics of the mid-Atlantic TC landfalls particularly at the highest atmospheric horizontal spectral resolution of T1279 (16-km grid spacing). Historical forecasts are further interrogated to identify regional and large-scale environmental conditions associated with these rare TC tracks to better quantify their predictability on synoptic time scales, and their dependence on model resolution. Evolution of the large-scale atmospheric flow patterns leading to mid-Atlantic TC landfalls is analyzed using local finite-amplitude wave activity (LWA). We have identified large-amplitude quasi-stationary features in the LWA and sea surface temperature (SST) anomaly distributions that persist up to about a week leading to these land-falling events. A statistical model utilizing indices based on the LWA and SST anomalies as predictors is developed that exhibits skill (mostly at T1279) in predicting mid-Atlantic TC landfalls several days in advance. Implications of these results for longer time-scale predictions of mid-Atlantic TC landfalls including climate change projections are discussed.

The North American Multi-Model Ensemble (NMME)-Phase II models are evaluated in terms of their retrospective seasonal forecast skill of the North Atlantic (NA) tropical cyclone (TC) activity, with a focus on TC frequency. The TC identification and tracking algorithm is modified to accommodate model data at daily resolution. It is also applied to three reanalysis products at the spatial and temporal resolution of the NMME-Phase II ensemble to allow for a more objective estimation of forecast skill. When used with the reanalysis data, the TC tracking generates realistic climatological distributions of the NA TC formation and tracks, and represents the interannual variability of the NA TC frequency quite well. Forecasts with the multi-model ensemble (MME) when initialized in April and later tend to have skill in predicting the NA seasonal TC counts (and TC days). At longer leads, the skill is low or marginal, although one of the models produces skillful forecasts when initialized as early as January and February. At short lead times, while demonstrating the highest skill levels the MME also tends to significantly outperform the individual models and attain skill comparable to the reanalysis. In addition, the short-lead MME forecasts are quite reliable. At regional scales, the skill is rather limited and mostly present in the western tropical NA and the Caribbean Sea. It is found that the overall MME forecast skill is limited by poor representation of the low-frequency variability in the predicted TC frequency, and large fluctuations in skill on decadal time scales. Addressing these deficiencies is thought to increase the value of the NMME ensemble in providing operational guidance.

Systematic error and forecast skill for temperature and precipitation in two regions of Southern Asia are investigated using hindcasts initialized May 1 from the North American Multi-Model Ensemble. We focus on two contiguous but geographically and dynamically diverse regions: the Extended Indian Monsoon Rainfall (70-100E, 10-30 N) and the nearby mountainous area of Pakistan and Afghanistan (60-75E, 23-39 N). Forecast skill is assessed using the Sign test framework, a rigorous statistical method that can be applied to non-Gaussian variables such as precipitation and to different ensemble sizes without introducing bias. We find that models show significant systematic error in both precipitation and temperature for both regions. The multi-model ensemble mean (MMEM) consistently yields the lowest systematic error and the highest forecast skill for both regions and variables. However, we also find that the MMEM consistently provides a statistically significant increase in skill over climatology only in the first month of the forecast. While the MMEM tends to provide higher overall skill than climatology later in the forecast, the differences are not significant at the 95% level. We also find that MMEMs constructed with a relatively small number of ensemble members per model can equal or outperform MMEMs constructed with more members in skill. This suggests some ensemble members either provide no contribution to overall skill or even detract from it.

An updated version of the Heated Condensation Framework (HCF) is implemented as a convective triggering criterion into the National Centers for Environmental Prediction (NCEP) Climate Forecast System version 2 (CFSv2). The new trigger replaces the original criteria in both the deep (Simplified Arakawa-Schubert - SAS) and shallow (SAS based) convective schemes. The performance of the original and new triggering criteria is first compared against radiosonde observations. Then, a series of hindcasts are performed to evaluate the influence of the triggering criterion in the CFSv2 representation of summer precipitation, the diurnal cycle of precipitation, and hurricanes that made landfall. The observational analysis shows that the HCF trigger better captures the frequency of convection, where the original SAS trigger initiates convection too often. When implemented in CFSv2, the HCF trigger improves the seasonal forecast of the Indian summer monsoon rainfall, including the representation of the onset dates of the rainy season over India. On the other hand, the HCF trigger increases error in the seasonal forecast of precipitation over the eastern United States. The HCF trigger also improves the representation of the intensity of hurricanes. Moreover, the simulation of hurricanes provides insights on the mechanism whereby the HCF trigger impacts the representation of convection.

The impact of sampling variability on the correlation between all-India rainfall (AIR) and the El Niño-Southern Oscillation is investigated in a large ensemble of seasonal climate simulations made using the European Centre for Medium-Range Weather Forecasting Ensemble Prediction System at T319 (64 km). The analyzed runs consist of 51 ensemble members initialized each May 1 for the period 1980–2011 and integrated for 7 months. 10,000 pairs of 32-year timeseries of June–September (JJAS) mean AIR and NINO3 indices are created from this database by randomly drawing one of the 51 ensemble members for each year. The correlation between each pair of AIR and NINO3 series is then calculated, generating a distribution of AIR–NINO3 correlation values. The model is reinitialized with observations each May 1 and thus all members are drawn from the same background state by construction and any differences in correlation are attributable to sampling variability. The spread in the calculated correlation values and the differences between 32-year segments are sufficient to explain the observed variations in AIR–NINO3 correlation since the beginning of the 1900s, including the sharp decrease in correlation strength since the late 1970s. Sampling variability thus represents a strong null hypothesis for the observed changes and one that cannot be rejected at the 95 % level based on our simulations.

Regional variations in seasonal mean Indian summer monsoon rainfall and circulation for the period 1979 2009 are investigated using multiple data products. The focus is on four separate regions: the Western Ghats (WG), the Ganges basin (GB), the Bay of Bengal (BB), and Bangladesh northeastern India (BD). Data reliability varies strongly by region, with particularly low correlations between different products for the BB and BD regions. Correlations between regions are generally not statistically significant, indicating rainfall varies independently in these four regions. The diagnosed associations between rainfall, circulation, and sea surface temperatures can be sensitive to the choice of rainfall product, and multiple precipitation products may need to be analyzed in this region to ensure that the results are robust.Enhanced precipitation in the BD region is associated with anomalous anticyclonic circulation at 850 mb and westerly anomalies along the foothills of the Tibetan Plateau, while precipitation in the other regions is associated with cyclonic flow and easterlies. These associations provide a dynamical explanation for previously reported weak, negative correlations between BD and the other regions.In addition to observed products, atmosphere-only simulations made using the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS) during Project Athena are analyzed. While the simulations do not reproduce the observed interannual variations in rainfall, the fidelity of the simulated precipitation and circulation structure is comparable to or even outperforms the different state-of-the-art reanalysis products considered. Accuracy in representing interannual variability and regional structure thus appears to be independent.

This study examines El Nino Southern Oscillation (ENSO) prediction in Project Minerva, a recent collaboration between the Center for Ocean-Land-Atmosphere Studies (COLA) and the European Centre for Medium-Range Weather Forecasts (ECMWF). The focus is primarily on the impact of the atmospheric horizontal resolution on ENSO prediction, but the effect from different ensemble sizes is also discussed. Particularly, three sets of 7-month hindcasts performed with ECMWF prediction system are compared, starting from 1 May (1 November) during 1982-2011 (1982-2010): spectral T319 atmospheric resolution with 15 ensembles, spectral T639 with 15 ensembles, and spectral T319 with 51 ensembles. The analysis herein shows that simply increasing either ensemble size from 15 to 51 or atmospheric horizontal resolution from T319 to T639 does not necessarily lead to major improvement in the ENSO prediction skill with current climate models. For deterministic prediction skill metrics, the three sets of predictions do not produce a significant difference in either anomaly correlation or root-mean-square error (RMSE). For probabilistic metrics, the increased atmospheric horizontal resolution generates larger ensemble spread, and thus increases the ratio between the intraensemble spread and RMSE. However, there is little change in the categorical distributions of predicted SST anomalies, and consequently there is little difference among the three sets of hindcasts in terms of probabilistic metrics or prediction reliability.

Seasonal forecast skill of the basinwide and regional tropical cyclone (TC) activity in an experimental coupled prediction system based on the ECMWF System 4 is assessed. As part of a collaboration between the Center for Ocean-Land-Atmosphere Studies (COLA) and the ECMWF called Project Minerva, the system is integrated at the atmospheric horizontal spectral resolutions of T319, T639, and T1279. Seven-month hindcasts starting from 1 May for the years 1980-2011 are produced at all three resolutions with at least 15 ensemble members. The Minerva system demonstrates statistically significant skill for retrospective forecasts of TC frequency and accumulated cyclone energy (ACE) in the North Atlantic (NA), eastern North Pacific (EP), and western North Pacific. While the highest scores overall are achieved in the North Pacific, the skill in the NA appears to be limited by an overly strong influence of the tropical Pacific variability. Higher model resolution improves skill scores for the ACE and, to a lesser extent, the TC frequency, even though the influence of large-scale climate variations on these TC activity measures is largely independent of resolution changes. The biggest gain occurs in transition from T319 to T639. Significant skill in regional TC forecasts is achieved over broad areas of the Northern Hemisphere. The highest-resolution hindcasts exhibit additional locations with skill in the NA and EP, including land-adjacent areas. The feasibility of regional intensity forecasts is assessed. In the presence of the coupled model biases, the benefits of high resolution for seasonal TC forecasting may be underestimated.

How tropical cyclone (TC) activity in the northwestern Pacific might change in a future climate is assessed using multidecadal Atmospheric Model Intercomparison Project (AMIP)-style and time-slice simulations with the ECMWF Integrated Forecast System (IFS) at 16-km and 125-km global resolution. Both models reproduce many aspects of the present-day TC climatology and variability well, although the 16-km IFS is far more skillful in simulating the full intensity distribution and genesis locations, including their changes in response to El Nino-Southern Oscillation. Both IFS models project a small change in TC frequency at the end of the twenty-first century related to distinct shifts in genesis locations. In the 16-km IFS, this shift is southward and is likely driven by the southeastward penetration of the monsoon trough/subtropical high circulation system and the southward shift in activity of the synoptic-scale tropical disturbances in response to the strengthening of deep convective activity over the central equatorial Pacific in a future climate. The 16-km IFS also projects about a 50% increase in the power dissipation index, mainly due to significant increases in the frequency of the more intense storms, which is comparable to the natural variability in the model. Based on composite analysis of large samples of supertyphoons, both the development rate and the peak intensities of these storms increase in a future climate, which is consistent with their tendency to develop more to the south, within an environment that is thermodynamically more favorable for faster development and higher intensities. Coherent changes in the vertical structure of supertyphoon composites show system-scale amplification of the primary and secondary circulations with signs of contraction, a deeper warm core, and an upward shift in the outflow layer and the frequency of the most intense updrafts. Considering the large differences in the projections of TC intensity change between the 16-km and 125-km IFS, this study further emphasizes the need for high-resolution modeling in assessing potential changes in TC activity.

The importance of using dedicated high-end computing resources to enable high spatial resolution in global climate models and advance knowledge of the climate system has been evaluated in an international collaboration called Project Athena. Inspired by the World Modeling Summit of 2008 and made possible by the availability of dedicated high-end computing resources provided by the National Science Foundation from October 2009 through March 2010, Project Athena demonstrated the sensitivity of climate simulations to spatial resolution and to the representation of subgrid-scale processes with horizontal resolutions up to 10 times higher than contemporary climate models. While many aspects of the mean climate were found to be reassuringly similar, beyond a suggested minimum resolution, the magnitudes and structure of regional effects can differ substantially. Project Athena served as a pilot project to demonstrate that an effective international collaboration can be formed to efficiently exploit dedicated supercomputing resources. The outcomes to date suggest that, in addition to substantial and dedicated computing resources, future climate modeling and prediction require a substantial research effort to efficiently explore the fidelity of climate models when explicitly resolving important atmospheric and oceanic processes.

Project Athena is an international collaboration testing the efficacy of high-resolution global climate models. We compare results from 7-km mesh experiments of the Nonhydrostatic Icosahedral Atmospheric Model (NICAM) and 10-km mesh experiments of the Integrated Forecast System (IFS), focusing on the Intra-Seasonal Oscillation (ISO) and its relationship with tropical cyclones (TC) among the boreal summer period (21 May–31 Aug) of 8 years (2001–2002, 2004–2009). In the first month of simulation, both models capture the intra-seasonal oscillatory behavior of the Indian monsoon similar to the observed boreal summer ISO in approximately half of the 8-year samples. The IFS simulates the NW–SE-oriented rainband and the westerly location better, while NICAM marginally reproduces mesoscale organized convective systems and better simulates the northward migration of the westerly peak and precipitation, particularly in 2006. The reproducibility of the evolution of MJO depends on the given year; IFS simulates the MJO signal well for 2002, while NICAM simulates it well for 2006. An empirical orthogonal function analysis shows that both models statistically reproduce MJO signals similar to observations, with slightly better phase speed reproduced by NICAM. Stronger TCs are simulated in NICAM than in IFS, and NICAM shows a wind-pressure relation for TCs closer to observations. TC cyclogenesis is active during MJO phases 3 and 4 in NICAM as in observations. The results show the potential of high-resolution global atmospheric models in reproducing some aspects of the relationship between MJO and TCs and the statistical behavior of TCs.

Global simulations have been conducted with the European Centre for Medium-Range Weather Forecasts operational model run at T1279 resolution for multiple decades representing climate from the late twentieth and late twenty-first centuries. Changes in key components of the water cycle are examined, focusing on variations at short time scales. Metrics of coupling and feedbacks between soil moisture and surface fluxes and between surface fluxes and properties of the planetary boundary layer (PBL) are inspected. Features of precipitation and other water cycle trends from coupled climate model consensus projections are well simulated. Extreme 6-hourly rainfall totals become more intense over much of the globe, suggesting an increased risk for flash floods. Seasonal-scale droughts are projected to escalate over much of the subtropics and midlatitudes during summer, while tropical and winter droughts become less likely. These changes are accompanied by an increase in the responsiveness of surface evapotranspiration to soil moisture variations. Even though daytime PBL depths increase over most locations in the next century, greater latent heat fluxes also occur over most land areas, contributing a larger energy effect per unit mass of air, except over some semiarid regions. This general increase in land atmosphere coupling is represented in a combined metric as a "land coupling index" that incorporates the terrestrial and atmospheric effects together. The enhanced feedbacks are consistent with the precipitation changes, but a causal connection cannot be made without further sensitivity studies. Nevertheless, this approach could be applied to the output of traditional climate change simulations to assess changes in land atmosphere feedbacks.

Northern Hemisphere tropical cyclone (TC) activity is investigated in multiyear global climate simulations with the ECMWF Integrated Forecast System (IFS) at 10-km resolution forced by the observed records of sea surface temperature and sea ice. The results are compared to analogous simulations with the 16-, 39-, and 125-km versions of the model as well as observations. In the North Atlantic, mean TC frequency in the 10-km model is comparable to the observed frequency, whereas it is too low in the other versions. While spatial distributions of the genesis and track densities improve systematically with increasing resolution, the 10-km model displays qualitatively more realistic simulation of the track density in the western subtropical North Atlantic. In the North Pacific, the TC count tends to be too high in the west and too low in the east for all resolutions. These model errors appear to be associated with the errors in the large-scale environmental conditions that are fairly similar in this region for all model versions. The largest benefits of the 10-km simulation are the dramatically more accurate representation of the TC intensity distribution and the structure of the most intense storms. The model can generate a supertyphoon with a maximum surface wind speed of 68.4 m s −1 . The life cycle of an intense TC comprises intensity fluctuations that occur in apparent connection with the variations of the eyewall/rainband structure. These findings suggest that a hydrostatic model with cumulus parameterization and of high enough resolution could be efficiently used to simulate the TC intensity response (and the associated structural changes) to future climate change.

The effects of horizontal resolution and the treatment of convection on simulation of the diurnal cycle of precipitation during boreal summer are analyzed in several innovative weather and climate model integrations. The simulations include: season-long integrations of the Non-hydrostatic Icosahedral Atmospheric Model (NICAM) with explicit clouds and convection; year-long integrations of the operational Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts at three resolutions (125, 39 and 16 km); seasonal simulations of the same model at 10 km resolution; and seasonal simulations of the National Center for Atmospheric Research (NCAR) low-resolution climate model with and without an embedded two-dimensional cloud-resolving model in each grid box. NICAM with explicit convection simulates best the phase of the diurnal cycle, as well as many regional features such as rainfall triggered by advancing sea breezes or high topography. However, NICAM greatly overestimates mean rainfall and the magnitude of the diurnal cycle. Introduction of an embedded cloud model within the NCAR model significantly improves global statistics of the seasonal mean and diurnal cycle of rainfall, as well as many regional features. However, errors often remain larger than for the other higher-resolution models. Increasing resolution alone has little impact on the timing of daily rainfall in IFS with parameterized convection, yet the amplitude of the diurnal cycle does improve along with the representation of mean rainfall. Variations during the day in atmospheric prognostic fields appear quite similar among models, suggesting that the distinctive treatments of model physics account for the differences in representing the diurnal cycle of precipitation.