Deciphering Changing Probabilities of Extreme Climate Events in Climate Models and Measurements

This project was developed during the 2023-2024 session of NASA GISS' Climate Change Research Initiative. Source code for this project is available in this webpage's associated Github repository.

This plot shows a histogram of precipitation in our present climate compared to a histogram of precipitation in our potential future climate under climate change where the future climate shows an increase in both intensity and frequency.

Figure 1: Future Climate Projection for Precipitation (Myhre et al. 2019).
Compared to the present climate, precipitation in the future climate may have an increase in both frequency and intensity, where there may be more precipitation in total and a tradeoff where light precipitation may decrease while extreme precipitation may increase leading to both more precipitation and more extreme precipitation.

Abstract

Mean changes in temperature and precipitation under a changing climate are well-understood, but future changes to extreme weather, particularly precipitation, are less known. In this project, we explored the implications of warmer and wetter weather conditions induced by climate change on extratropical cyclones (ETCs) in the Northern Hemisphere. ETCs, including Bomb Cyclones (a particularly intense form of ETC), are responsible for 80% of winter precipitation in mid-latitudes (Naud et al. 2022). However, the future frequency and intensity of ETCs under a changing climate are unknown. Additionally, atmospheric rivers will become more frequent and intense under a warming climate (Zhang et al. 2024), and ETCs that are concurrent with these ARs tend to be more intense (Zhang et al. 2019). We aimed to investigate how the frequency and intensity of bombing extratropical cyclones, including those that are concurrent with an atmospheric river, are changing under a warming climate, and if these changes are accurately reflected in climate models. We developed a method for identifying extratropical cyclones in gridded datasets which we applied to MERRA-2 and ERA-5 reanalysis data. We found that within the reanalysis data, extratropical cyclone frequency is higher during PNA+ and NAO- events, and bomb cyclone frequency is higher during PNA+ events. In future research, trend analysis and application of the described method on CMIP6 simulations should be done to further understand how ETCs and Bomb Cyclones are changing.

Plain Language Summary

Climate change is causing general changes in temperature and precipitation that scientists understand well, but how extreme weather events, especially heavy rainfall, will change is less clear. This project studied how warmer and wetter conditions from climate change affect extratropical cyclones (ETCs) in the Northern Hemisphere. ETCs, including a severe type called Bomb Cyclones, produce 80% of winter rainfall in mid-latitude regions. However, it’s uncertain how their frequency and intensity will change in the future. Additionally, atmospheric rivers (ARs), which will become more frequent and stronger with climate change, often make ETCs more intense. Our research aimed to find out how often and how intensely Bomb Cyclones occur, especially when combined with atmospheric rivers, under a warming climate. We used a method to identify these cyclones in specific datasets (MERRA-2 and ERA-5 reanalysis data). We discovered that extratropical cyclones are more common during certain weather patterns (PNA+ and NAO- events), and Bomb Cyclones happen more often during PNA+ events. Future research should use this method on climate model simulations (CMIP6) to better understand how ETCs and Bomb Cyclones are changing.

Team Members

Background

Due to anthropogenic activities since the Industrial Revolution, the earth’s climate has been experiencing substantial changes. From 2011 to 2020, global surface temperature has reached 1.1 °C above pre-industrial levels (GISTEMP v.4 2024) and is projected to rise. The increase in global temperature will throw off the equilibrium of the earth's systems, creating environmental, societal, and economic impacts.

Human-induced climate change is increasing the frequency and intensity of extreme climate events in every region of the world, including those related to precipitation. Under a changing climate, weather will be warmer and thus wetter. This phenomenon is supported by the Clausius-Clapeyron equation, which describes how the saturation vapor pressure of water increases exponentially with temperature, leading to more moisture in the atmosphere and potentially heavier precipitation events (Martinkova et al. 2020). Future climate is likely to have an increase in both frequency and intensity of precipitation events (Myhre et al. 2019).

Since the 1950s, most regions of the world have experienced increases in hot extremes, 21 out of 46 global regions experienced an increase in heavy precipitation, and West and Central Africa as well as parts of Europe and Asia experienced more occurrences of drought (IPCC AR6).

