Search

Cumulus (Cu) from shallow convection is one of the dominant cloud types over the Tibetan Plateau (TP) in the summer according to CloudSat- CALIPSO observations. Its thermodynamic effects on the atmospheric environment and impacts on the large-scale atmospheric circulation are studied in this paper using the Community Atmospheric Model, version 5.3 (CAM5.3). It is found that the model can reasonably simulate the unique distribution of diabatic heating and Cu over the TP. Shallow convection provides the dominant diabatic heating and drying to the lower and middle atmosphere over the TP. A sensitivity experiment indicates that without Cu over the TP, large-scale condensation and stratiform clouds would increase dramatically, which induces enhanced low-level wind and moisture convergence toward the TP, resulting in significantly enhanced monsoon circulation with remote impact on the areas far beyond the TP. Cu therefore acts as a safety valve to modulate the atmospheric environment that prevents the formation of superclusters of stratiform clouds and precipitation over the TP. [ABSTRACT FROM AUTHOR]

Deep convection frequently occurs on the eastern side of upper-level troughs, or potential vorticity (PV) anomalies. This is consistent with uplift ahead of a cyclonic PV anomaly, and consequent reduction in static stability and increase of convective available potential energy (CAPE). Nevertheless, the causal link between upper-level PV and deep convection has not been proven, and given that lift, moisture, and instability must all be present for deep convection to occur it is not clear that upper-level forcing is sufficient. In this paper a convective rainband that intensified ahead of a cyclonic PV anomaly in an environment with little CAPE (~10 J kg−1) is examined to determine the factors responsible for its intensification. The key feature was a low-level convergence line, arising from the remnants of an occluded front embedded in the low-level cyclonic flow. The rainband's intensity and morphology was influenced by the remnants of a tropopause fold that capped convection at midlevels in the southern part of the band, and by a reduction in upper-level static stability in the northern part of the band that allowed the convection to reach the tropopause. Ascent ahead of the trough appears to have played only a minor role in conditioning the atmosphere to convection: in most cases the ascending airstream had previously descended in the flow west of the trough axis. Thus, simple 'PV thinking' is not capable of describing the development of the rainband, and it is concluded that preexisting low-level wind and humidity features played the dominant role. [ABSTRACT FROM AUTHOR]

The relationship between n = 0 mixed Rossby-gravity waves (MRGs) and eastward inertio-gravity waves (EIGs) from Matsuno's shallow-water theory on an equatorial beta plane is studied using statistics of satellite brightness temperature T b and dynamical fields from ERA-Interim data. Unlike other observed convectively coupled equatorial waves, which have spectral signals well separated into eastward and westward modes, there is a continuum of MRG-EIG power standing above the background that peaks near wavenumber 0. This continuum is also present in the signals of dry stratospheric MRGs. While hundreds of papers have been written on MRGs, very little work on EIGs has appeared in the literature to date. The authors attribute this to the fact that EIG circulations are much weaker than those of MRGs for a given amount of divergence, making them more difficult to observe even though they strongly modulate convection. Empirical orthogonal function (EOF) and cross-spectral analysis of 2-6-day-filtered T b isolate zonally standing modes of synoptic-scale convection originally identified by Wallace in 1971. These display antisymmetric T b signals about the equator that propagate poleward with a period of around 4 days, along with westward-propagating MRG-like circulations that move through the T b patterns. Further analysis here and in Part II shows that these signatures are not artifacts of the EOF approach but result from a mixture of MRG or EIG modes occurring either in isolation or at the same time. [ABSTRACT FROM AUTHOR]

Atmospheric motion vectors (AMVs) have been produced for decades and remain an important source of wind information. Many studies have suggested that the traditional interpretation of AMVs as representative of the wind at cloud top is suboptimal and that they are more representative of the winds within the cloud. This paper investigates the vertical representativity of cloudy AMVs using both first-guess departure [observation − background ( O − B)] statistics and the simulation-study technique. A state-of-the-art convection-permitting mesoscale model ('UKV') is used in conjunction with a radiative transfer model and the Nowcasting Satellite Application Facility (NWCSAF) AMV package to produce synthetic AMVs over a 1-month period. The simulated upper-level AMVs suffered from large height-assignment errors uncharacteristic of those in reality; these issues were partially alleviated by using the model cloud top instead of the assigned height. In agreement with previous studies, both the simulated and real AMVs were found to have the closest fit to a layer mean of the model winds with the majority of the layer below the estimated cloud top. However, improvements in the fit between the AMVs and the model were also found by simply lowering the assigned height. A short NWP trial hinted that height reassignment might lead to short-range forecast improvements. The results of this study indicate that the simulation technique was able to match the usefulness of O − B statistics for AMVs associated with low- and medium-level clouds (albeit at a higher computational cost); however, challenges remain in the simulation of upper-level clouds. [ABSTRACT FROM AUTHOR]

Dynamical models of organized mesoscale convective systems have identified the important features that help maintain their overarching structure and longevity. The standard model is the trailing stratiform archetype, featuring a front-to-rear ascending circulation, a mesoscale downdraft circulation, and a cold pool/density current that affects the propagation speed and the maintenance of the system. However, this model does not represent all types of mesoscale convective systems, especially in moist environments where the evaporation-driven cold pools are weak and the convective inhibition is small. Moreover, questions remain about the role of gravity waves in creating and maintaining organized systems and affecting their propagation speed. This study presents simulations and dynamical models of self-organizing convection in a moist, low-convective inhibition environment and examines the long-lived convective regimes that emerge spontaneously. This paper, which is Part I of this study, specifically examines the structure, kinematics, and maintenance of long-lived, upshear-propagating convective systems that differ in important respects from the standard model of long-lived convective systems. Linear theory demonstrates the role of ducted gravity waves in maintaining the long-lived, upshear-propagating systems. A steady nonlinear model approximates the dynamics of upshear-propagating density currents that are key to the maintenance of the mesoscale convective system. [ABSTRACT FROM AUTHOR]

Atmospheric superrotation with prograde equatorial winds and an equatorial angular momentum maximum is ubiquitous in planetary atmospheres. It is clear that eddy fluxes of angular momentum toward the equator are necessary to generate it. But under what conditions superrotation arises has remained unclear. This paper presents simulations and a scaling theory that establish conditions under which superrotation occurs in terrestrial atmospheres. Whether superrotation arises depends on the relative importance of factors that favor or disfavor superrotation. Convection preferentially generates Rossby waves near the equator, where the Rossby number is O(1). Since the Rossby waves transport angular momentum toward their source regions, this favors superrotation. Meridional temperature gradients preferentially lead to baroclinic instability and wave generation away from the equator. Eddy transport of angular momentum toward the baroclinic source region implies transport out of low latitudes, which disfavors superrotation. Simulations with an idealized GCM show that superrotation tends to arise when the equatorial convective generation of wave activity and its associated eddy angular momentum flux convergence exceed the baroclinic eddy angular momentum flux divergence. Convective and baroclinic wave activity generation is related through scaling arguments to mean-flow properties, such as planetary rotation rates and meridional temperature gradients. The scaling arguments show, for example, that superrotation is favored when the off-equatorial baroclinicity and planetary rotation rates are low, as they are, for example, on Venus. Similarly, superrotation is favored when the convective heating strengthens, which may account for the superrotation seen in extreme global warming simulations. [ABSTRACT FROM AUTHOR]

