Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cloud parameterizations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

The Influence of Shallow Cloud Populations on Transitions to Deep Convection in the Amazon

Abstract In this study, a pair of convection-permitting (2-km grid spacing), month-long, wet-season Weather Research and Forecasting (WRF) Model simulations with and without the eddy-diffusivity mass-flux (EDMF) scheme are performed for a portion of the Green Ocean Amazon (GoAmazon) 2014/15 field campaign period. EDMF produces an ensemble of subgrid-scale convective plumes that evolve in response to the boundary layer meteorological conditions and can develop into shallow clouds. The objective of this study is to determine how different treatments of shallow cumulus clouds (i.e., with and without EDMF) impact the total cloud population and precipitation across the Amazonian rain forest, with emphasis on impacts on the likelihood of shallow-to-deep convection transitions. Results indicate that the large-scale synoptic conditions in the EDMF and control simulations are nearly identical; however, on the local scale their rainfall patterns diverge drastically and the biases decrease in EDMF. The EDMF scheme significantly increases the frequency of shallow clouds, but the frequencies of deep clouds are similar between the simulations. Deep convective clouds are tracked using a cloud-tracking algorithm to examine the impact of shallow cumulus on the surrounding ambient environment where deep convective clouds initiate. Results suggest that a rapid increase of low-level cloudiness acts to cool and moisten the low to midtroposphere during the day, favoring the transition to deep convection.

54 ENVIRONMENTAL SCIENCES↗

Sensitivity of a global climate model to the specification of convective updraft and downdraft mass fluxes

The response of the GISS global climate model to different parameterizations of moist convective mass flux is studied. A control run with arbitrarily specified updraft mass flux is compared to experiments predicting cumulus mass fulx on the basis of low-level convergence, convergence plus surface evaporation, or convergence and evaporation modified by varying boundary layer height. Also, an experiment that includes a simple parameterization of saturated convective-scale downdrafts is discussed. It is found that the model correctly simulates the correlation between deep convection strength and tropical sea surface temperature in each experiment with the parameterization of cumulus mass flux having little effect. The implications of the experiments for cloud effects on climate sensitivity are examined.

Del Genio, Anthony D.↗

Dynamical Conditions of Ice Supersaturation and Ice Nucleation in Convective Systems: A Comparative Analysis Between in Situ Aircraft Observations and WRF Simulations

