Spectral Retrieval of Latent Heating Profiles from TRMM PR Data. Part II: Algorithm Improvement and Heating Estimates over Tropical Ocean Regions

2007 ◽  
Vol 46 (7) ◽  
pp. 1098-1124 ◽  
Author(s):  
Shoichi Shige ◽  
Yukari N. Takayabu ◽  
Wei-Kuo Tao ◽  
Chung-Lin Shie

Abstract The spectral latent heating (SLH) algorithm was developed for the Tropical Rainfall Measuring Mission (TRMM) precipitation radar (PR) in Part I of this study. The method uses PR information [precipitation-top height (PTH), precipitation rates at the surface and melting level, and rain type] to select heating profiles from lookup tables. Heating-profile lookup tables for the three rain types—convective, shallow stratiform, and anvil rain (deep stratiform with a melting level)—were derived from numerical simulations of tropical cloud systems from the Tropical Ocean and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE) utilizing a cloud-resolving model (CRM). To assess its global application to TRMM PR data, the universality of the lookup tables from the TOGA COARE simulations is examined in this paper. Heating profiles are reconstructed from CRM-simulated parameters (i.e., PTH, precipitation rates at the surface and melting level, and rain type) and are compared with the true CRM-simulated heating profiles, which are computed directly by the model thermodynamic equation. CRM-simulated data from the Global Atmospheric Research Program Atlantic Tropical Experiment (GATE), South China Sea Monsoon Experiment (SCSMEX), and Kwajalein Experiment (KWAJEX) are used as a consistency check. The consistency check reveals discrepancies between the SLH-reconstructed and Goddard Cumulus Ensemble (GCE)-simulated heating above the melting level in the convective region and at the melting level in the stratiform region that are attributable to the TOGA COARE table. Discrepancies in the convective region are due to differences in the vertical distribution of deep convective heating due to the relative importance of liquid and ice water processes, which varies from case to case. Discrepancies in the stratiform region are due to differences in the level separating upper-level heating and lower-level cooling. Based on these results, improvements were made to the SLH algorithm. Convective heating retrieval is now separated into upper-level heating due to ice processes and lower-level heating due to liquid water processes. In the stratiform region, the heating profile is shifted up or down by matching the melting level in the TOGA COARE lookup table with the observed one. Consistency checks indicate the revised SLH algorithm performs much better for both the convective and stratiform components than does the original one. The revised SLH algorithm was applied to PR data, and the results were compared with heating profiles derived diagnostically from SCSMEX sounding data. Key features of the vertical profiles agree well—in particular, the level of maximum heating. The revised SLH algorithm was also applied to PR data for February 1998 and February 1999. The results are compared with heating profiles derived by the convective–stratiform heating (CSH) algorithm. Because observed information on precipitation depth is used in addition to precipitation type and intensity, differences between shallow and deep convection are more distinct in the SLH algorithm in comparison with the CSH algorithm.

2009 ◽  
Vol 22 (20) ◽  
pp. 5577-5594 ◽  
Author(s):  
Shoichi Shige ◽  
Yukari N. Takayabu ◽  
Satoshi Kida ◽  
Wei-Kuo Tao ◽  
Xiping Zeng ◽  
...  

Abstract The spectral latent heating (SLH) algorithm was developed to estimate latent heating profiles for the Tropical Rainfall Measuring Mission Precipitation Radar (TRMM PR). The method uses TRMM PR information (precipitation-top height, precipitation rates at the surface and melting level, and rain type) to select heating profiles from lookup tables (LUTs). LUTs for the three rain types—convective, shallow stratiform, and anvil rain (deep stratiform with a melting level)—were derived from numerical simulations of tropical cloud systems from the Tropical Ocean and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE) using a cloud-resolving model (CRM). The two-dimensional (2D) CRM was used in previous studies. The availability of exponentially increasing computer capabilities has resulted in three-dimensional (3D) CRM simulations for multiday periods becoming increasingly prevalent. In this study, LUTs from the 2D and 3D simulations are compared. Using the LUTs from 3D simulations results in less agreement between the SLH-retrieved heating and sounding-based heating for the South China Sea Monsoon Experiment (SCSMEX). The level of SLH-estimated maximum heating is lower than that of the sounding-derived maximum heating. This is explained by the fact that using the 3D LUTs results in stronger convective heating and weaker stratiform heating above the melting level than is the case if using the 2D LUTs. More condensate is generated in and carried from the convective region in the 3D model than in the 2D model, and less condensate is produced by the stratiform region’s own upward motion.