Parameterizing deep convection at model resolutions of a few kilometers or less requires diverging from several simplifying assumptions valid at coarser resolutions. The separation or complementarity between the deep-convection scheme and the model cloud scheme must be addressed properly to prevent a double counting of some phenomena, account for evolution in time, and keep consistent results while approaching resolutions where deep convection can be treated explicitly (without parameterization). In this paper, the author formulates and tests a perturbation approach of the bulk mass-flux representation of deep convective updrafts and downdrafts. The subgrid deep-convection scheme represents only the effect of the unresolved part of the real updrafts, complementing an explicit part associated with the mean gridbox vertical velocity. Special attention is paid to the ordering and interactions of the moist parameterizations, the formulation of the closure and of the triggering, and the accounting of time evolution aspects. The multiresolution behavior of the scheme is assessed in the operational numerical prediction model ALARO [a version of the Action de Recherche Petite Echelle Grande Echelle-Aire Limitée Adaptation Dynamique Développement International (ARPEGE-ALADIN) operational limited area model with a revised and modular structure of the physical parameterizations]. Unlike most mass-flux schemes, the parameterized part gradually fades out when the model resolution is increased, allowing the results to approach those of an explicit treatment of deep convection. [ABSTRACT FROM AUTHOR]

It is now acknowledged that representing model uncertainty in atmospheric simulators is essential for the production of reliable probabilistic forecasts, and a number of different techniques have been proposed for this purpose. This paper presents new perturbed parameter schemes for use in the European Centre for Medium-Range Weather Forecasts (ECMWF) convection scheme. Two types of scheme are developed and implemented. Both schemes represent the joint uncertainty in four of the parameters in the convection parameterization scheme, which was estimated using the Ensemble Prediction and Parameter Estimation System (EPPES). The first scheme developed is a fixed perturbed parameter scheme, where the values of uncertain parameters are varied between ensemble members, but held constant over the duration of the forecast. The second is a stochastically varying perturbed parameter scheme. The performance of these schemes was compared to the ECMWF operational stochastic scheme, stochastically perturbed parameterization tendencies (SPPT), and to a model that does not represent uncertainty in convection. The skill of probabilistic forecasts made using the different models was evaluated. While the perturbed parameter schemes improve on the stochastic parameterization in some regards, the SPPT scheme outperforms the perturbed parameter approaches when considering forecast variables that are particularly sensitive to convection. Overall, SPPT schemes are the most skillful representations of model uncertainty owing to convection parameterization. [ABSTRACT FROM AUTHOR]

This study examines spiral rainbands in a numerical simulation of Hurricane Bill (2009). This paper, the first part of the study, evaluates the structures of spiral rainbands and compares them to previous observations. Four types of spiral rainbands are identified: principal, secondary, distant, and inner rainbands. Principal rainbands tend to be stationary relative to the storm center, while secondary rainbands are more transient and move around the storm center. Both principal and secondary rainbands are tilted radially outward with height and have many of the commonly observed kinematic features, such as overturning secondary circulation and enhanced tangential velocity on their radially outward sides. Principal rainbands are bounded by very dry air on their radially outward sides. Distant rainbands are radially inward-tilting convective features that have dense cold pools near the surface. Inner rainbands are made of shallow convection that appears to have originated from near the eyewall. Differences in the structures of spiral rainbands between observations and the Hurricane Bill simulation are noted. The second part of the study investigates how inner rainbands propagate and makes comparison with previously proposed hypotheses such as vortex Rossby waves. [ABSTRACT FROM AUTHOR]

Tropospheric circulation shifts have strong potential to impact surface climate. However, the magnitude of these shifts in a changing climate and the attending regional hydrological changes are difficult to project. Part of this difficulty arises from the lack of understanding of the physical mechanisms behind the circulation shifts themselves. To better delineate circulation shifts and their respective causes the circulation response is decomposed into 1) the 'direct' response to radiative forcings themselves and 2) the 'indirect' response to changing sea surface temperatures. Using ensembles of 90-day climate model simulations with immediate switch-on forcings, including perturbed greenhouse gas concentrations, stratospheric ozone concentrations, and sea surface temperatures, this paper documents the direct and indirect transient responses of the zonal-mean general circulation, and investigates the roles of previously proposed mechanisms in shifting the midlatitude jet. It is found that both the direct and indirect wind responses often begin in the lower stratosphere. Changes in midlatitude eddies are ubiquitous and synchronous with the midlatitude zonal wind response. Shifts in the critical latitude of wave absorption on either flank of the jet are not indicted as primary factors for the poleward-shifting jet, although some evidence for increasing equatorward wave reflection over the Southern Hemisphere in response to sea surface warming is seen. Mechanisms for the Northern Hemisphere jet shift are less clear. [ABSTRACT FROM AUTHOR]

This paper proposes a new formulation of the deep convection triggering for general circulation model convective parameterizations. This triggering is driven by evolving properties of the strongest boundary layer thermals. To investigate this, a statistical analysis of large-eddy simulation cloud fields in a case of transition from shallow to deep convection over a semiarid land is carried out at different stages of the transition from shallow to deep convection. Based on the dynamical and geometrical properties at cloud base, a new computation of the triggering is first proposed. The analysis of the distribution law of the maximum size of the thermals suggests that, in addition to this necessary condition, another triggering condition is required, that is, that this maximum horizontal size should exceed a certain threshold. This is explicitly represented stochastically. Therefore, the new formulation integrates the whole transition process from the first cloud to the first deep convective cell and can be decomposed into three steps: (i) the appearance of clouds, (ii) crossing of the inhibition layer, and (iii) deep convection triggering. [ABSTRACT FROM AUTHOR]

