Search

The Kuroshio Extension (KE), an eastward-flowing jet located in the Pacific western boundary current system, exhibits prominent seasonal-to-decadal variability, which is crucial for understanding climate variations in the northern midlatitudes. We explore the representation and prediction skill for the KE in the GFDL SPEAR (Seamless System for Prediction and Earth System Research) coupled model. Two different approaches are used to generate coupled reanalyses and forecasts: 1) restoring the coupled model’s SST and atmospheric variables toward existing reanalyses, or 2) assimilating SST and subsurface observations into the coupled model without atmospheric assimilation. Both systems use an ocean model with 1° resolution and capture the largest sea surface height (SSH) variability over the KE region. Assimilating subsurface observations appears to be essential to reproduce the narrow front and related oceanic variability of the KE jet in the coupled reanalysis. We demonstrate skillful retrospective predictions of KE SSH variability in monthly (up to 1 year) and annual-mean (up to 5 years) KE forecasts in the seasonal and decadal prediction systems, respectively. The prediction skill varies seasonally, peaking for forecasts initialized in January and verifying in September due to the winter intensification of North Pacific atmospheric forcing. We show that strong large-scale atmospheric anomalies generate deterministic oceanic forcing (i.e., Rossby waves), leading to skillful long-lead KE forecasts. These atmospheric anomalies also drive Ekman convergence and divergence, which forms ocean memory, by sequestering thermal anomalies deep into the winter mixed layer that re-emerge in the subsequent autumn. The SPEAR forecasts capture the recent negative-to-positive transition of the KE phase in 2017, projecting a continued positive phase through 2022.

A subseasonal-to-seasonal (S2S) prediction system was recently developed using the GFDL Seamless System for Prediction and Earth System Research (SPEAR) global coupled model. Based on 20-yr hindcast results (2000–19), the boreal wintertime (November–April) Madden–Julian oscillation (MJO) prediction skill is revealed to reach 30 days measured before the anomaly correlation coefficient of the real-time multivariate (RMM) index drops to 0.5. However, when the MJO is partitioned into four distinct propagation patterns, the prediction range extends to 38, 31, and 31 days for the fast-propagating, slow-propagating, and jumping MJO patterns, respectively, but falls to 23 days for the standing MJO. A further improvement of MJO prediction requires attention to the standing MJO given its large gap with its potential predictability (38 days). The slow-propagating MJO detours southward when traversing the Maritime Continent (MC), and confronts the MC prediction barrier in the model, while the fast-propagating MJO moves across the central MC without this prediction barrier. The MJO diversity is modulated by stratospheric quasi-biennial oscillation (QBO): the standing (slow-propagating) MJO coincides with significant westerly (easterly) phases of QBO, partially explaining the contrasting MJO prediction skill between these two QBO phases. The SPEAR model shows its capability, beyond the propagation, in predicting their initiation for different types of MJO along with discrete precursory convection anomalies. The SPEAR model skillfully predicts the observed distinct teleconnections over the North Pacific and North America related to the standing, jumping, and fast-propagating MJO, but not the slow-propagating MJO. These findings highlight the complexities and challenges of incorporating MJO prediction into the operational prediction of meteorological variables.

This study investigates precipitation variability over the Arabian Peninsula (AP) during its wet season. The wet season is split into winter (November – February) and spring (March and April) seasons, and early (1950–1986) and late (1986–2021) periods to understand sub-seasonal characteristics of precipitation variability and long-term changes in global teleconnections. The first three Empirical Orthogonal Functions explain ~70% of the interannual wet season precipitation variance, which shows an increase (decrease) in the late period winter (spring). Linear regression of the sea surface temperatures and geopotential height onto associated principal components reveals many oceanic and atmospheric variability patterns, which exhibit significant differences between winter and spring and early and late periods. Further, linear regressions of AP precipitation onto 14 natural modes of climate variability reveal a complex network of global teleconnections. El Niño-Southern Oscillation (ENSO) is one of the key contributors to precipitation variability but considering ENSO diversity is crucial to fully understand its influence. While the direct ENSO influence only becomes robust after the 1980s, its indirect effect persists through projection onto atmospheric modes, such as East Atlantic West Russia Pattern and East Atlantic Mode, or inter-basin interaction (e.g., via the Indian Ocean). The Northern Hemisphere atmospheric modes also mediate influences of other natural modes in tropical Indian and Atlantic oceans and extra-tropical regions over the AP. Several precipitation teleconnections exhibit a shift in the 1980s. Some may be related to the introduction of satellite data, but further investigations are warranted to understand the causes of these shifts.

The excitation of the Pacific–North American (PNA) teleconnection pattern by the Madden–Julian oscillation (MJO) has been considered one of the most important predictability sources on subseasonal time scales over the extratropical Pacific and North America. However, until recently, the interactions between tropical heating and other extratropical modes and their relationships to subseasonal prediction have received comparatively little attention. In this study, a linear inverse model (LIM) is applied to examine the tropical–extratropical interactions. The LIM provides a means of calculating the response of a dynamical system to a small forcing by constructing a linear operator from the observed covariability statistics of the system. Given the linear assumptions, it is shown that the PNA is one of a few leading modes over the extratropical Pacific that can be strongly driven by tropical convection while other extratropical modes present at most a weak interaction with tropical convection. In the second part of this study, a two-step linear regression is introduced that leverages a LIM and large-scale climate variability to the prediction of hydrological extremes (e.g., atmospheric rivers) on subseasonal time scales. Consistent with the findings of the first part, most of the predictable signals on subseasonal time scales are determined by the dynamics of the MJO–PNA teleconnection while other extratropical modes are important only at the shortest forecast leads.

