skip to main content


Title: Tropical Cyclone Characteristics in the MERRA‐2 Reanalysis and AMIP Simulations
Abstract

This study evaluates the tropical cyclone (TC) activity in two high‐resolution data sets—MERRA‐2 Reanalysis (Modern‐Era Retrospective Analysis for Research and Applications, Version 2) and MERRA‐2 AMIP (Atmospheric Model Intercomparison Project). These data sets use the same atmospheric model, the Goddard Earth Observing System Model, Version 5 (GEOS‐5) during the same period. However, while MERRA‐2 AMIP is a free‐running atmospheric simulation forced only with sea surface temperature (SST), MERRA‐2 Reanalysis uses an advanced data assimilation system to include a large variety of data sets. Thus, we analyze (1) the sensitivity of TC activity to the model forcing, (2) how well the TCs in both data sets replicate observed TC characteristics, (3) the sensitivity of these results to tracking schemes and thresholds. Standard diagnostics such as the number of tropical cyclones and their intensity distribution are very similar in the AMIP model and the reanalysis. TCs in both data sets are weaker than observed, as is typical for the spatial resolution of these global models. Overall, the use of data assimilation in the MERRA‐2 Reanalysis does not lead to a significantly better TC climatology than in AMIP. Furthermore, comparison of the MERRA‐2 Reanalysis to two other reanalysis data sets shows that MERRA‐2 generates fewer, but more intense TCs, than those reanalysis products.

 
more » « less
NSF-PAR ID:
10374657
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Earth and Space Science
Volume:
8
Issue:
3
ISSN:
2333-5084
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Studies of projected changes in tropical cyclones under anthropogenic climate change, as well as their modulation by internal climate modes, make use of global climate models. To this end, tropical cyclones can be tracked in the output of higher resolution models. Using climate models to make future projections of tropical cyclones relies upon having a baseline of the characteristics of model storms under the current climate. This study focuses on two high-resolution datasets – the NASA GEOS-5 Model (Goddard Earth Observing System Model, Version 5) and the MERRA-2 Reanalysis (Modern-Era Retrospective analysis for Research and Applications, Version 2). Both of these datasets were created using exactly the same atmospheric model during the same period. However, while GEOS-5 is a free-running atmospheric model forced only with sea surface temperature, MERRA-2 is a reanalysis product, i.e. the model assimilates data from a large variety of data sources. Thus, by comparing tropical cyclones tracked in these datasets to each other and global best track datasets in the period 1980-1999, this project aims to evaluate 1) the sensitivity of this model to how it is forced and 2) how well the storms tracked in GEOS-5 and MERRA-2 replicate observed tropical cyclones’ characteristics. We used two different tracking schemes on both datasets and found no significant difference in the performance of the model and the reanalysis in simulating tropical cyclones. Standard diagnostics for tropical cyclones, such as the mean number, intensity distribution, as well as their interannual variability are very similar in the free-running model and the reanalysis. Both GEOS-5 and MERRA-2 show a bias towards weaker tropical cyclones than observed and GEOS-5 has storms that occur closer to the equator than in the observed record. Neither GEOS-5 nor MERRA-2 accurately reproduce tropical cyclone modulation by ENSO. Additionally, comparison of MERRA-2 to other reanalysis datasets shows that MERRA-2 on average generates fewer total but also more intense storms than the European Centre for Medium-Range Weather Forecasts Interim Reanalysis (ERA-Interim) and Japanese 55-Year Reanalysis (JRA-55). Further research must be performed to understand why this data assimilation is failing to provide a positive impact on the tropical cyclone simulation in this model. 
    more » « less
  2. Abstract This study compares the spread in climatological tropical cyclone (TC) precipitation across eight different reanalysis datasets: NCEP-CFSR, ERA-20C, ERA-40, ERA5, ERA-Interim, JRA-55, MERRA-2, and NOAA-20C. TC precipitation is assigned using manual tracking via a fixed 500-km radius from each TC center. The reanalyses capture similar general spatial patterns of TC precipitation and TC precipitation fraction, defined as the fraction of annual precipitation assigned to TCs, and the spread in TC precipitation is larger than the spread in total precipitation across reanalyses. The spread in TC precipitation relative to the inter-reanalysis mean TC precipitation, or relative spread, is larger in the east Pacific than in the west Pacific. Partitioned by reanalysis intensity, the largest relative spread across reanalyses in TC precipitation is from high-intensity TCs. In comparison with satellite observations, reanalyses show lower climatological mean annual TC precipitation over most areas. A comparison of area-averaged precipitation rate in TCs composited over reanalysis intensity shows the spread across reanalyses is larger for higher intensity TCs. Testing the sensitivity of TC precipitation assignment to tracking method shows that climatological mean annual TC precipitation is systematically larger when assigned via manual tracking versus objective tracking. However, this tendency is minimized when TC precipitation is normalized by TC density. Overall, TC precipitation in reanalyses is affected by not only horizontal output resolution or any TC preprocessing, but also data assimilation and parameterization schemes. The results indicate that improvements in the representation of TCs and their precipitation in reanalyses are needed to improve overall precipitation. 
    more » « less
  3. Abstract

    Tropical weather phenomena—including tropical cyclones (TCs) and equatorial waves—are influenced by planetary‐to‐convective‐scale processes; yet, existing data sets and tools can only capture a subset of those processes. This study introduces a convection‐permitting aquaplanet simulation that can be used as a laboratory to study TCs, equatorial waves, and their interactions. The simulation was produced with the Model for Prediction Across Scales‐Atmosphere (MPAS‐A) using a variable resolution mesh with convection‐permitting resolution (i.e., 3‐km cell spacing) between 10°S and 30°N. The underlying sea‐surface temperature is given by a zonally symmetric profile with a peak at 10°N, which allows for the formation of TCs. A comparison between the simulation and satellite, reanalysis, and airborne dropsonde data is presented to determine the realism of the simulated phenomena. The simulation captures a realistic TC intensity distribution, including major hurricanes, but their lifetime maximum intensities may be limited by the stronger vertical wind shear in the simulation compared to the observed tropical Pacific region. The simulation also captures convectively coupled equatorial waves, including Kelvin waves and easterly waves. Despite the idealization of the aquaplanet setup, the simulated three‐dimensional structure of both groups of waves is consistent with their observed structure as deduced from satellite and reanalysis data. Easterly waves, however, have peak rotation and meridional winds at a slightly higher altitude than in the reanalysis. Future studies may use this simulation to understand how convectively coupled equatorial waves influence the multi‐scale processes leading to tropical cyclogenesis.

     
    more » « less
  4. Abstract

    Tropical cyclones (TCs) are one of the greatest threats to coastal communities along the US Atlantic and Gulf coasts due to their extreme wind, rainfall and storm surge. Analyzing historical TC climatology and modeling TC hazards can provide valuable insight to planners and decision makers. However, detailed TC size information is typically only available from 1988 onward, preventing accurate wind, rainfall, and storm surge modeling for TCs occurring earlier in the historical record. To overcome temporally limited TC size data, we develop a database of size estimates that are based on reanalysis data and a physics‐based model. Specifically, we utilize ERA5 reanalysis data to estimate the TC outer size, and a physics‐based TC wind model to estimate the radius of maximum wind. We evaluate our TC size estimates using two high‐resolution wind data sets as well as Best Track information for a wide variety of TCs. Using the estimated size information plus the TC track and intensity, we reconstruct historical storm tides from 1950 to 2020 using a basin‐scale hydrodynamic model and show that our reconstructions agree well with observed peak storm tide and storm surge. Finally, we demonstrate that incorporating an expanded set of historical modeled storm tides beginning in 1950 can enhance our understanding of US coastal hazard. Our newly developed database of TC sizes and associated storm tides/surges can aid in understanding North Atlantic TC climatology and modeling TC wind, storm surge, and rainfall hazard along the US Atlantic and Gulf coasts.

     
    more » « less
  5. Abstract

    The internal atmospheric variability (IAV) of the net surface heat flux (NHF) in the observed 20th/21st century atmosphere is estimated as the residual after removing the sea surface temperature (SST) and externally forced atmospheric response derived from the four atmospheric models of the Atmospheric Model Intercomparison Project (AMIP) simulations under phase 5 of the Coupled Model Intercomparison Project (CMIP5). The mean NHF of four atmospheric reanalysis datasets is an estimate of the observed NHF. Although the AMIP models are forced with the same SST and external forcing, the forced responses differ significantly among AMIP models, suggestive of uncertainty in the models. Besides, the uncertainty of IAV in the reanalyses could also arise from the uncertainty in reanalyses as observations contain errors and reanalysis includes interpolation by models. It is concluded that: (a) The SST/NHF and SST/forced NHF correlations are significantly negative over most of world ocean in the AMIP models, indicating damping of the SST anomalies by the NHF. (b) The IAV of the AMIP models is not correlated with SST, while the positive IAV/SST correlations in the reanalyses suggests the role of IAV in forcing the SST variability in the extra‐tropics. (c) The standard deviation (STD) of the IAV of AMIP models is indistinguishable from that of the mean reanalysis over a majority of world ocean, and the STD of the NHF of the AMIP models is larger than that of the mean reanalysis in the subtropics and midlatitudes. (d) The IAV in the mean reanalysis plays a role in forcing the SST variability in the extra‐tropics (e.g., Atlantic Multidecadal Variability), while it may not be an important forcing in the tropical oceans (e.g., ENSO).

     
    more » « less