Understanding the shifting probabilities of extreme climate events is incredibly important to prepare the world for the escalating occurrences of severe weather phenomena. As precipitation patterns shift, it is paramount for us to be ready for both extreme precipitation events (which can have consequences such as severe flooding, landslides, infrastructure damage, and proliferation of waterborne diseases) and droughts (which can have consequences such as water supply shortages, agricultural impacts, and increased risk of wildfires).

A High Level View of Precipitation Patterns

How Precipitation Patterns are Changing in the United States

As global temperature increases, there will be an increase in annual hottest-days and wettest-days (IPCC AR6). The frequency of extreme precipitation events is projected to increase significantly with global warming, with total precipitation from these events nearly doubling per degree of warming, primarily due to increased frequency rather than intensity. This trend, observed in both climate models and actual observations, suggests that the most intense precipitation events will become almost twice as common with each additional degree of global warming (Myhre et al. 2019).

A precipitation return value is a statistical estimate representing the magnitude of a precipitation event that is expected to occur once in a given period. For example, a 100-year return value (sometimes called the '100 year storm' or 'once in a 100 year storm') refers to an extreme rainfall event that has a 1% chance of happening in any given year. Since the 1950s, precipitation return values have generally increased in the Eastern United States and decreased in the Western United States (Wu 2015).

An animation showing precipitation 10, 20, 50, and 100 year return values plotted in the contiguous USA from 1955-2015.

Figure 2: Precipitation Return Values across the contiguous United States based on Weather Station Data over Several Decades (1955-2015) (Risser et al. 2019)

Findings from the 2023 National Climate Assessment (NCA-5):

Climate Models and Precipitation

Climate models are essential for assessing climate variability and designing infrastructure, but they may not accurately predict precipitation extremes, which are expected to increase in frequency and intensity due to climate change (Abdelmoaty et al. 2021). Climate models often struggle to accurately predict precipitation extremes due to several inherent biases and limitations. One key issue is the "drizzling bias," where models tend to simulate precipitation events that are too frequent and prolonged but with lower intensity compared to observations. This results in an overestimation of light rain events and an underestimation of heavy precipitation events. Coarser models exacerbate this problem as they rely heavily on parameterized convective precipitation, which is more drizzle-like (Chen et al. 2021).

Working with NASA's GISS ModelE2.1, a coarser resolution general circulation model (GCM) amongst the CMIP6 group (Miller et al. 2020), we found differences between the model's representation of precipitation distribution and observational data. We used IMERG v.7, a remote sensing product distributed by the NASA Global Precipitation Measurement Mission, as our observational comparison (GPM_3IMERGHHE).

probability density functions for precipitation in GISS vs IMERG in NYC 2011-2020, GISS' average and 95th percentiles are slightly higher than IMERG but IMERG's 99th percentile is much higher than GISS

Figure 3: Comparative Analysis of Precipitation Distribution: The NASA GISS Climate Model vs. IMERG v.7 Remote Sensing Data
Compared to sattelite retrievals from IMERG v.7, the GISS climate model (ModelE2.1) does not capture the most extreme of precipitation in New York City 2011-2020. While the GISS model's precipitation values had a slightly higher average and 95th percentile than IMERG for this region, the observed IMERG precipitation values had a 25% higher 99th percentile value than GISS. The IMERG data used in this plot was regridded to the coarser GISS grid for comparison at the same resolution.
The source code for the regridding, preprocessing, and plotting used in this figure is available in this project's Github repository in the 'Probability Density Functions' folder.

Extreme Precipitation Events

Having looked at extreme precipitation from a total distribution standpoint, we moved forward in the project examining individual precipition events. To start, we examined IMERG v.7 precipitation data for the New York City area from 2001 to 2022 to find precipitation events and rank them by total precipitation, precipitation rate, and precipitation events with the highest individual measurements. The source code used for this analysis is available in this project's Github repository in the 'Top Precipitation Events' folder.

After we identified 100+ top precipitation events in New York City within IMERG v.7, we looked at the events individually by identifying if the event was in a National Weather Service (NWS) or NOAA database, if the event was declared a disaster by FEMA, and if the event had substantial news coverage. During this process, we discovered that many of the events we identified were Nor'easters (including this one on December 5, 2020) and Bomb Cyclones (including this one on October 17, 2019). We then shifted our research goal to focus on Bomb Cyclones (described below) and winter precipitation.

