Engineering Papers⌕ Search

Engineering topics

Tao, Wei-Kuo

Publications and source records attributed to Tao, Wei-Kuo.

At least 181 records · Page 10

Goddard Cumulus Ensemble (GCE) Model: Application for Understanding Precipitation Processes

One of the most promising methods to test the representation of cloud processes used in climate models is to use observations together with Cloud Resolving Models (CRMs). The CRMs use more sophisticated and realistic representations of cloud microphysical processes, and they can reasonably well resolve the time evolution, structure, and life cycles of clouds and cloud systems (size about 2-200 km). The CRMs also allow explicit interaction between out-going longwave (cooling) and incoming solar (heating) radiation with clouds. Observations can provide the initial conditions and validation for CRM results. The Goddard Cumulus Ensemble (GCE) Model, a cloud-resolving model, has been developed and improved at NASA/Goddard Space Flight Center over the past two decades. Dr. Joanne Simpson played a central role in GCE modeling developments and applications. She was the lead author or co-author on more than forty GCE modeling papers. In this paper, a brief discussion and review of the application of the GCE model to (1) cloud interactions and mergers, (2) convective and stratiform interaction, (3) mechanisms of cloud-radiation interaction, (4) latent heating profiles and TRMM, and (5) responses of cloud systems to large-scale processes are provided. Comparisons between the GCE model's results, other cloud-resolving model results and observations are also examined.

Tao, Wei-Kuo↗

Parameterizations of Depositional Growth of Cloud Ice in a Bulk Microphysical Scheme

The conversion of cloud ice to snow by depositional growth, designated P(sub SFI), in the Goddard Cumulus Ensemble Model cloud physics parameterization is examined. The original formulation of P(sub SFI) is shown to produce excessive conversion of cloud ice to snow because of an implicit assumption that the relative humidity is 100% with respect to water even though the air may actually be quite less humid. Two possible corrections to this problem are proposed, the first involving application of a relative humidity dependent correction factor to the original formulation of P(sub SFI), and the second involving a new formulation of P(sub SFI) based on the equation for depositional growth of cloud ice.

Braun, Scott A.↗

Microphysics, Radiation and Surface Processes in the Goddard Cumulus Ensemble (GCE) Model

The response of cloud systems to their environment is an important link in a chain of processes responsible for monsoons, frontal depression, El Nino Southern Oscillation (ENSO) episodes and other climate variations (e.g., 30-60 day intra-seasonal oscillations). Numerical models of cloud properties provide essential insights into the interactions of clouds with each other, with their surroundings, and with land and ocean surfaces. Significant advances are currently being made in the modeling of rainfall and rain-related cloud processes, ranging in scales from the very small up to the simulation of an extensive population of raining cumulus clouds in a tropical- or midlatitude-storm environment. The Goddard Cumulus Ensemble (GCE) model is a multi-dimensional nonhydrostatic dynamic/microphysical cloud resolving model. It has been used to simulate many different mesoscale convective systems that occurred in various geographic locations. In this paper, recent GCE model improvements (microphysics, radiation and surface processes) will be described as well as their impact on the development of precipitation events from various geographic locations. The performance of these new physical processes will be examined by comparing the model results with observations. In addition, the explicit interactive processes between cloud, radiation and surface processes will be discussed.

Tao, Wei-Kuo↗

A Numerical Study of Tropical Sea-Air Interactions Using a Cloud Resolving Model Coupled with an Ocean Mixed-Layer Model

Coupling a cloud resolving model (CRM) with an ocean mixed layer (OML) model can provide a powerful tool for better understanding impacts of atmospheric precipitation on sea surface temperature (SST) and salinity. The objective of this study is twofold. First, by using the three dimensional (3-D) CRM-simulated (the Goddard Cumulus Ensemble model, GCE) diabatic source terms, radiation (longwave and shortwave), surface fluxes (sensible and latent heat, and wind stress), and precipitation as input for the OML model, the respective impact of individual component on upper ocean heat and salt budgets are investigated. Secondly, a two-way air-sea interaction between tropical atmospheric climates (involving atmospheric radiative-convective processes) and upper ocean boundary layer is also examined using a coupled two dimensional (2-D) GCE and OML model. Results presented here, however, only involve the first aspect. Complete results will be presented at the conference.