Occurrence frequency and dynamical conditions of ice supersaturation (ISS, where relative humidity with respect to ice (RHi) greater than 100%) are examined in the upper troposphere around convective activity. Comparisons are conducted between in situ airborne observations and the Weather Research and Forecasting model simulations using four double-moment microphysical schemes at temperatures less than or or equal to -40degdegC. All four schemes capture both clear-sky and in-cloud ISS conditions. However, the clear-sky (in-cloud) ISS conditions are completely (significantly) limited to the RHi thresholds of the Cooper parameterization. In all of the simulations, ISS occurrence frequencies are higher by approximately 3-4 orders of magnitude at higher updraft speeds (greater than 1 m s(exp -1) than those at the lower updraft speeds when ice water content (IWC) greater than 0.01 gm(exp -3), while observations show smaller differences up to approximately 1-2 orders of magnitude. The simulated ISS also occurs less frequently at weaker updrafts and downdrafts than observed. These results indicate that the simulations have a greater dependence on stronger updrafts to maintain/generate ISS at higher IWC. At lower IWC (less than or equal or 0.01 gm(exp -3), simulations unexpectedly show lower ISS frequencies at stronger updrafts. Overall, the Thompson aerosol-aware scheme has the closest magnitudes and frequencies of ISS greater than 20% to the observations, and the modified Morrison has the closest correlations between ISS frequencies and vertical velocity at higher IWC and number density. The Cooper parameterization often generates excessive ice crystals and therefore suppresses the frequency and magnitude of ISS, indicating that it should be initiated at higher ISS (e.g.,lees than or equal to 25%).

D’Alessandro, John J.↗

Understanding the Microphysical Control and Spatial-Temporal Variability of Warm Rain Probability Using CloudSat and MODIS Observations

By combining measurements from MODIS and the CloudSat radar, we develop a parameterization scheme to quantify the combined microphysical controls by liquid water path (LWP) and cloud droplet number concentration (CDNC) of the probability of precipitation (PoP) in marine low cloud over tropical oceans. We demonstrate that the spatial-temporal variation of grid-mean in-cloud can be largely explained by the variation of the joint probability density function of LWP and CDNC in the phase space specified by the bivariate PoP (LWP and CDNC) function. Through a series of sensitivity tests guided by this understanding, we find that in the Southeastern Pacific and Atlantic the stratocumulus to cumulus transition of the is mainly due to the variation of CDNC while the annual cycle is mainly due to the variation of LWP. The results of this study provide a viable way to diagnose the root cause of warm rain problems in global climate models.

Zhibo Zhang↗

Electrification Within Wintertime Stratiform Regions Sampled During the 2020/2022 NASA IMPACTS Field Campaign

Two nor'easter events—sampled during the NASA Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) field campaign—were examined to characterize the microphysics in relation to the underlying electrification processes within wintertime stratiform regions. A theoretical model was developed to determine whether accretion or diffusion growth regimes were preferential during periods of greatest electrification. Model simulation with electrification parameterization was used to provide supplemental context to the physical processes of in-cloud microphysics and electrification. The strongest electric fields (i.e., ∼80 V m −1 at 20 km) during the 2020 NASA IMPACTS deployment was associated with large non-rimed ice crystals colliding with each other. During the 29–30 January 2022 science flight, the NASA P-3 microphysical probe data demonstrated that non-inductive charging was possible off the coastline of Cape Cod, Massachusetts. Later in the science flight, when the NASA P-3 and ER-2 were coordinating with each other, measured electric fields consistently were less than 8 V m −1 and electrification was subdued owing to reduced concentrations of graupel and large ice hydrometeors. Altogether, the in-situ observations provide evidence for the non-riming collisional charging mechanism and demonstrates that graupel and supercooled liquid water may not be necessary for weak electrification within wintertime stratiform regions. Model output from simulation of both events suggested that the main synoptic snowbands were associated with elevated hydrometeor snow charge density and electric fields.

snowfall↗

Interconnection of Aerosol‐Cloud Interactions and Cloud Feedback Through Warm Rain Process

Recent research has revealed a correlation within the Coupled Model Intercomparison Project Phase 6 (CMIP6) simulations: models exhibiting more pronounced warming due to positive cloud feedback simultaneously show enhanced cooling from aerosol‐cloud interactions (ACI). However, the interplay between ACI and cloud feedback is not well understood in these models. Our study addresses this gap by modifying the autoconversion rate in two Earth system models (ESMs), elucidating how ACI could connect with cloud feedback through the warm rain process. We find that higher autoconversion rates, which are associated with stronger ACI, suppress the increase in cloud liquid water path (LWP) and cloud optical depth with warming by enhancing the precipitation efficiency over the extratropical regions, leading to a larger positive cloud feedback. This study offers new insights into the compensatory mechanism between ACI and cloud feedback through the warm rain process and highlights the importance of constraining the autoconversion parameterization in models.

54 ENVIRONMENTAL SCIENCES↗

Improving the Simulation of Tropical Convective Cloud-Top Heights in CAM5 with CloudSat Observations

Using 4 years of CloudSat data, the simulation of tropical convective cloud-top heights (CCTH) above 6 km simulated by the convection scheme in the Community Atmosphere Model, version 5 (CAM5), is evaluated. Compared to CloudSat observations, CAM5 underestimates CCTH by more than 2 km on average. Further analysis of model results suggests that the dilute CAPE calculation, which has been incorporated into the convective parameterization since CAM4, is a main factor restricting CCTH to much lower levels. After removing this restriction, more convective clouds develop into higher altitudes, although convective clouds with tops above 12 km are still underestimated significantly. The environmental conditions under which convection develops in CAM5 are compared with CloudSat observations for convection with similar CCTHs. It is shown that the model atmosphere is much more unstable compared to CloudSat observations, and there is too much entrainment in CAM5. Since CCTHs are closely associated with cloud radiative forcing, the impacts of CCTH on model simulation are further investigated. Results show that the change of CCTH has important impacts on cloud radiative forcing and precipitation. With increased CCTHs, there is more cloud radiative forcing in tropical Africa and the eastern Pacific, but less cloud radiative forcing in the western Pacific. The contribution to total convective precipitation from convection with cloud tops above 9 km is also increased substantially.

54 ENVIRONMENTAL SCIENCES↗

A riming‐dependent parameterization of scattering by snowflakes using the self‐similar Rayleigh–Gans approximation

Abstract Riming is a key process of precipitation formation in ice‐containing clouds, but quantifying riming from observations is challenging, limiting our ability to evaluate the riming process in numerical weather models. One challenge for radar observations is that riming changes both the physical properties (mass, area cross‐section) and scattering properties of ice particles. These changes need to be implemented consistently as a function of riming in radar forward operators, which are required for retrievals and model evaluation in observation space. In this study, mass–size, cross‐section area–size, and backscattering cross‐section relations are developed as a function of the normalized rime mass for aggregates composed of various monomer types (columns, dendrites, needles, plates, and rosettes). The proposed framework allows us to simulate scattering properties of aggregated ice particles consistently as a function of riming in retrievals and radar forward operators. The parameterizations are developed from a large data set of simulated rimed aggregates of different sizes and monomer crystal types. The backscattering cross‐section parameterization (the “riming‐dependent parameterization”) is evaluated for radar frequencies of 35.6 and 94.0 GHz and is based on the Self‐Similar Rayleigh–Gans approximation (SSRGA), which is increasingly used to calculate microwave scattering of ice crystals and snowflakes. Compared with parameterizations from the literature that do not consider riming, the riming‐dependent parameterization leads to significantly smaller biases in terms of backscattering cross‐section. When using the particle masses and scattering properties of the individual particles simulated by the aggregation and riming model as a reference, the bias of our parameterization is below 1 dB when integrating over an exponential particle size distribution with sizes from 0.1–10 mm.

54 ENVIRONMENTAL SCIENCES↗

Sensitivity of Simulated Deep Convection to a Stochastic Ice Microphysics Framework

Ice microphysics parameterizations in models must make major simplifications relative to observations, typically employing empirical relationships to represent average functional properties of particles. However, previous studies have established that ice particle properties vary even in similar cloud types and thermodynamic environments, and it remains unclear how this so-called “natural variability” impacts simulated deep convection. This uncertainty is addressed by implementing a stochastic framework into the Predicted Particle Properties microphysics scheme in the Weather Research and Forecasting model. The approach stochastically varies the coefficients of the mass-size (m-D) relationship (m=aD b ) for unrimed and partially rimed ice. Using guidance from aircraft in situ measurements obtained during the Midlatitude Continental Convective Clouds Experiment (MC3E), the scheme samples from distributions of the prefactor (a) and the exponent (b) of the m-D relationship. Simulations of two MC3E deep convective cases indicate that the stochastic m-D scheme produces considerable variability of anvil cirrus cloud optical depth (τ) distributions, even for the same ice water path (IWP). Thus, the stochastic scheme produces variable cloud radiative forcing that is independent of IWP. This τ-IWP relationship variability is nonexistent using the deterministic m-D ensemble. Additional sensitivity tests are performed in which the fallspeed-size relationship (V=cD d ) is stochastically varied, resulting in variable precipitation amounts and rain rate distributions. Results are presented in the context of satellite and precipitation observations and include comparison with other ensemble configurations using perturbed initial and lateral boundary conditions and small-amplitude noise added to the potential temperature field.

54 ENVIRONMENTAL SCIENCES↗

Separating physical impacts from natural variability using piggybacking technique

Abstract. In a chaotic system, like moist convection, it is difficult to separate the impact of a physical process from effects of natural variability. This is because modifying even a small element of the system physics typically leads to a different system evolution and it is difficult to tell whether the difference comes from the physical impact or it merely represents a different flow realization. This paper discusses a relatively simple and computationally efficient modelling methodology that allows separation of the two. The methodology is referred to as the piggybacking approach. The idea is to use two sets of thermodynamic variables (the temperature, water vapor, and all aerosol, cloud, and precipitation variables) in a single cloud simulation. The two sets differ in a specific element of the physics, such as aerosol properties, microphysics parameterization, large-scale forcing, environmental profiles, etc. One thermodynamic set is coupled to the dynamics and drives the simulated flow, and the other set piggybacks the flow, that is, thermodynamic variables are carried by the flow but they do not affect it. By switching the two sets (i.e. the set driving the simulation becomes the piggybacking one, and vice versa), the impact on the cloud dynamics can be evaluated. This paper provides details of the method and reviews results of its application to such problems as the postulated deep convection invigoration in polluted environments, the impact of changes in environmental profiles (e.g., due to climate change) on convective dynamics, and the role of cloud-layer heterogeneities for shallow convective cloud field evolution. Prospects for applying piggybacking technique to other areas of atmospheric simulation (e.g., weather prediction or geoengineering) are also mentioned.

Grabowski, Wojciech W.↗

Quantification of marine aerosol subgrid variability and its correlation with clouds based on high‐resolution regional modeling

Abstract One limitation of most global climate models (GCMs) is that with the horizontal resolutions they typically employ, they cannot resolve the subgrid variability (SGV) of clouds and aerosols, adding extra uncertainties to the aerosol radiative forcing estimation. To inform the development of an aerosol subgrid variability parameterization, here we analyze the aerosol SGV over the southern Pacific Ocean simulated by the high‐resolution Weather Research and Forecasting model coupled to Chemistry. We find that within a typical GCM grid, the aerosol mass subgrid standard deviation is 15% of the grid‐box mean mass near the surface on a 1 month mean basis. The fraction can increase to 50% in the free troposphere. The relationships between the sea‐salt mass concentration, meteorological variables, and sea‐salt emission rate are investigated in both the clear and cloudy portion. Under clear‐sky conditions, marine aerosol subgrid standard deviation is highly correlated with the standard deviations of vertical velocity, cloud water mixing ratio, and sea‐salt emission rates near the surface. It is also strongly connected to the grid box mean aerosol in the free troposphere (between 2 km and 4 km). In the cloudy area, interstitial sea‐salt aerosol mass concentrations are smaller, but higher correlation is found between the subgrid standard deviations of aerosol mass and vertical velocity. Additionally, we find that decreasing the model grid resolution can reduce the marine aerosol SGV but strengthen the correlations between the aerosol SGV and the total water mixing ratio (sum of water vapor, cloud liquid, and cloud ice mixing ratios).

Lin, Guangxing↗

Comparisons of Ice Water Path in Deep Convective Systems Among Ground-Based, GOES, and CERES-MODIS Retrievals

Retrievals of convective cloud microphysical properties based on passive satellite imagery are difficult. To help quantify their uncertainties, ice water paths (IWPs) retrieved from the NASA Clouds and the Earth's Radiant Energy System project using Geostationary Operational Environmental Satellite (GOES) and Terra/Aqua MODerate-resolution Imaging Spectroradiometer observations are compared with IWPs retrieved from Next-Generation Radar (NEXRAD) observations over a large domain (32°N to 40°N and 105°W to 91°W) during the 2011 Midlatitude Continental Convective Clouds Experiment field campaign. Based on comparisons of pixel-level (4 km × 4 km) daytime IWP retrievals from NEXRAD and GOES, it is found that NEXRAD- and GOES-retrieved mean IWPs are 2.03 and 1.83 kg m -2 , respectively, for ice-phase cloud in thick anvil area. Their mean difference of 0.20 kg m -2 (with 95% confidence interval: 0.14-0.26 kg m -2 ) is within the uncertainty of NEXRAD retrievals. However, the low correlation between pixel-to-pixel comparisons indicates a large variation in GOES-retrieved IWP. For mixed-phase clouds in thick anvil areas, in addition to IWPs, total water paths (TWPs, sum of ice and liquid water path) are estimated with aid of aircraft measurements for NEXRAD retrievals and corrected using a TWP parameterization for GOES retrievals. The mean values of estimated TWPs from NEXRAD (corrected using aircraft in situ measurements) and GOES are similar. GOES and Clouds and the Earth's Radiant Energy System-MODerate-resolution Imaging Spectroradiometer-retrieved IWPs/TWPs generally do not exceed 5 kg m - 2. Large differences and low correlations exist between satellite and NEXRAD retrievals in stratiform rain areas. Possible reasons for the differences between retrievals are discussed.

54 ENVIRONMENTAL SCIENCES↗

The Impacts of Microphysics and Planetary Boundary Layer Physics on Model Simulations of U. S. Deep South Summer Convection

Inspection of output from various configurations of high-resolution, explicit convection forecast models such as the Weather Research and Forecasting (WRF) model indicates significant sensitivity to the choices of model physics parameterizations employed. Some of the largest apparent sensitivities are related to the specifications of the cloud microphysics and planetary boundary layer physics packages. In addition, these sensitivities appear to be especially pronounced for the weakly-sheared, multicell modes of deep convection characteristic of the Deep South of the United States during the boreal summer. Possible ocean-land sensitivities also argue for further examination of the impacts of using unique ocean-land surface initialization datasets provided by the NASA Short-term Prediction Research and Transition (SPoRT Center to select NOAA/NWS weather forecast offices. To obtain better quantitative understanding of these sensitivities and also to determine the utility of the ocean-land initialization data, we have executed matrices of regional WRF forecasts for selected convective events near Mobile, AL (MOB), and Houston, TX (HGX). The matrices consist of identically initialized WRF 24-h forecasts using any of eight microphysics choices and any of three planetary boundary layer choices. The resulting 24 simulations performed for each event within either the MOB or HGX regions are then compared to identify the sensitivities of various convective storm metrics to the physics choices. Particular emphasis is placed on sensitivities of precipitation timing, intensity, and coverage, as well as amount and coverage of lightning activity diagnosed from storm kinematics and graupel in the mixed phase layer. The results confirm impressions gleaned from study of the behavior of variously configured WRF runs contained in the ensembles produced each spring at the Center for the Analysis and Prediction of Storms, but with the benefit of more straightforward control of the physics package choices. The design of the experiments thus allows for more direct interpretation of the sensitivities to each possible physics combination. The results should assist forecasters in their efforts to anticipate and correct for possible biases in simulated WRF convection patterns, and help the modeling community refine their model parameterizations.

McCaul, E. W., Jr.↗

Horizontal Resolution Sensitivity of the Simple Convection-Permitting E3SM Atmosphere Model in a Doubly-Periodic Configuration

We develop a doubly periodic version of the Simple Convection-Permitting E3SM Atmosphere Model (SCREAM) to provide an efficient configuration for this global convection permitting model (GCPM), akin to a single column model often found in conventional general circulation models. The design details are explained, in addition to the extensive case library associated with the doubly periodic SCREAM (DP-SCREAM) configuration. We demonstrate that doubly periodic cloud resolving models are useful tools to explore the horizontal resolution sensitivity of GCPMs, in addition to replicating biases seen in the global models. Using DP-SCREAM, we show that SCREAM exhibits behaviors of a scale aware model as it is able to naturally partition between sub-grid scale (SGS) and resolved vertical transport across the gray zone of turbulence. We show that SCREAM is reasonably scale insensitive when run at resolutions from 1 to 5 km, but can exhibit sensitivity, particularly for the shallow convective regime, when run at resolutions approaching that of large eddy simulations. We conclude that SGS parameterization improvements are likely needed to reduce this scale sensitivity.

54 ENVIRONMENTAL SCIENCES↗

Vertical Transport Processes for Inert and Scavenged Species: TRACE-A Measurements

The TRACE-A mission of the NASA DC-8 aircraft made a large-scale survey of the tropical and subtropical atmosphere in September and October of 1992. Both In-situ measurements of CO (G. Sachsen NASA Langley) and aerosol size (J. Browell group, NASA Langley) provide excellent data sets with which to constrain vertical transport by planetary boundary layer mixing and deep-cloud cumulus convection. Lidar profiles of aerosol-induced scattering and ozone (also by Bremen) are somewhat require more subtle interpretation as tracers, but the vertical information on layering largely compensates for these complexities. The reason this DC-8 dataset is so useful is that very large areas of biomass burning over Africa and South America provide surface sources of appropriate sizes with which to characterize vertical and horizontal motions; the major limitation of our source description is that biomass burning patterns move considerably every few days, and daily burning inventories are a matter of concurrent, intensive research. We use the Penn State / NCAR MM5 model in an assimilation mode on the synoptic and intercontinental scale, and assess the success it shows in vertical transport descriptions. We find that the general level of emissions suggested by the climatological approach (Will. Has, U. of Montana) appears to be approximately correct, possibly a bit low, for this October, 1992, time period. Vertical transport in planetary boundary layer mixing to 5.5 kin was observed and reproduced in our simulations. Furthermore we find evidence that Blackader "transilient" or matrix-transport scheme is needed, but may require some adaptation in our tracer model: CO seems to exhibit very high values at the top of the planetary boundary layer, a process that stretches the eddy-diffusion parameterization. We will report on progress in improving the deep convective transport of carbon monoxide: the Grail scheme as we used it at 100 kin resolution did not transport enough material to the upper troposphere. We expect to be able to attribute this to either parameterization reasons (inadequacy of this parameterization at the large 100km scale) or other reasons. Nevertheless, the qualitative nature of deep transport by clouds shows up well in the simulations. As for scavengable species, the simulations predict tens of micrograms per standard cubic meter of smoke aerosol in the boundary layer. In a straightforward illustration of our simple bulk-mass scavenging parameterization, to one or two micrograms per standard cubic meter of smoke aerosol in the free troposphere just above the source regions: very high concentrations for the free troposphere. We expect to report on comparisons of these predictions to a variety of observations.

Chatfield, Robert B.↗

Quantifying the Observational Requirements of a Space-borne LiDAR Snow Mission

This study quantifies the level of observational accuracy required from a spaceborne light detection and ranging (LiDAR) snow depth retrieval mission for enabling beneficial impacts for snow estimation. The study is conducted over a region in Western Colorado using a suite of observing system simulation experiments(OSSEs).The Joint UK Land Environment Simulator, version 5.0 (JULES v5.0) is employed to simulate a suite of idealized LiDAR observations, considering a range of LiDAR snow depth retrieval errors, different hypothetical sensor swath widths, and the impact of cloud cover on observability. These simulated observations are then assimilated into the Noah land surface modelwith multi-parameterization options, version 3.6 (Noah-MP v3.6) model. This data assimilation setup is used to systematically evaluate the potential utility of LiDAR observations for improving modeled snow water equivalent (SWE)estimates and water budget variables such as runoff. Results from the OSSE runsshow that, in general, assimilation of synthetic LiDAR observations provide beneficial impacts when theLiDAR snow depth retrieval error standard deviation (σerror) is below 60 cm.Based on comparisons between the realistic (i.e., swath-limited and cloud-attenuated) case and the idealized (i.e., infinite swath width in the absence of cloud cover) case,this study concludes that observations with a conservative error standard deviation threshold of 40 cm (i.e., upper limit of the snow depth retrieval error that adds value to the SWE estimates viaassimilation) are needed for improving modeled snow estimates. More than a 33% reduction in SWEroot mean square errors and more than a 15% increase in correlation coefficientsare achieved when σ error ≤ 40 cmusing a 170-km sensor swath width in the presence of cloud attenuation effects. Further, the integrated hydrologic response, as represented by total (surface and subsurface) runoff estimates during the snow ablation season, are also enhanced when assimilating synthetic LiDAR snow depth retrievals with errors below this level.

Yonghwan Kwon↗

Structurally flexible cloud microphysics, observationally constrained at all scales via ML-accelerated Bayesian inference

We discuss the challenge of developing observationally informed parameterizations of microphysics for use at a hierarchy of modeling scales. Our proposed approach is applicable to any domain that suffers from a two-fold parameterization problem, where physical processes are not resolved at the model scale (the first problem), and those processes are uncertain at any scale (the second problem). For such problems, a physical approach facilitates modeling across scales, as well as systematic observational inference accelerated by machine learning (ML) surrogate models, and quantification of physical uncertainties.

54 ENVIRONMENTAL SCIENCES↗

A statistical analysis of the dependency of closure assumptions in cumulus parameterization on the horizontal resolution

Simulated data from the UCLA cumulus ensemble model are used to investigate the quasi-universal validity of closure assumptions used in existing cumulus parameterizations. A closure assumption is quasi-universally valid if it is sensitive neither to convective cloud regimes nor to horizontal resolutions of large-scale/mesoscale models. The dependency of three types of closure assumptions, as classified by Arakawa and Chen, on the horizontal resolution is addressed in this study. Type I is the constraint on the coupling of the time tendencies of large-scale temperature and water vapor mixing ratio. Type II is the constraint on the coupling of cumulus heating and cumulus drying. Type III is a direct constraint on the intensity of a cumulus ensemble. The macroscopic behavior of simulated cumulus convection is first compared with the observed behavior in view of Type I and Type II closure assumptions using 'quick-look' and canonical correlation analyses. It is found that they are statistically similar to each other. The three types of closure assumptions are further examined with simulated data averaged over selected subdomain sizes ranging from 64 to 512 km. It is found that the dependency of Type I and Type II closure assumptions on the horizontal resolution is very weak and that Type III closure assumption is somewhat dependent upon the horizontal resolution. The influences of convective and mesoscale processes on the closure assumptions are also addressed by comparing the structures of canonical components with the corresponding vertical profiles in the convective and stratiform regions of cumulus ensembles analyzed directly from simulated data. The implication of these results for cumulus parameterization is discussed.

Xu, Kuan-Man↗