Author(s):  
Yuya Hamaguchi ◽  
Yukari N. Takayabu

AbstractIn this study, the statistical relationship between tropical upper-tropospheric troughs (TUTTs) and the initiation of summertime tropical-depression type disturbances (TDDs) over the western and central North Pacific is investigated. By applying a spatiotemporal filter to the 34-year record of brightness temperature and using JRA-55 reanalysis products, TDD-event initiations are detected and classified as trough-related (TR) or non-trough-related (non-TR). The conventional understanding is that TDDs originate primarily in the lower-troposphere; our results refine this view by revealing that approximately 30% of TDDs in the 10°N-20°N latitude ranges are generated under the influence of TUTTs. Lead-lag composite analysis of both TR- and non-TR-TDDs clarifies that TR-TDDs occur under relatively dry and less convergent large-scale conditions in the lower-troposphere. This result suggests that TR-TDDs can form in a relatively unfavorable low-level environment. The three-dimensional structure of the wave activity flux reveals southward and downward propagation of wave energy in the upper troposphere that converges at the mid-troposphere around the region where TR-TDDs occur, suggesting the existence of extratropical forcing. Further, the role of dynamic forcing associated with the TUTT on the TR-TDD-initiation is analyzed using the quasi-geostrophic omega equation. The result reveals that moistening in the mid-to-upper troposphere takes place in association with the sustained dynamical ascent at the southeast side of the TUTT, which precedes the occurrence of deep convective heating. Along with a higher convective available potential energy due to the destabilizing effect of TUTTs, the moistening in the mid-to-upper troposphere also helps to prepare the environment favorable to TDDs initiation.


2008 ◽  
Vol 47 (7) ◽  
pp. 1929-1939 ◽  
Author(s):  
Carlton W. Ulbrich ◽  
David Atlas

Abstract Raindrop size distributions (DSDs) for tropical convective storms are used to examine the relationships between the parameters of a gamma DSD, with special emphasis on their variation with the stage of the storm. Such a distinction has rarely been made before. Several storms from a variety of tropical locations are divided into storm stages according to the temporal dependence of their reflectivity factor Z, rainfall rate R, and median volume diameter D0. In most cases it is found that the DSD parameter D0 is approximately constant in time during the convective, or C, stage, which leads to a Z–R relation of the form Z = AR, that is, a linear relationship between Z and R. This finding implies the existence of equilibrium DSDs during the C stage. The convective stage is sometimes marked by pulsations in draft strength so that D0, R, and Z and associated values of the shape parameter μ decrease in a quasi-transition stage before increasing once more. Theoretical relations between the differential reflectivity ZDR and the ratio Z/R as functions of the DSD parameter μ are derived by assuming a gamma DSD and an accurate raindrop fall speed law. It is found that data derived from disdrometer observations lie along a μ = 5 isopleth for tropical continental C stages (Puerto Rico and Brazil) and along a μ = 12 isopleth for tropical maritime C stages [Tropical Ocean and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE)]. Small values of μ that occur in the weak updraft intervals do not impact the rainfall measurements because they correspond to relatively small R. The latter features imply that the measurement of rainfall for the convective stages can be performed with standard polarimetry involving only two measurables rather than three, provided knowledge of μ is available a priori. A new rain parameter diagram is presented in which isopleths of the generalized number concentration and D0 are superimposed on the Z–R plot. It is proposed that it is possible to estimate D0 from climatological and observable storm structural features, which, with Z, provide estimates of R. Such an approach is necessary for use with conventional radars until polarimetric radars are more widely available.


2009 ◽  
Vol 66 (1) ◽  
pp. 3-21 ◽  
Author(s):  
Xiaowen Li ◽  
Wei-Kuo Tao ◽  
Alexander P. Khain ◽  
Joanne Simpson ◽  
Daniel E. Johnson

Abstract A two-dimensional cloud-resolving model is used to study the sensitivities of two microphysical schemes, a bulk scheme and an explicit spectral bin scheme, in simulating a midlatitude summertime squall line [Preliminary Regional Experiment for Storm-Scale Operational and Research Meteorology (PRE-STORM), 10–11 June 1985]. In this first part of a two-part paper, the developing and mature stages of simulated storms are compared in detail. Some variables observed during the field campaign are also presented for validation. It is found that both schemes agree well with each other, and also with published observations and retrievals, in terms of storm structures and evolution, average storm flow patterns, pressure and temperature perturbations, and total heating profiles. The bin scheme is able to produce a much more extensive and homogeneous stratiform region, which compares better with observations. However, instantaneous fields and high temporal resolution analyses show distinct characteristics in the two simulations. During the mature stage, the bulk simulation produces a multicell storm with convective cells embedded in its stratiform region. Its leading convection also shows a distinct life cycle (strong evolution). In contrast, the bin simulation produces a unicell storm with little temporal variation in its leading cell regeneration (weak evolution). More detailed, high-resolution observations are needed to validate and, perhaps, generalize these model results. Interactions between the cloud microphysics and storm dynamics that produce the sensitivities described here are discussed in detail in Part II of this paper.