Shie, Chung-Lin↗

Convective Systems Observed and Simulated During TRMM Field Campaigns

Recently completed TRMM field campaigns (TEFLUN1998, SCSMEX-1998, TRMM.LBA-1999, and KWAJEX 1999) have obtained direct measurements of microphysical data associated with convective systems from various geographical locations. These TRMM field experiments were designed to contribute to fundamental understanding of cloud dynamics and microphysics, as well as for validation,, testing assumptions and error estimates of cloud-resolving models, forward radiative transfer models, algorithms used to estimate rainfall statistics and vertical structure of precipitation and latent heating from both surface-based radar and satellites.

Tao, Wei-Kuo↗

The Goddard Cumulus Ensemble Model: Model Description and Its Application for Studying the TOGA COARE and GATE Convective Systems

The Goddard Cumulus Ensemble (GCE) model was utilized in two and three dimensions in order to examine the behavior and response of simulated deep tropical cloud systems occurred in west Pacific warm pool region and Atlantic ocean. The periods chosen for simulation were convectively active period over the TOGA-COARE IFA (19-27 December 1992) and GATE (September 1 to 7, 1974). The TOGA COARE IFA period was also in the framework of the GEWEX Cloud System Study (GCSS) WG4 case 2. We will examine the differences between the microphysics (warm rain and ice processes, evaporation/sublimation and condensation/deposition), Q1 (Temperature) and Q2 (Water vapor) budgets between these two convective events occurred in different large-scale environments. The contribution of stratiform precipitation and its relationship to the vertical shear of the large-scale horizontal wind will also be examined. The results from GCSS model intercomparsion will be presented. The new improvements (i.e., microphysics, cloud radiation interaction, surface processes and numerical advection scheme) of the GCE model as well as their sensitivity to the model results will be discussed.

Tao, Wei-Kuo↗

The Impact of TRMM Data on Mesoscale Numerical Simulation of Super Typhoon Paka

This paper assesses the impact of TRMM Microwave Imager (TMI) derived surface rainfall data on the numerical simulation of Super Typhoon Paka (1997). A series of mesoscale numerical simulations is conducted for simulating Super Typhoon Paka during its mature stage by using the Penn State/NCAR MM5 model with initial conditions derived from the Goddard Earth Observing System (GEOS) global analyses with and without assimilation of the TMI surface rainfall data. Simulation results clearly demonstrate that the GEOS analysis with TMI rainfall data leads to a improved MM5 simulation of Typhoon Paka in terms of its intensity and kinematical and precipitation structures, since inclusion of the rainfall data into initial and boundary condition improves the storm environmental conditions (e.g., moisture and circulation). Further improvement can also be obtained by incorporating TMI data with mesoscale bogus vortex information into the initial conditions.

Pu, Zhao-Xia↗

Cloud Modeling

Numerical cloud models have been developed and applied extensively to study cloud-scale and mesoscale processes during the past four decades. The distinctive aspect of these cloud models is their ability to treat explicitly (or resolve) cloud-scale dynamics. This requires the cloud models to be formulated from the non-hydrostatic equations of motion that explicitly include the vertical acceleration terms since the vertical and horizontal scales of convection are similar. Such models are also necessary in order to allow gravity waves, such as those triggered by clouds, to be resolved explicitly. In contrast, the hydrostatic approximation, usually applied in global or regional models, does allow the presence of gravity waves. In addition, the availability of exponentially increasing computer capabilities has resulted in time integrations increasing from hours to days, domain grids boxes (points) increasing from less than 2000 to more than 2,500,000 grid points with 500 to 1000 m resolution, and 3-D models becoming increasingly prevalent. The cloud resolving model is now at a stage where it can provide reasonably accurate statistical information of the sub-grid, cloud-resolving processes poorly parameterized in climate models and numerical prediction models.

Tao, Wei-Kuo↗

The Impact of TRMM Data on Mesoscale Numerical Simulation of Super Typhoon Paka

Accurate measurement of the spatial and temporal variations of tropical rainfall around the globe had remained as a critical problem in meteorology until the recent launch of the Tropical Rainfall Measuring Mission (TRMM). TRMM offers a unique opportunity to improve understanding of tropical meteorology and also offers a great opportunity to evaluate the impact of rainfall data on tropical weather forecasts. This study assesses the impact of TRMM Microwave Imager (TMI) derived surface rainfall data on the numerical simulation of Super Typhoon Paka (1997).

