Large-eddy simulation of stratocumulus-topped boundary layer with an explicit and a new bulk microphysics scheme.

The realism of the model is evaluated by a direct comparison of the model predictions with the aircraft observations of the STBL. The first case study is based on the UKMRF flight 526 measurements collected over the North Sea on 22 July 1982; the second case study corresponds to the ASTEX flight A209 flown on 12-13 June 1992. The model is able to reproduce reasonably well most of the observed boundary layer parameters, including turbulent fluxes and variances of various fields, intensity and vertical distribution of the turbulent kinetic energy, upward and downward radiation fluxes, and the cloud drop spectra.I designed a new bulk microphysical parameterization using the explicit model as a benchmark for comparison. The liquid water is divided into two categories--non-precipitable cloud water and drizzle, similar to traditional Kessler-type parameterizations. The water content and drop concentration are predicted for each category. The source/sink terms such as autoconversion of cloud water into drizzle are deduced directly from the drop size spectra predicted by the explicit microphysical model. The predictions of the LES model using the new bulk microphysics are compared with the predictions using explicit microphysics for two cases: non-drizzling and heavy-drizzling STBL. The results show that the new bulk microphysical model satisfactorily reproduces many characteristics of the STBL as simulated by explicit microphysical model.A case of stratocumulus-to-cumulus transition triggered by the depletion of CCN is simulated. It is shown that the response of the STBL to the increase in drizzle due to CCN depletion is the reduction of its cloud fractional cover and change of the character of circulation toward the cumulus convection. The boundary layer after the Sc-to-Cu transition consists of two layers: the well-mixed cloud free surface layer driven by surface heat fluxes and shear, and the conditionally unstable upper layer capped by the inversion with embedded cumulus clouds connected to the moisture and CCN supply in the surface layer.A new LES dynamical framework coupled with an explicit microphysical module has been developed. It is verified against analytical solution (linear mountain wave test) and against predictions from the other LES models. The results of the tests of the microphysical module convincingly show that the drop spectrum resolution in our model is adequate to accurately predict the cloud microphysics parameters

Simulation of the intraseasonal variability over the Eastern Pacific ITCZ in climate models

During boreal summer, convective activity over the eastern Pacific (EPAC) inter-tropical convergence zone (ITCZ) exhibits vigorous intraseasonal variability (ISV). Previous observational studies identified two dominant ISV modes over the EPAC, i.e., a 40-day mode and a quasi-biweekly mode (QBM). The 40-day ISV mode is generally considered a local expression of the Madden-Julian Oscillation. However, in addition to the eastward propagation, northward propagation of the 40-day mode is also evident. The QBM mode bears a smaller spatial scale than the 40-day mode, and is largely characterized by northward propagation. While the ISV over the EPAC exerts significant influences on regional climate/weather systems, investigation of contemporary model capabilities in representing these ISV modes over the EPAC is limited. In this study, the model fidelity in representing these two dominant ISV modes over the EPAC is assessed by analyzing six atmospheric and three coupled general circulation models (GCMs), including one super-parameterized GCM (SPCAM) and one recently developed high-resolution GCM (GFDL HIRAM) with horizontal resolution of about 50 km. While it remains challenging for GCMs to faithfully represent these two ISV modes including their amplitude, evolution patterns, and periodicities, encouraging simulations are also noted. In general, SPCAM and HIRAM exhibit relatively superior skill in representing the two ISV modes over the EPAC. While the advantage of SPCAM is achieved through explicit representation of the cumulus process by the embedded 2-D cloud resolving models, the improved representation in HIRAM could be ascribed to the employment of a strongly entraining plume cumulus scheme, which inhibits the deep convection, and thus effectively enhances the stratiform rainfall. The sensitivity tests based on HIRAM also suggest that fine horizontal resolution could also be conducive to realistically capture the ISV over the EPAC, particularly for the QBM mode. Further analysis illustrates that the observed 40-day ISV mode over the EPAC is closely linked to the eastward propagating ISV signals from the Indian Ocean/Western Pacific, which is in agreement with the general impression that the 40-day ISV mode over the EPAC could be a local expression of the global Madden-Julian Oscillation (MJO). In contrast, the convective signals associated with the 40-day mode over the EPAC in most of the GCM simulations tend to originate between 150°E and 150°W, suggesting the 40-day ISV mode over the EPAC might be sustained without the forcing by the eastward propagating MJO. Further investigation is warranted towards improved understanding of the origin of the ISV over the EPAC