Compared to the Arctic, seasonal predictions of Antarctic sea ice have received relatively little attention. In this work, we utilize three coupled dynamical prediction systems developed at the Geophysical Fluid Dynamics Laboratory to assess the seasonal prediction skill and predictability of Antarctic sea ice. These systems, based on the FLOR, SPEAR_LO, and SPEAR_MED dynamical models, differ in their coupled model components, initialization techniques, atmospheric resolution, and model biases. Using suites of retrospective initialized seasonal predictions spanning 1992–2018, we investigate the role of these factors in determining Antarctic sea ice prediction skill and examine the mechanisms of regional sea ice predictability. We find that each system is capable of skillfully predicting regional Antarctic sea ice extent (SIE) with skill that exceeds a persistence forecast. Winter SIE is skillfully predicted 11 months in advance in the Weddell, Amundsen and Bellingshausen, Indian, and West Pacific sectors, whereas winter skill is notably lower in the Ross sector. Zonally advected upper ocean heat content anomalies are found to provide the crucial source of prediction skill for the winter sea ice edge position. The recently-developed SPEAR systems are more skillful than FLOR for summer sea ice predictions, owing to improvements in sea ice concentration and sea ice thickness initialization. Summer Weddell SIE is skillfully predicted up to 9 months in advance in SPEAR_MED, due to the persistence and drift of initialized sea ice thickness anomalies from the previous winter. Overall, these results suggest a promising potential for providing operational Antarctic sea ice predictions on seasonal timescales.

Positive precipitation biases over western North America have remained a pervasive problem in the current generation of coupled global climate models. These biases are substantially reduced, however, in a version of the Geophysical Fluid Dynamics Laboratory Forecast-Oriented Low Ocean Resolution (FLOR) coupled climate model with systematic sea surface temperature (SST) biases artificially corrected through flux adjustment. This study examines how the SST biases in the Atlantic and Pacific Oceans contribute to the North American precipitation biases. Experiments with the FLOR model in which SST biases are removed in the Atlantic and Pacific are carried out to determine the contribution of SST errors in each basin to precipitation statistics over North America. Tropical and North Pacific SST biases have a strong impact on northern North American precipitation, while tropical Atlantic SST biases have a dominant impact on precipitation biases in southern North America, including the western United States. Most notably, negative SST biases in the tropical Atlantic in boreal winter induce an anomalously strong Aleutian low and a southward bias in the North Pacific storm track. In boreal summer, the negative SST biases induce a strengthened North Atlantic subtropical high and Great Plains low-level jet. Each of these impacts contributes to positive annual mean precipitation biases over western North America. Both North Pacific and North Atlantic SST biases induce SST biases in remote basins through dynamical pathways, so a complete attribution of the effects of SST biases on precipitation must account for both the local and remote impacts.

In mountain terrain, well-configured high-resolution atmospheric models are able to simulate total annual rain and snowfall better than spatial estimates derived from in situ observational networks of precipitation gauges, and significantly better than radar or satellite-derived estimates. This conclusion is primarily based on comparisons with streamflow and snow in basins across the western United States and in Iceland, Europe, and Asia. Even though they outperform gridded datasets based on gauge networks, atmospheric models still disagree with each other on annual average precipitation and often disagree more on their representation of individual storms. Research to address these difficulties must make use of a wide range of observations (snow, streamflow, ecology, radar, satellite) and bring together scientists from different disciplines and a wide range of communities.

Midlatitude baroclinic waves drive extratropical weather and climate variations, but their predictability beyond 2 weeks has been deemed low. Here we analyze a large ensemble of climate simulations forced by observed sea surface temperatures (SSTs) and demonstrate that seasonal variations of baroclinic wave activity (BWA) are potentially predictable. This potential seasonal predictability is denoted by robust BWA responses to SST forcings. To probe regional sources of the potential predictability, a regression analysis is applied to the SST-forced large ensemble simulations. By filtering out variability internal to the atmosphere and land, this analysis identifies both well-known and unfamiliar BWA responses to SST forcings across latitudes. Finally, we confirm the model-indicated predictability by showing that an operational seasonal prediction system can leverage some of the identified SST-BWA relationships to achieve skillful predictions of BWA. Our findings help to extend long-range predictions of the statistics of extratropical weather events and their impacts.

Climate-influenced changes in hydrology affect water-food-energy security that may impact up to two billion people downstream of the High Mountain Asia (HMA) region. Changes in water supply affect energy, industry, transportation, and ecosystems (agriculture, fisheries) and as a result, also affect the region's social, environmental, and economic fabrics. Sustaining the highly interconnected food-energy-water nexus (FEWN) will be a fundamental and increasing challenge under a changing climate regime. High variability in topography and distribution of glaciated and snow-covered areas in the HMA region, and scarcity of high resolution (in-situ) data make it difficult to model and project climate change impacts on individual watersheds. We lack basic understanding of the spatial and temporal variations in climate, surface impurities in snow and ice such as black carbon and dust that alter surface albedo, and glacier mass balance and dynamics. These knowledge gaps create challenges in predicting where and when the impact of changes in river flow will be the most significant economically and ecologically. In response to these challenges, the United States National Aeronautics and Space Administration (NASA) established the High Mountain Asia Team (HiMAT) in 2016 to conduct research to address knowledge gaps. This paper summarizes some of the advances HiMAT made over the past 5 years, highlights the scientific challenges in improving our understanding of the hydrology of the HMA region, and introduces an integrated assessment framework to assess the impacts of climate changes on the FEWN for the HMA region. The framework, developed under a NASA HMA project, links climate models, hydrology, hydropower, fish biology, and economic analysis. The framework could be applied to develop scientific understanding of spatio-temporal variability in water availability and the resultant downstream impacts on the FEWN to support water resource management under a changing climate regime.