Pu, Zhao-Xia↗

Latent Heating Retrievals Using the TRMM Precipitation Radar: A Multi-Seasonal Study

Rainfall is a key link in the hydrologic cycle and is a primary heat source for the atmosphere. The vertical distribution of latent heat release, which is accompanied by rainfall, modulates the large-scale circulations of the tropics and in turn can impact midlatitude weather. This latent heat release is a consequence of phase changes between vapor, liquid, and solid water. Present largescale weather and climate models can simulate latent heat release only crudely, thus reducing their confidence in predictions on both global and regional scales. This paper represents the first attempt to use NASA Tropical Rainfall Measuring Mission (TRMM) rainfall information to estimate the four-dimensional structure of global monthly latent heating profiles over the global tropics from December 1997 to October 2000. The Goddard Convective-Stratiform. Heating (CSH) algorithm and TRMM precipitation radar data are used for this study. We will examine and compare the latent heating structures between 1997-1998 (winter) ENSO and 1998-2000 (non-ENSO). We will also examine over the tropics. The seasonal variation of heating over various geographic locations (i.e., oceanic vs continental; Indian oceans vs west Pacific; Africa vs S. America) will be also examined and compared. In addition, we will examine the relationship between latent heating (max heating level) and SST. The period of interest also coincides with several TRMM field campaigns that recently occurred over the South China Sea in 1998 (SCSMEX), Brazil in 1999 (TRMM-LBA), and in the central Pacific in 1999 (KWAJEX). Sounding diagnosed Q1 budgets from these experiments could provide a means of validating the retrieved profiles of latent heating from the CSH algorithm.

Tao, Wei-Kuo↗

The Sensitivity of West African Squall Line Water Budgets to Land Cover

This study used a two-dimensional coupled land/atmosphere (cloud-resolving) model to investigate the influence of land cover on the water budgets of squall lines in the Sahel. Study simulations used the same initial sounding and one of three different land covers, a sparsely vegetated semi-desert, a grassy savanna, and a dense evergreen broadleaf forest. All simulations began at midnight and ran for 24 hours to capture a full diurnal cycle. In the morning, the latent heat flux, boundary layer mixing ratio, and moist static energy in the boundary layer exhibited notable variations among the three land covers. The broadleaf forest had the highest latent heat flux, the shallowest, moistest, slowest growing boundary layer, and significantly more moist static energy per unit area than the savanna and semi-desert. Although all simulations produced squall lines by early afternoon, the broadleaf forest had the most intense, longest-lived squall lines with 29% more rainfall than the savanna and 37% more than the semi-desert. The sensitivity of the results to vegetation density, initial sounding humidity, and grid resolution was also assessed. There were greater differences in rainfall among land cover types than among simulations of the same land cover with varying amounts of vegetation. Small changes in humidity were equivalent in effect to large changes in land cover, producing large changes in the condensate and rainfall. Decreasing the humidity had a greater effect on rainfall volume than increasing the humidity. Reducing the grid resolution from 1.5 km to 0.5 km decreased the temperature and humidity of the cold pools and increased the rain volume.

Mohr, Karen I.↗

On the Numerical Study of Heavy Rainfall in Taiwan

Heavy rainfall events are frequently observed over the western side of the CMR (central mountain range), which runs through Taiwan in a north-south orientation, in a southwesterly flow regime and over the northeastern side of the CMR in a northeasterly flow regime. Previous studies have revealed the mechanisms by which the heavy rainfall events are formed. Some of them have examined characteristics of the heavy rainfall via numerical simulations. In this paper, some of the previous numerical studies on heavy rainfall events around Taiwan during the Mei-Yu season (May and June), summer (non-typhoon cases) and autumn will be reviewed. Associated mechanisms proposed from observational studies will be reviewed first, and then characteristics of numerically simulated heavy rainfall events will be presented. The formation mechanisms of heavy rainfall from simulated results and from observational analysis are then compared and discussed. Based on these previous modeling studies, we will also discuss what are the major observations and modeling processes which will be needed for understanding the heavy precipitation in the future.

Tao, Wei-Kuo↗

Momentum Transport: 2D and 3D Cloud Resolving Model Simulations