2007 ◽  
Vol 135 (8) ◽  
pp. 2854-2868 ◽  
Author(s):  
Changhai Liu ◽  
Mitchell W. Moncrieff

Abstract This paper investigates the effects of cloud microphysics parameterizations on simulations of warm-season precipitation at convection-permitting grid spacing. The objective is to assess the sensitivity of summertime convection predictions to the bulk microphysics parameterizations (BMPs) at fine-grid spacings applicable to the next generation of operational numerical weather prediction models. Four microphysical parameterization schemes are compared: simple ice (Dudhia), four-class mixed phase (Reisner et al.), Goddard five-class mixed phase (Tao and Simpson), and five-class mixed phase with graupel (Reisner et al.). The experimentation involves a 7-day episode (3–9 July 2003) of U.S. midsummer convection under moderate large-scale forcing. Overall, the precipitation coherency manifested as eastward-moving organized convection in the lee of the Rockies is insensitive to the choice of the microphysics schemes, and the latent heating profiles are also largely comparable among the BMPs. The upper-level condensate and cloudiness, upper-level radiative cooling/heating, and rainfall spectrum are the most sensitive, whereas the domain-mean rainfall rate and areal coverage display moderate sensitivity. Overall, the three mixed-phase schemes outperform the simple ice scheme, but a general conclusion about the degree of sophistication in the microphysics treatment and the performance is not achievable.


2014 ◽  
Vol 142 (4) ◽  
pp. 1385-1402 ◽  
Author(s):  
Nick Guy ◽  
David P. Jorgensen

Abstract This study presents characteristics of convective systems observed during the Dynamics of the Madden–Julian oscillation (DYNAMO) experiment by the instrumented NOAA WP-3D aircraft. Nine separate missions, with a focus on observing mesoscale convective systems (MCSs), were executed to obtain data in the active and inactive phase of a Madden–Julian oscillation (MJO) in the Indian Ocean. Doppler radar and in situ thermodynamic data are used to contrast the convective system characteristics during the evolution of the MJO. Isolated convection was prominent during the inactive phases of the MJO, with deepening convection during the onset of the MJO. During the MJO peak, convection and stratiform precipitation became more widespread. A larger population of deep convective elements led to a larger area of stratiform precipitation. As the MJO decayed, convective system top heights increased, though the number of convective systems decreased, eventually transitioning back to isolated convection. A distinct shift of echo top heights and contoured frequency-by-altitude diagram distributions of radar reflectivity and vertical wind speed indicated that some mesoscale characteristics were coupled to the MJO phase. Convective characteristics in the climatological initiation region (Indian Ocean) were also apparent. Comparison to results from the Tropical Ocean and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE) in the western Pacific indicated that DYNAMO MCSs were linearly organized more parallel to the low-level shear and without strong cold pools than in TOGA COARE. Three-dimensional MCS airflow also showed a different dynamical structure, with a lack of the descending rear inflow present in shear perpendicularly organized TOGA COARE MCSs. Weaker, but deeper updrafts were observed in DYNAMO.


2020 ◽  
Vol 77 (10) ◽  
pp. 3461-3477
Author(s):  
Rebecca D. Adams-Selin

AbstractThe sensitivity of low-frequency gravity waves generated during the development and mature stages of an MCS to variations in the characteristics of the rimed ice parameterization were tested through idealized numerical simulations over a range of environment shears and instabilities. Latent cooling in the simulations with less dense, graupel-like rimed ice was more concentrated aloft near the melting level, while cooling in simulations with denser, hail-like rimed ice extended from the melting level to the surface. However, the cooling profiles still had significant internal variability across different environments and over each simulation’s duration. Initial wave production during the MCS developing stage was fairly similar in the hail and graupel simulations. During the mature stages, graupel simulations showed stronger perturbations in CAPE due to the cooling and associated wave vertical motion being farther aloft; hail simulations showed stronger perturbations in LFC due to cooling and wave vertical motion being concentrated at lower levels. The differences in the cooling profiles were not uniform enough to produce consistently different higher-order wave modes. However, the initiation of discrete cells ahead of the convective line was found to be highly sensitive to the nature of the prior destabilizing wave. Individual events of discrete propagation were suppressed in some of the graupel simulations due to the higher location of both peak cooling and vertical wave motion. Such results underscore the need to fully characterize MCS microphysical heating profiles and their low-frequency gravity waves to understand their structure and development.