This paper presents a stochastic triggering parameterization for deep convection and its implementation in the latest standard version of the Laboratoire de Météorologie Dynamique-Zoom (LMDZ) general circulation model: LMDZ5B. The derivation of the formulation of this parameterization and the justification, based on large-eddy simulation results, for the main hypothesis was proposed in Part I of this study. Whereas the standard triggering formulation in LMDZ5B relies on the maximum vertical velocity within a mean bulk thermal, the new formulation presented here (i) considers a thermal size distribution instead of a bulk thermal, (ii) provides a statistical lifting energy at cloud base, (iii) proposes a three-step trigger (appearance of clouds, inhibition crossing, and exceeding of a cross-section threshold), and (iv) includes a stochastic component. Here the complete implementation is presented, with its coupling to the thermal model used to treat shallow convection in LMDZ5B. The parameterization is tested over various cases in a single-column model framework. A sensitivity study to each parameter introduced is also carried out. The impact of the new triggering is then evaluated in the single-column version of LMDZ on several case studies and in full 3D simulations. It is found that the new triggering (i) delays deep convection triggering, (ii) suppresses it over oceanic trade wind cumulus zones, (iii) increases the low-level cloudiness, and (iv) increases the convective variability. The scale-aware nature of this parameterization is also discussed. [ABSTRACT FROM AUTHOR]

The purpose of this paper is to validate the representation of topographic flows and moist convection over the European Alps in a convection-parameterizing simulation (CPM; Ax = 6.6 km) and two cloud-resolving simulations (CRM; Ax = 1.1 and 2.2 km). All simulations and further sensitivity experiments are validated against a large set of observations for an 18-day fair-weather summer period. The episode considered is characterized by pronounced plain-valley pressure gradients, strong daytime upvalley flows, and weak nighttime down-valley flows. In addition, convective precipitation is recorded during the late afternoon and is preceded by a phase of shallow convection. The observed transition from shallow to deep convection occurs within a 3-h period. The results indicate good agreement between both CRMs and the observed diurnal evolution in terms of near-surface winds, cloud formation, and precipitation. The differences between the two CRMs are surprisingly small. In contrast, the CPM produces too-early peaks of cloud cover and precipitation that are due to a too-early activation of deep convection. Detailed sensitivity experiments show that the convection scheme, rather than the underresolved small-scale topography, is responsible for the poor performance of the CPM. In addition, observations and simulations show that late-morning mass convergence does not correlate with afternoon precipitation. Rather, it is found that enhanced convective activity is related to increased conditional instability. [ABSTRACT FROM AUTHOR]

Time indices of the Madden-Julian oscillation (MJO) are often used to generate empirical forecasts of the global atmospheric circulation. Moist deep convection associated with the MJO initiates eastward-propagating Rossby waves that disperse into the midlatitudes. The background circulation then guides extratropical waves back into the tropics of the eastern Pacific Ocean. Previous works have shown that equatorial Rossby (ER) waves occur following intrusion of extratropical Rossby waves into the tropics. Westward-propagating ER waves and the MJO modulate the total convection. This convection modulates the zonal wind, which influences the location and existence of westerly wind ducts. These wind ducts, in turn, guide extratropical waves into the tropics. This paper demonstrates through a simple composite analysis that a simultaneous assessment of MJO and ER waves yields more information about the extratropical circulation during boreal winter than can be obtained based on either type of disturbance alone, or from a sum of the signals associated with the MJO and ER waves composited separately. This analysis, together with previous results, suggests a feedback loop between the MJO, these waves, and the extratropical circulation. Thus, assessment of the ER wave state during a particular phase of the MJO might yield better empirical prediction of the global atmospheric circulation that follows. [ABSTRACT FROM AUTHOR]

Persistent superrotation is seen in the atmospheres of other terrestrial bodies (Venus, Titan) but not in that of present Earth, which is distinguished by equatorial easterlies. Nevertheless, superrotation has appeared in numerical simulations of Earth's atmosphere, from two-layer models to multilevel comprehensive GCMs. Simulations of warm climates that generate enhanced tropical convective variability seem particularly prone to superrotation, which has led to hypotheses that the warmer atmospheres of the early Pliocene and Eocene may have been superrotating, and that the phenomenon may be relevant to future climate projections. This paper considers a positive feedback leading to superrotation based on an equatorial wave resonance that occurs in a westerly background flow. The authors present simulations with an idealized multilevel GCM forced with a zonally varying equatorial heating, which show abrupt transitions to strongly superrotating states. Linear shallow water theory is used to show that these transitions occur as the superrotating jet velocity approaches the phase speed of free equatorial Rossby wave modes, leading to a resonant amplification of the response to eddy heating and its associated equatorward momentum flux. The resonance and transition are most prominent in simulations where the meridional temperature gradient has been reduced, and hysteresis behavior is seen when the gradient is eliminated completely. No evidence is found in these simulations for the midlatitude wave feedback believed to drive abrupt transitions in two-layer models, and there is only a minor role for the axisymmetric feedback based on vertical advection by the Hadley circulation. [ABSTRACT FROM AUTHOR]

This paper explores the response of the tropical hydrologic cycle to surface warming through the lens of large-domain cloud-system-resolving model experiments run in a radiative-convective equilibrium framework. Simulations are run for 55 days and are driven with fixed insolation and constant sea surface temparatures (SSTs) of 298 K, 300 K, and 302 K. In each experiment, convection organizes into coherent regions of large-scale ascent separated by areas with relatively clear air and troposphere-deep descent. Aspects of the simulations correspond to observed features of the tropical climate system, including the transition to large precipitation rates above a critical value of total column water vapor, and an increase in convective intensity with SST amidst weakening of the large-scale overturning circulation. However, the authors also find notable changes to the interaction between convection and the environment as the surface warms. In particular, organized convection in simulations with SSTs of 298 and 300 K is inhibited by the presence of a strong midtropospheric stable layer and dry upper troposphere. As a result, there is a decrease in the vigor of deep convection and an increase in stratiform precipitation fraction with an increase in SST from 298 to 300 K. With an increase in SST to 302 K, moistening of the middletroposphere and increase in lower-tropospheric buoyancy serve to overcome these limitations, leading to an overall increase in convective intensity and larger increase in upper-tropospheric relative humidity. The authors conclude that, while convective intensity increases with SST, the aggregate nature of deep convection is strongly affected by the details of the thermodynamic environment in which it develops. In particular, the positive feedback between increasing SST and a moistening upper troposphere found in the simulations, operates as a nonmonotonic function of SST and is modulated by a complex interaction between deep convection and the environmental relative humidity and static stability profile. The results suggest that projected changes in convection that assume a monotonic dependence on SST may constitute an oversimplification. [ABSTRACT FROM AUTHOR]

Gross moist stability (GMS), a measure of how efficiently divergent flow exports column-integrated moist static energy (MSE), is a widely used quantity in current simplified models of the tropical mean circulation and intraseasonal variabilities such as the Madden--Julian oscillation (MJO), where it is often assumed to be constant. In this paper, it is shown, with cloud-system-resolving model experiments that incorporate feedback from the large-scale flow, that the GMS is smaller at longer wavelengths. The reason for this wavelength dependence is that temperature anomalies required to maintain a given divergent flow increase with wavelength. At long wavelengths, the required temperature anomalies become sufficiently strong to affect the shape of convective heating. As a consequence, the divergent flow is forced to be less top heavy in order to maintain the balance of momentum, heat, and moisture, as well as consistency with the behavior of cumulus convection. A simple model is constructed to illustrate this behavior. Given the ongoing theoretical efforts that view the MJO as resulting from instability in column-integrated MSE, the results presented here provide a planetary-scale selection for such instability, which is absent in current theoretical models that assume a constant GMS. [ABSTRACT FROM AUTHOR]

