Search

Several heavy precipitation episodes occurred over Taiwan from August 10 to 13, 1994. Precipitation patterns and characteristics are quite different between the precipitation events that occurred from August 10 and I I and from August 12 and 13. In Part I (Chen et al. 2001), the environmental situation and precipitation characteristics are analyzed using the EC/TOGA data, ground-based radar data, surface rainfall patterns, surface wind data, and upper air soundings. In this study (Part II), the Penn State/NCAR Mesoscale Model (MM5) is used to study the precipitation characteristics of these heavy precipitation events. Various physical processes (schemes) developed at NASA Goddard Space Flight Center (i.e., cloud microphysics scheme, radiative transfer model, and land-soil-vegetation surface model) have recently implemented into the MM5. These physical packages are described in the paper, Two way interactive nested grids are used with horizontal resolutions of 45, 15 and 5 km. The model results indicated that Cloud physics, land surface and radiation processes generally do not change the location (horizontal distribution) of heavy precipitation. The Goddard 3-class ice scheme produced more rainfall than the 2-class scheme. The Goddard multi-broad-band radiative transfer model reduced precipitation compared to a one-broad band (emissivity) radiation model. The Goddard land-soil-vegetation surface model also reduce the rainfall compared to a simple surface model in which the surface temperature is computed from a Surface energy budget following the "force-re store" method. However, model runs including all Goddard physical processes enhanced precipitation significantly for both cases. The results from these runs are in better agreement with observations. Despite improved simulations using different physical schemes, there are still some deficiencies in the model simulations. Some potential problems are discussed. Sensitivity tests (removing either terrain or radiative processes) are performed to identify the physical processes that determine the precipitation patterns and characteristics for heavy rainfall events. These sensitivity tests indicated that terrain can play a major role in determining the exact location for both precipitation events. The terrain can also play a major role in determining the intensity of precipitation for both events. However, it has a large impact on one event but a smaller one on the other. The radiative processes are also important for determining, the precipitation patterns for one case but. not the other. The radiative processes can also effect the total rainfall for both cases to different extents.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

The global hydrological cycle is central to climate system interactions and the key to understanding their behavior. Rainfall and its associated precipitation processes are a key link in the hydrologic cycle. Fresh water provided by tropical rainfall and its variability can exert a large impact upon the structure of the upper ocean layer. In addition, approximately two-thirds of the global rain falls in the Tropics, while the associated latent heat release accounts for about three-fourths of the total heat energy for the Earth's atmosphere. Precipitation from convective cloud systems comprises a large portion of tropical heating and rainfall. Furthermore, the vertical distribution of convective latent-heat releases modulates large-scale tropical circulations (e.g., the 30-60-day intraseasonal oscillation), which, in turn, impacts midlatitude weather through teleconnection patterns such as those associated with El Nino. Shifts in these global circulations can result in prolonged periods of droughts and floods, thereby exerting a tremendous impact upon the biosphere and human habitation. And yet, monthly rainfall over the tropical oceans is still not known within a factor of two over large (5 degrees latitude by 5 degrees longitude) areas. Hence, the Tropical Rainfall Measuring Mission (TRMM), a joint U.S./Japan space project, can provide a more accurate measurement of rainfall as well as estimate the four-dimensional structure of diabatic heating over the global tropics. The distributions of rainfall and inferred heating can be used to advance our understanding of the global energy and water cycle. In addition, this information can be used for global circulation and climate models for testing and improving their parameterizations.

