Engineering topics
Xu, Kuan-Man
Publications and source records attributed to Xu, Kuan-Man.
Changes in Tropical Clouds and Atmospheric Circulation Associated with Rapid Adjustment Induced by Increased Atmospheric CO2 – A Multiscale Modeling Framework Study
The radiative heating increase due to increased CO2 concentration is the primary source for the rapid adjustment of atmospheric circulation and clouds. In this study, we investigate the rapid adjustment resulting from doubling of CO2 and its physical mechanism using a multiscale modeling framework (MMF). The MMF includes an advanced higher-order turbulence closure in its cloud-resolving model component and simulates realistic shallow and deep cloud climatology and boundary layer turbulence. The rapid adjustment over the tropics is characterized by 1) reduced ascent and descent strengths over the ocean, 2) increased lower tropospheric stability (LTS) over the subsidence region, 3) shoaling of planetary boundary layers over the ocean, 4) increased deep convection over lands and shift of cloud coverage from the ocean to lands, and 5) reduced sensible (SH) and latent heat (LH) fluxes over the oceanic deep convective regions. Unlike conventional general circulation models and another MMF, a reduction in the global-mean shortwave cloud radiative cooling is not simulated, due to the increase in low clouds at lower altitudes over the ocean, resulting from reduced cloud-top entrainment due to strengthened inversion. Changes in regional circulation play a key role in cloud changes and shift of cloud coverage to lands. Weaker energy transport resulting from water vapor and cloud CO2 masking effects reduces the upward motion and convective clouds in the oceanic regions. The ocean-land transports are linked to the partitioning of surface SH and LH fluxes that increases humidity over lands and enhances deep convection over the tropical lands.
Small vs. Large Convective Cloud Objects from CERES Aqua Observations: Where are the Intraseasonal Variation Signals?
During inactive phases of Madden-Julian oscillation (MJO), there are plenty of deep but small convective systems and far fewer deep and large ones. During active phases of MJO, a manifestation of an increase in the occurrence of large and deep cloud clusters results from an amplification of large-scale motions by stronger convective heating. This study is designed to quantitatively examine the roles of small and large cloud clusters during the MJO life cycle. We analyze the cloud object data from Aqua CERES observations for tropical deep convective (DC) and cirrostratus (CS) cloud object types according to the real-time multivariate MJO index. The cloud object is a contiguous region of the earth with a single dominant cloud-system type. The size distributions, defined as the footprint numbers as a function of cloud object diameters, for particular MJO phases depart greatly from the combined (8-phase) distribution at large cloud-object diameters due to the reduced/increased numbers of cloud objects related to changes in the large-scale environments. The medium diameter corresponding to the combined distribution is determined and used to partition all cloud objects into "small" and "large" groups of a particular phase. The two groups corresponding to the combined distribution have nearly equal numbers of footprints. The medium diameters are 502 km for DC and 310 km for cirrostratus. The range of the variation between two extreme phases (typically, the most active and depressed phases) for the small group is 6-11% in terms of the numbers of cloud objects and the total footprint numbers. The corresponding range for the large group is 19-44%. In terms of the probability density functions of radiative and cloud physical properties, there are virtually no differences between the MJO phases for the small group, but there are significant differences for the large groups for both DC and CS types. These results suggest that the intreseasonal variation signals reside at the large cloud clusters while the small cloud clusters represent the background noises resulting from various types of the tropical waves with different wavenumbers and propagation directions/speeds.
Application of the CERES Flux-by-Cloud Type Simulator to GCM Output
The CERES Flux By CloudType data product produces CERES top-of-atmosphere (TOA) fluxes by region and cloud type. Here, the cloud types are defined by cloud optical depth (t) and cloud top pressure (pc), with bins similar to those used by ISCCP (International Satellite Cloud Climatology Project). This data product has the potential to be a powerful tool for the evaluation of the clouds produced by climate models by helping to identify which physical parameterizations have problems (e.g., boundary-layer parameterizations, convective clouds, processes that affect surface albedo). Also, when the flux-by-cloud type and frequency of cloud types are simultaneously used to evaluate a model, the results can determine whether an unrealistically large or small occurrence of a given cloud type has an important radiative impact for a given region. A simulator of the flux-by-cloud type product has been applied to three-hourly data from the year 2008 from the UK Met Office HadGEM2-A model using the Langley Fu-Lour radiative transfer model to obtain TOA SW and LW fluxes.
Covariance Between Arctic Sea Ice and Clouds Within Atmospheric State Regimes at the Satellite Footprint Level
Understanding the cloud response to sea ice change is necessary for modeling Arctic climate. Previous work has primarily addressed this problem from the interannual variability perspective. This paper provides a refined perspective of sea ice-cloud relationship in the Arctic using a satellite footprint-level quantification of the covariance between sea ice and Arctic low cloud properties from NASA A-Train active remote sensing data. The covariances between Arctic low cloud properties and sea ice concentration are quantified by first partitioning each footprint into four atmospheric regimes defined using thresholds of lower tropospheric stability and mid-tropospheric vertical velocity. Significant regional variability in the cloud properties is found within the atmospheric regimes indicating that the regimes do not completely account for the influence of meteorology. Regional anomalies are used to account for the remaining meteorological influence on clouds. After accounting for meteorological regime and regional influences, a statistically significant but weak covariance between cloud properties and sea ice is found in each season for at least one atmospheric regime. Smaller average cloud fraction and liquid water are found within footprints with more sea ice. The largest-magnitude cloud-sea ice covariance occurs between 500m and 1.2 km when the lower tropospheric stability is between 16 and 24 K. The covariance between low cloud properties and sea ice is found to be largest in fall and is accompanied by significant changes in boundary layer temperature structure where larger average near-surface static stability is found at larger sea ice concentrations.
How Difficult is it to Reduce Low-Level Cloud Biases With the Higher-Order Turbulence Closure Approach in Climate Models?
Low-level clouds cover nearly half of the Earth and play a critical role in regulating the energy and hydrological cycle. Despite the fact that a great effort has been put to advance the modeling and observational capability in recent years, low-level clouds remains one of the largest uncertainties in the projection of future climate change. Low-level cloud feedbacks dominate the uncertainty in the total cloud feedback in climate sensitivity and projection studies. These clouds are notoriously difficult to simulate in climate models due to its complicated interactions with aerosols, cloud microphysics, boundary-layer turbulence and cloud dynamics. The biases in both low cloud coverage/water content and cloud radiative effects (CREs) remain large. A simultaneous reduction in both cloud and CRE biases remains elusive. This presentation first reviews the effort of implementing the higher-order turbulence closure (HOC) approach to representing subgrid-scale turbulence and low-level cloud processes in climate models. There are two HOCs that have been implemented in climate models. They differ in how many three-order moments are used. The CLUBB are implemented in both CAM5 and GDFL models, which are compared with IPHOC that is implemented in CAM5 by our group. IPHOC uses three third-order moments while CLUBB only uses one third-order moment while both use a joint double-Gaussian distribution to represent the subgrid-scale variability. Despite that HOC is more physically consistent and produces more realistic low-cloud geographic distributions and transitions between cumulus and stratocumulus regimes, GCMs with traditional cloud parameterizations outperform in CREs because tuning of this type of models is more extensively performed than those with HOCs. We perform several tuning experiments with CAM5 implemented with IPHOC in an attempt to produce the nearly balanced global radiative budgets without deteriorating the low-cloud simulation. One of the issues in CAM5-IPHOC is that cloud water content is much higher than in CAM5, which is combined with higher low-cloud coverage to produce larger shortwave CREs in some low-cloud prevailing regions. Thus, the cloud-radiative feedbacks are exaggerated there. The turning exercise is focused on microphysical parameters, which are also commonly used for tuning in climate models. The results will be discussed in this presentation.
Background Noises Versus Intraseasonal Variation Signals: Small vs. Large Convective Cloud Objects From CERES Aqua Observations
During inactive phases of Madden-Julian Oscillation (MJO), there are plenty of deep but small convective systems and far fewer deep and large ones. During active phases of MJO, a manifestation of an increase in the occurrence of large and deep cloud clusters results from an amplification of large-scale motions by stronger convective heating. This study is designed to quantitatively examine the roles of small and large cloud clusters during the MJO life cycle. We analyze the cloud object data from Aqua CERES (Clouds and the Earth's Radiant Energy System) observations between July 2006 and June 2010 for tropical deep convective (DC) and cirrostratus (CS) cloud object types according to the real-time multivariate MJO index, which assigns the tropics to one of the eight MJO phases each day. The cloud object is a contiguous region of the earth with a single dominant cloud-system type. The criteria for defining these cloud types are overcast footprints and cloud top pressures less than 400 hPa, but DC has higher cloud optical depths (=10) than those of CS (<10). The size distributions, defined as the footprint numbers as a function of cloud object diameters, for particular MJO phases depart greatly from the combined (8-phase) distribution at large cloud-object diameters due to the reduced/increased numbers of cloud objects related to changes in the large-scale environments. The medium diameter corresponding to the combined distribution is determined and used to partition all cloud objects into "small" and "large" groups of a particular phase. The two groups corresponding to the combined distribution have nearly equal numbers of footprints. The medium diameters are 502 km for DC and 310 km for cirrostratus. The range of the variation between two extreme phases (typically, the most active and depressed phases) for the small group is 6-11% in terms of the numbers of cloud objects and the total footprint numbers. The corresponding range for the large group is 19-44%. In terms of the probability density functions of radiative and cloud physical properties, there are virtually no differences between the MJO phases for the small group, but there are significant differences for the large groups for both DC and CS types. These results suggest that the intreseasonal variation signals reside at the large cloud clusters while the small cloud clusters represent the background noises resulting from various types of the tropical waves with different wavenumbers and propagation speeds/directions.
Cloud Feedbacks on Greenhouse Warming in a Multi-Scale Modeling Framework with a Higher-Order Turbulence Closure
Five-year simulation experiments with a multi-scale modeling Framework (MMF) with a advanced intermediately prognostic higher-order turbulence closure (IPHOC) in its cloud resolving model (CRM) component, also known as SPCAM-IPHOC (super parameterized Community Atmospheric Model), are performed to understand the fast tropical (30S-30N) cloud response to an instantaneous doubling of CO2 concentration with SST held fixed at present-day values. SPCAM-IPHOC has substantially improved the low-level representation compared with SPCAM. It is expected that the cloud responses to greenhouse warming in SPCAM-IPHOC is more realistic. The change of rising motion, surface precipitation, cloud cover, and shortwave and longwave cloud radiative forcing in SPCAM-IPHOC from the greenhouse warming will be presented in the presentation.
Recent Trends of the Tropical Hydrological Cycle Inferred from Global Precipitation Climatology Project and International Satellite Cloud Climatology Project data
Scores of modeling studies have shown that increasing greenhouse gases in the atmosphere impact the global hydrologic cycle; however, disagreements on regional scales are large, and thus the simulated trends of such impacts, even for regions as large as the tropics, remain uncertain. The present investigation attempts to examine such trends in the observations using satellite data products comprising Global Precipitation Climatology Project precipitation and International Satellite Cloud Climatology Project cloud and radiation. Specifically, evolving trends of the tropical hydrological cycle over the last 20-30 years were identified and analyzed. The results show (1) intensification of tropical precipitation in the rising regions of the Walker and Hadley circulations and weakening over the sinking regions of the associated overturning circulation; (2) poleward shift of the subtropical dry zones (up to 2deg/decade in June-July-August (JJA) in the Northern Hemisphere and 0.3-0.7deg/decade in June-July-August and September-October-November in the Southern Hemisphere) consistent with an overall broadening of the Hadley circulation; and (3) significant poleward migration (0.9-1.7deg/decade) of cloud boundaries of Hadley cell and plausible narrowing of the high cloudiness in the Intertropical Convergence Zone region in some seasons. These results support findings of some of the previous studies that showed strengthening of the tropical hydrological cycle and expansion of the Hadley cell that are potentially related to the recent global warming trends.
An Estimate of Low-Cloud Feedbacks from Variations of Cloud Radiative and Physical Properties with Sea Surface Temperature on Interannual Time Scales
Simulations of climate change have yet to reach a consensus on the sign and magnitude of the changes in physical properties of marine boundary layer clouds. In this study, the authors analyze how cloud and radiative properties vary with SST anomaly in low-cloud regions, based on five years (March 2000 - February 2005) of Clouds and the Earth s Radiant Energy System (CERES) -- Terra monthly gridded data and matched European Centre for Medium-Range Weather Forecasts (ECMWF) meteorological reanalaysis data. In particular, this study focuses on the changes in cloud radiative effect, cloud fraction, and cloud optical depth with SST anomaly. The major findings are as follows. First, the low-cloud amount (-1.9% to -3.4% /K) and the logarithm of low-cloud optical depth (-0.085 to -0.100/K) tend to decrease while the net cloud radiative effect (3.86 W/m(exp 2)/ K) becomes less negative as SST anomalies increase. These results are broadly consistent with previous observational studies. Second, after the changes in cloud and radiative properties with SST anomaly are separated into dynamic, thermodynamic, and residual components, changes in the dynamic component (taken as the vertical velocity at 700 hPa) have relatively little effect on cloud and radiative properties. However, the estimated inversion strength decreases with increasing SST, accounting for a large portion of the measured decreases in cloud fraction and cloud optical depth. The residual positive change in net cloud radiative effect (1.48 W/m(exp 2)/ K) and small changes in low-cloud amount (-0.81% to 0.22% /K) and decrease in the logarithm of optical depth (-0.035 to -0.046/ K) with SST are interpreted as a positive cloud feedback, with cloud optical depth feedback being the dominant contributor. Last, the magnitudes of the residual changes differ greatly among the six low-cloud regions examined in this study, with the largest positive feedbacks (approximately 4 W/m(exp 2)/ K) in the southeast and northeast Atlantic regions and a slightly negative feedback (-0.2 W/m(exp 2)/ K) in the south-central Pacific region. Because the retrievals of cloud optical depth and/or cloud fraction are difficult in the presence of aerosols, the transport of heavy African continental aerosols may contribute to the large magnitudes of estimated cloud feedback in the two Atlantic regions.
Variations of Cloud and Radiative Properties of Boundary-layer and Deep Convective Systems with Sea Surface Temperature Anomalies
Gridded monthly-mean satellite data contain compositing information from different cloud system types and clear-sky environments. To isolate the variations of cloud physical properties of an individual cloud system type with its environment, orbital data are needed. In this study, we will analyze the variations of cloud and radiative properties of boundary-layer clouds and deep convective cloud systems with sea surface temperature (SST) anomalies. We use Terra-CERES (Clouds and the Earth s Radiant Energy System) Level 2 data to classify distinct cloud objects defined by cloud-system types (deep convection, boundary-layer cumulus, stratocumulus and overcast clouds), sizes, geographic locations, and matched large-scale environments. This analysis method identifies a cloud object as a contiguous region of the Earth with a single dominant cloud-system type. It determines the shape and size of the cloud object from the satellite data and the cloud-system selection criteria. The statistical properties of the identified cloud objects are analyzed in terms of probability density functions (PDFs) of a single property or joint PDFs between two properties. The SST anomalies are defined as the differences from five-year annual-cycle means. Individual cloud objects are sorted into one of five equal size subsets, with the matched SST anomalies ranging from the most negative to the most positive values, for a given size category of deep convective cloud objects, boundary-layer cumulus, stratocumulus and overcast cloud objects. The PDFs of cloud and radiative properties for deep convective cloud objects (between 30 S and 30 N) are found to largely similar among the five SST anomaly subsets except for the lowest SST anomaly subset. The different characteristics from this SST anomaly subset may be related to some cloud objects resulting from equatorward movement of extratropical cloud systems. This result holds true for all three different size categories (measured by equivalent diameters of 100-150, 150-300 and > 300 km). The PDFs of cloud and radiative properties for boundary-layer cloud types are more variable for macrophysical properties than for microphysical properties, due probably to different large-scale environments where these cloud objects are formed. Implications of these results for the cloud feedback will also be discussed.
Toward Realistic Simulation of low-Level Clouds Using a Multiscale Modeling Framework With a Third-Order Turbulence Closure in its Cloud-Resolving Model Component
This study presents preliminary results from a multiscale modeling framework (MMF) with an advanced third-order turbulence closure in its cloud-resolving model (CRM) component. In the original MMF, the Community Atmosphere Model (CAM3.5) is used as the host general circulation model (GCM), and the System for Atmospheric Modeling with a first-order turbulence closure is used as the CRM for representing cloud processes in each grid box of the GCM. The results of annual and seasonal means and diurnal variability are compared between the modified and original MMFs and the CAM3.5. The global distributions of low-level cloud amounts and precipitation and the amounts of low-level clouds in the subtropics and middle-level clouds in mid-latitude storm track regions in the modified MMF show substantial improvement relative to the original MMF when both are compared to observations. Some improvements can also be seen in the diurnal variability of precipitation.
Evaluation of Cloud Physical Properties of ECMWF Analysis and Re-Analysis (ERA-40 and ERA Interim) against CERES Tropical Deep Convective Cloud Object Observations
This study presents an approach that converts the vertical profiles of grid-averaged cloud properties from large-scale models to probability density functions (pdfs) of subgrid-cell cloud physical properties measured at satellite footprints. Cloud physical and radiative properties, rather than just cloud and precipitation occurrences, of assimilated cloud systems by the European Center for Medium-range Weather Forecasts (ECMWF) operational analysis (EOA) and ECMWF Re-Analyses (ERA-40 and ERA Interim) are validated against those obtained from Earth Observing System satellite cloud object data for January-August 1998 and March 2000 periods. These properties include ice water path (IWP), cloud-top height and temperature, cloud optical depth and solar and infrared radiative fluxes. Each cloud object, a contiguous region with similar cloud physical properties, is temporally and spatially matched with EOA and ERA-40 data. Results indicate that most pdfs of EOA and ERA-40 cloud physical and radiative properties agree with those of satellite observations of the tropical deep convective cloud-object type for the January-August 1998 period. There are, however, significant discrepancies in selected ranges of the cloud property pdfs such as the upper range of EOA cloud top height. A major discrepancy is that the dependence of the pdfs on the cloud object size for both EOA and ERA-40 is not as strong as in the observations. Modifications to the cloud parameterization in ECMWF that occurred in October 1999 eliminate the clouds near the tropopause but shift power of the pdf to lower cloud-top heights and greatly reduce the ranges of IWP and cloud optical depth pdfs. These features persist in ERA-40 due to the use of the same cloud parameterizations. The downgrade of data assimilation technique and the lack of snow water content information in ERA-40, not the coarser horizontal grid resolution, are also responsible for the disagreements with observed pdfs of cloud physical properties although the detection rates of cloud object occurrence are improved for small size categories. A possible improvement to the convective parameterization is to introduce a stronger dependence of updraft penetration heights with grid-cell dynamics. These conclusions will be rechecked using the ERA Interim data, due to recent changes in the ECMWF convective parameterization (Bechtold et al. 2004, 2008). Results from the ERA Interim will be presented at the meeting.
Arctic Mixed-phase Clouds Simulated by a Cloud-Resolving Model: Comparison with ARM Observations and Sensitivity to Microphysics Parameterizations
Single-layer mixed-phase stratiform (MPS) Arctic clouds, which formed under conditions of large surface heat flux combined with general subsidence during a subperiod of the Atmospheric Radiation Measurement (ARM) Program Mixed-Phase Arctic Cloud Experiment (M-PACE), are simulated with a cloud resolving model (CRM). The CRM is implemented with either an advanced two-moment (M05) or a commonly used one-moment (L83) bulk microphysics scheme and a state-of-the-art radiative transfer scheme. The CONTROL simulation, that uses the M05 scheme and observed aerosol size distribution and ice nulei (IN) number concentration, reproduces the magnitudes and vertical structures of cloud liquid water content (LWC), total ice water content (IWC), number concentration and effective radius of cloud droplets as suggested by the M-PACE observations. It underestimates ice crystal number concentrations by an order of magnitude and overestimates effective radius of ice crystals by a factor of 2-3. The OneM experiment, that uses the L83 scheme, produces values of liquid water path (LWP) and ice plus snow water path (ISWP) that were about 30% and 4 times, respectively, of those produced by the CONTROL. Its vertical profile of IWC exhibits a bimodal distribution in contrast to the constant distribution of IWC produced in the CONTROL and observations.
Multi-Layer Arctic Mixed-Phase Clouds Simulated by a Cloud-Resolving Model: Comparison with ARM Observations and Sensitivity Experiments
A cloud-resolving model (CRM) is used to simulate the multiple-layer mixed-phase stratiform (MPS) clouds that occurred during a three-and-a-half day subperiod of the Department of Energy-Atmospheric Radiation Measurement Program s Mixed-Phase Arctic Cloud Experiment (M-PACE). The CRM is implemented with an advanced two-moment microphysics scheme, a state-of-the-art radiative transfer scheme, and a complicated third-order turbulence closure. Concurrent meteorological, aerosol, and ice nucleus measurements are used to initialize the CRM. The CRM is prescribed by time-varying large-scale advective tendencies of temperature and moisture and surface turbulent fluxes of sensible and latent heat. The CRM reproduces the occurrences of the single- and double-layer MPS clouds as revealed by the M-PACE observations. However, the simulated first cloud layer is lower and the second cloud layer thicker compared to observations. The magnitude of the simulated liquid water path agrees with that observed, but its temporal variation is more pronounced than that observed. As in an earlier study of single-layer cloud, the CRM also captures the major characteristics in the vertical distributions and temporal variations of liquid water content (LWC), total ice water content (IWC), droplet number concentration and ice crystal number concentration (nis) as suggested by the aircraft observations. However, the simulated mean values differ significantly from the observed. The magnitude of nis is especially underestimated by one order of magnitude. Sensitivity experiments suggest that the lower cloud layer is closely related to the surface fluxes of sensible and latent heat; the upper cloud layer is probably initialized by the large-scale advective cooling/moistening and maintained through the strong longwave (LW) radiative cooling near the cloud top which enhances the dynamical circulation; artificially turning off all ice-phase microphysical processes results in an increase in LWP by a factor of 3 due to interactions between the excessive LW radiative cooling and extra cloud water; heating caused by phase change of hydrometeors could affect the LWC and cloud top height by partially canceling out the LW radiative cooling. It is further shown that the resolved dynamical circulation appears to contribute more greatly to the evolution of the MPS cloud layers than the parameterized subgrid-scale circulation.
Simulation of Boundary-Layer Cumulus and Stratocumulus Clouds using a Cloud-Resolving Model With Low- and Third-Order Turbulence Closures
The effects of subgrid-scale condensation and transport become more important as the grid spacings increase from those typically used in large-eddy simulation (LES) to those typically used in cloud-resolving models (CRMs). Incorporation of these effects can be achieved by a joint probability density function approach that utilizes higher-order moments of thermodynamic and dynamic variables. This study examines how well shallow cumulus and stratocumulus clouds are simulated by two versions of a CRM that is implemented with low-order and third-order turbulence closures (LOC and TOC) when a typical CRM horizontal resolution is used and what roles the subgrid-scale and resolved-scale processes play as the horizontal grid spacing of the CRM becomes finer. Cumulus clouds were mostly produced through subgrid-scale transport processes while stratocumulus clouds were produced through both subgrid-scale and resolved-scale processes in the TOC version of the CRM when a typical CRM grid spacing is used. The LOC version of the CRM relied upon resolved-scale circulations to produce both cumulus and stratocumulus clouds, due to small subgrid-scale transports. The mean profiles of thermodynamic variables, cloud fraction and liquid water content exhibit significant differences between the two versions of the CRM, with the TOC results agreeing better with the LES than the LOC results. The characteristics, temporal evolution and mean profiles of shallow cumulus and stratocumulus clouds are weakly dependent upon the horizontal grid spacing used in the TOC CRM. However, the ratio of the subgrid-scale to resolved-scale fluxes becomes smaller as the horizontal grid spacing decreases. The subcloud-layer fluxes are mostly due to the resolved scales when a grid spacing less than or equal to 1 km is used. The overall results of the TOC simulations suggest that a 1-km grid spacing is a good choice for CRM simulation of shallow cumulus and stratocumulus.
Coincident Occurrences of Tropical Individual Cirrus Clouds and Deep Convective Systems Derived from TRMM Observations
Measurements of cloud properties and atmospheric radiation taken between January and August 1998 by the Tropical Rainfall Measuring Mission (TRMM) satellite were used to investigate the effect of spatial and temporal scales on the coincident occurrences of tropical individual cirrus clouds (ICCs) and deep convective systems (DCSs). It is found that there is little or even negative correlation between instantaneous occurrences of ICC and DCS in small areas, in which both types of clouds cannot grow and expand simultaneously. When spatial and temporal domains are increased, ICCs become more dependent on DCSs due to the origination of many ICCs from DCSs and moisture supply from the DCS in the upper troposphere for the ICCs to grow, resulting in significant positive correlation between the two types of tropical high clouds in large spatial and long temporal scales. This result may suggest that the decrease of tropical high clouds with SST from model simulations is likely caused by restricted spatial domains and limited temporal periods. Finally, the radiative feedback due to the change in tropical high cloud area coverage with sea surface temperature appears small and about -0.14 W/sq m per degree Kelvin.
Statistical Analyses of Satellite Cloud Object Data from CERES: Comparison with Cloud-Resolving Model Simulations of Tropical Convective Clouds - Part III
The present study evaluates the ability of a cloud-resolving model (CRM) to simulate the physical properties of tropical deep convective cloud objects identified from a Clouds and the Earth s Radiant Energy System (CERES) data product. The emphasis of this study is the comparisons among the small-, medium- and large-size categories of cloud objects observed during March 1998 and between the large-size categories of cloud objects observed during March 1998 (strong El Ni o) and March 2000 (weak La Ni a). Results from the CRM simulations are analyzed in a way that is consistent with the CERES retrieval algorithm and they are averaged to match the scale of the CERES satellite footprints. Cloud physical properties are analyzed in terms of their summary histograms for each category. It is found that there is a general agreement in the overall shapes of all cloud physical properties between the simulated and observed distributions. Each cloud physical property produced by the CRM also exhibits different degrees of disagreement with observations over different ranges of the property. The simulated cloud tops are generally too high and cloud top temperatures are too low except for the large-size category of March 1998. The probability densities of the simulated top-of-the-atmosphere (TOA) albedos for all four categories are underestimated for high albedos, while those of cloud optical depth are overestimated at its lowest bin. These disagreements are mainly related to uncertainties in the cloud microphysics parameterization and inputs such as cloud ice effective size to the radiation calculation. Summary histograms of cloud optical depth and TOA albedo from the CRM simulations of the large-size category of cloud objects do not differ significantly between the March 1998 and 2000 periods, consistent with the CERES observations. However, the CRM is unable to reproduce the significant differences in the observed cloud top height while it overestimates the differences in the observed outgoing longwave radiation and cloud top temperature between the two periods. Comparisons between the CRM results and the observations for most parameters in March 1998 consistently show that both the simulations and observations have larger differences between the large- and small-size categories than between the large- and medium-size, or between the medium- and small-size categories. However, the simulated cloud properties do not change as much with size as observed. These disagreements are likely related to the spatial averaging of the forcing data and the mismatch in time and in space between the numerical weather prediction model from which the forcing data are produced and the CERES observed cloud systems.