2000 ◽  
Vol 39 (5) ◽  
pp. 680-685 ◽  
Author(s):  
Qihang Li ◽  
Ralph Ferraro ◽  
Norman Grody

Abstract Until recently, monthly rainfall products using the National Oceanic and Atmospheric Administration National Environmental Satellite, Data, and Information Service Office of Research and Applications Special Sensor Microwave Imager (SSM/I) rainfall algorithm have been generated on a global 2.5° × 2.5° grid. The rainfall estimates are based on a subsampled set of the full-resolution SSM/I data, with a resulting spatial density of about one-third of what is possible at SSM/I’s highest spatial resolution. The reduction in the spatial resolution was introduced in 1992 as a compromise dictated by data processing capabilities. Currently, daily SSM/I data processing at full resolution has been established and is being operated in parallel with the subsampled set. Reprocessing of the entire SSM/I time series based on the full-resolution data is plausible but requires the reprocessing of over 10 yr of retrospective data. Because the Global Precipitation Climatology Project is considering the generation of a daily 1° × 1° rainfall product, it is important that the effects of using the reduced spatial resolution be reexamined. In this study, error due to using the reduced-resolution versus the full-resolution SSM/I data in the gridded products at 2.5° and 1° grid sizes is examined. The estimates are based on statistics from radar-derived rain data and from SSM/I data taken over the Tropical Ocean and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE) radar site. SSM/I data at full resolution were assumed to provide rain estimates with 12.5-km spacing. Subsampling with spacings of 25, 37.5 (which corresponds to the present situation of ⅓° latitude–longitude spatial resolution), and 50 km were considered. For the instantaneous 2.5° × 2.5° product, the error due to subsampling, expressed as a percentage of the gridbox mean, was estimated using radar-derived data and was 6%, 10%, and 15% at these successively poorer sampling densities. For monthly averaged products on a 2.5° × 2.5° grid, it was substantially lower: 3%, 4%, and 7%, respectively. Subsampling errors for monthly averages on a 1° × 1° grid were 8%, 16%, and 23%, respectively. Estimates based on SSM/I data at full resolution gave errors that were somewhat larger than those from radar-based estimates. It was concluded that a rain product of monthly averages on a 1° × 1° grid must use the full-resolution SSM/I data. More work is needed to determine how applicable these estimates are to other areas of the globe with substantially different rain statistics.


2005 ◽  
Vol 18 (8) ◽  
pp. 1190-1202 ◽  
Author(s):  
D. J. Bernie ◽  
S. J. Woolnough ◽  
J. M. Slingo ◽  
E. Guilyardi

Abstract The intraseasonal variability of SST associated with the passage of the Madden–Julian oscillation (MJO) is well documented; yet coupled model integrations generally underpredict the magnitude of this SST variability. Observations from the Improved Meteorological Instrument (IMET) mooring in the western Pacific during the intensive observing period (IOP) of the Tropical Ocean Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE) showed a large diurnal signal in SST that is modulated by the passage of the MJO. In this study, observations from the IOP of the TOGA COARE and a one-dimensional (1D) ocean mixed layer model incorporating the K-Profile Parameterization (KPP) vertical mixing scheme have been used to investigate the rectification of the intraseasonal variability of SST by the diurnal cycle and the implied impact of the absence of a representation of this process on the modeled intraseasonal variability in coupled GCMs. Analysis of the SST observations has shown that the increase of the daily mean SST by the diurnal cycle of SST accounts for about one-third of the magnitude of intraseasonal variability of SST associated with the Madden–Julian oscillation in the western Pacific warm pool. Experiments from the 1D model forced with fluxes at a range of temporal resolutions and with differing vertical resolution of the model have shown that to capture 90% of the diurnal variability of SST, and hence 95% of the intraseasonal variability of SST, requires a 3-h or better temporal resolution of the fluxes and a vertical grid with an upper-layer thickness of the order of 1 m. In addition to the impact of the representation of the diurnal cycle on the intraseasonal variability of SST, the strength of the mixing across the thermocline was found to be enhanced by the proper representation of the nighttime deep mixing in the ocean, implying a possible impact of the diurnal cycle onto the mean climate of the tropical ocean.


Sign in / Sign up

Export Citation Format

Share Document