The 1993 US Midwest Flood produced record levels of flooding, in the Mississippi River Basin. This flooding resulted from repeated frontal passages and mesoscale convective complexes (MCCs) during the months of June and July. A better understanding of processes that influenced MCC development during the 1993 Flood may lead to improved forecasts of heavy precipitation and flooding. Here, we consider the impact of soil moisture on MCC development during a two-day period (June 23-24) of the 1993 US Midwest Flood. The purpose of this study is to assess the importance of soil moisture distribution on the timing, intensity, and location of heavy precipitation. In this study, the MM5-PLACE Atmosphere/Land-Surface Model is utilized. The atmospheric component consists of the Penn State/NCAR MM5 mesoscale model, and the land-surface component consists of the Goddard Parameterization for Land Atmosphere-Cloud Exchange (PLACE). Initial soil moisture is provided from two sources: 1) NCEP reanalysis, and 2) Antecedent Precipitation Index (API) using NOAA rain gauge measurements as a proxy for soil moisture. NCEP reanalysis provides coarse resolution initial soil moisture (2.5 degree), while API provides high resolution initial soil moisture (10-200 km depending on NOAA rain gauge spacing). Initial results indicate that the distribution of soil moisture has a significant impact on the timing and location of heavy precipitation during this two-day flood event. Precipitation in simulations with high resolution initial soil moisture agrees more closely with observed precipitation. These results suggest that high resolution soil moisture observations are necessary to accurately predict severe storm development, heavy precipitation, and subsequent flooding.

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.

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.

Two aspects of the cloud ice parameterization in the Goddard Cumulus Ensemble Model cloud physics parameterization are examined: the conversion of cloud ice to snow by depositional growth, designated PSFI, and the saturation adjustment scheme. The original formulation of PSFI 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 PSFI, and the second involving a new formulation of PSFI based on the equation for depositional growth of cloud ice. The sensitivity of these formulations of PSFI to the assumed masses of the ice particles is examined. Possible problems associated with using a saturation adjustment scheme for cloud ice are discussed and simulations of a squall line with and without application of the adjustment scheme for ice are compared.

A sophisticated land-surface model, PLACE, the Parameterization for Land Atmospheric Convective Exchange, has been coupled to a 1.5-order turbulent kinetic energy (TKE) turbulence sub-model. Both have been incorporated into the Penn State/National Center for Atmospheric Research (PSU/NCAR) mesoscale model MM5. Such model improvements should have their greatest effect in conditions where surface contrasts dominate over dynamic processes, such as the simulation of warm-season, convective events. A validation study used the newly coupled model, MM5 TKE-PLACE, to simulate the evolution of Florida sea-breeze moist convection during the Convection and Precipitation Electrification Experiment (CaPE). Overall, eight simulations tested the sensitivity of the MM5 model to combinations of the new and default model physics, and initialization of soil moisture and temperature. The TKE-PLACE model produced more realistic surface sensible heat flux, lower biases for surface variables, more realistic rainfall, and cloud cover than the default model. Of the 8 simulations with different factors (i.e., model physics or initialization), TKE-PLACE compared very well when each simulation was ranked in terms of biases of the surface variables and rainfall, and percent and root mean square of cloud cover. A factor separation analysis showed that a successful simulation required the inclusion of a multi-layered, land surface soil vegetation model, realistic initial soil moisture, and higher order closure of the planetary boundary layer (PBL). These were needed to realistically model the effect of individual, joint, and synergistic contributions from the land surface and PBL on the CAPE sea-breeze, Lake Okeechobee lake breeze, and moist convection.

Latent heating profiles associated with three (TOGA COARE) Tropical Ocean and Global Atmosphere Coupled Ocean Atmosphere Response Experiment active convective episodes (December 10-17 1992; December 19-27 1992; and February 9-13 1993) are examined using the Goddard Cumulus Ensemble (GCE) Model and retrieved by using the Goddard Convective and Stratiform Heating (CSH) algorithm . The following sources of rainfall information are input into the CSH algorithm: Special Sensor Microwave Imager (SSM/1), Radar and the GCE model. Diagnostically determined latent heating profiles calculated using 6 hourly soundings are used for validation. The GCE model simulated rainfall and latent heating profiles are in excellent agreement with those estimated by soundings. In addition, the typical convective and stratiform heating structures (or shapes) are well captured by the GCE model. Radar measured rainfall is smaller than that both estimated by the GCE model and SSM/I in all three different COARE IFA periods. SSM/I derived rainfall is more than the GCE model simulated for the December 19-27 and February 9-13 periods, but is in excellent agreement with the GCE model for the December 10-17 period. The GCE model estimated stratiform amount is about 50% for December 19-27, 42% for December 11-17 and 56% for the February 9-13 case. These results are consistent with large-scale analyses. The accurate estimates of stratiform amount is needed for good latent heating retrieval. A higher (lower) percentage of stratiform rain can imply a maximum heating rate at a higher (lower) altitude. The GCE model always simulates more stratiform rain (10 to 20%) than the radar for all three convective episodes. SSM/I derived stratiform amount is about 37% for December 19-27, 48% for December 11-17 and 41% for the February 9-13 case. Temporal variability of CSH algorithm retrieved latent heating profiles using either GCE model simulated or radar estimated rainfall and stratiform amount is in good agreement with that diagnostically determined for all three periods. However, less rainfall and a smaller stratiform percentage estimated by radar resulted in a weaker (underestimated) latent heating profile and a lower maximum latent heating level compared to those determined diagnostically. Rainfall information from SSM/I can not retrieve individual convective events due to poor temporal sampling. Nevertheless, this study suggests that a good 4r, rainfall retrieval from SSM/I for a convective event always leads to a good latent heating retrieval. Sensitivity testing has been performed and the results indicate that the SSM/I derived time averaged stratiform amount may be underestimated for December 19-27. Time averaged heating profiles derived from SSM/I, however, are not in bad agreement with those derived by soundings for the December 10-17 convective period. The heating retrievals may be more accurate for longer time scales provided there is no bias in the sampling.