The Boundary Layer Evolution (BLE) missions of the International H2O Project (IHOP_2002) were designed to provide comprehensive observations of the distribution of water vapor in the quiescent boundary layer and its evolution during the early morning. The case study discussed in this paper presents detailed observations of the development of the boundary layer from before sunrise through to the period of growth of the mature convective boundary layer (CBL) during the 14 June 2002 BLE mission. The large number of remote sensing platforms, including the multiple instruments collocated at the Homestead Profiling Site, provided a detailed set of measurements of the growth and structure of the CBL. The observations describe the classic evolution of a daytime CBL, beginning with a shallow nocturnal boundary layer (NBL) below the remnants of the previous day’s mixed layer, or residual layer. The vertical distribution of humidity in these layers during the early morning was affected by advection of dry air and by gravity waves. About an hour after sunrise a CBL developed, and gradually deepened with time as it mixed out the NBL and residual layer. The growth of the top of the CBL was particularly well observed because of the strong vertical gradients in temperature, humidity, and aerosol concentration. As the CBL deepened and the average CBL wind speed decreased, the mode of convective organization evolved from horizontal convective rolls to open-celled convection. A unique set of detailed measurements of the structure of the open cells was obtained from multiple instruments including the Doppler-on-Wheels radar, the Mobile Integrated Profiling System wind profiler, and the Scanning Raman lidar. They showed the relationship between open cells, thermals, mantle echoes, and the CBL top. [ABSTRACT FROM AUTHOR]

Wildfires are capable of inducing atmospheric circulations that result predominantly from large temperature anomalies produced by the fire. The fundamental dynamics through which a forest fire and the atmosphere interact to yield different convective regimes is still not well understood. This study uses the Advanced Regional Prediction System (ARPS) model to investigate the impact of the environmental (i.e., far upstream, undisturbed by fire) wind profile on dry convection above a prescribed heat source of an intensity and spatial scale comparable to a wildfire. Dimensional analysis of the fire–atmosphere problem provides two relevant parameters: a surface buoyancy parameter that addresses the amount of heat a parcel of air receives in transiting above the fire and an advection parameter that addresses the degree to which the environmental wind advects updrafts away from the fire. Two-dimensional simulations are performed in which the upstream surface wind speed and mixed layer mean wind speed are varied independently to better understand the fundamental processes governing the organizational mode and updraft strength. The result of these experiments is the identification of two primary classes of dry convection: plume and multicell. Simulated plume cases exhibit weak advection by the mean wind and are subdivided into intense plume and hybrid classes based on the degree of steadiness within the convection column. Hybrid cases contain columns of largely discrete updrafts versus the more continuous updraft column associated with the intense plume mode. Multicell cases develop with strong mixed layer advection and are subdivided into strong and weak classes based on the depth of convection. Intense plume and strong multicell (hybrid and weak multicell) cases occur when the surface buoyancy is large (small). Parcel analyses are performed to more closely examine the forcing of convection within different areas of the parameter space. The multicell (strong and weak) and intense plume modes are forced by a combination of buoyancy and dynamic pressure gradient forcing associated with the perturbation wind field, whereas the hybrid mode is forced by a combination of buoyancy and dynamic pressure gradient forcing associated with the strong background shear. The paper concludes with a discussion of the degree of nonlinearity that is likely to exist at the fire front for each of the convective modes; nonlinear fire behavior is most likely for the hybrid mode and least likely for the weak multicell mode. Knowledge of the sensitivity of the convective mode to upstream conditions can provide information about the degree of nonlinear or erratic fire behavior expected for a given wind profile upstreamof the fire. [ABSTRACT FROM AUTHOR]

With many operational centers moving toward order 1-km-gridlength models for routine weather forecasting, this paper presents a systematic investigation of the properties of high-resolution versions of the Met Office Unified Model for short-range forecasting of convective rainfall events. The authors describe a suite of configurations of the Met Office Unified Model running with grid lengths of 12, 4, and 1 km and analyze results from these models for a number of convective cases from the summers of 2003, 2004, and 2005. The analysis includes subjective evaluation of the rainfall fields and comparisons of rainfall amounts, initiation, cell statistics, and a scale-selective verification technique. It is shown that the 4- and 1-km-gridlength models often give more realistic-looking precipitation fields because convection is represented explicitly rather than parameterized. However, the 4-km model representation suffers from large convective cells and delayed initiation because the grid length is too long to correctly reproduce the convection explicitly. These problems are not as evident in the 1-km model, although it does suffer from too numerous small cells in some situations. Both the 4- and 1-km models suffer from poor representation at the start of the forecast in the period when the high-resolution detail is spinning up from the lower-resolution (12 km) starting data used. A scale-selective precipitation verification technique implies that for later times in the forecasts (after the spinup period) the 1-km model performs better than the 12- and 4-km models for lower rainfall thresholds. For higher thresholds the 4-km model scores almost as well as the 1-km model, and both do better than the 12-km model. [ABSTRACT FROM AUTHOR]

In the Midwest U.S. Corn Belt, the 1999 and 2000 summer seasons (15 June–15 September) expressed contrasting spatial patterns and magnitudes of precipitation (1999: dry; 2000: normal to moist). Distinct from the numerical modeling approach often used in studies of land surface–climate interactions, a “synoptic climatological” (i.e., stratified composite) approach is applied to observation data (e.g., precipitation, radar, and atmospheric reanalyses) to determine the relative influences of “top-down” synoptic atmospheric circulation (Part I, this paper) and “bottom-up” land surface mesoscale conditions (Part II) on the predominantly convective precipitation variations. Because mesoscale modeling suggests that the free-atmosphere wind speed (“background wind”) regulates the land surface–atmosphere mesoscale interaction, each day’s spatial range of wind speed at 500 hPa [V(500)] over the Central Corn Belt (CCB) is classified into one of five categories ranging from “weak flow” to “jet maximum.” Deep convective activity (i.e., presence/absence and morphological signature type) is determined for each afternoon and early evening period from the Next Generation Weather Radar (NEXRAD) imagery. Frequencies of the resulting background wind–convection joint occurrence types for the 1999 and 2000 summer seasons are examined in the context of the statistics determined for summers in the longer period of 1996–2001, and also compose categories for which NCEP–NCAR reanalysis (NNR) fields are averaged to yield synoptic composite environments for the two study seasons. The latter composites are compared visually with high-resolution (spatial) composites of precipitation to help identify the influence of top-down climate controls. The analysis confirms that reduced (increased) organization of radar-indicated deep convection tends to occur with weaker (stronger) background flow. The summers of 1999 and 2000 differ from one another in terms of background flow and convective activity, but more so with respect to the six-summer averages, indicating that a fuller explanation of the precipitation differences in the two summers must be sought in the analysis of additional synoptic meteorological variables. The composite synoptic conditions on convection (CV) days (no convection (NC) days) in 1999 and 2000 are generalized as follows: low pressure incoming from the west (high pressure or ridging), southerly (northerly) lower-tropospheric winds, positive (negative) anomalies of moisture in the lower troposphere, rising (sinking) air in the midtroposphere, and a location south of the upper-tropospheric jet maximum (absence of an upper-tropospheric jet or one located just south of the area). Features resembling the “northerly low-level jets” identified in previous studies for the Great Plains are present on some NC-day composites. On CV days the spatial synchronization of synoptic features implying baroclinity increases with increasing background wind speed. The CV and NC composites differ least on days of weaker flow, and there are small areas within the CCB having no obvious association between precipitation elevated amounts and synoptic circulation features favoring the upward motion of air. These spatial incongruities imply a contributory influence of “stationary” (i.e., climatic) land surface mesoscale processes in convective activity, which are examined in Part II. [ABSTRACT FROM AUTHOR]