Extratropical Cyclones (ETCs)

An extratropical cyclone is a large-scale weather system that forms outside the tropics, typically featuring a low-pressure center, cold and warm fronts, and associated weather patterns such as strong winds, rain, and sometimes snow. Extratropical cyclones are a major component of mid-latitude precipitation in the Northern Hemisphere, accounting for approximately 80% of winter precipitation in this region (Naud et al. 2022). ETCs which undergo bombogenesis, or ‘bomb cyclones’, account for particularly heavy precipitation (Seiler et al. 2016). There is a lack of consensus on how extratropical cyclones (ETCs) will change under a warming climate (Naud et al. 2022), and the majority of CMIP5 climate models struggle to accurately predict the frequency of bomb cyclone occurrence (Seiler et al. 2016). Additionally, ETCs have stronger effects when they are concurrent with Atmospheric Rivers, or ARs, (Zhang et al. 2019) and the frequency and intensity of ARs is predicted to increase under the changing climate (Zhang et al. 2024).

Screenshot of news article about Nor'easter Screenshot of news article about bomb cyclone Screenshot of news article about atmospheric river Screenshot of news article about bomb cyclone

Figure 4: Recent News Stories about Extratropical Cyclones and Bomb Cyclones impacting the U.S. East Coast. Click the images to read more.

Bomb Cyclones

Bomb Cyclones are extratropical cyclones which undergo rapid intensification due to a rapid drop in pressure at their storm center. Bomb Cyclones are generally characterized by a pressure drop of at least 24 millibars within 24 hours, but this pressure drop threshold is latitude dependent and may be less than 24 millibars in some cases (Sanders and Gyakum 1980). Bomb Cyclones can produce severe weather conditions, including heavy snowfall, strong winds, and coastal flooding, making them particularly dangerous and impactful. They often originate from the western North Atlantic in the winter months.

This plot is an animation of a bomb cyclone showing air temperature and pressure.

Figure 5: An Example of a Bomb Cyclone (source: weather.com).
This animation shows a cold air mass (high pressure system) moving southeast while a warmer air mass (low pressure system) moves northward. The steep gradient in temperature and air pressure created a dangerous bomb cyclone on this date.

Atmospheric Rivers