Using the Penn State/NCAR MM5 mesoscale model, a westerly wind burst (WWB) that occurred during the period from 19 to 30 December 1992 over the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE) has been simulated and compared with observation. This event is characterized by the presence of super cloud clusters and the occurrence of a major WWB that extended over the western and central Pacific Ocean. Although several of the observed convective systems were not precisely simulated by MM5, the model did capture many other observed characteristics, such as the explosive development of convection, the cyclonic circulation and the WWB. The WWB resulted from the coalescence of three types of tropical disturbances. The first type was a low-level westerly jet (LWJ) that developed at the equator and may be associated with the eastward propagation of an ISO (Intraseasonal Oscillation). The second type featured an easterly wave-like disturbance that originated in the south central Pacific Ocean and propagated westward. Finally, the third type involved a cross-equatorial flow that deflected Northern Hemispheric easterlies into the Southern Hemisphere and may be caused by inertial instability. These disturbances worked in concert, resulting in intense convection over the TOGA COARE region. Once intense convection developed, a large-scale circulation was produced over the western Pacific warm pool, propagated eastward, and initiated a WWB.

The MM5 mesoscale model is used to simulate Hurricane Bob (1991) using grids nested to high resolution (4 km). Tests are conducted to determine the sensitivity of the simulation to the available planetary boundary layer parameterizations, including the bulk-aerodynamic, Blackadar, Medium-RanGe Forecast (MRF) model, and Burk-Thompson boundary-layer schemes. Significant sensitivity is seen, with minimum central pressures varying by up to 17 mb. The Burk-Thompson and bulk-aerodynamic boundary-layer schemes produced the strongest storms while the MRF scheme produced the weakest storm. Precipitation structure of the simulated hurricanes also varied substantially with the boundary layer parameterizations. Diagnostics of boundary-layer variables indicated that the intensity of the simulated hurricanes generally increased as the ratio of the surface exchange coefficients for heat and momentum, C(sub h)/C(sub M), although the manner in which the vertical mixing takes place was also important. Findings specific to the boundary-layer schemes include: 1) the MRF scheme produces mixing that is too deep and causes drying of the lower boundary layer in the inner-core region of the hurricane; 2) the bulk-aerodynamic scheme produces mixing that is probably too shallow, but results in a strong hurricane because of a large value of C(sub h)/C(sub M) (approximately 1.3); 3) the MRF and Blackadar schemes are weak partly because of smaller surface moisture fluxes that result in a reduced value of C(sub h)/C(sub M) (approximately 0.7); 4) the Burk-Thompson scheme produces a strong storm with C(sub h)/C(sub M) approximately 1; and 5) the formulation of the wind-speed dependence of the surface roughness parameter, z(sub 0), is important for getting appropriate values of the surface exchange coefficients in hurricanes based upon current estimates of these parameters.