This paper deals with the dynamical aspect of variational data assimilation in the tropics and the role of the background-error covariances in the observing system simulation experiments for the tropics. The study uses a model that describes the horizontal structure of the potential temperature and wind fields in regions of deep tropical convection. The assimilation method is three- and four-dimensional variational data assimilation. The background-error covariance model for the assimilation is a multivariate model that includes the mass–wind couplings representative of equatorial inertio-gravity modes and equatorial Kelvin and mixed Rossby–gravity modes in addition to those representative of balanced equatorial Rossby waves. Spectra of the background errors based on these waves are derived from the tropical forecast errors of the European Centre for Medium-Range Weather Forecasts (ECMWF) model. Tropical mass–wind (im)balances are illustrated by studying the potential impact of the spaceborne Doppler wind lidar (DWL) Atmospheric Dynamic Mission (ADM)-Aeolus, which measures horizontal line-of-sight (LOS) wind components. Several scenarios with two DWLs of ADM-Aeolus type are compared under different flow conditions and using different assumptions about the quality of the background-error covariances. Results of three-dimensional variational data assimilation (3DVAR) illustrate the inefficiency of multivariate assimilation in the tropics. The consequence for the assimilation of LOS winds is that the missing part of the wind vector can hardly be reconstructed from the mass-field observations and applied balances as in the case of the midlatitudes. Results of four-dimensional variational data assimilation (4DVAR) show that for large-scale tropical conditions and using reliable background-error statistics, differences among various DWL scenarios are not large. As the background-error covariances becomes less reliable, horizontal scales become smaller and the flow becomes less zonal, the importance of obtaining information about the wind vector increases. The added value of another DWL satellite increases as the quality of the background-error covariances deteriorates and it can be more than twice as large as in the case of reliable covariances. [ABSTRACT FROM AUTHOR]

The multiple equilibria of the thermohaline circulation (THC: used here in the sense of the meridional overturning circulation) as function of the surface freshwater flux has been studied intensively following a Stommel paper from 1961. It is shown here that multistability and hysteresis of the THC also exist when the wind stress amplitude is varied as a control parameter. Both the Massachusetts Institute of Technology ocean general circulation model (MITgcm) and a simple three-box model are used to study and explain different dynamical regimes of the THC and THC variability as a function of the wind stress amplitude. Starting with active winds and a thermally dominant thermohaline circulation state, the wind stress amplitude is slowly reduced to zero over a time period of ∼40 000 yr (40 kyr) and then increased again to its initial value over another ∼40 kyr. It is found that during the decreasing wind stress phase, the THC remains thermally dominant until very low wind stress amplitude at which pronounced Dansgaard–Oeschger-like THC relaxation oscillations are initiated. However, while the wind stress amplitude is increased, these relaxation oscillations are present up to significantly larger wind stress amplitude. The results of this study thus suggest that under the same wind stress amplitude, the THC can be either in a stable thermally dominant state or in a pronounced relaxation oscillations state. The simple box model analysis suggests that the observed hysteresis is due to the combination of the Stommel hysteresis and the Winton and Sarachik “deep decoupling” oscillations. [ABSTRACT FROM AUTHOR]

The generally accepted view of the meridional circulation in the tropical east Pacific is that of a single deep overturning cell driven by deep convective heating in the intertropical convergence zone (ITCZ), similar to the zonal mean Hadley circulation. However, recent observations of the atmosphere from the tropical eastern Pacific have called this view into question. In several independent datasets, significant meridional return flows out of the ITCZ region were observed, not only at high altitudes, but also at low altitudes, just above the atmospheric boundary layer. This paper presents a theory and idealized simulations to understand the causes and dynamics of this shallow meridional circulation (SMC). Fundamentally, the SMC can be seen as a large-scale sea-breeze circulation driven by sea surface temperature gradients when deep convection is absent in the ITCZ region. A simple model of this circulation is presented. Using observed values, the sea-breeze model shows that the pressure gradient above the boundary can indeed reverse, leading to the pressure force that drives the shallow return flow out of the ITCZ. The Weather Research and Forecast Model (WRF) is used to simulate an idealized Hadley circulation driven by moist convection in a tropical channel. The SMC is reproduced, with reasonable similarity to the circulation observed in the east Pacific. The simulations confirm that the SMC is driven by a reversal of the pressure gradient above the boundary layer, and that the return flow is strongest when deep convection is absent in the ITCZ, and weakest when deep convection is active. The model also shows that moisture transport out of the ITCZ region is far greater in the low-level shallow return flow than in the high-altitude return flow associated with the deep overturning, and that a budget for water transport in and out of the ITCZ region is grossly incomplete without it. Much of the moisture carried in the shallow return flow is recycled into the boundary layer, but does not appear to contribute to enhanced cloudiness in the subtropical stratocumulus poleward of the ITCZ. [ABSTRACT FROM AUTHOR]

The accurate parameterization of moist convection presents a major challenge for the accurate prediction of weather and climate through numerical models. Superparameterization is a promising recent alternative strategy for including the effects of moist convection through explicit turbulent fluxes calculated from a cloud-resolving model. Basic scales for cloud-resolving modeling are the microscales on the order of 10 km in space on time scales on the order of 15 min, where vertical and horizontal motions are comparable and moist processes are strongly nonlinear (meso-gamma scale). In this paper, systematic multiscale asymptotic analysis is utilized to develop simplified microscale mesoscale dynamic (MMD) models for interaction between the microscales and spatiotemporal mesoscales on the order of 100 km and 2.5 h (meso-beta scale). The new MMD models lead to a systematic framework for superparameterization for numerical weather prediction (NWP) generalizing the traditional column modeling framework. The MMD formulation also provides a flexible systematic framework for devising new parameterization strategies for NWP intermediate between the two extremes of column modeling and detailed cloud-resolving modeling. It is also established here that these MMD models fit crudely into the recent systematic multiscale framework developed to explain the observed larger-scale statistical self-similarity of tropical convection, and therefore provide a systematic framework for superparameterization. Finally, it is shown that the new MMD models have the structure of a heterogeneous multiscale method so that many numerical techniques recently developed in the applied mathematics literature can be applied to this formulation. [ABSTRACT FROM AUTHOR]

