A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w _*. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).     

A GCSS Boundary-Layer Cloud Model Intercomparison Study Of The First Astex Lagrangian Experiment

Three single-column models (all with an explicit liquid water budget and compara-tively high vertical resolution) and three two-dimensional eddy-resolving models (including one with bin-resolved microphysics) are compared with observations from the first ASTEX Lagrangian experiment. This intercomparison was a part of the second GCSS boundary-layer cloud modelling workshop in August 1995.In the air column tracked during the first ASTEX Lagrangian experiment, a shallow subtropical drizzling stratocumulus-capped marine boundary layer deepens after two days into a cumulus capped boundary layer with patchy stratocumulus. The models are forced with time varying boundary conditions at the sea-surface and the capping inversion to simulate the changing environment of the air column.The models all predict the observed deepening and decoupling of the boundary layer quite well, with cumulus cloud evolution and thinning of the overlying stratocumulus. Thus these models all appear capable of predicting transitions between cloud and boundary-layer types with some skill. The models also produce realistic drizzle rates, but there are substantial quantitative differences in the cloud cover and liquid water path between models. The differences between the eddy-resolving model results are nearly as large as between the single column model results. The eddy resolving models give a more detailed picture of the boundary-layer evolution than the single-column models, but are still sensitive to the choice of microphysical and radiative parameterizations, sub-grid-scale turbulence models, and probably model resolution and dimensionality. One important example of the differences seen in these parameterizations is the absorption of solar radiation in a specified cloud layer, which varied by a factor of four between the model radiation parameterizations.     

THE STUDY ON THE WINTERTIME STRATOCUMULUS CLOUDS OVER ÜRÜMQI,XINJIANG OF CHINA——NUMERICAL MODELING OF STRATOCUMULUS CLOUDS

On the basis of the conceptual model,the development of wintertime stratocumulus cloud has been simulated by using a planetary boundary layer/cloud physics model.The main characteristics of the cloud structure and evolution obtained from the simulation are consistent with those from observations.By analyzing the modeled results,some details of the cloud structures and microphysics processes are discussed.     

中尺度降水集合预报随机参数扰动方法敏感性试验

中尺度降水模式预报具有很大的不确定性,为更好地描述与模式降水预报密切相关的物理过程关键参数的不确定性,基于中国气象局GRAPES（Global/Regional Assimilation and Prediction System）中尺度区域集合预报模式,从对模式降水预报不确定性有较大影响的积云对流、云微物理、边界层及近地面层等4个参数化方案中选取了18个关键参数,设计了一种随机参数扰动方案（Stochastically Perturbed Parameterization,SPP）,并通过2015年6—7月总计10 d的随机扰动集合预报试验,对比分析了SPP方案对不同物理过程参数扰动敏感性、随机场时、空尺度敏感性、能量变化特征及其集合预报效果。结果显示,对所选择的任一物理过程参数化方案增加SPP扰动后,降水及等压面要素的概率预报技巧优于无SPP扰动的预报,而扰动积云对流和边界层过程中的参数较扰动云微物理过程中的参数影响更显著,且同时扰动积云对流、云微物理、边界层及近地面层参数化方案中的18个参数的集合预报效果优于扰动任何单一物理过程中的部分参数,表明SPP方案能够有效地提高中尺度降水概率预报技巧;从能量变化特征可知,不同物理过程的参数扰动对动能、内能和总能量的影响层次和特征有所不同,但总体而言,扰动前后各项能量基本相同;随机场时、空尺度敏感性试验发现,SPP扰动随机场时间、空间相关尺度对集合预报效果有明显影响,当扰动随机场选用12 h抗相关时间及截断波数20时,集合预报结果最优。上述结果表明,SPP随机参数扰动方案不仅能够有效提高集合概率预报效果,还能够提高集合降水概率预报技巧,具有良好的业务应用与发展前景。      

A Sensitivity Study Of The Subtropical Ocean Surface Energy Balance To The Parameterization Of Precipitation From Stratocumulus Clouds