Extreme heat occurred over Alaska in June–July 2019, posing risks to infrastructure, ecosystem, and human health. It is vital to improve our understanding of the causes of such events and the extent to which anthropogenic forcing may alter their likelihood and magnitude. Here, we use multiple large ensembles of climate models, comprising thousands of simulated years, to investigate these issues. Our results suggest that the presence of anthropogenic radiative forcing increased the likelihood of the 2019 extreme heat event by as much as 6%. Further we show the rate of occurrence of such an extreme heat event is likely to substantially increase in the future with increasing levels of atmospheric greenhouse gases. While uncertainty in projected climate risk from model choice leads to a broad range of future extreme heat event probabilities, some models project that with rapidly increasing levels of greenhouse gases the likelihood of such events would exceed 75% by 2090.

Responses of tropical cyclones (TCs) to CO2 doubling are explored using coupled global climate models (GCMs) with increasingly refined atmospheric/land horizontal grids (~ 200 km, ~ 50 km and ~ 25 km). The three models exhibit similar changes in background climate fields thought to regulate TC activity, such as relative sea surface temperature (SST), potential intensity, and wind shear. However, global TC frequency decreases substantially in the 50 km model, while the 25 km model shows no significant change. The ~ 25 km model also has a substantial and spatially-ubiquitous increase of Category 3–4–5 hurricanes. Idealized perturbation experiments are performed to understand the TC response. Each model’s transient fully-coupled 2 × CO2 TC activity response is largely recovered by “time-slice” experiments using time-invariant SST perturbations added to each model’s own SST climatology. The TC response to SST forcing depends on each model’s background climatological SST biases: removing these biases leads to a global TC intensity increase in the ~ 50 km model, and a global TC frequency increase in the ~ 25 km model, in response to CO2-induced warming patterns and CO2 doubling. Isolated CO2 doubling leads to a significant TC frequency decrease, while isolated uniform SST warming leads to a significant global TC frequency increase; the ~ 25 km model has a greater tendency for frequency increase. Global TC frequency responds to both (1) changes in TC “seeds”, which increase due to warming (more so in the ~ 25 km model) and decrease due to higher CO2 concentrations, and (2) less efficient development of these“seeds” into TCs, largely due to the nonlinear relation between temperature and saturation specific humidity.

High-impact extratropical cyclones (ETCs) cause considerable damage along the northeast coast of the United States through strong winds and inundation, but these relatively rare events are difficult to analyze owing to limited historical records. Using a 1505-yr simulation from the GFDL FLOR coupled model, statistical analyses of extreme events are performed including exceedance probability computations to compare estimates from shorter segments to estimates that could be obtained from a record of considerable length. The most extreme events possess characteristics including exceptionally low central pressure, hurricane-force winds, and a large surge potential, which would greatly impact nearby regions. Return level estimates of metrics of ETC intensity using shorter, historical-length segments of the FLOR simulation are underestimated compared to levels determined using the full simulation. This indicates that if the underlying distributions of observed ETC metrics are similar to those of the 1505-yr FLOR distributions, the actual frequency of extreme ETC events could also be underestimated. Comparisons between FLOR and reanalysis products suggest that not all features of simulated high-impact ETCs are representative of observations. Spatial track densities are similar, but FLOR exhibits a negative bias in central pressure and a positive bias in wind speed, particularly for more intense events. Although the existence of these model biases precludes the quantitative use of model-derived return statistics as a substitute for those derived from shorter observational records, this work suggests that statistics from future models of higher fidelity could be used to better constrain the probability of extreme ETC events and their impacts.

This study shows that the frequency of North American summer hot days are skillfully predictable months in advance in the newly-developed GFDL (Geophysical Fluid Dynamics Laboratory) SPEAR (Seamless System for Prediction and EArth System Research) seasonal forecast system. We also demonstrate that climate change, the Pacific Decadal Oscillation (PDO), the Atlantic Multidecadal Oscillation (AMO) and atmosphere-land feedback all contribute to the seasonal predictive skill of the frequency of North American summer hot days. Using a statistical optimization method (Average Predictability Time) we identify two large-scale components of the frequency of North American summer hot days that are predictable with significant correlation skill. One component shows an increase in the frequency of summer hot days everywhere over North America and is highly predictable at least 9 months in advance. This component is related to a secular warming trend. Another predictable component shows largest loadings over the central U.S., and is significantly predictable 6 months ahead. This second component is related to the PDO and the AMO, and is significantly correlated with the central U.S. soil water. These findings have potential implications for predictions of North American summer hot days on seasonal time scales.