This paper describes a technique for photogrammetric analysis of stereo pairs of images that is applied to the study of orographic convection. The technique is designed for use with digital images and assumes detailed knowledge of the camera properties (focal length and imaging chip) and that the position and orientation are known as a first guess. An iterative scheme using known landmarks on the frame is used to determine the camera orientation. The scheme is accurate to 10–100 m at a distance of 15 km from the camera pair. The transition from shallow to deep convection over the Santa Catalina Mountains in southern Arizona on 26 July 2005 is presented. The three-dimensional structure of the visible portion of the cloud is determined and compared with the composite reflectivity from the National Weather Service Weather Surveillance Radar-1988 Doppler radar and the tropopause height from the 1200 UTC sounding in Tucson, Arizona, providing additional validation of the scheme. The shallow to deep transition is characterized by tracking individual turrets and determining the maximum height of the cloud top. The cloud tops were limited to beneath 6000 m MSL for the first 1.5 h followed by the development of deep convection. The motion of the turrets and location of the eventual deep convection were consistent with the idea that moistening by shallow convection conditions the atmosphere for further development. [ABSTRACT FROM AUTHOR]

This paper describes an analysis of large-scale [O(1000 km)] convectively coupled gravity waves simulated using a two-dimensional cloud-resolving model. The waves develop spontaneously under uniform radiative cooling and approximately zero-mean-flow conditions, with wavenumber 2 of the domain appearing most prominently and right-moving components dominating over left-moving components for random reasons. The analysis discretizes the model output in two ways. First, a vertical-mode transform projects profiles of winds, temperature, and heating onto the vertical modes of the model’s base-state atmosphere. Second, a cloud-partitioning algorithm sorts sufficiently cloudy grid columns into three categories: shallow convective, deep convective, and stratiform anvil. Results show that much of the tilted structures of the waves can be captured by just two main vertical spectral “bands,” each consisting of a pair of vertical modes. The “slow” modes have propagation speeds of 16 and 18 m s-1 (and roughly a full-wavelength vertical structure through the troposphere), while the “fast” modes have speeds of 35 and 45 m s-1 (and roughly a half-wavelength structure). Deep convection anomalies in the waves are more or less in phase with the low-level cold temperature anomalies of the slow modes and in quadrature with those of the fast modes. Owing to the characteristic life cycle of deep convective cloud systems, shallow convective heating peaks ∼2 h prior to maximum deep convective heating, while stratiform heating peaks ∼3 h after. The onset of deep convection in the waves is preceded by a gradual deepening of shallow convection lasting a period of many hours. Results of this study are in broad agreement with simple two-mode models of unstable large-scale wave growth, under the name “stratiform instability.” Differences here are that 1) the key dynamical modes have speeds in the range 16–18 m s-1, rather than 23–25 m s-1 (owing to a shallower depth of imposed radiative cooling), and 2) deep convective heating, as well as stratiform heating, is essential for the generation and maintenance of the slow modes. [ABSTRACT FROM AUTHOR]

A better conceptual understanding and more realistic parameterizations of convective boundary layers in climate and weather prediction models have been major challenges in meteorological research. In particular, parameterizations of the dry convective boundary layer, in spite of the absence of water phase-changes and its consequent simplicity as compared to moist convection, typically suffer from problems in attempting to represent realistically the boundary layer growth and what is often referred to as countergradient fluxes. The eddy-diffusivity (ED) approach has been relatively successful in representing some characteristics of neutral boundary layers and surface layers in general. The mass-flux (MF) approach, on the other hand, has been used for the parameterization of shallow and deep moist convection. In this paper, a new approach that relies on a combination of the ED and MF parameterizations (EDMF) is proposed for the dry convective boundary layer. It is shown that the EDMF approach follows naturally from a decomposition of the turbulent fluxes into 1) a part that includes strong organized updrafts, and 2) a remaining turbulent field. At the basis of the EDMF approach is the concept that nonlocal subgrid transport due to the strong updrafts is taken into account by the MF approach, while the remaining transport is taken into account by an ED closure. Large-eddy simulation (LES) results of the dry convective boundary layer are used to support the theoretical framework of this new approach and to determine the parameters of the EDMF model. The performance of the new formulation is evaluated against LES results, and it is shown that the EDMF closure is able to reproduce the main properties of dry convective boundary layers in a realistic manner. Furthermore, it will be shown that this approach has strong advantages over the more traditional countergradient approach, especially in the entrainment layer. As a result, this EDMF approach opens the way to parameterize the clear and cumulus-topped boundary layer in a simple and unified way. [ABSTRACT FROM AUTHOR]

This paper presents the revision and evaluation of the interface between the convective parameterization by Emanuel and Živković-Rothman and the Lagrangian particle dispersion model “FLEXPART” based on meteorological data from the European Centre for Medium-Range Weather Forecasts (ECMWF). The convection scheme relies on the ECMWF grid-scale temperature and humidity and provides a matrix necessary for the vertical convective particle displacement. The benefits of the revised interface relative to its previous version are presented. It is shown that, apart from minor fluctuations caused by the stochastic convective redistribution of the particles, the well-mixed criterion is fulfilled in simulations that include convection. Although for technical reasons the calculation of the displacement matrix differs somewhat between the forward and the backward simulations in time, the mean relative difference between the convective mass fluxes in forward and backward simulations is below 3% and can therefore be tolerated. A comparison of the convective mass fluxes and precipitation rates with those archived in the 40-yr ECMWF Reanalysis (ERA-40) data reveals that the convection scheme in FLEXPART produces upward mass fluxes and precipitation rates that are generally smaller by about 25% than those from ERA-40. This result is interpreted as positive, because precipitation is known to be overestimated by the ECMWF model. Tracer transport simulations with and without convection are compared with surface and aircraft measurements from two tracer experiments and to 222Rn measurements from two aircraft campaigns. At the surface no substantial differences between the model runs with and without convection are found, but at higher altitudes the model runs with convection produced better agreement with the measurements in most of the cases and indifferent results in the others. However, for the tracer experiments only few measurements at higher altitudes are available, and for the aircraft campaigns the 222Rn emissions are highly uncertain. Other datasets better suitable for the validation of convective transport in models are not available. Thus, there is a clear need for reliable datasets suitable to validate vertical transport in models. [ABSTRACT FROM AUTHOR]

