skip to main content

Title: Predictive Statistical Representations of Observed and Simulated Rainfall Using Generalized Linear Models

This study explores the feasibility of predicting subdaily variations and the climatological spatial patterns of rain in the tropical Pacific from atmospheric profiles using a set of generalized linear models: logistic regression for rain occurrence and gamma regression for rain amount. The prediction is separated into different rain types from TRMM satellite radar observations (stratiform, deep convective, and shallow convective) and CAM5 simulations (large-scale and convective). Environmental variables from MERRA-2 and CAM5 are used as predictors for TRMM and CAM5 rainfall, respectively. The statistical models are trained using environmental fields at 0000 UTC and rainfall from 0000 to 0600 UTC during 2003. The results are used to predict 2004 rain occurrence and rate for MERRA-2/TRMM and CAM5 separately. The first EOF profile of humidity and the second EOF profile of temperature contribute most to the prediction for both statistical models in each case. The logistic regression generally performs well for all rain types, but does better in the east Pacific compared to the west Pacific. The gamma regression produces reasonable geographical rain amount distributions but rain rate probability distributions are not predicted as well, suggesting the need for a different, higher-order model to predict rain rates. The results of this study suggest that statistical models applied to TRMM radar observations and MERRA-2 environmental parameters can predict the spatial patterns and amplitudes of tropical rainfall in the time-averaged sense. Comparing the observationally trained models to models that are trained using CAM5 simulations points to possible deficiencies in the convection parameterization used in CAM5.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Journal of Climate
Page Range / eLocation ID:
p. 3409-3427
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Predicting rain from large-scale environmental variables remains a challenging problem for climate models and it is unclear how well numerical methods can predict the true characteristics of rainfall without smaller (storm) scale information. This study explores the ability of three statistical and machine learning methods to predict 3-hourly rain occurrence and intensity at 0.5° resolution over the tropical Pacific Ocean using rain observations the Global Precipitation Measurement (GPM) satellite radar and large-scale environmental profiles of temperature and moisture from the MERRA-2 reanalysis. We also separated the rain into different types (deep convective, stratiform, and shallow convective) because of their varying kinematic and thermodynamic structures that might respond to the large-scale environment in different ways. Our expectation was that the popular machine learning methods (i.e., the neural network and random forest) would outperform a standard statistical method (a generalized linear model) because of their more flexible structures, especially in predicting the highly skewed distribution of rain rates for each rain type. However, none of the methods obviously distinguish themselves from one another and each method still has issues with predicting rain too often and not fully capturing the high end of the rain rate distributions, both of which are common problems in climate models. One implication of this study is that machine learning tools must be carefully assessed and are not necessarily applicable to solving all big data problems. Another implication is that traditional climate model approaches are not sufficient to predict extreme rain events and that other avenues need to be pursued.

    more » « less
  2. Abstract

    Convective organization has a large impact on precipitation and feeds back on larger‐scale circulations in the tropics. The degree of this convective organization changes with modes of climate variability like the El Niño–Southern Oscillation (ENSO), but because organization is not represented in current climate models, a quantitative assessment of these shifts has not been possible. Here, we construct multidecade satellite climatologies of occurrence of tropical convective organization and its properties and assess changes with ENSO phase. The occurrence of organized deep convection becomes more concentrated, increasing threefold in the eastern and central Pacific during El Niño and decreasing twofold outside of these regions. Both horizontal extent of the cold cloud shield and convective depth increase in regions of positive sea surface temperature anomaly (SSTa); however, the regions of greatest convective deepening are those of large‐scale ascent, rather than those of warmest SSTa. Extent decreases with SSTa at a rate of about 20 km/K, while the SSTa dependence of depth is only about 0.2 K/K. We introduce two values to describe convective changes with ENSO more succinctly: (1) an information entropy metric to quantify the clustering of convective system occurrences and (2) a growth metric to quantify deepening relative to spreading over the system lifetime. Finally, with collocated precipitation data, we see that rainfall attributable to convective organization jumps up to 5% with warming. Rain intensity and amount increase for a given system size during El Niño, but a given rain amount may actually fall with higher intensity during La Niña.

    more » « less
  3. Abstract

    Spatial aggregation of deep convection and its possible role in larger-scale atmospheric behavior have received growing attention. Here we seek aggregation-correlated statistical properties of convective events in 5° × 5° boxes over the tropical Indian Ocean. Events are identified by box-averaged rainfall exceeding 5 mm day−1at the center of a 4-day time window, and aggregation is estimated by an index [simple convective aggregation index (SCAI)] based on contiguous cold cloud areas and their geometrical distances in infrared imagery. A physical framework using gross moist stability (GMS) helps to interpret relationships between aggregation, box-scale ascent profiles, moist static energy budgets, and time evolution both within composite events and on longer time scales. For a given precipitation rate, more-aggregated events (with fewer and larger cloud objects on average) exhibit a drier area mean, greater horizontal gradient of moisture, more bottom-heavy ascent profile, and a greater prevalence of low-altitude cloud tops, especially for lower rain rates. In the GMS budget, this bottom-heavy ascent implies net energy import into the atmospheric column during the 4-day event composite. Consistently, net energy variations filtered to reveal longer time scales do indeed exhibit more-aggregated rain events in their growth phase than in their flat and decaying phases. More-aggregated scenes also have more drying by analysis than less-aggregated scenes in MERRA-2’s assimilation budgets. This suggests that parameterized convection (lacking any organization effect) is raining out less water than nature’s real, aggregated convection in such scenes.

    more » « less
  4. Abstract The parametric hurricane rainfall model (PHRaM), firstly introduced in 2007, has been widely used to forecast and quantify tropical-cyclone-induced rainfall (TC rainfall). The PHRaM is much more computationally efficient than global climate models, but PHRaM cannot be effectively utilized in the context of climate change because it does not have any parameters to capture the increase of tropospheric water vapor under the warming world. This study develops a new model that incorporates tropospheric water vapor to the PHRaM framework, named as the PHRaM with moisture (PHRaMM). The PHRaMM is trained to best fit the TC rainfall over the western North Pacific (WNP) unlike the PHRaM trained with the TCs over the continental US. The PHRaMM reliably simulates radial profile of TC rainfall and spatial distribution of accumulated rainfall during landfall in the present climate with the better prediction skills than existing statistical and operational numerical models. Using the PHRaMM, we evaluated the impacts of TC intensity and environmental moisture increase on TC rainfall change in a future climate. An increased TC intensity causes TC rainfall to increase in the inner-core region but to decrease in the outer region, whereas an increased environmental moisture causes the TC rainfall to increase over the entire TC area. According to the both effects of increased TC intensity and environmental moisture, the PHRaMM projected that the WNP TC rainfall could increase by 4.61–8.51% in the inner-core region and by 17.96–20.91% over the entire TC area under the 2-K warming scenario. 
    more » « less
  5. Abstract

    Organized deep convective activity has been routinely monitored by satellite precipitation radar from the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Mission (GPM). Organized deep convective activity is found to increase not only with sea surface temperature (SST) above 27°C, but also with low-level wind shear. Precipitation shows a similar increasing relationship with both SST and low-level wind shear, except for the highest low-level wind shear. These observations suggest that the threshold for organized deep convection and precipitation in the tropics should consider not only SST, but also vertical wind shear. The longwave cloud radiative feedback, measured as the tropospheric longwave cloud radiative heating per amount of precipitation, is found to generally increase with stronger organized deep convective activity as SST and low-level wind shear increase. Organized deep convective activity, the longwave cloud radiative feedback, and cirrus ice cloud cover per amount of precipitation also appear to be controlled more strongly by SST than by the deviation of SST from its tropical mean. This study hints at the importance of non-thermodynamic factors such as vertical wind shear for impacting tropical convective structure, cloud properties, and associated radiative energy budget of the tropics.

    Significance Statement

    This study uses tropical satellite observations to demonstrate that vertical wind shear affects the relationship between sea surface temperature and tropical organized deep convection and precipitation. Shear also affects associated cloud properties and how clouds affect the flow of radiation in the atmosphere. Although how vertical wind shear affects convective organization has long been studied in the mesoscale community, the study attempts to apply mesoscale theory to explain the large-scale mean organization of tropical deep convection, cloud properties, and radiative feedbacks. The study also provides a quantitative observational baseline of how vertical wind shear modifies cloud radiative effects and convective organization, which can be compared to numerical simulations.

    more » « less