The recent multi-year 2015–2019 drought after a multi-decadal drying trend over Central America raises the question of whether anthropogenic climate change (ACC) played a role in exacerbating these events. While the occurrence of the 2015–2019 drought in Central America has been asserted to be associated with ACC, we lack an assessment of natural vs anthropogenic contributions. Here, we use five different large ensembles—including high-resolution ensembles (i.e., 0.5◦ horizontally)—to estimate the contribution of ACC to the probability of occurrence of the 2015–2019 event and the recent multi-decadal trend. The comparison of ensembles forced with natural and natural plus anthropogenic forcing suggests that the recent 40-year trend is likely associated with internal climate variability. However, the 2015–2019 rainfall deficit has been made more likely by ACC. The synthesis of the results from model ensembles supports the notion of a significant increase, by a factor of four, over the last century for the 2015–2019 meteorological drought to occur because of ACC. All the model results further suggest that, under intermediate and high emission scenarios, the likelihood of similar drought events will continue to increase substantially over the next decades. Alma Mater Studiorum - Universit`a di Bologna/[]/UNIBO/Italia UCR::Vicerrectoría de Investigación::Unidades de Investigación::Ciencias Básicas::Centro de Investigaciones Geofísicas (CIGEFI)

Extreme precipitation from the South-Asian monsoon season combines with significant topographic relief within the Himalayan region to cause landslides that result in hundreds to thousands of fatalities each year. While there are few consistent and publicly available in-situ estimates of rainfall across this region, satellite products and global climate models provide insight into the extreme precipitation patterns that may impact the frequency of landsliding. In this work, we analyzed several extreme precipitation indices using data from a global climate model and the satellite-based Tropical Rainfall Measuring Mission Multi-satellite Precipitation Analysis product to represent extreme precipitation over High Mountain Asia. We then compared the temporal distribution of extreme precipitation to a global database of landslides to better understand the spatiotemporal distribution of potential landslide triggering factors. We found that these indices successfully model the seasonality of landslide activity across the region, but other aspects of spatiotemporal variability require additional information and analysis before they can be applied more broadly.

High Mountain Asia is impacted by extreme monsoonal rainfall that triggers landslides in large proportions relative to global distributions, resulting in substantial human impacts and damage to infrastructure each year. Previous landslide research has qualitatively estimated how patterns in landslide activity may change based on climate change scenarios. We present the first quantitative view of potential modulation in future landslide activity over the High Mountain Asia region leveraging a new landslide hazard model and precipitation data from satellite and Global Climate Model sources. In doing so, we find that the rate of increase in landslide activity at the end of the century is expected to be greatest over areas covered by current glaciers and glacial lakes, potentially exacerbating the impacts of cascading hazards on populations downstream. This work demonstrates the potential of Global Climate Models and satellite-based precipitation estimates to characterize landslide hazards at time scales affected by climate change.

We document the development and simulation characteristics of the next generation modeling system for seasonal to decadal prediction and projection at the Geophysical Fluid Dynamics Laboratory (GFDL). SPEAR (Seamless System for Prediction and EArth System Research) is built from component models recently developed at GFDL—the AM4 atmosphere model, MOM6 ocean code, LM4 land model, and SIS2 sea ice model. The SPEAR models are specifically designed with attributes needed for a prediction model for seasonal to decadal time scales, including the ability to run large ensembles of simulations with available computational resources. For computational speed SPEAR uses a coarse ocean resolution of approximately 1.0° (with tropical refinement). SPEAR can use differing atmospheric horizontal resolutions ranging from 1° to 0.25°. The higher atmospheric resolution facilitates improved simulation of regional climate and extremes. SPEAR is built from the same components as the GFDL CM4 and ESM4 models but with design choices geared toward seasonal to multidecadal physical climate prediction and projection. We document simulation characteristics for the time mean climate, aspects of internal variability, and the response to both idealized and realistic radiative forcing change. We describe in greater detail one focus of the model development process that was motivated by the importance of the Southern Ocean to the global climate system. We present sensitivity tests that document the influence of the Antarctic surface heat budget on Southern Ocean ventilation and deep global ocean circulation. These findings were also useful in the development processes for the GFDL CM4 and ESM4 models.

Three consecutive dry winters (2015-2017) in southwestern South Africa (SSA) resulted in the Cape Town "Day Zero" drought in early 2018. The contribution of anthropogenic global warming to this prolonged rainfall deficit has previously been evaluated through observations and climate models. However, model adequacy and insufficient horizontal resolution make it difficult to precisely quantify the changing likelihood of extreme droughts, given the small regional scale. Here, we use a high-resolution large ensemble to estimate the contribution of anthropogenic climate change to the probability of occurrence of multiyear SSA rainfall deficits in past and future decades. We find that anthropogenic climate change increased the likelihood of the 2015-2017 rainfall deficit by a factor of five to six. The probability of such an event will increase from 0.7 to 25% by the year 2100 under an intermediate-emission scenario (Shared Socioeconomic Pathway 2-4.5 [SSP2-4.5]) and to 80% under a high-emission scenario (SSP5-8.5). These results highlight the strong sensitivity of the drought risk in SSA to future anthropogenic emissions.

Eastern North America contains densely populated, highly developed areas, making winter storms with strong winds and high snowfall among the costliest storm types. For this reason, it is important to determine how the frequency of high-impact winter storms, specifically, those combining significant snowfall and winds, will change in this region under increasing greenhouse gas concentrations. This study uses a high-resolution coupled global climate model to simulate the changes in extreme winter conditions from the present climate to a future scenario with doubled CO2 concentrations (2XC). In particular, this study focuses on changes in high-snowfall, extreme-wind (HSEW) events, which are defined as the occurrence of 2-day snowfall and high winds exceeding thresholds based on extreme values from the control simulation, where greenhouse gas concentrations remain fixed. Mean snowfall consistently decreases across the entire region, but extreme snowfall shows a more inconsistent pattern, with some areas experiencing increases in the frequency of extreme-snowfall events. Extreme-wind events show relatively small changes in frequency with 2XC, with the exception of high-elevation areas where there are large decreases in frequency. As a result of combined changes in wind and snowfall, HSEW events decrease in frequency in the 2XC simulation for much of eastern North America. Changes in the number of HSEW events in the 2XC environment are driven mainly by changes in the frequency of extreme-snowfall events, with most of the region experiencing decreases in event frequency, except for certain inland areas at higher latitudes.

Floods in the Mississippi basin can have large negative societal, natural, and economic impacts. Understanding the drivers of floods, now and in the future, is relevant for risk management and infrastructure-planning purposes. We investigate the drivers of 100-yr-return lower Mississippi River floods using a global coupled climate model with an integrated surface water module. The model provides 3400 years of physically consistent data from a static climate, in contrast to available observational data (relatively short records, incomplete land surface data, transient climate). In the months preceding the model’s 100-yr floods, as indicated by extreme monthly discharge, above-average rain and snowfall lead to moist subsurface conditions and the buildup of snowpack, making the river system prone to these major flooding events. The meltwater from snowpack in the northern Missouri and upper Mississippi catchments primes the river system, sensitizing it to subsequent above-average precipitation in the Ohio and Tennessee catchments. An ensemble of transient forcing experiments is used to investigate the impacts of past and projected anthropogenic climate change on extreme floods. There is no statistically significant projected trend in the occurrence of 100-yr floods in the model ensemble, despite significant increases in extreme precipitation, significant decreases in extreme snowmelt, and significant decreases in less extreme floods. The results emphasize the importance of considering the fully coupled land–atmosphere system for extreme floods. This initial analysis provides avenues for further investigation, including comparison to characteristics of less extreme floods, the sensitivity to model configuration, the role of human water management, and implications for future flood-risk management.

Unprecedented high-intensity flooding induced by extreme precipitation was reported over Chennai in India during November–December of 2015, which led to extensive damage to human life and property. It is of utmost importance to determine the odds of occurrence of such extreme floods in the future, and the related climate phenomena, for planning and mitigation purposes. Here, a suite of simulations from GFDL high-resolution coupled climate models are used to investigate the odds of occurrence of extreme floods induced by extreme precipitation over Chennai and the role of radiative forcing and/or large-scale SST forcing in enhancing the probability of such events in the future. The climate of twentieth-century experiments with large ensembles suggest that the radiative forcing may not enhance the probability of extreme floods over Chennai. Doubling of CO2 experiments also fails to show evidence for an increase of such events in a global warming scenario. Further, this study explores the role of SST forcing from the Indian and Pacific Oceans on the odds of occurrence of Chennai-like floods. Neither El Niño nor La Niña enhances the probability of extreme floods over Chennai. However, a warm Bay of Bengal tends to increase the odds of occurrence of extreme Chennai-like floods. In order to trigger a Chennai like-flood, a conducive weather event, such as a tropical depression over the Bay of Bengal with strong transport of moisture from a moist atmosphere over the warm Bay, is necessary for the intense precipitation.

The Mediterranean region, located in the transition zone between the dry subtropical and wet European mid-latitude climate, is very sensitive to changes in the global mean climate state. Projecting future changes of the Mediterranean hydroclimate under global warming therefore requires dynamic climate models to reproduce the main mechanisms controlling regional hydroclimate with sufficiently high resolution to realistically simulate climate extremes. To assess future winter precipitation changes in the Mediterranean region we use the Geophysical Fluid Dynamics Laboratory high-resolution general circulation model for control simulations with pre-industrial greenhouse gas and aerosol concentrations which are compared to future scenario simulations. Here we show that the coupled model is able to reliably simulate the large-scale winter circulation, including the North Atlantic Oscillation and Eastern Atlantic patterns of variability, and its associated impacts on the mean Mediterranean hydroclimate. The model also realistically reproduces the regional features of daily heavy rainfall, which are absent in lower-resolution simulations. A five-member future projection ensemble, which assumes comparatively high greenhouse gas emissions (RCP8.5) until 2100, indicates a strong winter decline in Mediterranean precipitation for the coming decades. Consistent with dynamical and thermodynamical consequences of a warming atmosphere, derived changes feature a distinct bipolar behavior, i.e. wetting in the north—and drying in the south. Changes are most pronounced over the northwest African coast, where the projected winter precipitation decline reaches 40% of present values. Despite a decrease in mean precipitation, heavy rainfall indices show drastic increases across most of the Mediterranean, except the North African coast, which is under the strong influence of the cold Canary Current.

A stationary low pressure system and elevated levels of precipitable water provided a nearly continuous source of precipitation over Louisiana, United States (US), starting around 10 August 2016. Precipitation was heaviest in the region broadly encompassing the city of Baton Rouge, with a 3-day maximum found at a station in Livingston, LA (east of Baton Rouge), from 12 to 14 August 2016 (648.3 mm, 25.5 inches). The intense precipitation was followed by inland flash flooding and river flooding and in subsequent days produced additional backwater flooding. On 16 August, Louisiana officials reported that 30 000 people had been rescued, nearly 10 600 people had slept in shelters on the night of 14 August and at least 60 600 homes had been impacted to varying degrees. As of 17 August, the floods were reported to have killed at least 13 people. As the disaster was unfolding, the Red Cross called the flooding the worst natural disaster in the US since Super Storm Sandy made landfall in New Jersey on 24 October 2012. Before the floodwaters had receded, the media began questioning whether this extreme event was caused by anthropogenic climate change. To provide the necessary analysis to understand the potential role of anthropogenic climate change, a rapid attribution analysis was launched in real time using the best readily available observational data and high-resolution global climate model simulations. The objective of this study is to show the possibility of performing rapid attribution studies when both observational and model data and analysis methods are readily available upon the start. It is the authors' aspiration that the results be used to guide further studies of the devastating precipitation and flooding event. Here, we present a first estimate of how anthropogenic climate change has affected the likelihood of a comparable extreme precipitation event in the central US Gulf Coast. While the flooding event of interest triggering this study occurred in south Louisiana, for the purposes of our analysis, we have defined an extreme precipitation event by taking the spatial maximum of annual 3-day inland maximum precipitation over the region of 29–31° N, 85–95° W, which we refer to as the central US Gulf Coast. Using observational data, we find that the observed local return time of the 12–14 August precipitation event in 2016 is about 550 years (95 % confidence interval (CI): 450–1450). The probability for an event like this to happen anywhere in the region is presently 1 in 30 years (CI 11–110). We estimate that these probabilities and the intensity of extreme precipitation events of this return time have increased since 1900. A central US Gulf Coast extreme precipitation event has effectively become more likely in 2016 than it was in 1900. The global climate models tell a similar story; in the most accurate analyses, the regional probability of 3-day extreme precipitation increases by more than a factor of 1.4 due to anthropogenic climate change. The magnitude of the shift in probabilities is greater in the 25 km (higher-resolution) climate model than in the 50 km model. The evidence for a relation to El Niño half a year earlier is equivocal, with some analyses showing a positive connection and others none.

The Angola low is a summertime low pressure system that affects the convergence of low-level moisture fluxes into southern Africa. Interannual variations of the Angola low reduce the seasonal prediction skills for this region that arise from coupled atmosphere–ocean variability. Despite its importance, the interannual dynamics of the Angola low, and its relationship with El Niño–Southern Oscillation (ENSO) and other coupled modes of variability, are still poorly understood, mostly because of the scarcity of atmospheric data and short-term duration of atmospheric reanalyses in the region. To bypass this issue, we use a long-term (3500 year) run from a 50-km-resolution global coupled model capable of simulating the summertime southern African large-scale circulation and teleconnections. We find that the meridional displacement and strength of the Angola low are moderately modulated by local sea surface temperature anomalies, especially those in proximity of the southeastern African coast, and to a lesser extent by ENSO and the subtropical Indian Ocean dipole. Comparison of the coupled run with a 1000-yr run driven by climatological sea surface temperatures reveals that the interannual excursions of the Angola low are in both cases associated with geopotential height anomalies over the southern Atlantic and Indian Ocean related to extratropical atmospheric variability. Midlatitude atmospheric variability explains almost 60% of the variance of the Angola low variability in the uncoupled run, but only 20% in the coupled run. Therefore, while the Angola low appears to be intrinsically controlled by atmospheric extratropical variability, the interference of the atmospheric response forced by sea surface temperature anomalies weakens this influence.

Precipitation extremes have a widespread impact on societies and ecosystems; it is therefore important to understand current and future patterns of extreme precipitation. Here, a set of new global coupled climate models with varying atmospheric resolution has been used to investigate the ability of these models to reproduce observed patterns of precipitation extremes and to investigate changes in these extremes in response to increased atmospheric CO2 concentrations. The atmospheric resolution was increased from 2° × 2° grid cells (typical resolution in the CMIP5 archive) to 0.25° × 0.25° (tropical cyclone permitting). Analysis has been confined to the contiguous United States (CONUS). It is shown that, for these models, integrating at higher atmospheric resolution improves all aspects of simulated extreme precipitation: spatial patterns, intensities, and seasonal timing. In response to 2 × CO2 concentrations, all models show a mean intensification of precipitation rates during extreme events of approximately 3%–4% K−1. However, projected regional patterns of changes in extremes are dependent on model resolution. For example, the highest-resolution models show increased precipitation rates during extreme events in the hurricane season in the U.S. Southeast; this increase is not found in the low-resolution model. These results emphasize that, for the study of extreme precipitation there is a minimum model resolution that is needed to capture the weather phenomena generating the extremes. Finally, the observed record and historical model experiments were used to investigate changes in the recent past. In part because of large intrinsic variability, no evidence was found for changes in extreme precipitation attributable to climate change in the available observed record.

Western US snowpack—snow that accumulates on the ground in the mountains—plays a critical role in regional hydroclimate and water supply, with 80% of snowmelt runoff being used for agriculture. While climate projections provide estimates of snowpack loss by the end of the century and weather forecasts provide predictions of weather conditions out to 2 weeks, less progress has been made for snow predictions at seasonal timescales (months to 2 years), crucial for regional agricultural decisions (e.g., plant choice and quantity). Seasonal predictions with climate models first took the form of El Niño predictions 3 decades ago, with hydroclimate predictions emerging more recently. While the field has been focused on single-season predictions (3 months or less), we are now poised to advance our predictions beyond this timeframe. Utilizing observations, climate indices, and a suite of global climate models, we demonstrate the feasibility of seasonal snowpack predictions and quantify the limits of predictive skill 8 months in advance. This physically based dynamic system outperforms observation-based statistical predictions made on July 1 for March snowpack everywhere except the southern Sierra Nevada, a region where prediction skill is nonexistent for every predictor presently tested. Additionally, in the absence of externally forced negative trends in snowpack, narrow maritime mountain ranges with high hydroclimate variability pose a challenge for seasonal prediction in our present system; natural snowpack variability may inherently be unpredictable at this timescale. This work highlights present prediction system successes and gives cause for optimism for developing seasonal predictions for societal needs.

A newly developed modelling approach reveals how future global climate change might severely dampen economic growth in poorer countries, while increasing the variability of growth in both poorer and richer countries.

The seasonal predictability of extratropical storm tracks in the Geophysical Fluid Dynamics Laboratory’s (GFDL)’s high-resolution climate model has been investigated using an average predictability time analysis. The leading predictable components of extratropical storm tracks are the ENSO-related spatial patterns for both boreal winter and summer, and the second predictable components are mostly due to changes in external radiative forcing and multidecadal oceanic variability. These two predictable components for both seasons show significant correlation skill for all leads from 0 to 9 months, while the skill of predicting the boreal winter storm track is consistently higher than that of the austral winter. The predictable components of extratropical storm tracks are dynamically consistent with the predictable components of the upper troposphere jet flow for both seasons. Over the region with strong storm-track signals in North America, the model is able to predict the changes in statistics of extremes connected to storm-track changes (e.g., extreme low and high sea level pressure and extreme 2-m air temperature) in response to different ENSO phases. These results point toward the possibility of providing skillful seasonal predictions of the statistics of extratropical extremes over land using high-resolution coupled models.

This study demonstrates skillful seasonal prediction of 2-m air temperature and precipitation over land in a new high-resolution climate model developed by the Geophysical Fluid Dynamics Laboratory and explores the possible sources of the skill. The authors employ a statistical optimization approach to identify the most predictable components of seasonal mean temperature and precipitation over land and demonstrate the predictive skill of these components. First, the improved skill of the high-resolution model over the previous lower-resolution model in seasonal prediction of the Niño-3.4 index and other aspects of interest is shown. Then, the skill of temperature and precipitation in the high-resolution model for boreal winter and summer is measured, and the sources of the skill are diagnosed. Last, predictions are reconstructed using a few of the most predictable components to yield more skillful predictions than the raw model predictions. Over three decades of hindcasts, the two most predictable components of temperature are characterized by a component that is likely due to changes in external radiative forcing in boreal winter and summer and an ENSO-related pattern in boreal winter. The most predictable components of precipitation in both seasons are very likely ENSO-related. These components of temperature and precipitation can be predicted with significant correlation skill at least 9 months in advance. The reconstructed predictions using only the first few predictable components from the model show considerably better skill relative to observations than raw model predictions. This study shows that the use of refined statistical analysis and a high-resolution dynamical model leads to significant skill in seasonal predictions of 2-m air temperature and precipitation over land.

Recent developments in global dynamical climate prediction systems have allowed for skillful predictions of climate variables relevant to living marine resources (LMRs) at a scale useful to understanding and managing LMRs. Such predictions present opportunities for improved LMR management and industry operations, as well as new research avenues in fisheries science. LMRs respond to climate variability via changes in physiology and behavior. For species and systems where climate-fisheries links are well established, forecasted LMR responses can lead to anticipatory and more effective decisions, benefitting both managers and stakeholders. Here, we provide an overview of climate prediction systems and advances in seasonal to decadal prediction of marine-resource relevant environmental variables. We then describe a range of climate-sensitive LMR decisions that can be taken at lead-times of months to decades, before highlighting a range of pioneering case studies using climate predictions to inform LMR decisions. The success of these case studies suggests that many additional applications are possible. Progress, however, is limited by observational and modeling challenges. Priority developments include strengthening of the mechanistic linkages between climate and marine resource responses, development of LMR models able to explicitly represent such responses, integration of climate driven LMR dynamics in the multi-driver context within which marine resources exist, and improved prediction of ecosystem-relevant variables at the fine regional scales at which most marine resource decisions are made. While there are fundamental limits to predictability, continued advances in these areas have considerable potential to make LMR managers and industry decision more resilient to climate variability and help sustain valuable resources. Concerted dialog between scientists, LMR managers and industry is essential to realizing this potential.

Future changes in the North American monsoon, a circulation system that brings abundant summer rains to vast areas of the North American Southwest, could have significant consequences for regional water resources. How this monsoon will change with increasing greenhouse gases, however, remains unclear, not least because coarse horizontal resolution and systematic sea-surface temperature biases limit the reliability of its numerical model simulations. Here we investigate the monsoon response to increased atmospheric carbon dioxide (CO 2) concentrations using a 50-km-resolution global climate model which features a realistic representation of the monsoon climatology and its synoptic-scale variability. It is found that the monsoon response to CO 2 doubling is sensitive to sea-surface temperature biases. When minimizing these biases, the model projects a robust reduction in monsoonal precipitation over the southwestern United States, contrasting with previous multi-model assessments. Most of this precipitation decline can be attributed to increased atmospheric stability, and hence weakened convection, caused by uniform sea-surface warming. These results suggest improved adaptation measures, particularly water resource planning, will be required to cope with projected reductions in monsoon rainfall in the American Southwest.

Tropical cyclones (TCs) are a hazard to life and property and a prominent element of the global climate system; therefore, understanding and predicting TC location, intensity, and frequency is of both societal and scientific significance. Methodologies exist to predict basinwide, seasonally aggregated TC activity months, seasons, and even years in advance. It is shown that a newly developed high-resolution global climate model can produce skillful forecasts of seasonal TC activity on spatial scales finer than basinwide, from months and seasons in advance of the TC season. The climate model used here is targeted at predicting regional climate and the statistics of weather extremes on seasonal to decadal time scales, and comprises high-resolution (50 km × 50 km) atmosphere and land components as well as more moderate-resolution (~100 km) sea ice and ocean components. The simulation of TC climatology and interannual variations in this climate model is substantially improved by correcting systematic ocean biases through “flux adjustment.” A suite of 12-month duration retrospective forecasts is performed over the 1981–2012 period, after initializing the climate model to observationally constrained conditions at the start of each forecast period, using both the standard and flux-adjusted versions of the model. The standard and flux-adjusted forecasts exhibit equivalent skill at predicting Northern Hemisphere TC season sea surface temperature, but the flux-adjusted model exhibits substantially improved basinwide and regional TC activity forecasts, highlighting the role of systematic biases in limiting the quality of TC forecasts. These results suggest that dynamical forecasts of seasonally aggregated regional TC activity months in advance are feasible.

Glaciers in the Karakoram mountains have been stable in mass, whereas in nearby regions, mass loss has prevailed. Climate model simulations reveal a unique seasonal cycle in Karakoram snowfall that contributes to this pattern.

Mathematical and Physical Sciences: 2nd Place (The Ohio State University Edward F. Hayes Graduate Research Forum)

This study assesses the ability of a newly developed high-resolution coupled model from the Geophysical Fluid Dynamics Laboratory to simulate the cold-season hydroclimate in the present climate and examines its response to climate change forcing. Output is assessed from a 280-yr control simulation that is based on 1990 atmospheric composition and an idealized 140-yr future simulation in which atmospheric carbon dioxide increases at 1% yr−1 until doubling in year 70 and then remains constant. When compared with a low-resolution model, the high-resolution model is found to better represent the geographic distribution of snow variables in the present climate. In response to idealized radiative forcing changes, both models produce similar global-scale responses in which global-mean temperature and total precipitation increase while snowfall decreases. Zonally, snowfall tends to decrease in the low to midlatitudes and increase in the mid- to high latitudes. At the regional scale, the high- and low-resolution models sometimes diverge in the sign of projected snowfall changes; the high-resolution model exhibits future increases in a few select high-altitude regions, notably the northwestern Himalaya region and small regions in the Andes and southwestern Yukon, Canada. Despite such local signals, there is an almost universal reduction in snowfall as a percent of total precipitation in both models. By using a simple multivariate model, temperature is shown to drive these trends by decreasing snowfall almost everywhere while precipitation increases snowfall in the high altitudes and mid- to high latitudes. Mountainous regions of snowfall increases in the high-resolution model exhibit a unique dominance of the positive contribution from precipitation over temperature.

A stationary low pressure system and elevated levels of precipitable water provided a nearly continuous source of precipitation over Louisiana, United States (U.S.) starting around 10 August, 2016. Precipitation was heaviest in the region broadly encompassing the city of Baton Rouge, with a three-day maximum found at a station in Livingston, LA (east of Baton Rouge) from 12–14 August, 2016 (648.3 mm, 25.5 inches). The intense precipitation was followed by inland flash flooding and river flooding and in subsequent days produced additional backwater flooding. On 16 August, Louisiana officials reported that 30,000 people had been rescued, nearly 10,600 people had slept in shelters on the night of 14 August, and at least 60,600 homes had been impacted to varying degrees. As of 17 August, the floods were reported to have killed at least thirteen people. As the disaster was unfolding, the Red Cross called the flooding the worst natural disaster in the U.S. since Super Storm Sandy made landfall in New Jersey on 24 October, 2012. Before the floodwaters had receded, the media began questioning whether this extreme event was caused by anthropogenic climate change. To provide the necessary analysis to understand the potential role of anthropogenic climate change, a rapid attribution analysis was launched in real-time using the best readily available observational data and high-resolution global climate model simulations. The objective of this study is to show the possibility of performing rapid attribution studies when both observational and model data, and analysis methods are readily available upon the start. It is the authors aspiration that the results be used to guide further studies of the devastating precipitation and flooding event. Here we present a first estimate of how anthropogenic climate change has affected the likelihood of a comparable extreme precipitation event in the Central U.S. Gulf Coast. While the flooding event of interest triggering this study occurred in south Louisiana, for the purposes of our analysis, we have defined an extreme precipitation event by taking the spatial maximum of annual 3-day inland maximum precipitation over the region: 29–31º N, 85–95º W, which we refer to as the Central U.S. Gulf Coast. Using observational data, we find that the observed local return time of the 12–14 August precipitation event in 2016 is about 550 years (95 % confidence interval (C.I.): 450–1450). The probability for an event like this to happen anywhere in the region is presently 1 in 30 years (C.I. 11–110). We estimate that these probabilities and the intensity of extreme precipitation events of this return time have increased since 1900. A Central U.S. Gulf Coast extreme precipitation event has effectively become more likely in 2016 than it was in 1900. The global climate models tell a similar story, with the regional probability of 3-day extreme precipitation increasing due to anthropogenic climate change by a factor of more than a factor 1.4 in the most accurate analyses. The magnitude of the shift in probabilities is greater in the 25 km (higher resolution) climate model than in the 50 km model. The evidence for a relation to El Niño half a year earlier is equivocal, with some analyses showing a positive connection and others none.

The impact of atmosphere and ocean horizontal resolution on the climatology of North American monsoon Gulf of California (GoC) moisture surges is examined in a suite of global circulation models (CM2.1, FLOR, CM2.5, CM2.6, and HiFLOR) developed at the Geophysical Fluid Dynamics Laboratory (GFDL). These models feature essentially the same physical parameterizations but differ in horizontal resolution in either the atmosphere (≃200, 50, and 25 km) or the ocean (≃1°, 0.25°, and 0.1°). Increasing horizontal atmospheric resolution from 200 to 50 km results in a drastic improvement in the model’s capability of accurately simulating surge events. The climatological near-surface flow and moisture and precipitation anomalies associated with GoC surges are overall satisfactorily simulated in all higher-resolution models. The number of surge events agrees well with reanalyses, but models tend to underestimate July–August surge-related precipitation and overestimate September surge-related rainfall in the southwestern United States. Large-scale controls supporting the development of GoC surges, such as tropical easterly waves (TEWs), tropical cyclones (TCs), and trans-Pacific Rossby wave trains (RWTs), are also well captured, although models tend to underestimate the TEW and TC magnitude and number. Near-surface GoC surge features and their large-scale forcings (TEWs, TCs, and RWTs) do not appear to be substantially affected by a finer representation of the GoC at higher ocean resolution. However, the substantial reduction of the eastern Pacific warm sea surface temperature bias through flux adjustment in the Forecast-Oriented Low Ocean Resolution (FLOR) model leads to an overall improvement of tropical–extratropical controls on GoC moisture surges and the seasonal cycle of precipitation in the southwestern United States.