The reported study examines the dynamics of entrainment and its effects on the evolution of the dry atmospheric convective boundary layer (CBL) when wind shear is present. The sheared CBL can be studied by means of direct measurements in the atmosphere, laboratory studies, and numerical techniques. The advantages and disadvantages of each technique are discussed in the present paper, which also describes the methodological background for studying the dynamics of entrainment in sheared CBLs. For the reported study, large-eddy simulation (LES) was chosen as the primary method of convective entrainment investigation. Twenty-four LES runs were conducted for CBLs growing under varying conditions of surface buoyancy flux, free-atmospheric stratification, and wind shear. The simulations were divided into three categories: CBL with no mean wind (NS), CBL with a height-constant geostrophic wind of 20 m s-1 (GC), and CBL with geostrophic wind shear (GS). In the simulated cases, the sheared CBLs grew fastest, relative to the NS CBLs, when the surface buoyancy flux was weak and the atmospheric stratification was moderate or weak. Three fundamental findings resulted from the investigated CBL cases: (i) the entrainment zone shear is much more important than the surface shear in enhancing CBL entrainment, although entrainment zone shear is indirectly affected by surface shear; (ii) the sheared entrainment zone features a sublayer of nearly constant flux Richardson number, which points to a balance between shear production and buoyancy consumption of turbulence kinetic energy (TKE) that regulates entrainment; and (iii) the fraction of entrainment zone shear-generated TKE spent on the entrainment is lower than suggested by earlier studies. [ABSTRACT FROM AUTHOR]

The detailed analysis of the three-dimensional structure of a dryline observed over the Oklahoma panhandle during the International H2O Project (IHOP_2002) on 11 June 2002 is presented. High-resolution observations obtained from the National Center for Atmospheric Research Electra Doppler Radar (ELDORA), S-band dual-polarization Doppler radar (S-Pol), water vapor differential absorption lidar (DIAL) Lidar pour l'Etude des Interactions Aérosols Nuages Dynamique Rayonnement et du Cycle de l'Eau (LEANDRE II; translated as Lidar for the Study of Aerosol–Cloud–Dynamics–Radiation Interactions and of the Water Cycle) as well as Learjet dropsondes are used to reveal the evolution of the dryline structure during late afternoon hours when the dryline was retreating to the northwest. The dryline reflectivity shows significant variability in the along-line direction. Dry air was observed to overrun the moist air in vertical cross sections similar to a density current. The updrafts associated with the dryline were 2–3 m s-1 and were able to initiate boundary-layer-based clouds along the dryline. The formation of this dryline was caused by high equivalent potential temperature air pushing northwestward toward a stationary front in the warm sector. Middle-level clouds with radar reflectivity greater than 18 dBZe near the dryline were detected by ELDORA. A roll boundary, which was associated with larger convergence and moisture content, was evident in the S-Pol data. It is found that the instability parameters most favorable for convection initiation were actually associated with the roll boundary, not the dryline. A storm was initiated near the roll boundary probably as a result of the combination of the favorable instability parameters and stronger upward forcing. It is noted that both the 11 June 2002 dryline and the roll boundary presented in this paper would not be identified if the special datasets from IHOP_2002 were not available. Although all model runs [fifth-generation Pennsylvania State University–NCAR Mesoscale Model (MM5), Meso Eta, and Rapid Update Cycle (RUC)] suggested deep convection over the Oklahoma panhandle and several cloud lines were observed near the dryline, the dryline itself did not initiate any storms. The reasons why the dryline failed to produce any storm inside the IHOP_2002 intensive observation region are discussed. Both synoptic-scale and mesoscale conditions that were detrimental to convection initiation in this case are investigated in great detail. [ABSTRACT FROM AUTHOR]

Results from a high-resolution numerical simulation of the 24 May 2002 dryline convective initiation (CI) case are presented. The simulation uses a 400 km × 700 km domain with a 1-km horizontal resolution grid nested inside a 3-km domain and starts from an assimilated initial condition at 1800 UTC. Routine as well as special upper-air and surface observations collected during the International H2O Project (IHOP_2002) are assimilated into the initial condition. The initiation of convective storms at around 2015 UTC along a section of the dryline south of the Texas panhandle is correctly predicted, as is the noninitiation of convection at a cold-front–dryline intersection (triple point) located farther north. The timing and location of predicted CI are accurate to within 20 min and 25 km, respectively. The general evolution of the predicted convective line up to 6 h of model time also verifies well. Mesoscale convergence associated with the confluent flow around the dryline is shown to produce an upward moisture bulge, while surface heating and boundary layer mixing are responsible for the general deepening of the boundary layer. These processes produce favorable conditions for convection but the actual triggering of deep moist convection at specific locations along the dryline depends on localized forcing. Interaction of the primary dryline convergence boundary with horizontal convective rolls on its west side provides such localized forcing, while convective eddies on the immediate east side are suppressed by a downward mesoscale dryline circulation. A companion paper analyzes in detail the exact processes of convective initiation along this dryline. [ABSTRACT FROM AUTHOR]

This paper introduces the use of very large ensembles for detailed sensitivity analysis and applies this technique to study the sensitivity of model forecast rainfall to initial boundary layer and soil moisture fields for a particular case from the International H2O Project (IHOP_2002) field program. In total, an aggregate ensemble of over 12 000 mesoscale model forecasts are made, with each forecast having different perturbations of boundary layer moisture, boundary layer wind, or soil moisture. Sensitivity fields are constructed from this ensemble, producing detailed sensitivity fields of defined forecast functions to initial perturbations. This work is based on mesoscale model forecasts of convection of the 24 May 2002 IHOP case, which saw convective initiation along a dryline in western Texas as well as precipitation along a cold front. The large ensemble technique proves to be highly sensitive, and both significant and subtle connections between model state variables are revealed. A number of interesting sensitivity results are obtained. It is found that soil moisture and ABL moisture have opposite effects on the amount of precipitation along the dryline; that moisture on both the dry and moist sides of the dryline was equally important; and that some small perturbations were alone responsible for entire convective storm cells near the cold front, a result implying a high level of nonlinearity. These sensitivity analyses strongly indicate the importance of accurate low-level water vapor characterization for quantitative precipitation forecasting. [ABSTRACT FROM AUTHOR]