An atmospheric river is a narrow corridor of concentrated moisture in the atmosphere that can transport vast amounts of water vapor from tropical regions to mid-latitudes, and they are often associated with extratropical cyclones (LeGrande et al. 2024). ARs can stretch for thousands of kilometers and are often responsible for significant rainfall when they make landfall, potentially leading to flooding. Atmospheric rivers play a crucial role in the global water cycle but can also cause severe weather impacts, such as heavy precipitation and strong winds (NASA GHRC). Atmospheric Rivers (ARs) are most often associated with western coastal regions, but they also have impacts beyond those regions, including across the U.S. East Coast (O'Brien et al. 2024). Both extratropical cyclones and atmospheric rivers tend to be stronger when they are combined (Zhang et al. 2019). Approximately 80% of cyclones in North Atlantic and North Pacific basins that undergo bombogenesis (i.e. bomb cyclones) are linked to the presence of an AR; in contrast, less than half of cyclones that do not undergo bombogenesis (i.e. non-bombing cyclones) are associated with an atmospheric river (Eiras-Barca et al. 2018). There is projected to be a 50% increase in AR occurrence globally by the late 21st century (Payne et al. 2020).

This plot is an animation of atmospheric rivers and associated precipitation impacting the U.S. West Coast

Figure 6: An Atmospheric River moves Water Vapor and Precipitation from the Tropics to the U.S. West Coast (source: JPL).
This visualization uses satellite data to show the movement of water vapor and precipitation as an atmospheric river slams into California.

Research Question

With a focus on the U.S. East Coast, we aimed to investigate how the frequency and intensity of extratropical cyclones, and specifically bomb cyclones, are represented in various datasets (including reanalysis and climate models) and how their characteristics are changing under climate change. Additionally, because ETCs have stronger effects when they are concurrent with atmospheric rivers, (Zhang et al. 2019) and the frequency and intensity of ARs is predicted to increase under the changing climate (Zhang et al. 2024), we also aimed to investigate atmospheric concurrent ETCs and BCs.

Back to the table of contents

Methods

Datasets and Tools

Method for Identifying Extratropical Cyclones, including Bombogenesis and Atmospheric River Concurrency, in Gridded Datasets

Identifying Extratropical Cyclones

Our first method involved capturing extratropical cyclones identified in MERRA-2 and ERA-5 reanalysis using the MAP Climatology of Mid-latitude Storminess (MCMS) method (Booth et al. 2013) which passed through our U.S. East and West Coast region masks in order to organize ETCs by relevant region. The global datasets of ETCs identified with MCMS were provided to us by Catherine M. Naud (MERRA-2) and James Booth (ERA-5) (Booth et al. 2013). U.S. East and West Coast ETCs were captured from 2007-2021 in MERRA-2 and 1950-2023 in ERA-5.

a plot of the United States showing region masks for the west and east coast with boundaries for land and ocean

Figure 7: U.S. East and West Coast Region Masks. Coastal boundaries are based on RECCAP2 (Laruelle et al. 2017) In this plot, the white grid boxes represent land, and the purple grid boxes represent coastal ocean. These mask files are available in this project's Github repository.

Identifying Bombogenesis

After regionally relevant ETCs were identified, each of the ETCs’ storm paths were searched for pressure drops associated with bombogenesis as defined by the latitude dependent equations set forth by Sanders and Gyakum 1980. Pressure drops were considered at the storm’s center and within 4 degrees orthogonal and diagonal to the storm’s center, all along the storm’s path. If the pressure drop threshold was met during a 24 hour time period at any center point along the storm's track, the ETC was marked as “bomb”; if no pressure drops that met the threshold were found, the storm was marked as “non-bomb”.

Identifying Atmospheric River Concurrency

The area around each ETC was checked for Atmospheric Rivers to determine AR-concurrency. Identified atmospheric rivers were sourced from GISS ModelE2.1 nudged to MERRA-2 horizontal winds (LeGrande et al. 2024), but in the future should be sourced from ERA5 or MERRA-2 (Guan and Waliser 2024, database of ARs available on UCLA Dataverse). A composite of an AR-concurrent ETC from Zhang et al. 2019 was used to determine a proximity box in which ARs were associated with the ETC in question. If an atmospheric river existed within the proximity box identified at any timestep of the storm, the storm was marked as 'AR-concurrent'.

a plot showing a composite of an extratropical cyclone including the ETC's storm center and an associated atmospheric river

Figure 8: Atmospheric River Concurrency - Proximity Box (Zhang et al. 2019). The probability (colors, %) of atmospheric river (AR) occurrence around the composite extratropical cyclone (EC) center (black dot). The red square shows the proximity box used to identify Atmospheric River concurrency in our project.

Identifying Storm-Specific Precipitation

We developed a method to identify storm-specific precipitation by finding the Sea Level Pressure (SLP) closed loop farthest away from the storm center which contains the storm center. To achieve this, we first loaded the SLP data and used the Haversine formula to calculate the distances between the storm center and each grid point in the dataset. We then identified the grid point closest to the storm center.

Next, we identified closed SLP loops by labeling contiguous areas of relatively low SLP values. Specifically, we considered regions where SLP values were below the median (50th percentile) of all SLP values in the dataset. This threshold allowed us to identify a manageable number of closed loops while ensuring that these loops represented significant low-pressure systems. We then calculated the centroid of each closed loop to determine their central coordinates. For each loop, we checked if the storm center was within the loop's boundaries.

Among the loops containing the storm center, we selected the one with the centroid farthest from the storm center by calculating the distances from the storm center to the centroids of these loops. We then calculated the area of this selected loop and created a mask corresponding to this loop, which we used for subsequent precipitation analysis. This mask allowed us to isolate and study the precipitation associated with the identified closed SLP loop.

the left side plot shows the closed SLP loop plotted on a map and the right side plot shows the loop surrounding a masked set of precipitation data. The data is right off the east coast of the USA, and the precipitation plot labels a Bomb Cyclone center and an associated atmospheric river

Figure 9: Identifying Storm-Specific Precipitation Example of December 26, 2010 AR-Concurrent Bomb Cyclone
The figures show Sea Level Pressure (SLP) and precipitation data from MERRA-2 reanalysis for a single timestep of an identified East Coast ETC, which was also identified as an AR-concurrent Bomb Cyclone. Precipitation is masked to only include data within the chosen SLP closed loop.

This method was not applied to the entire ETC database, but future work on the project should include this. The source code for this method is available on this project's Github repository in the "Segmenting Storm Precipitation" folder.

an animation of 5 timesteps of a storm with each timestep's respective closed SLP loop mask outlined

Figure 10: Segmented Precipitation for an Example Extratropical Cyclone (Dec. 26, 2010)
This ETC was identified by our previously described methods as an atmospheric river concurrent bomb cyclone. This animation shows the identified SLP loop mask (outlined in red) found for each timestep of a portion of the storm. For each timestep, the precipitation (from ERA-5 reanalysis) was masked with the timestep's SLP loop mask so that only precipitation data within the SLP loop was included. The animation shows this masked precipitation for five timesteps of the storm.

Output: a Database of ETCs

Extratropical Cyclones (ETCs) were identified using the above methods for the U.S. East Coast and West Coast in MERRA-2 and ERA5 reanalysis data and stored in databases used in the analysis described below. The database containing ETCs within MERRA-2 included bomb cyclone identification and atmospheric river concurrency identification based on atmospheric rivers identified in GISS ModelE2.1 nudged to MERRA-2 horizontal winds (LeGrande et al. 2024) due to availability at the time of creating the database. The MERRA-2 database included data from 2007-2021.

The database containing ETCs within ERA5 included bomb cyclone identification but did not include atmospheric river concurrency metrics due to data availability and resolution processing issues. The MERRA-2 database included data from 1950-2023. In the future, the algorithms used to create the ETC databases should be rerun using atmospheric river data in MERRA-2 and ERA5 identified by Guan and Waliser 2024, of which the global database of ARs in both MERRA-2 and ERA5 is available on UCLA Dataverse).