The major objective of this study is to investigate the momentum budgets associated with several convective systems that developed during the TOGA COARE IOP (west Pacific warm pool region) and GATE (east Atlantic region). The tool for this study is the improved Goddard Cumulas Ensemble (GCE) model which includes a 3-class ice-phase microphysical scheme, explicit cloud radiative interactive processes and air-sea interactive surface processes. The model domain contains 256 x 256 grid points (with 2 km resolution) in the horizontal and 38 grid points (to a depth of 22 km) in the vertical. The 2D domain has 1024 grid points. The simulations were performed over a 7-day time period (December 19-26, 1992, for TOGA COARE and September 1-7, 1994 for GATE). Cyclic literal boundary conditions are required for this type of long-term integration. Two well organized squall systems (TOGA, COARE February 22, 1993, and GATE September 12, 1994) were also simulated using the 3D GCE model. Only 9 h simulations were required to cover the life time of the squall systems. the lateral boundary conditions were open for these two squall systems simulations. the following will be examined: (1) the momentum budgets in the convective and stratiform regions, (2) the relationship between momentum transport and cloud organization (i.e., well organized squall lines versus less organized convective), (3) the differences and similarities in momentum transport between 2D and 3D simulated convective systems, and (4) the differences and similarities in momentum budgets between cloud systems simulated with open and cyclic lateral boundary conditions. Preliminary results indicate that there are only small differences between 2D and 3D simulated momentum budgets. Major differences occur, however, between momentum budgets associated with squall systems simulated using different lateral boundary conditions.

Tao, Wei-Kuo↗

A Melting Layer Model for Passive/Active Microwave Remote Sensing Applications: Simulation of TRMM Observations - Part 2

The one-dimensional, steady-state melting layer model developed in Part I of this study is used to calculate both the microphysical and radiative properties of melting precipitation, based upon the computed concentrations of snow and graupel just above the freezing level at applicable horizontal gridpoints of 3-dimensional cloud resolving model simulations. The modified 3-dimensional distributions of precipitation properties serve as input to radiative transfer calculations of upwelling radiances and radar extinction/reflectivities at the TRMM Microwave Imager (TMI) and Precipitation Radar (PR) frequencies, respectively. At the resolution of the cloud resolving model grids (approx. 1 km), upwelling radiances generally increase if mixed-phase precipitation is included in the model atmosphere. The magnitude of the increase depends upon the optical thickness of the cloud and precipitation, as well as the scattering characteristics of ice-phase precipitation aloft. Over the set of cloud resolving model simulations utilized in this study, maximum radiance increases of 43, 28, 18, and 10 K are simulated at 10.65, 19.35 GHz, 37.0, and 85.5 GHz, respectively. The impact of melting on TMI-measured radiances is determined not only by the physics of the melting particles but also by the horizontal extent of the melting precipitation, since the lower-frequency channels have footprints that extend over 10''s of kilometers. At TMI resolution, the maximum radiance increases are 16, 15, 12, and 9 K at the same frequencies. Simulated PR extinction and reflectivities in the melting layer can increase dramatically if mixed-phase precipitation is included, a result consistent with previous studies. Maximum increases of 0.46 (-2 dB) in extinction optical depth and 5 dBZ in reflectivity are simulated based upon the set of cloud resolving model simulations.

Olson, William S.↗

A Melting Layer Model for Passive/Active Microwave Remote Sensing Applications: Model Formulation and Comparison with Observations - Part 1