This study presents the simulation of the Madden–Julian oscillation (MJO) in the NCAR CCM3 using a modified Zhang–McFarlane convection parameterization scheme. It is shown that, with the modified scheme, the intraseasonal (20–80 day) variability in precipitation, zonal wind, and outgoing longwave radiation (OLR) is enhanced substantially compared to the standard CCM3 simulation. Using a composite technique based on the empirical orthogonal function (EOF) analysis, the paper demonstrates that the simulated MJOs are in better agreement with the observations than the standard model in many important aspects. The amplitudes of the MJOs in 850-mb zonal wind, precipitation, and OLR are comparable to those of the observations, and the MJOs show clearly eastward propagation from the Indian Ocean to the Pacific. In contrast, the simulated MJOs in the standard CCM3 simulation are weak and have a tendency to propagate westward in the Indian Ocean. Nevertheless, there remain several deficiencies that are yet to be addressed. The time period of the MJOs is shorter, about 30 days, compared to the observed time period of 40 days. The spatial scale of the precipitation signal is smaller than observed. Examination of convective heating from both deep and shallow convection and its relationship with moisture anomalies indicates that near the mature phase of the MJO, regions of shallow convection developing ahead of the deep convection coincide with regions of positive moisture anomalies in the lower troposphere. This is consistent with the recent observations and theoretical development that shallow convection helps to precondition the atmosphere for MJO by moistening the lower troposphere. Sensitivity tests are performed on the individual changes in the modified convection scheme. They show that both change of closure and use of a relative humidity threshold for the convection trigger play important roles in improving the MJO simulation. Use of the new closure leads to the eastward propagation of the MJO and increases the intensity of the MJO signal in the wind field, while imposing a relative humidity threshold enhances the MJO variability in precipitation. [ABSTRACT FROM AUTHOR]

Algorithms are now currently used for the retrieval of cloud optical thickness and droplet effective radius from multispectral radiance measurements. This paper extends their application to the retrieval of cloud droplet number concentration, cloud geometrical thickness, and liquid water path in shallow convective clouds, using an algorithm that was previously tested with airborne measurements of cloud radiances and validated against in situ measurements of the same clouds. The retrieval is based on a stratified cloud model of liquid water content and droplet spectrum. Radiance measurements in visible and near-infrared channels of the Moderate Resolution Imaging Spectroradiometer (MODIS), which is operated from the NASA platforms Terra and Aqua, are analyzed. Because of uncertainties in the simulation of the continental surface reflectance, the algorithm is presently limited to the monitoring of the microphysical structure of boundary layer clouds over the ocean. Two MODIS scenes of extended cloud fields over the North Atlantic Ocean trade wind region are processed. A transport and dispersion model (the Hybrid Single-Particle Lagrangian Integrated Trajectory Model, HYSPLIT4) is also used to characterize the origin of the air masses and hence their aerosol regimes. One cloud field formed in an air mass that was advected from southern Europe and North Africa. It shows high values of the droplet concentration when compared with the second cloud system, which developed in a more pristine environment. The more pristine case also exhibits a higher geometrical thickness and, thus, liquid water path, which counterbalances the expected cloud albedo increase of the polluted case. Estimates of cloud liquid water path are then compared with retrievals from the Special Sensor Microwave Imager (SSM/I). SSM/I-derived liquid water paths are in good agreement with the MODIS-derived values. [ABSTRACT FROM AUTHOR]

A realistic representation of the evolution of the dry convective boundary layer in mesoscale and large-scale atmospheric models has been an elusive goal for many years. In this paper the performance of a new mixing-length formulation for the dry convective boundary layer is evaluated in the context of the Coupled Ocean– Atmosphere Mesoscale Prediction System (COAMPS). In this new formulation, the mixing length is proportional to a time scale and to the square root of the turbulent kinetic energy. The model results are tested against observations from the Climate Impact of Changes in Land Use (CICLUS) field experiment in the south of Portugal. It is shown that COAMPS with the new formulation produces a more realistic simulation of the boundary layer growth. A data assimilation experiment performed with COAMPS shows that the improvements provided by the new formulation are significant, particularly in terms of the humidity vertical distribution. Finally, one-dimensional simulations are used to confirm that the new formulation provides more accurate results because of a more realistic representation of the entrainment and of the vertical mixing in general. [ABSTRACT FROM AUTHOR]

A new parameterization of shallow cumulus convection is presented. The parameterization consists of a mass flux scheme based on a buoyancy-sorting, entrainment–detrainment plume model. The mass flux scheme is coupled to a 1.5-order turbulence closure model with an entrainment closure for convective boundary layers. Model performance is verified using single-column-model simulations at relatively high vertical resolution of pure trade-cumulus convection and a cumulus to stratocumulus transition. Mixing rates, cloud cover, and vertical flux profiles, as deduced from previously published large-eddy simulation studies, are well reproduced by the parameterization. The model is used to demonstrate that height variations of lateral mixing rates can be successfully captured by a simple implementation of a buoyancy-sorting mechanism in the updraft cloud model. A companion paper describes the implementation of the scheme in a mesoscale model. [ABSTRACT FROM AUTHOR]

The impact of physical parameterizations on simulations of cloud-topped marine boundary layers is investigated using the fifth-generation Pennsylvania State University–NCAR Mesoscale Model (MM5). Three-month MM5 simulations of the northeast and southeast Pacific during June–August 1987 are boundary forced with time- varying ECMWF analyses. Runs with four planetary boundary layer (PBL) parameterizations already implemented in MM5 are compared with runs using new parameterizations of boundary layer turbulence and shallow cumulus convection (ShCu) described in a companion paper. Numerous modifications to the MM5 that allow it to be used as a regional climate model are described. The simulated 3-month mean shortwave cloud radiative forcing (SWCF) and vertical structure of cloud-topped boundary layers in the northeast Pacific are sensitive to the PBL/shallow convection schemes. All four current MM5 PBL schemes [the Blackadar, Medium-Range Forecast (MRF), Burk–Thompson, and Gayno–Seaman schemes] produce overly shallow boundary layers with excessive SWCF throughout this region, especially in the transition from stratocumulus to trade cumulus where their SWCF errors range from 130 (Gayno–Seaman) to 200 W m[sup -2] (MRF). These errors likely reflect inadequate vertical mixing by parameterized turbulence and shallow convection. The only shallow convection scheme available for MM5, the Grell scheme, was almost totally inactive in this region, so no shallow convection scheme was used for the above simulations. The Grenier– Bretherton scheme, which entrains more aggressively above stratocumulus-capped convective layers, has much better regional SWCF and vertical structure. Without shallow cumulus convection, the scheme still produces excessive cloud in the transition regions; the main impact of the ShCu parameterization is to remove this bias. With all schemes, the near-surface air has a cool, dry bias, and surface turbulent fluxes are somewhat larger than observed. Sensitivity studies show that the SWCF is sensitive to a halving of the cloud droplet concentration, to plausible uncertainties in parameterized penetrative mixing at cumulus cloud tops and in stratocumulus entrainment, and (in the coastal zone) to horizontal resolution. Southeast Pacific simulations show that the GB01+ShCu scheme can also accurately represent another subtropical boundary layer cloud regime. [ABSTRACT FROM AUTHOR]