Back to the table of contents

Extratropical Cyclone Identification Analysis

Storm Track Density Plots

Once the database of ETCs was completed, we began by plotting the ETC storm tracks. Because of the large number of ETCs in the dataset, these storm track plots were not very readable past the monthly time scale, so we developed storm track density plots to translate this information. Below is an example of a storm track plot of which the density plots are based on.

bomb cyclone tracks plotted as line, the majority are red (AR-concurrent)

Figure 11: Bombogenesis Tracks of Bomb Cyclones in MERRA-2, East Coast USA 2007-2021
The lines represent the path of bombogenesis of each storm where the west most point is the 'start time' of bombogenesis and the east most point is the 'end time' of bombogenesis where the pressure drop threshold was met. The red lines represent Bomb Cyclones that were identified as AR-concurrent and the blue lines represent Bomb Cyclones that were identified as not being AR-concurrent.

Note that the density plots are smoothed using a Gaussian filter (sigma=1) for easier viewing. The code used to create these plots is available in this project's Github repository in the "Storm Track and Density Plots" folder.

ETCs Identified in MERRA-2 Reanalysis (2007-2021)

U.S. East Coast

storm track density plot of ETCs on the East Coast, most etcs are concentrated between 30 and 50 N latitude just off the east coast

Figure 12: Density Plot of ETC Storm Tracks identified in MERRA-2, East Coast USA 2007-2021

U.S. West Coast

storm track density plot of ETCs on the West Coast, most etcs are concentrated between 30 and 50 N latitude, primarily around 40N latitude

Figure 13: Density Plot of ETC Storm Tracks identified in MERRA-2, West Coast USA 2007-2021

ETCs Identified in ERA5 Reanalysis

U.S. East Coast

storm track density plot of ETCs on the East Coast, most etcs are concentrated between 30 and 50 N latitude

Figure 14: Density Plot of ETC Storm Tracks identified in ERA5, East Coast USA 1950-2023