In this study, a 1-D steady-state microphysical model which describes the vertical distribution of melting precipitation particles is developed. The model is driven by the ice-phase precipitation distributions just above the freezing level at applicable gridpoints of "parent" 3-D cloud-resolving model (CRM) simulations. It extends these simulations by providing the number density and meltwater fraction of each particle in finely separated size categories through the melting layer. The depth of the modeled melting layer is primarily determined by the initial material density of the ice-phase precipitation. The radiative properties of melting precipitation at microwave frequencies are calculated based upon different methods for describing the dielectric properties of mixed phase particles. Particle absorption and scattering efficiencies at the Tropical Rainfall Measuring Mission Microwave Imager frequencies (10.65 to 85.5 GHz) are enhanced greatly for relatively small (approx. 0.1) meltwater fractions. The relatively large number of partially-melted particles just below the freezing level in stratiform regions leads to significant microwave absorption, well-exceeding the absorption by rain at the base of the melting layer. Calculated precipitation backscatter efficiencies at the Precipitation Radar frequency (13.8 GHz) increase in proportion to the particle meltwater fraction, leading to a "bright-band" of enhanced radar reflectivities in agreement with previous studies. The radiative properties of the melting layer are determined by the choice of dielectric models and the initial water contents and material densities of the "seeding" ice-phase precipitation particles. Simulated melting layer profiles based upon snow described by the Fabry-Szyrmer core-shell dielectric model and graupel described by the Maxwell-Garnett water matrix dielectric model lead to reasonable agreement with radar-derived melting layer optical depth distributions. Moreover, control profiles that do not contain mixed-phase precipitation particles yield optical depths that are systematically lower than those observed. Therefore, the use of the melting layer model to extend 3-D CRM simulations appears justified, at least until more realistic spectral methods for describing melting precipitation in high-resolution, 3-D CRM's are implemented.

Olson, William S.↗

The Influence of Soil Moisture, Coastline Curvature, and Land-Breeze Circulations on Sea-Breeze Initiated Precipitation

Idealized numerical simulations are performed with a coupled atmosphere/land-surface model to identify the roles of initial soil moisture, coastline curvature, and land breeze circulations on sea breeze initiated precipitation. Data collected on 27 July 1991 during the Convection and Precipitation Electrification Experiment (CAPE) in central Florida are used. The 3D Goddard Cumulus Ensemble (GCE) cloud resolving model is coupled with the Goddard Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) land surface model, thus providing a tool to simulate more realistically land-surface/atmosphere interaction and convective initiation. Eight simulations are conducted with either straight or curved coast-lines, initially homogeneous soil moisture or initially variable soil moisture, and initially homogeneous horizontal winds or initially variable horizontal winds (land breezes). All model simulations capture the diurnal evolution and general distribution of sea-breeze initiated precipitation over central Florida. The distribution of initial soil moisture influences the timing, intensity and location of subsequent precipitation. Soil moisture acts as a moisture source for the atmosphere, increases the connectively available potential energy, and thus preferentially focuses heavy precipitation over existing wet soil. Strong soil moisture-induced mesoscale circulations are not evident in these simulations. Coastline curvature has a major impact on the timing and location of precipitation. Earlier low-level convergence occurs inland of convex coastlines, and subsequent precipitation occurs earlier in simulations with curved coastlines. The presence of initial land breezes alone has little impact on subsequent precipitation. however, simulations with both coastline curvature and initial land breezes produce significantly larger peak rain rates due to nonlinear interactions.

Baker, David R.↗

Kinematic and Microphysical Structures of Hurricane Bob (1991) Determined from a 1.3-km-Resolution Numerical Simulation

Hurricane Bob (1991) is simulated using the Penn State/NCAR mesoscale model MM5. The simulation is conducted for a 24-h period at 4-km resolution and for a 6-h period at 1.3-km resolution. The 4-km simulation is able to fairly realistically capture the intensity and structure of the storm. The 1.3-km simulation depicts very small-scale convective structures and produces a convective band outside of the eye wall that did not occur in the 4-km simulation. The 1.3-km results are used to characterize several kinematic and cloud microphysical structures in the storm. Characteristics of air parcels flowing into and rising within the eye wall will be examined through trajectory calculations.

Braun, Scott A.↗

Kinematic and Thermodynamic Structures of Hurricane Bob (1991) Determined From A 1.3-KM Resolution Numerical Simulation

Hurricane Bob (1991) is simulated using the Penn State/NCAR mesoscale model MM5. The simulation is conducted for a 24-h period at 4-km resolution and for a 6-h period at 1.3-km resolution. The 4-km simulation is able to fairly realistically capture the intensity and structure of the storm. The 1.3-km simulation depicts very small-scale convective structures and produces a convective rain band outside of the eye wall that did not occur in the 4-km simulation. The 1.3-km results are used to characterize several kinematic structures in the storm, including low-level outflow in the eye wall, multiple outflows at upper levels, and the convective rain band structure outside of the eye wall. Thermodynamic characteristics of air parcels flowing into and rising within the eye wall will also be examined through trajectory calculations.

Braun, Scott A.↗