Intraseasonal Variability in a Cloud-Permitting Near-Global Equatorial Aquaplanet Model

Recent studies have suggested that the Madden-Julian oscillation is a result of an instability driven mainly by cloud-radiation feedbacks, similar in character to self-aggregation of convection in nonrotating, cloud-permitting simulations of radiative-convective equilibrium (RCE). Here we bolster that inference by simulating radiative-convective equilibrium states on a rotating sphere with constant sea surface temperature, using the cloud-permitting System for Atmospheric Modeling (SAM) with 20-km grid spacing and extending to walls at 46° latitude in each hemisphere. Mechanism-denial experiments reveal that cloud-radiation interaction is the quintessential driving mechanism of the simulated MJO-like disturbances, but wind-induced surface heat exchange (WISHE) feedbacks are the primary driver of its eastward propagation. WISHE may also explain the faster Kelvin-like modes in the simulations. These conclusions are supported by a linear stability analysis of RCE states on an equatorial beta plane. ©2018 American Meteorological Society.NSF Grant (AGS1418508)NSF Grant (AGS1418309

Large-eddy simulation of maritime deep tropical convection

This study represents an attempt to apply Large-Eddy Simulation (LES) resolution to simulate deep tropical convection in near equilibrium for 24 hours over an area of about 205 x 205 km2, which is comparable to that of a typical horizontal grid cell in a global climate model. The simulation is driven by large-scale thermodynamic tendencies derived from mean conditions during the GATE Phase III field experiment. The LES uses 2048 x 2048 x 256 grid points with horizontal grid spacing of 100 m and vertical grid spacing ranging from 50 m in the boundary layer to 100 m in the free troposphere. The simulation reaches a near equilibrium deep convection regime in 12 hours. The simulated vertical cloud distribution exhibits a trimodal vertical distribution of deep, middle and shallow clouds similar to that often observed in Tropics. A sensitivity experiment in which cold pools are suppressed by switching off the evaporation of precipitation results in much lower amounts of shallow and congestus clouds. Unlike the benchmark LES where the new deep clouds tend to appear along the edges of spreading cold pools, the deep clouds in the no-cold-pool experiment tend to reappear at the sites of the previous deep clouds and tend to be surrounded by extensive areas of sporadic shallow clouds. The vertical velocity statistics of updraft and downdraft cores below 6 km height are compared to aircraft observations made during GATE. The comparison shows generally good agreement, and strongly suggests that the LES simulation can be used as a benchmark to represent the dynamics of tropical deep convection on scales ranging from large turbulent eddies to mesoscale convective systems. The effect of horizontal grid resolution is examined by running the same case with progressively larger grid sizes of 200, 400, 800, and 1600 m. These runs show a reasonable agreement with the benchmark LES in statistics such as convective available potential energy, convective inhibition, cloud fraction, precipitation rates, and surface latent and sensible fluxes. All runs reveal a tri-model cloud distribution in the vertical. However, there are differences in the updraft-core cloud statistics, and convergence of statistical properties is found only between the LES benchmark and the run with 200 m grid size. The effect of vertical grid resolution is also investigated with another run that uses a typical cloud-resolving model (CRM) horizontal grid size on the order of 1 km and only 64 vertical levels. A comparison to the run with 256 vertical levels shows different vertical cloud distributions. It is concluded that representation of the often observed trimodal vertical distribution of clouds requires a vertical grid spacing in the range of 50-100 m in mid-to-low troposphere

An intercomparison of large-eddy simulations of the stable boundary layer

Results are presented from the first intercomparison of Large-eddy simulation (LES) models for the stable boundary layer (SBL), as part of the GABLS (Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study) initiative. A moderately stable case is used, based on Arctic observations. All models produce successful simulations, inasmuch as they reflect many of the results from local scaling theory and observations. Simulations performed at 1 m and 2 m resolution show only small changes in the mean profiles compared to coarser resolutions. Also, sensitivity to sub-grid models for individual models highlights their importance in SBL simulation at moderate resolution (6.25 m). Stability functions are derived from the LES using typical mixing lengths used in Numerical Weather Prediction (NWP) and climate models. The functions have smaller values than those used in NWP. There is also support for the use of K-profile similarity in parametrizations. Thus, the results provide improved understanding and motivate future developments of the parametrization of the SBL

Simulation of the intraseasonal variability over the Eastern Pacific ITCZ in climate models

During boreal summer, convective activity over the eastern Pacific (EPAC) inter-tropical convergence zone (ITCZ) exhibits vigorous intraseasonal variability (ISV). Previous observational studies identified two dominant ISV modes over the EPAC, i.e., a 40-day mode and a quasi-biweekly mode (QBM). The 40-day ISV mode is generally considered a local expression of the Madden-Julian Oscillation. However, in addition to the eastward propagation, northward propagation of the 40-day mode is also evident. The QBM mode bears a smaller spatial scale than the 40-day mode, and is largely characterized by northward propagation. While the ISV over the EPAC exerts significant influences on regional climate/weather systems, investigation of contemporary model capabilities in representing these ISV modes over the EPAC is limited. In this study, the model fidelity in representing these two dominant ISV modes over the EPAC is assessed by analyzing six atmospheric and three coupled general circulation models (GCMs), including one super-parameterized GCM (SPCAM) and one recently developed high-resolution GCM (GFDL HIRAM) with horizontal resolution of about 50 km. While it remains challenging for GCMs to faithfully represent these two ISV modes including their amplitude, evolution patterns, and periodicities, encouraging simulations are also noted. In general, SPCAM and HIRAM exhibit relatively superior skill in representing the two ISV modes over the EPAC. While the advantage of SPCAM is achieved through explicit representation of the cumulus process by the embedded 2-D cloud resolving models, the improved representation in HIRAM could be ascribed to the employment of a strongly entraining plume cumulus scheme, which inhibits the deep convection, and thus effectively enhances the stratiform rainfall. The sensitivity tests based on HIRAM also suggest that fine horizontal resolution could also be conducive to realistically capture the ISV over the EPAC, particularly for the QBM mode. Further analysis illustrates that the observed 40-day ISV mode over the EPAC is closely linked to the eastward propagating ISV signals from the Indian Ocean/Western Pacific, which is in agreement with the general impression that the 40-day ISV mode over the EPAC could be a local expression of the global Madden-Julian Oscillation (MJO). In contrast, the convective signals associated with the 40-day mode over the EPAC in most of the GCM simulations tend to originate between 150A degrees E and 150A degrees W, suggesting the 40-day ISV mode over the EPAC might be sustained without the forcing by the eastward propagating MJO. Further investigation is warranted towards improved understanding of the origin of the ISV over the EPAC.close5

Intercomparison of model simulations of mixed-phase clouds observed during the ARM Mixed-Phase Arctic Cloud Experiment. I: single-layer cloud

Results are presented from an intercomparison of single-column and cloud-resolving model simulations of a cold-air outbreak mixed-phase stratocumulus cloud observed during the Atmospheric Radiation Measurement (ARM) program's Mixed-Phase Arctic Cloud Experiment. The observed cloud occurred in a well-mixed boundary layer with a cloud top temperature of -15 C. The observed average liquid water path of around 160 g m{sup -2} was about two-thirds of the adiabatic value and much greater than the average mass of ice crystal precipitation which when integrated from the surface to cloud top was around 15 g m{sup -2}. The simulations were performed by seventeen single-column models (SCMs) and nine cloud-resolving models (CRMs). While the simulated ice water path is generally consistent with the observed values, the median SCM and CRM liquid water path is a factor of three smaller than observed. Results from a sensitivity study in which models removed ice microphysics suggest that in many models the interaction between liquid and ice-phase microphysics is responsible for the large model underestimate of liquid water path. Despite this general underestimate, the simulated liquid and ice water paths of several models are consistent with the observed values. Furthermore, there is evidence that models with more sophisticated microphysics simulate liquid and ice water paths that are in better agreement with the observed values, although considerable scatter is also present. Although no single factor guarantees a good simulation, these results emphasize the need for improvement in the model representation of mixed-phase microphysics