In the `First Lagrangian' of the Atlantic Stratocumulus Experiment(ASTEX), a cloudy air mass was tracked as it was advected by thetrade winds toward higher sea surface temperatures. In this study,a full diurnal cycle observed during this experiment is simulated andthe impact of the precipitation parameterization is examined. The modelwe use is the one dimensional version of the hydrostatic primitiveequation model MAR (Modéle Atmosphérique Régional) developed at the Université catholique de Louvain (UCL).It includes an E- turbulence closure, a wide-band formulationof the radiative transfer, and a parameterized microphysical schemeallowing partial condensation. The model realistically reproducesthe diurnal clearing of the cloud layer as well as the formation ofcumulus clouds under the stratocumulus deck. Nevertheless, as thesurface warms and the boundary layer becomes more convective,the simulation progressively differs from the observed evolution.Further experiments are carried out with different precipitationparameterizations commonly used in mesoscale modelsand general circulation models (GCMs).A strong sensitivity of the simulated liquid water path evolution isfound. The impact on the surface energy flux and the solar fluxreflected by the cloud is also examined. For both fluxes averagedover 24 hours, differences as large as 20 W m^-2 are obtainedbetween the various simulations. Low cloudiness covers large areasover the ocean and such errors on the reflected solar flux may stronglyaffect the Earth's radiative budget in GCM simulations. We estimatethat the impact on the globally averaged outgoing solar flux could beas large as 5 W m^-2. Furthermore, when atmospheric models arecoupled to ocean models, errors in the surface energy exchanges mayinduce significant drift in the simulated climate.     

基于CWRF的中国近海热带气旋集合模拟及其对物理参数化方案的敏感性研究

为评估CWRF模式的降尺度能力和其热带气旋模拟对物理参数化方案的敏感性,本文利用ERI再分析资料驱动CWRF在30km网格上对1982-2016年中国近海热带气旋开展了一次集合模拟.结果表明:CWRF与ERI均能模拟出热带气旋的季节变化和年际变化形势且均存在低估,但相较ERI,CWRF的降尺度技术和集合模拟可以再现更多...     

WRF4.0模式陆气耦合对不同物理参数化方案的敏感性初探

模式的不同参数化方案会对模拟的天气和气候产生影响,但如何影响陆气耦合还不是很清楚。利用WRF 4.0模式设计了一组16种不同参数化方案组合的集合试验,对华北地区2013年夏季气候进行模拟。结果表明,不同积云对流参数化方案对模拟的降水平均态影响较大,其次是短波辐射方案,而微物理方案影响较小。进一步利用潜热通量(LH,Latent Heat flux)和土壤湿度(SM,Soil Moisture)、抬升凝结高度(LCL,Lifting Condensation Level)的相关系数R_(SM,LH)和R_(LH,LCL)表征陆气耦合过程,评估了陆气耦合对不同参数化方案的敏感性。模式模拟的R_(SM,LH)在不同区域随着平均降水、相对湿度的增加以及平均短波辐射的减少而减小,而R_(LH,LCL)的变化正好相反。它们内在的物理机制类似,与气候态的干湿变化密切相关,即随着土壤湿度趋于饱和,蒸散发过程逐渐受太阳辐射的限制,陆面对大气的影响减弱。另外,陆气耦合强度受积云对流方案、微物理方案和短波辐射方案的影响而产生差别,这与不同参数化方案模拟的气候态的不同密切相关,因此本研究对于WRF 4.0模式中参数化方案特别是积云对流方案的选择同样有着重要的指导意义。     

Parametric representation of heat and moisture fluxes in cloud-topped mixed layers

Aircraft measurements were made from the NCAR Electra in stratus and stratocumulus clouds off the coast of California in June 1976. Several types of cloud conditions were observed, including (1) a broken layer less than 100 m thick, capped by an inversion at ~1000 m, (2) a broken stratocumulus layer ~300 m thick with an inversion at ~500 m, and (3) a solid stratocumulus layer ~250 m thick with an inversion at ~500 m. Although these observations indicate that a variety of cloud conditions may exist in mixed layers, simple one-dimensional mixed-layer models implicitly assume a solid cloud layer with no unsaturated region within the cloud. In order to generalize these simple models, a parametric representation of the heat and moisture fluxes is considered. In this scheme, the fluxes are parameterized in terms of the product of a cloud mass flux and the characteristic difference between the thermodynamic properties of an updraft-downdraft circulation. This representation allows for an explicit representation of the buoyancy flux when the downdraft has no liquid water.Data collected during these flights were used to calculate heat and moisture fluxes and to obtain the mean difference in the thermodynamic properties of the updrafts and downdrafts at a given level. The mass flux was calculated using updraft-downdraft differences and the fluxes. The mass fluxes obtained using various thermodynamic quantities are examined for consistency. The vertical distribution of the mass flux is determined. These results indicate that a mass flux formulation could prove to be useful in modeling applications where cloud conditions may vary between solid and broken.     

利用TWP-ICE试验资料对比两种边界层参数化方案

利用高分辨率WRF单气柱模式,选取了两种边界层参数化方案(YSU,MYJ),对TWP-ICE(Tropical Warm Pool International Cloud Experiment)试验期间的个例进行数值模拟,比较了两种方案对边界层结构、云和降水模拟的影响。结果表明:季风活跃期,YSU方案模拟的湍流交换系数较小,湍流混合偏弱,边界层内热通量偏小,使地表热量和水汽不易向上输送,水汽含量在近地表明显偏多,而在边界层及其以上大气层具有显著的干偏差,因此该方案模拟的云中液态水和固态水含量偏低,云量偏少,降水率偏小;MYJ方案对于季风活跃期的边界层结构具有较好的模拟能力,其模拟的云和降水更为准确。季风抑制期,MYJ方案模拟的夜间边界层结构存在较大误差,这是因为该方案模拟的夜间湍流交换系数较大,湍流混合偏强,边界层内热通量偏大,模拟的位温和水汽混合比在边界层内随高度变化较小,而观测廓线在边界层内存在较大梯度。季风抑制期两种方案模拟的云和降水均比观测值偏多,方案之间的差异较小。     

Current problems in the stratocumulus-topped atmospheric boundary layer

Extended sheets of stratocumulus (Sc) in the upper part of the atmospheric boundary layer (ABL) often occur under appropriate meteorological conditions. These cloud decks are important both in climate studies and in weather forecasting. We review the current knowledge of the turbulent structure of the ABL capped by a cloud deck, in the light of recent observations and model studies. The most important physical processes determining this structure are longwave radiative cooling at cloud top, shortwave radiative wanning by absorption in the cloud, surface buoyancy flux, and wind shear in the ABL. As a result, turbulence can cause entrainment against the buoyancy jump at cloud top. In cases where only longwave radiative fluxes and surface buoyancy fluxes are important, the turbulent structure is relatively well understood. When shortwave radiative fluxes and/or wind shear are also important, the resulting turbulent structure may change considerably. A decoupling of the cloud from the sub-cloud layer or of the top of the cloud from the rest of the ABL is then regularly observed. In no cases are the details of the entrainment at cloud top understood well enough to derive a relatively simple formulation that is consistent with observations. Cloud-top entrainment instability may lead to the break-up of a cloud deck (but also to cloud deepening). The role of mesoscale circulations in determining fractional cloudiness is not yet well understood.     

Coamps<Superscript>®</Superscript>-Les: Model Evaluation and Analysis of Second-and Third-Moment Vertical Velocity Budgets

The Naval Research Laboratory Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS^®) has been extended to perform as a large-eddy simulation (LES) model. It has been validated with a series of boundary-layer experiments spanning a range of cloud nighttime, and includes a nighttime stratocumulus case, a trade wind cumulus layer, shallow cumulus convection over land, and a mixed regime consisting of cumulus clouds under broken stratocumulus. COAMPS-LES results are in good agreement with other models for all the cases simulated. Exact numerical budgets for the vertical velocity second\((\overline{w^{'2}})\) and third moment\((\overline{w^{'3}})\) have been derived for the stratocumulus and trade wind cumulus cases. For the\(\overline{w^{'3}}\) budget in the stratocumulus, the buoyancy contribution from the updraughts and downdraughts largely cancel each other due to their similar magnitudes but opposite signs. In contrast, for the cumulus layer, the negative buoyancy contribution from the environmental downdraughts is negligible and the positive contribution from the updraughts completely dominates due to the conditional instability in the environment. As a result,\(\overline{w^{'3}}\) is significantly larger in the cumulus than in the stratocumulus layer.     

CWRF-based ensemble simulation of tropical cyclone activity near China and its sensitivity to the model physical parameterization schemes: 基于CWRF的中国近海热带气旋集合模拟及其对物理参数化方案的敏感性研究

To evaluate the downscaling ability with respect to tropical cyclones (TCs) near China and its sensitivity to the model physics representation, the authors performed a multi-physics ensemble simulation with the regional Climate–Weather Research and Forecasting (CWRF) model at a 30 km resolution driven by ERA-Interim reanalysis data. The ensemble consisted of 28 integrations during 1979–2016 with varying CWRF physics configurations. Both CWRF and ERA-Interim can generally capture the seasonal cycle and interannual variation of the TC number near China, but evidently underestimate them. The CWRF downscaling and its multi-physics ensemble can notably reduce the underestimation and significantly improve the simulation of the TC occurrences. The skill enhancement is especially large in terms of the interannual variation, which is most sensitive to the cumulus scheme, followed by the boundary layer, surface and radiation schemes, but weakly sensitive to the cloud and microphysics schemes. Generally, the Noah surface scheme, CAML(CAM radiation scheme as implemented by Liang together with the diagnostic cloud cover scheme of Xu and Randall(1996)) radiation scheme, prognostic cloud scheme, and Thompson microphysics scheme stand out for their better performance in simulating the interannual variation of TC number. However, the Emanuel cumulus and MYNN boundary layer schemes produce severe interannual biases. Our study provides a valuable reference for CWRF application to improve the understanding and prediction of TC activity.摘要为评估CWRF模式的降尺度能力和其热带气旋模拟对物理参数化方案的敏感性, 本文利用ERI再分析资料驱动CWRF在30km网格上对1982-2016年中国近海热带气旋开展了一次集合模拟.结果表明:CWRF与ERI均能模拟出热带气旋的季节变化和年际变化形势且均存在低估, 但相较ERI, CWRF的降尺度技术和集合模拟可以再现更多的热带气旋, 显著减少低估.年际变化结果提升最为明显, 它对积云方案最为敏感, 其次是边界层, 陆面和辐射方案, 对云和微物理方案较弱.该研究为应用CWRF理解和预报热带气旋提供了参考.     

Spatial patterns and trends of daily rainfall regime in Peninsular Malaysia during the southwest and northeast monsoons: 1975–2004

Model precipitation can be produced implicitly through convective parameterization schemes or explicitly through cloud microphysics schemes. These two precipitation production schemes control the spatial and temporal distribution of precipitation and consequently can yield distinct vertical profiles of heating and moistening in the atmosphere. The partition between implicit and explicit precipitation can be different as the model changes resolutions. Within the range of mesoscale resolutions (about 20 km) and cumulus scale, hybrid solutions are suggested, in which cumulus convection parameterization is acting together with the explicit form of representation. In this work, it is proposed that, as resolution increases, the convective scheme should convert less condensed water into precipitation. Part of the condensed water is made available to the cloud microphysics scheme and another part evaporates. At grid sizes smaller than 3 km, the convective scheme is still active in removing convective instability, but precipitation is produced by cloud microphysics. The Eta model version using KF cumulus parameterization was applied in this study. To evaluate the quantitative precipitation forecast, the Eta model with the KF scheme was used to simulate precipitation associated with the South Atlantic Convergence Zone (SACZ) and Cold Front (CF) events. Integrations with increasing horizontal resolutions were carried out for up to 5 days for the SACZ cases and up to 2 days for the CF cases. The precipitation partition showed that most of precipitation was generated by the implicit scheme. As the grid size decreased, the implicit precipitation increased and the explicit decreased. However, as model horizontal resolution increases, it is expected that precipitation be represented more explicitly. In the KF scheme, the fraction of liquid water or ice, generated by the scheme, which is converted into rain or snow is controlled by a parameter S ₁. An additional parameter was introduced into KF scheme and the parameter acts to evaporate a fraction of liquid water or ice left in the model grid by S ₁ and return moisture to the resolved scale. An F parameter was introduced to combine the effects of S ₁ and S ₂ parameters. The F parameter gives a measure of the conversion of cloud liquid water or ice to convective precipitation. A function dependent on the horizontal resolution was introduced into the KF scheme to influence the implicit and explicit precipitation partition. The explicit precipitation increased with model resolution. This function reduced the positive precipitation bias at all thresholds and for the studied weather systems. With increased horizontal resolution, the maximum precipitation area was better positioned and the total precipitation became closer to observations. Skill scores for all events at different forecast ranges showed precipitation forecast improvement with the inclusion of the function F.     

The impact of boundary forcing on RegCM4.2 surface energy budget

The surface energy budget components from two simulations of the regional climate model RegCM4.2 over the European/North African domain during the period 1989–2005 are analysed. The simulations differ in specified boundary forcings which were obtained from ERA-Interim reanalysis and the HadGEM2-ES Earth system model. Surface radiative and turbulent fluxes are compared against ERA-Interim. Errors in surface radiative fluxes are derived with respect to the Global Energy and Water Cycle Experiment/Surface Radiation Budget satellite-based products. In both space and time, we find a high degree of realism in the RegCM surface energy budget components, but some substantial errors and differences between the two simulations are also present. The most prominent error is an overestimation of the net surface shortwave radiation flux of more than 50 W/m² over central and southeastern Europe during summer months. This error strongly correlates with errors in the representation of total cloud cover, and less strongly with errors in surface albedo. During other seasons, the amplitude of the surface energy budget components is more in line with reference datasets. The errors may limit the usefulness of RegCM simulations in applications (e.g. high-quality simulation-driven impact studies). However, by using a simple diagnostic model for error interpretation, we suggest potential sensitivity studies aiming to reduce the underestimation of cloud cover and overestimation of shortwave radiation flux.     

A treatment for the stratocumulus-to-cumulus transition in GCMs

Numerical models of climate have great difficulties with the simulation of marine low clouds in the subtropical Pacific and Atlantic Oceans. It has been especially difficult to reproduce the observed geographical distributions of the different cloud regimes in those regions. The present study discusses mechanisms proposed in previous works for changing one regime into another. One criterion is based on the theory of stratocumulus destruction through cloud top entrainment instability due to buoyancy reversal—situations in which the mixture of two air parcels becomes denser than either of the original parcels due to evaporation of cloud water. Another criterion is based on the existence of decoupling in the boundary layer. When decoupled, the stratocumulus regime changes to another in which these clouds can still exist together with cumulus. In a LES study, the authors have suggested that a combination of those two criteria can be used to diagnose whether, at a location, the cloud regime corresponds to a well-mixed stratocumulus regime, a shallow cumulus regime, or to a transitional regime where the boundary layer is decoupled. The concept is tested in the framework of an atmospheric general circulation model (GCM). It is found that several outstanding features of disagreement between simulation and observation can be interpreted as misrepresentations of the cloud regimes by the GCM. A novel criterion for switching among regimes is proposed to alleviate the effects of these misrepresentations.     

层积云覆盖的海洋边界层云详细微物理过程的数值模拟

文中建立了一个含显式分档的云微物理模式和辐射传输模式的一维 3阶湍流闭合模式 ,该模式可用于研究海洋边界层云中气溶胶和云的相互作用过程 ,同时提出了一种新的动力模式和微物理模式耦合方法 ,该方法可使动力模式中液态水相关项可以直接由微物理模式变量计算得到。作为模式的初步应用模拟了 2 0 0 1年APEX/ACE Asia在西太平洋上所观测到的一个个例。模拟结果和观测资料比较表明该模式基本上模拟出层积云覆盖的海洋边界层的基本结构     

WRF模式物理参数化方案对一次局地大暴雨预报的影响研究

基于WRF(Weather Research and Forecasting)模式及其3Dvar(3-Dimentional Variational)资料同化系统,采用36、12、4 km嵌套网格进行快速更新循环同化和不同的微物理及积云对流参数化方案对比试验,对2011年5月8日鲁中一次局地大暴雨过程进行了研究。结果表明,快速更新循环同化地面观测资料是影响模式降水落区预报准确性的关键因素,不同的微物理和积云对流参数化方案主要影响降水强度预报。采用不同的微物理参数化方案和积云对流参数化方案进行降水预报对比试验表明,LIN方案和WSM6(WRF Single-Moment 6-class)微物理参数化方案对降水预报均较好,LIN方案降水预报较WSM6方案略强。4 km网格预报使用K-F (Kain-Fritsch)积云对流参数化方案或不使用积云对流参数化方案,预报的降水均较好。4 km网格使用旧的K-F积云对流参数化方案,预报的近地层大气风场偏弱,导致大气动力抬升作用偏弱,从而造成模式降水预报偏弱。     

Simulation of surface energy fluxes using high-resolution non-hydrostatic simulations and comparisons with measurements for the LITFASS-2003 experiment

Land-surface heterogeneity effects on the subgrid scale of regional climate and numerical weather prediction models are of vital interest for the energy and mass exchange between the surface and the atmospheric boundary layer. High-resolution numerical model simulations can be used to quantify these effects, and are a tool used to obtain area-averaged surface fluxes over heterogeneous land surfaces. We present high-resolution model simulations for the LITFASS area near Berlin during the LITFASS-2003 experiment, which were carried out using the non-hydrostatic model FOOT3DK of the University of Köln with horizontal resolutions of 1 km and 250 m. The LITFASS-2003 experimental dataset is used for comparison. The screen level quantities show good quality for the simulated pressure, temperature, humidity and wind speed and direction. Averaged over the four week experimental period, simulated surface energy fluxes at land stations show a small bias for the turbulent heat fluxes and an underestimation of the net radiation caused by excessive cloudiness in the simulations. For eight selected days with low cloud amounts, the net radiation bias is close to zero, but the sensible heat flux shows a strong positive bias. Large differences are found for latent heat fluxes over a lake, which are partly due to local effects on the measurements, but an additional problem seems to be the overestimation of the turbulent exchange under stable conditions in the daytime internal boundary layer over the lake. In the area average over the LITFASS area of 20 ×  20 km², again a strong positive bias of 70 W m^?2 for the sensible heat is present. For the low soil moisture conditions during June 2003, the simulation of the turbulent heat fluxes is sensitive to variations in the soil type and its hydrological properties. Under these conditions, the supply of ground water to the lowest soil layer should be accounted for. Different area-averaging methods are tested. The experimental set-up of the LITFASS-2003 experiment is found to be well suited for the computation of area-averaged turbulent heat fluxes.     

Evaluation of significance of spray evaporation in the atmospheric surface layer for computing the trajectory and intensity of tropical cyclones

The results are presented ofcomputation oftrajectory and intensity oftwo typhoons in 2012 for two groups of parameterization schemes of ocean-atmosphere energy exchange in tropical cyclones (TCs), boundary and surface layers of TCs, radiation fluxes, microphysics, and convection. The surface layer computations were carried out with and without the relationships that take into account the surface layer cool ing due to spray evaporation. It is demonstrated that the accounting of this effect exerts positive influence on the results of TC trajectory modeling at different stages. Ifthere are specific features in TC movement, for example, the loop formation, this improvement may be crucial. The accuracy ofthe results of computation of TC intensity increases as well.     

Air mass modification over Lake Michigan

The modification of a relatively cold air mass over the warm water of Lake Michigan is studied by using a two-dimensional nonlinear mesoscale model. Considerable amounts of heat and water vapor are supplied from the water surface to the lower atmosphere by turbulent eddies. A convective mixed layer develops and grows toward the downwind region with stratocumulus clouds over the lake.The model simulates the warming and moistening of the mixed layer, the development of a boundary layer, the divergence and convergence of wind near the coastlines, and the turbulent fluxes.The model warming of the mixed layer across the lake was about 2.2 °K and the moistening of the mixed layer was about 0.8 g kg^–1, which are comparable to 2.7 °K and 0.8 g kg^–1 observed by Lenschow (1973). The convective boundary layer, which includes the cloud layer, subcloud layer, and superadiabatic layer near the water surface, is well simulated. The tilt of the inversion which coincides with the cloud top is also well reproduced. When a prescribed cooling rate is applied at the cloud top, stronger turbulence and a deeper cloud layer are generated. Without the cooling, the cloud is shallow and the shape of the cloud base is determined by surface conditions. The rise of the inversion is due to upward vertical motion, and deepening of the convective layer in the downwind region.