storm track density plot of Bomb Cyclones on the East Coast, most BCs are concentrated between 30 and 50 N latitude just off the east coast

Figure 15: Density Plot of Bomb Cyclone Storm Tracks identified in ERA-5, East Coast USA 1950-2023

U.S. West Coast

storm track density plot of ETCs on the West Coast, most etcs are concentrated between 40 and 50 N latitude near the Seattle/Vancouver gulf

Figure 16: Density Plot of ETC Storm Tracks identified in ERA5, West Coast USA 1950-2023

Key Takeaways

Based on ERA5 reanalysis 1950-2023, there were nearly 390% more extratropical cyclones on the U.S. East Coast than on the U.S. West Coast, and it was ~6 times more likely for an East Coast ETC to undergo bombogenesis than a West Coast ETC. Based on MERRA-2 reanalysis 2007-2021, there were nearly 200% more extratropical cyclones on the U.S. East Coast than on the U.S. West Coast. It was ~3.5 times more likely for an East Coast ETC to undergo bombogenesis than a West Coast ETC. Within MERRA-2 reanalysis, 84% of Bomb Cyclones on the East Coast were concurrent with an Atmospheric River. In both ERA5 and MERRA-2 reanalysis, extratropical cyclones tend to occur in the Northern Hemisphere between 30N and 50N latitude.

Back to the table of contents

Atmospheric River-Concurrent Bomb Cyclone Case Study

Geopotential Heights during December 26 2010 and January 5 Bomb Cyclones

Figure 17: MERRA-2 sea level pressure output for the northeast U.S. during the December 26-27, 2010 and January 3-4, 2018 bomb cyclones. Figures made using NASA Giovanni.
The MERRA-2 output above depicts a drop in sea level pressure as two separate bomb cyclones impacted the Northeastern U.S. We can see the significant drop in sea level pressure, predominantly along and off the Northeast coast as the cyclones tracked through the region. Note that the MERRA-2 reanalysis depicts that the greatest drop in sea level pressure in both cyclones occurred near or around the 40N-70W latitude, longitude lines. This seems to be a consistently common location for developing cyclones to bomb out as they impact the Northeast.

Geopotential Heights during December 26 2010 Bomb Cyclone

Figure 18: 500 millibar geopotential heights reanalysis (NOAA) during the December 26, 2010 bomb cyclone
NOAA daily weather maps show a seesaw-like pattern across North America, with far lower heights (trough) associated with the bomb cyclone over the eastern half and far higher heights (ridge) over the western half. This ridge sends cold air southward from Canada/arctic regions, which then clashes with the warmer tropical air mass situated over the southern central U.S. This clashing then leads to the rapid development and intensification of a low pressure system as it tracks Northeast toward and up the east coast. (Husain et al., 2018)

Back to the table of contents

ETC Teleconnections Analysis

We studied several atmospheric teleconnections which have some influence on the Eastern United States, to see if they affect extratropical or bomb cyclone formation frequency.

North Atlantic Oscillation (NAO)

500 mb geopotential height association with positive NAO

Figure 19: 500mb geopotential height association with NAO+.



NAO association with cyclone frequency, by month
NAO association with extratropical cyclone formation, up to 28 days before formation NAO association with bomb cyclone formation, up to 28 days before formation

Figure 20: NAO association with Extratropical and Bomb Cyclone frequency. Comparisons are made by month (top) and by number of days before formation (bottom).

The source code for oscillation/teleconnection analysis is available in this project's Github repository in the 'teleconnection-analysis' folder.

North Atlantic Oscillation (NAO) tracks the pressure difference between the Azores High and the Icelandic Low in the North Atlantic Ocean; an abnormally large pressure difference indicates an NAO+, and an abnormally low pressure difference indicates an NAO-. NAO- is associated with colder East Coast winters, while NAO+ is associated with milder winters there. (El Niño & Other Oscillations)

NAO was compared with extratropical and bomb cyclone formation rates by month, and up to 28 days before storm formation. NAO- within seven days of storm formation was found to be associated with more extratropical cyclones, and NAO- from three to fourteen days before storm formation was found to be associated with more bomb cyclones. Upon storm formation, NAO- is associated with more extratropical cyclones, but not with more bomb cyclones.