Mesoscale research and forecast models are increasingly being used at horizontal resolutions of 1-8 km to simulate a variety of precipitating systems. When the model is used to simulate convective systems, it is uncertain to what extent the dynamics and microphysics of convective updrafts can be resolved with grids larger than 1 km. In this study, two- and three-dimensional versions of the Goddard Cumulus Ensemble model are used to determine the impact of horizontal grid resolution on the behavior of the simulated storms and on the characteristics of the cloud microphysical fields. It will be shown that as resolution decreases from about 1 km to greater than 3 km, there is a fairly rapid degradation of the storm structure in the form of reduced convective mass fluxes, updraft tilts, and cloud microphysics. A high-resolution simulation of hurricane outer rainbands using the MM5 mesoscale model shows also that there can be a substantial modification of the key microphysical processes that contribute to rainfall as a result of reducing the horizontal resolution.

Land surface-atmosphere interaction plays a key role in the development of summertime convection and precipitation over the Florida peninsula. Land-ocean temperature contrasts induce sea-breeze circulations along both coasts. Clouds develop along sea-breeze fronts, and significant precipitation can occur during the summer months. However, other factors such as soil moisture distribution and coastline curvature may modulate the timing, location, and intensity of sea breeze initiated precipitation. Here, we investigate the role of soil moisture and coastline curvature on Florida precipitation using the 3-D Goddard Cumulus Ensemble (GCE) cloud model coupled with the Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) land surface model. This study utilizes data from the Convection and Precipitation Electrification Experiment (CaPE) collected on 27 July 1991. Our numerical simulations suggest that a realistic distribution of soil moisture influences the location and intensity of precipitation but not the timing of precipitation. In contrast, coastline curvature affects the timing and location of precipitation but has little influence on peak rainfall rates. However, both factors (soil moisture and coastline curvature) are required to fully account for observed rainfall amounts.

A set of relatively high resolution three-dimensional (3D) simulations were produced to investigate the triggering of moist convection by landscape generated mesoscale circulations. The local accumulated rainfall varied monotonically (linearly) with the size of individual landscape patches, demonstrating the need to develop a trigger function that is sensitive to the size of individual patches. A new triggering function that includes the effect of landscapes generated mesoscale circulations over patches of different sizes consists of a parcel's perturbation in vertical velocity (nu(sub 0)), temperature (theta(sub 0)), and moisture (q(sub 0)). Each variable in the triggering function was also sensitive to soil moisture gradients, atmospheric initial conditions, and moist processes. The parcel's vertical velocity, temperature, and moisture perturbation were partitioned into mesoscale and turbulent components. Budget equations were derived for theta(sub 0) and q(sub 0). Of the many terms in this set of budget equations, the turbulent, vertical flux of the mesoscale temperature and moisture contributed most to the triggering of moist convection through the impact of these fluxes on the parcel's temperature and moisture profile. These fluxes needed to be parameterized to obtain theta(sub 0) and q(sub 0). The mesoscale vertical velocity also affected the profile of nu(sub 0). We used similarity theory to parameterize these fluxes as well as the parcel's mesoscale vertical velocity.

Land surface processes play a key role in water and energy budgets of the hydrological cycle. For example, the distribution of soil moisture will affect sensible and latent heat fluxes, which in turn may dramatically influence the location and intensity of precipitation. However, mean wind conditions also strongly influence the distribution of precipitation. The relative importance of soil moisture and wind on rainfall location and intensity remains uncertain. Here, we examine the influence of soil moisture distribution and wind distribution on precipitation in the Florida peninsula using the 3-D Goddard Cumulus Ensemble (GCE) cloud model Coupled with the Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) land surface model. This study utilizes data collected on 27 July 1991 in central Florida during the Convection and Precipitation Electrification Experiment (CaPE). The idealized numerical experiments consider a block of land (the Florida peninsula) bordered on the east and on the west by ocean. The initial soil moisture distribution is derived from an offline PLACE simulation, and the initial environmental wind profile is determined from the CaPE sounding network. Using the factor separation technique, the precise contribution of soil moisture and wind to rainfall distribution and intensity is determined.

