Engineering topics
Publications and source records attributed to Tao, Wei-Kuo.
Conflicting views on role of the SAL pre- and post-genesis (Karyampudi and Carlson 1988, Dunion and Velden 2004, Braun 2010, among others). Early dust impact studies claimed negative impacts, but had unrealistic dust distributions. (Zhang et al. 2007, 2009). More recent work with more realistic dust suggest possible positive impacts in some cases (Herbener et al. 2014).
Numerous cloud microphysical schemes designed for cloud and mesoscale models are currently in use, ranging from simple bulk to multi-moment, multi-class to explicit bin schemes. This study details the benefits of adding a 4th ice class (hail) to an already improved 3-class ice bulk microphysics scheme developed for the Goddard Cumulus Ensemble model based on Rutledge and Hobbs (1983,1984). Besides the addition and modification of several hail processes from Lin et al. (1983), further modifications were made to the 3-ice processes, including allowing greater ice super saturation and mitigating spurious evaporationsublimation in the saturation adjustment scheme, allowing graupelhail to become snow via vapor growth and hail to become graupel via riming, and the inclusion of a rain evaporation correction and vapor diffusivity factor. The improved 3-ice snowgraupel size-mapping schemes were adjusted to be more stable at higher mixing rations and to increase the aggregation effect for snow. A snow density mapping was also added. The new scheme was applied to an intense continental squall line and a weaker, loosely-organized continental case using three different hail intercepts. Peak simulated reflectivities agree well with radar for both the intense and weaker case and were better than earlier 3-ice versions when using a moderate and large intercept for hail, respectively. Simulated reflectivity distributions versus height were also improved versus radar in both cases compared to earlier 3-ice versions. The bin-based rain evaporation correction affected the squall line case more but did not change the overall agreement in reflectivity distributions.
The Goddard microphysics scheme was recently improved by adding a 4th ice class (frozen dropshail). This new 4ICE scheme was implemented and tested in the Goddard Cumulus Ensemble model (GCE) for an intense continental squall line and a moderate,less-organized continental case. Simulated peak radar reflectivity profiles were improved both in intensity and shape for both cases as were the overall reflectivity probability distributions versus observations. In this study, the new Goddard 4ICE scheme is implemented into the regional-scale NASA Unified - Weather Research and Forecasting model (NU-WRF) and tested on an intense mesoscale convective system that occurred during the Midlatitude Continental Convective Clouds Experiment (MC3E). The NU42WRF simulated radar reflectivities, rainfall intensities, and vertical and horizontal structure using the new 4ICE scheme agree as well as or significantly better with observations than when using previous versions of the Goddard 3ICE (graupel or hail) schemes. In the 4ICE scheme, the bin microphysics-based rain evaporation correction produces more erect convective cores, while modification of the unrealistic collection of ice by dry hail produces narrow and intense cores, allowing more slow-falling snow to be transported rearward. Together with a revised snow size mapping, the 4ICE scheme produces a more horizontally stratified trailing stratiform region with a broad, more coherent light rain area. In addition, the NU-WRF 4ICE simulated radar reflectivity distributions are consistent with and generally superior to those using the GCE due to the less restrictive open lateral boundaries
Cloud-top temperature (CTT) is an important parameter for convective clouds and is usually different from the 11-micrometers brightness temperature due to non-blackbody effects. This paper presents an algorithm for estimating convective CTT by using simultaneous passive [Moderate Resolution Imaging Spectroradiometer (MODIS)] and active [CloudSat 1 Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO)] measurements of clouds to correct for the non-blackbody effect. To do this, a weighting function of the MODIS 11-micrometers band is explicitly calculated by feeding cloud hydrometer profiles from CloudSat and CALIPSO retrievals and temperature and humidity profiles based on ECMWF analyses into a radiation transfer model.Among 16 837 tropical deep convective clouds observed by CloudSat in 2008, the averaged effective emission level (EEL) of the 11-mm channel is located at optical depth; approximately 0.72, with a standard deviation of 0.3. The distance between the EEL and cloud-top height determined by CloudSat is shown to be related to a parameter called cloud-top fuzziness (CTF), defined as the vertical separation between 230 and 10 dBZ of CloudSat radar reflectivity. On the basis of these findings a relationship is then developed between the CTF and the difference between MODIS 11-micrometers brightness temperature and physical CTT, the latter being the non-blackbody correction of CTT. Correction of the non-blackbody effect of CTT is applied to analyze convective cloud-top buoyancy. With this correction, about 70% of the convective cores observed by CloudSat in the height range of 6-10 km have positive buoyancy near cloud top, meaning clouds are still growing vertically, although their final fate cannot be determined by snapshot observations.
Convection is the primary transport process in the Earth's atmosphere. About two-thirds of the Earth's rainfall and severe floods derive from convection. In addition, two-thirds of the global rain falls in the tropics, while the associated latent heat release accounts for three-fourths of the total heat energy for the Earth's atmosphere. Cloud-resolving models (CRMs) have been used to improve our understanding of cloud and precipitation processes and phenomena from micro-scale to cloud-scale and mesoscale as well as their interactions with radiation and surface processes. CRMs use sophisticated and realistic representations of cloud microphysical processes and can reasonably well resolve the time evolution, structure, and life cycles of clouds and cloud systems. CRMs also allow for explicit interaction between clouds, outgoing longwave (cooling) and incoming solar (heating) radiation, and ocean and land surface processes. Observations are required to initialize CRMs and to validate their results. The Goddard Cumulus Ensemble model (GCE) has been developed and improved at NASA/Goddard Space Flight Center over the past three decades. It is amulti-dimensional non-hydrostatic CRM that can simulate clouds and cloud systems in different environments. Early improvements and testing were presented in Tao and Simpson (1993) and Tao et al. (2003a). A review on the application of the GCE to the understanding of precipitation processes can be found in Simpson and Tao (1993) and Tao (2003). In this paper, recent model improvements (microphysics, radiation and land surface processes) are described along with their impact and performance on cloud and precipitation events in different geographic locations via comparisons with observations. In addition, recent advanced applications of the GCE are presented that include understanding the physical processes responsible for diurnal variation, examining the impact of aerosols (cloud condensation nuclei or CCN and ice nuclei or IN) on precipitation processes, utilizing a satellite simulator to improve the microphysics, providing better simulations for satellite-derived latent heating retrieval, and coupling with a general circulation model to improve the representation of precipitation processes.
One of the current challenges in tropical cyclone (TC) research is how to improve our understanding of TC interannual variability and the impact of climate change on TCs. Recent advances in global modeling, visualization, and supercomputing technologies at NASA show potential for such studies. In this article, the authors discuss recent scalability improvement to the multiscale modeling framework (MMF) that makes it feasible to perform long-term TC-resolving simulations. The MMF consists of the finite-volume general circulation model (fvGCM), supplemented by a copy of the Goddard cumulus ensemble model (GCE) at each of the fvGCM grid points, giving 13,104 GCE copies. The original fvGCM implementation has a 1D data decomposition; the revised MMF implementation retains the 1D decomposition for most of the code, but uses a 2D decomposition for the massive copies of GCEs. Because the vast majority of computation time in the MMF is spent computing the GCEs, this approach can achieve excellent speedup without incurring the cost of modifying the entire code. Intelligent process mapping allows differing numbers of processes to be assigned to each domain for load balancing. The revised parallel implementation shows highly promising scalability, obtaining a nearly 80-fold speedup by increasing the number of cores from 30 to 3,335.
This study investigates the impact of snow, graupel, and hail processes on simulated squall lines over the Southern Great Plains in the United States. The Weather Research and Forecasting (WRF) model is used to simulate two squall line events in Oklahoma during May 2007, and the simulations are validated against radar and surface observations. Several microphysics schemes are tested in this study, including the WRF 5-Class Microphysics (WSM5), WRF 6-Class Microphysics (WSM6), Goddard Cumulus Ensemble (GCE) Three Ice (3-ice) with graupel, Goddard Two Ice (2-ice), and Goddard 3-ice hail schemes. Simulated surface precipitation is sensitive to the microphysics scheme when the graupel or hail categories are included. All of the 3-ice schemes overestimate the total precipitation with WSM6 having the largest bias. The 2-ice schemes, without a graupel/hail category, produce less total precipitation than the 3-ice schemes. By applying a radar-based convective/stratiform partitioning algorithm, we find that including graupel/hail processes increases the convective areal coverage, precipitation intensity, updraft, and downdraft intensities, and reduces the stratiform areal coverage and precipitation intensity. For vertical structures, simulations have higher reflectivity values distributed aloft than the observed values in both the convective and stratiform regions. Three-ice schemes produce more high reflectivity values in convective regions, while 2-ice schemes produce more high reflectivity values in stratiform regions. In addition, this study has demonstrated that the radar-based convective/stratiform partitioning algorithm can reasonably identify WRF-simulated precipitation, wind, and microphysical fields in both convective and stratiform regions.
The present generation of general circulation models (GCM) use parameterized cumulus schemes and run at hydrostatic grid resolutions. To improve the representation of cloud-scale moist processes and landeatmosphere interactions, a global, Multi-scale Modeling Framework (MMF) coupled to the Land Information System (LIS) has been developed at NASA-Goddard Space Flight Center. The MMFeLIS has three components, a finite-volume (fv) GCM (Goddard Earth Observing System Ver. 4, GEOS-4), a 2D cloud-resolving model (Goddard Cumulus Ensemble, GCE), and the LIS, representing the large-scale atmospheric circulation, cloud processes, and land surface processes, respectively. The non-hydrostatic GCE model replaces the single-column cumulus parameterization of fvGCM. The model grid is composed of an array of fvGCM gridcells each with a series of embedded GCE models. A horizontal coupling strategy, GCE4fvGCM4Coupler4LIS, offered significant computational efficiency, with the scalability and I/O capabilities of LIS permitting landeatmosphere interactions at cloud-scale. Global simulations of 2007e2008 and comparisons to observations and reanalysis products were conducted. Using two different versions of the same land surface model but the same initial conditions, divergence in regional, synoptic-scale surface pressure patterns emerged within two weeks. The sensitivity of largescale circulations to land surface model physics revealed significant functional value to using a scalable, multi-model land surface modeling system in global weather and climate prediction.
A unique microphysical structure of rainfall is observed by the surface laser optical Particle Size and Velocity (Parsivel) disdrometers on 25 April 2011 during Midlatitude Continental Convective Clouds Experiment (MC3E). According to the systematic differences in rainfall rate and bulk effective droplet radius, the sampling data can be divided into two groups; the rainfall mostly from the deep convective clouds has relatively high rainfall rate and large bulk effective droplet radius, whereas the reverse is true for the rainfall from the shallow wrm clouds. The Weather Research and Forecasting model coupled with spectral bin microphysics (WRF-SBM) successfully reproduces the two distinct modes in the observed rainfall microphysical structure. The results show that the up-to-date model can demonstrate how the cloud physics and the weather condition on the day are involved in forming the unique rainfall characteristic.
Two distinct snowfall events are observed over the region near the Great Lakes during 19-23 January 2007 under the intensive measurement campaign of the Canadian CloudSat/CALIPSO validation project (C3VP). These events are numerically investigated using the Weather Research and Forecasting model coupled with a spectral bin microphysics (WRF-SBM) scheme that allows a smooth calculation of riming process by predicting the rimed mass fraction on snow aggregates. The fundamental structures of the observed two snowfall systems are distinctly characterized by a localized intense lake-effect snowstorm in one case and a widely distributed moderate snowfall by the synoptic-scale system in another case. Furthermore, the observed microphysical structures are distinguished by differences in bulk density of solid-phase particles, which are probably linked to the presence or absence of supercooled droplets. The WRF-SBM coupled with Goddard Satellite Data Simulator Unit (G-SDSU) has successfully simulated these distinctive structures in the three-dimensional weather prediction run with a horizontal resolution of 1 km. In particular, riming on snow aggregates by supercooled droplets is considered to be of importance in reproducing the specialized microphysical structures in the case studies. Additional sensitivity tests for the lake-effect snowstorm case are conducted utilizing different planetary boundary layer (PBL) models or the same SBM but without the riming process. The PBL process has a large impact on determining the cloud microphysical structure of the lake-effect snowstorm as well as the surface precipitation pattern, whereas the riming process has little influence on the surface precipitation because of the small height of the system.
During the past decade, numerical weather and global non-hydrostatic models have started using more complex microphysical schemes originally developed for high resolution cloud resolving models (CRMs) with 1-2 km or less horizontal resolutions. These microphysical schemes affect the dynamic through the release of latent heat (buoyancy loading and pressure gradient) the radiation through the cloud coverage (vertical distribution of cloud species), and surface processes through rainfall (both amount and intensity). Recently, several major improvements of ice microphysical processes (or schemes) have been developed for cloud-resolving model (Goddard Cumulus Ensemble, GCE, model) and regional scale (Weather Research and Forecast, WRF) model. These improvements include an improved 3-ICE (cloud ice, snow and graupel) scheme (Lang et al. 2010); a 4-ICE (cloud ice, snow, graupel and hail) scheme and a spectral bin microphysics scheme and two different two-moment microphysics schemes. The performance of these schemes has been evaluated by using observational data from TRMM and other major field campaigns. In this talk, we will present the high-resolution (1 km) GeE and WRF model simulations and compared the simulated model results with observation from recent field campaigns [i.e., midlatitude continental spring season (MC3E; 2010), high latitude cold-season (C3VP, 2007; GCPEx, 2012), and tropical oceanic (TWP-ICE, 2006)].
Numerical weather prediction (NWP) simulations using the Japan Meteorological Agency NonhydrostaticModel (JMA-NHM) are conducted for three precipitation events observed by shipborne or spaceborneW-band cloud radars. Spectral bin and single-moment bulk cloud microphysics schemes are employed separatelyfor an intercomparative study. A radar product simulator that is compatible with both microphysicsschemes is developed to enable a direct comparison between simulation and observation with respect to theequivalent radar reflectivity factor Ze, Doppler velocity (DV), and path-integrated attenuation (PIA). Ingeneral, the bin model simulation shows better agreement with the observed data than the bulk modelsimulation. The correction of the terminal fall velocities of snowflakes using those of hail further improves theresult of the bin model simulation. The results indicate that there are substantial uncertainties in the masssizeand sizeterminal fall velocity relations of snowflakes or in the calculation of terminal fall velocity of snowaloft. For the bulk microphysics, the overestimation of Ze is observed as a result of a significant predominanceof snow over cloud ice due to substantial deposition growth directly to snow. The DV comparison shows thata correction for the fall velocity of hydrometeors considering a change of particle size should be introducedeven in single-moment bulk cloud microphysics.
The diurnal variation of precipitation processes in the United States (US) is well recognized but incompletely understood (Cabone et al. 2002). The diurnal cycle of precipitation has been studied using surface rainfall data, radar reflectivity data, and satellite-derived cloudiness and precipitation (Wallace 1975; Dai et al. 1999; Carbone et al. 2002; Carbone and Tuttle, 2008; Parker and Ahijevych, 2007; Matsui et al. 2010 and others). These observations indicate that the summer-time precipitation most of the North America and typically feature late-afternoon precipitation maxima. These diurnal variation of precipitation can also be generally categorized into three different types: 1) afternoon rainfall maxima due to mesoscale and local circulations over the south and east of the Mississippi and Ohio valleys, 2) nocturnal rainfall maxima from eastward-propagating mesoscale convective systems (MCSs) over the Lee side of Rocky Mountain regions and 3) afternoon rainfall maxima in the Appalachian Mountains, and then propagate eastward toward the coast. The main objective of this paper is to use a regional cloud-scale model with very high-resolution (i.e., WRF) to examine the WRF ability to simulate diurnal variation of precipitation. Specifically, the study will (1) identify the physical processes responsible to diurnal variation of precipitation, (2) examine the sensitivity of resolution (2, 6, 18, and 30 km) to model simulated diurnal variation of precipitation and (3) identify the relationships between microphysics and cumulus parameterization schemes.
Aerosols are a critical.factor in the atmospheric hydrological cycle and radiation budget. As a major agent for clouds to form and a significant attenuator of solar radiation, aerosols affect climate in several ways. Current research suggests that aerosols have a major impact on the dynamics, microphysics, and electrification properties of continental mixed-phase convective clouds. In addition, high aerosol concentrations in urban environments could affect precipitation variability by providing a significant source of cloud condensation nuclei (CCN). Such pollution . effects on precipitation potentially have enormous climatic consequences both in terms of feedbacks involving the land surface via rainfall as well as the surface energy budget and changes in latent heat input to the atmosphere. Basically, aerosol concentrations can influence cloud droplet size distributions, the warm-rain process, the cold-rain process, cloud-top heights, the depth of the mixed-phase region, and the occurrence of lightning. Recently, many cloud resolution models (CRMs) have been used to examine the role of aerosols on mixed-phase convective clouds. These modeling studies have many differences in terms of model configuration (two- or three-dimensional), domain size, grid spacing (150-3000 m), microphysics (two-moment bulk, simple or sophisticated spectral-bin), turbulence (1st or 1.5 order turbulent kinetic energy (TKE)), radiation, lateral boundary conditions (i.e., closed, radiative open or cyclic), cases (isolated convection, tropical or midlatitude squall lines) and model integration time (e.g., 2.5 to 48 hours). Among these modeling studies, the most striking difference is that cumulative precipitation can either increase or decrease in response to higher concentrations of CCN. In this presentation, we review past efforts and summarize our current understanding of the effect of aerosols on convective precipitation processes. Specifically, this paper addresses the following topics: observational evidence of the effect of aerosols on precipitation processes, and results from (CRM) simulations. Note that this presentation is mainly based on a recent paper published in Geophy. Rev. (Tao et al. 2012).
Recently, a detailed spectral-bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (GCE) model. Atmospheric aerosols are also described using number density size-distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and summertime convection over a mid-latitude continent with different concentrations of CCN: a low clean concentration and a high dirty concentration. The impact of atmospheric aerosol concentration on cloud and precipitation will be investigated.
In this study, it is proposed that twin tropical cyclones (TCs), Kesiny and 01A, in May 2002 formed in association with the scale interactions of three gyres that appeared as a convectively coupled mixed Rossby gravity (ccMRG) wave during an active phase of the Madden-Julian Oscillation (MJO). This is shown by analyzing observational data, including NCEP reanalysis data and METEOSAT 7 IR satellite imagery, and performing numerical simulations using a global mesoscale model. A 10-day control run is initialized at 0000 UTC 1 May 2002 with grid-scale condensation but no sub-grid cumulus parameterizations. The ccMRG wave was identified as encompassing two developing and one non-developing gyres, the first two of which intensified and evolved into the twin TCs. The control run is able to reproduce the evolution of the ccMRG wave and thus the formation of the twin TCs about two and five days in advance as well as their subsequent intensity evolution and movement within an 8-10 day period. Five additional 10-day sensitivity experiments with different model configurations are conducted to help understand the interaction of the three gyres, leading to the formation of the TCs. These experiments suggest the improved lead time in the control run may be attributed to the realistic simulation of the ccMRG wave with the following processes: (1) wave deepening (intensification) associated with a reduction in wavelength and/or the intensification of individual gyres, (2) poleward movement of gyres that may be associated with boundary layer processes, (3) realistic simulation of moist processes at regional scales in association with each of the gyres, and (4) the vertical phasing of low- and mid-level cyclonic circulations associated with a specific gyre.
In this study, it is proposed that twin tropical cyclones (TCs), Kesiny and 01A, in May 2002 formed in association with the scale interactions of three gyres that appeared as a convectively-coupled mixed Rossby gravity (ccMRG) wave during an active phase of the Madden-Julian Oscillation (MJO). This is shown by analyzing observational data and performing simulations using a global mesoscale model. A 10-day control run is initialized at 0000 UTC 1 May 2002 with grid-scale condensation but no cumulus parameterizations. The ccMRG wave was identified as encompassing two developing and one non-developing gyres, the first two of which intensified and evolved into the twin TCs. The control run is able to reproduce the evolution of the ccMRG wave and the formation of the twin TCs about two and five days in advance as well as their subsequent intensity evolution and movement within an 8-10 day period. Five additional 10-day sensitivity experiments with different model configurations are conducted to help understand the interaction of the three gyres. These experiments suggest the improved lead time in the control run may be attributed to the realistic simulation of the ccMRG wave with the following processes: (I) wave deepening associated with wave shortening and/or the intensification of individual gyres, (2) poleward movement of gyres that may be associated with bOlll1dary layer processes, (3) realistic simulation of moist processes at regional scales in association with each of the gyres, and (4) the vertical phasing of low- and mid-level cyclonic circulations associated with a specific gyre.
One of major CRM approaches to studying precipitation processes is sometimes referred to as "cloud ensemble modeling". This approach allows many clouds of various sizes and stages of their lifecycles to be present at any given simulation time. Large-scale effects derived from observations are imposed into CRMs as forcing, and cyclic lateral boundaries are used. The advantage of this approach is that model results in terms of rainfall and QI and Q2 usually are in good agreement with observations. In addition, the model results provide cloud statistics that represent different types of clouds/cloud systems during their lifetime (life cycle). The large-scale forcing derived from MC3EI will be used to drive GCE model simulations. The model-simulated results will be compared with observations from MC3E. These GCE model-simulated datasets are especially valuable for LH algorithm developers. In addition, the regional scale model with very high-resolution, NASA Unified WRF is also used to real time forecast during the MC3E campaign to ensure that the precipitation and other meteorological forecasts are available to the flight planning team and to interpret the forecast results in terms of proposed flight scenarios. Post Mission simulations are conducted to examine the sensitivity of initial and lateral boundary conditions to cloud and precipitation processes and rainfall. We will compare model results in terms of precipitation and surface rainfall using GCE model and NU-WRF