Arctic Oscillation (AO)

AO association with cyclone frequency, by month
AO association with extratropical cyclone formation, up to 28 days before formation AO association with bomb cyclone formation, up to 28 days before formation

Figure 21: AO association with Extratropical and Bomb Cyclone frequency. Comparisons are made by month (top) and by number of days before formation (bottom).

Arctic Oscillation (AO) tracks the pressure difference between the Arctic and the Mid-latitudes, and the strength of the westerly winds surrounding the Arctic. An AO- indicates weak westerly winds, allowing Arctic air to escape equatorward, while an AO+ indicates strong westerly winds, leading to colder temperatures in the Arctic and milder temperatures in the mid-latitudes. (El Niño & Other Oscillations)

AO was compared with extratropical and bomb cyclone formation rates by month, and up to 28 days before storm formation. Between October and March, An AO- is associated with higher extratropical cyclone formation rates from zero to 28 days before storm formation, and is associated with higher bomb cyclone formation rates from three to 28 days before storm formation. Upon storm formation, AO- is associated with more extratropical cyclones, but not with more bomb cyclones.

Note the large degree of noise in the two middle columns. This is because AO has a relatively high standard deviation, close to 1.6. there are relatively few days with an AO amplitude of less than 0.5.

Pacific-North American Oscillation (PNA)

PNA association with cyclone frequency, by month
PNA association with extratropical cyclone formation, up to 28 days before formation PNA association with bomb cyclone formation, up to 28 days before formation

Figure 22: PNA association with Extratropical and Bomb Cyclone frequency. Comparisons are made by month (top) and by number of days before formation (bottom).

Pacific-North American Oscillation (PNA) tracks atmospheric patterns in the North Pacific and over the North American continent. Throughout most of the winter, a PNA+ is associated with a cooler East Coast and a warmer West Coast; the reverse is associated with a -PNA. (El Niño & Other Oscillations)

PNA was compared with extratropical and bomb cyclone formation rates by month, and up to 28 days before storm formation. Between October and March, a PNA+ is associated with higher extratropical and bomb cyclone formation rates from zero to 21 days before storm formation. Upon storm formation, PNA+ is associated with more extratropical and bomb cyclones.

El Niño / La Niña (ENSO)

ENSO association with cyclone frequency, by month

Figure 23: ENSO association with Extratropical and Bomb Cyclone frequency, by month.

El Niño / La Niña state tracks the sea surface temperature anomaly in the Eastern Equatorial Pacific Ocean; El Niño indicates positive anomalies in this region, while La Niña indicates negative anomalies. This oscillation has an irregular period of two to ten years. (El Niño & Other Oscillations)

Extratropical and bomb cyclone formation rates were also compared to ENSO state across all months. However, no statistically significant difference in ETC or BC formation rates was observed based on ENSO state, in extended winter (ONDJFM) or winter (DJF).

Key Takeaways

Extratropical cyclone frequency is higher during PNA+ and NAO- events, and 3 to 28 days after AO- events, and bomb cyclone frequency is higher during PNA+ events, and 3 to 28 days after AO- events. We observed little measured influence on extratropical or bomb cyclone frequency from ENSO state.

El Niño/La Niña (ENSO) Case Study

Sea level pressures during winter of 2011-2012 Precipitation during winter of 2011-2012 Sea level pressures during winter of 2016-2017 Precipitation during winter of 2016-2017 Sea level pressures during winter of 2022-2023 Precipitation during winter of 2022-2023

Figure 24: Time averaged maps of monthly SLP and satellite precip estimates during 3 La Niña winters (2011-2012), (2016-2017), and (2022-2023) respectively (top to bottom)

The La Niña ENSO state during the winter months is strongly associated with higher mean sea level pressures over the southeast U.S. This area of sinking air unsurprisingly cuts off the subtropical jet and weakens/prevents the development and cyclonic organization of many atmospheric impulses that may move through the region. We can see a lack of significant precipitation near and off of the southeast U.S. coast during the three La Niña winters, with the core of the highest precip amounts falling far to the Northeast out in the open atlantic. This implies that ETC development is delayed until the impulses move further out to sea (away from the influence of the dominant southeast ridge.) It is this delay in intensification that hampers the frequency of bomb cyclones along and off the east coast during La Niña winters.