Test data are required to test algorithms for the TRMM Precipitation Radar. These data are needed to test the design of the computer codes under development for the operational phase of the mission, and also to test and evaluate alternative or improved precipitation retrieval algorithms. Over a number of years we have developed and used a 3-dimensional radar model for simulating spaceborne precipitation radars. We have adapted this code to produce data files as close as possible to the TRMM file specifications. In this paper, we will describe the model as it is currently implemented, and show some samples of the synthetic data sets.

A preliminary comparison of the GEOS-1 (Goddard Earth Observing System) data assimilation system convective cloud mass fluxes with fluxes from a cloud-resolving model (the Goddard Cumulus Ensemble Model, GCE) is reported. A squall line case study (10-11 June 1985 Oklahoma PRESTORM episode) is the basis of the comparison. Regional (central U. S.) monthly total convective mass flux for June 1985 from GEOS-1 compares favorably with estimates from a statistical/dynamical approach using GCE simulations and satellite-derived cloud observations. The GEOS-1 convective mass fluxes produce reasonable estimates of monthly-averaged regional convective venting of CO from the boundary layer at least in an urban-influenced continental region, suggesting that they can be used in tracer transport simulations.

We have constructed a regional budget for boundary layer carbon monoxide over the central United States (32.5 deg - 50 deg N, 90 deg - 105 deg W), emphasizing a detailed evaluation of deep convective vertical fluxes appropriate for the month of June. Deep convective venting of the boundary layer (upward) dominates other components of the CO budget, e.g., downward convective transport, loss of CO by oxidation, anthropogenic emissions, and CO produced from oxidation of methane, isoprene, and anthropogenic nonmethane hydrocarbons (NMHCs). Calculations of deep convective venting are based on the method pf Pickering et al.(1992a) which uses a satellite-derived deep convective cloud climatology along with transport statistics from convective cloud model simulations of observed prototype squall line events. This study uses analyses of convective episodes in 1985 and 1989 and CO measurements taken during several midwestern field campaigns. Deep convective venting of the boundary layer over this moderately polluted region provides a net (upward minus downward) flux of 18.1 x 10(exp 8) kg CO/month to the free troposphere during early summer. Shallow cumulus and synoptic-scale weather systems together make a comparable contribution (total net flux 16.2 x 10(exp 8) kg CO/month). Boundary layer venting of CO with other O3 precursors leads to efficient free troposheric O3 formation. We estimate that deep convective transport of CO and other precursors over the central United States in early summer leads to a gross production of 0.66 - 1.1 Gmol O3/d in good agreement with estimates of O3 production from boundary layer venting in a continental-scale model (Jacob et al., 1993a, b). On this respect the central U.S. region acts as s `chimney' for the country, and presumably this O3 contributes to high background levels of O3 in the eastern United States and O3 export to the North Atlantic.

It is found from model simulations of trace gas and meteorological data from aircraft campaigns that deep convection may enhance the potential for photochemical ozone production in the middle and upper troposphere by up to a factor of 60. Examination of half a dozen individual convective episodes show that the degree of enhancement is highly variable. Factors affecting enhancement include boundary layer NO(x) mixing ratios, differences in the strength and structure of convective cells, as well as variation in the amount of background pollution already in the free troposphere.

The topics covered include the following: Regional Acid Deposition Model (RADM) -- a coupled chemistry/mesoscale model; convection in RADM; unresolved issues for mesoscale modeling with chemistry -- nonprecipitating clouds; unresolved issues for mesoscale modeling with chemistry -- aerosols; tracer studies with Goddard Cumulus Ensemble Model (GCEM); field observations of trace gas transport in convection; and photochemical consequences of convection.

The chemistry and dynamics of a convective system observed on February 2, 1987 in the EMEX and STEP campaigns are analyzed. Chemical and thermodynamic profiles in undisturbed air near the EMEX 9 system indicate that the troposphere was well mixed by previous convection. As a consequence there was little direct transport from the boundary layer to the upper troposphere and air transported upward was detrained throughout the middle and upper troposphere. There was minimal effect on O3 production. Other convective complexes located 800-900 km upstream produced greater perturbations on trace gas profile immediately below the tropopause. Ozone production was reduced by about 0.25 ppbv/d at these altitudes, representing a reduction in P(O3) of 15-20 percent over the column from 14.5 to 17 km. P(O3) from 12 to 17 km in a region distant from active convection and subject to lightning was 2-3 times higher than it would have been without lightning.

It is shown that rapid vertical transport of air from urban plumes through deep convective clouds can cause substantial enhancement of the rate of O3 production in the free troposphere. Simulation of convective redistribution and subsequent photochemistry of an urban plume from Oklahoma City during the 1985 PRESTORM campaign shows enhancement of O3 production in the free tropospheric cloud outflow layer by a factor of almost 4. In contrast, simulation of convective transport of an urban plume from Manaus, Brazil, into a prestine free troposphere during GTE/ABLE 2B (1987), followed by a photochemical simulation, showed enhancement of O3 production by a factor of 35. The reasons for the different enhancements are (1) intensity of cloud vertical motion; (2) initial boundary layer O3 precursor concentrations; and (3) initial amount of background free tropospheric NO(x). Convective transport of ozone precursors to the middle and upper troposphere allows the resulting O3 to spread over large geographic regions, rather than being confined to the lower troposphere where loss processes are much more rapid. Conversely, as air with lower NO descends and replaces more polluted air, there is greater O3 production efficiency per molecule of NO in the boundary layer following convective transport. As a result, over 30 percent more ozone could be produced in the entire tropospheric column in the first 24 hours following convective transport of urban plumes.

A regional-scale estimate of the fraction of biomass burning emissions that are transported to the free troposphere by deep convection is presented. The focus is on CO and the study region is a part of Brazil that underwent intensive deforestation in the 1980s. The method of calculation is stepwise, scaling up from a prototype convective event, the dynamics of which are well-characterized, to the vertical mass flux of carbon monoxide over the region. Given uncertainties in CO emissions from biomass burning and the representativeness of the prototype event, it is estimated that 10-40 percent of CO emissions from the burning region may be rapidly transported to the free troposphere over the burning region. These relatively fresh emissions will produce O3 efficiently in the free troposphere where O3 has a longer lifetime than in the boundary layer.

The present 2D nonhydrostatic, time-dependent numerical cloud model, which is applied to the structure and development of a subtropical squall line observed during the Taiwan Area Mesoscale Experiment, includes (1) a parameterized ice-phase microphysical scheme, (2) longwave and shortwave radiative-transfer processes, and (3) heat and moisture fluxes from the ocean surface. A comparative study is made of the vertical transport of horizontal momentum and of the latent heat release by the simulated subtropical squall system, as well as the heat release of squall systems occurring in other geographic locations.

The performance of several ice parameterizations has been evaluated through a numerical cloud model. Ice effects using different schemes are contrasted with each other and with an ice-free control by incorporating them into the cloud model and by applying them to simulations of tropical squall systems. The latter are simulated in 2D so that a large domain can be used to incorporate a complete anvil. Nonsquall-type convective lines are simulated in 3D owing to their smaller horizontal scale. It is concluded that inclusion of ice microphysics in the cloud model enhanced the agreement of the simulated convection with some features of observed convection, including the proportion of surface rainfall in the anvil region and the intensity and structure of the radar brightband near the melting level in the anvil. In the experiments with bulk microphysics, three ice categories produced much better results than two ice categories, which in turn was better than no ice. For the tropical squall-type and nonsquall-type systems the optimal mix was ice, snow, and graupel.

A three-dinensional cloud model-microwave radiative transfer model combination is used to study the relations among the precipitation and other microphysical characteristics of a tropical oceanic squall line and the upwelling radiance at pertinent microwave frequencies. Complex brightness temperature-rain rate relations are evident at the full horizontal resolution (1.5 km) of the models, with spatial averaging producing smoother, shifter relations in most cases. Nonprecipitating cloud water is shown to be important in understanding the resulting distribution of brightness temperature. At the mature stage, convective portions of the cloud system are shown to produce different brightness temperature relations than the stratiform portion, primarily related to the distribution of cloud water. The evolution of the convective system from a small convective complex through its mature stage and the beginning of its dissipation also is shown to result in a variation of brightness temperature-rain relations, related to the distribution of cloud water and the evolution of ice in the precipitating system. The results of the study point to the need to take into account the evolution of nonprecipitating cloud water and precipitation-sized ice in the retrieval of rain from passive microwave space observations. This effect is evident for both the life cycle of individual convective elements and the life cycle of the convective system as a whole.

The role of convection was examined in trace gas transport and ozone production in a tropical dry season squall line sampled on August 3, 1985, during NASA Global Tropospheric Experiment/Amazon Boundary Layer Experiment 2A (NASA GTE/ABLE 2A) in Amazonia, Brazil. Two types of analyses were performed. Transient effects within the cloud are examined with a combination of two-dimensional cloud and one-dimensional photochemical modeling. Tracer analyses using the cloud model wind fields yield a series of cross sections of NO(x), CO, and O3 distribution during the lifetime of the cloud; these fields are used in the photochemical model to compute the net rate of O3 production. At noon, when the cloud was mature, the instantaneous ozone production potential in the cloud is between 50 and 60 percent less than in no-cloud conditions due to reduced photolysis and cloud scavenging of radicals. Analysis of cloud inflows and outflows is used to differentiate between air that is undisturbed and air that has been modified by the storm. These profiles are used in the photochemical model to examine the aftereffects of convective redistribution in the 24-hour period following the storm. Total tropospheric column O3 production changed little due to convection because so little NO(x) was available in the lower troposphere. However, the integrated O3 production potential in the 5- to 13-km layer changed from net destruction to net production as a result of the convection. The conditions of the August 3, 1985, event may be typical of the early part of the dry season in Amazonia, when only minimal amounts of pollution from biomass burning have been transported into the region.

The basic features of a new and improved bulk-microphysical parameterization capable of simulating the hydrometeor structure of convective systems in all types of large-scale environments (with minimal adjustment of coefficients) are studied. Reflectivities simulated from the model are compared with radar observations of an intense midlatitude convective system. Simulated reflectivities using the novel four-class ice scheme with a microphysical parameterization rain distribution at 105 min are illustrated. Preliminary results indicate that this new ice scheme works efficiently in simulating midlatitude continental storms.

A simple algorithm to estimate the latent heating of cloud systems from their vertical hydrometeor profiles is proposed. The derivation as well as the validation of the algorithm is based on output generated by a nonhydrostatic cloud model with parameterized microphysical processes. Mature and decaying stages of a GATE squall-type convective system have been tested. The algorithm-derived heating budget is in reasonable agreement with the budget predicted by the cloud model. The input to the proposed algoritm can be obtained from either a rain retrieval technique based on information from multichannel passive microwave signals or a kinematic cloud model based on information from Doppler radar wind fields and radar reflectivity patterns. Such an application would have significant implications for spaceborne remote sensing and the large-scale weather prediction data assimilation problem.

During the second Amazon Boundary Layer Experiment (ABLE 2B), meteorological observations, chemical measurements, and model simulations are utilized in order to interpret convective cloud draft structure and to analyze its role in transport and vertical distribution of trace gases. One-dimensional photochemical model results suggest that the observed poststorm changes in ozone concentration can be attributed to convective transports rather than photochemical production and the results of a two-dimensional time-dependent cloud model simulation are presented for the May 6, 1987 squall system. The mesoscale convective system exhibited evidence of significant midlevel detrainment in addition to transports to anvil heights. Chemical measurements of O3 and CO obtained in the convective environment are used to predict photochemical production within the troposphere and to corroborate the cloud model results.

The relationship between the rain rate and the brightness temperature (Tb) was investigated using a cloud model/microwave radiative transfer model combination to obtain the rain-rate/Tb relations for four different frequencies: 10, 19, 37, and 86 GHz. The results at 19, 37, and 86 GHz were found to be significantly affected by ice in the modeled convective system, while the results at 10 GHz showed very little effect. Nonprecipitating cloud water was found to affect Tb in two ways. First, at low rain rates, the presence of significant cloud water produced higher Tb values than in cases with little cloud water. The second effects occurs at 19, 37, and 86 GHz at higher rainrates associated with significant ice formation; the scattering by ice lowered the Tb.