Sea level pressures during winter of 2009-2010 Precipitation during winter of 2009-2010

Figure 25: Time averaged map of monthly SLP and satellite precip estimates during Non-La Niña winter (2009-2010)

Here, a randomly selected Non-La Niña winter shows the characteristics of a far more active subtropical jet. We can see that the mean sea level pressures are significantly lower when compared to the other three La Niña winters. This implies a far less hostile environment for the rapid development of ETCs near and off the southeast coast. It is also evident that the highest precip amounts are located far Southwest relative to the other three La Niña winter case studies, which indicates that ETCs and bomb cyclones develop far quicker, and closer to the coast during non-La Niña winters when compared to La Niña winters.

Back to the table of contents

Research Conclusions and Next Steps

Of the identified extratropical cyclones and the eastern and western coasts of the United States, there were markedly less extratropical cyclones and bomb cyclones on the west coast than the east coast. Extratropical Cyclones on the U.S. East Coast are more likely to occur and experience bombogenesis, and to be concurrent with an atmospheric river, than on the west coast.

We found that within the reanalysis data, extratropical cyclone frequency is higher during PNA+ and NAO- events, and after AO- events, and bomb cyclone frequency is higher during PNA+ events, and after AO- events. We observed little measured influence on extratropical or bomb cyclone frequency from ENSO or AMO state.

Next Steps

Back to the table of contents

Education and Community Outreach

The education component of the CCRI experience centers on creating a unit plan designed to make complex scientific concepts about extreme precipitation accessible and engaging for high school students. The unit plan leverages ArcGIS StoryMaps, allowing students to explore and present data about extreme weather events in their own neighborhoods or hometowns. This approach ensures that learning is relevant and impactful by connecting it to students' personal experiences and communities. Additionally, active involvement in community STEM engagement events fosters a love for science and technology among diverse groups of people including students, teachers, school and district administrators, and broader audience who are interested in STEM education. These events range from workshops and presentations to interactive activities, all aimed at making STEM education accessible and exciting for everyone, regardless of their background. By involving the community in these educational experiences, the goal is to inspire the next generation of scientists and engineers and to demonstrate the practical importance of scientific research in everyday life.

Educator's Unit Plan

an ArcGIS storymap about extreme rain

Example of Student Capstone Project

Community STEM Engagement Events

Educator’s takeaway from CCRI experience (from Sangmin Pak)

Back to the table of contents

More about CCRI and Alignment with NASA Objectives

The Climate Change Research Initiative (CCRI) program is an interdisciplinary, collaborative, year-long STEM engagement, and experiential learning opportunity for educators and graduate students to work directly with NASA scientists and lead research teams in a NASA research project hosted at NASA's Goddard Institute for Space Studies; Columbia University's Lamont Doherty Earth Observatory; CUNY City College of Technology in New York City; or NASA's Goddard Space Flight Center in Greenbelt, MD. The summer component of each CCRI project also includes undergraduate and high school interns.

Project Alignment with NASA Objectives:

NASA Decadal Survey for Earth and Environmental Science
Contributes to several “most important” categories
Coupling of the water and energy cycles
Extending and improving weather and air quality forecasts
Reducing climate uncertainty and informing societal response
Science 2020-2024: A Vision for Science Excellence
STRATEGY 1.3: Advance discovery in emerging fields by identifying and exploiting cross-disciplinary opportunities between traditional science disciplines
NASA Strategic Plan
Strategic Objective 1.1: Understand the Earth system and its climate
NASA’s Advancing Climate Strategy 2023
Priority 1.1: Advance climate and Earth science through novel observations, research, and modeling
OSTEM Mission and Vision Objectives
Build a diverse future STEM workforce by engaging students in authentic learning experiences with NASA’s people, content and facilities

Back to the table of contents

Project Codebase

The analysis and visualiation code used in this project is open-sourcea and available in the project's Github repository.

Back to the table of contents

References


Note: this site was last updated August 15, 2024.