scholarly journals Modeling the Interaction between Cumulus Convection and Linear Gravity Waves Using a Limited-Domain Cloud System–Resolving Model

2008 ◽  
Vol 65 (2) ◽  
pp. 576-591 ◽  
Author(s):  
Zhiming Kuang

Abstract A limited-domain cloud system–resolving model (CSRM) is used to simulate the interaction between cumulus convection and two-dimensional linear gravity waves, a single horizontal wavenumber at a time. With a single horizontal wavenumber, soundings obtained from horizontal averages of the CSRM domain allow the large-scale wave equation to be evolved, and thereby its interaction with cumulus convection is modeled. It is shown that convectively coupled waves with phase speeds of 8–13 m s−1 can develop spontaneously in such simulations. The wave development is weaker at long wavelengths (>∼10 000 km). Waves at short wavelengths (∼2000 km) also appear weaker, but the evidence is less clear because of stronger influences from random perturbations. The simulated wave structures are found to change systematically with horizontal wavelength, and at horizontal wavelengths of 2000–3000 km they exhibit many of the basic features of the observed 2-day waves. The simulated convectively coupled waves develop without feedback from radiative processes, surface fluxes, or wave radiation into the stratosphere, but vanish when moisture advection by the large-scale waves is disabled. A similar degree of vertical tilt is found in the simulated convective heating at all wavelengths considered, consistent with observational results. Implications of these results to conceptual models of convectively coupled waves are discussed. In addition to being a useful tool for studying wave–convection interaction, the present approach also represents a useful framework for testing the ability of coarse-resolution CSRMs and single-column models in simulating convectively coupled waves.

2004 ◽  
Vol 132 (11) ◽  
pp. 2570-2594 ◽  
Author(s):  
T. N. Krishnamurti ◽  
J. Sanjay ◽  
A. K. Mitra ◽  
T. S. V. Vijaya Kumar

Abstract This paper addresses a procedure to extract error estimates for the physical and dynamical components of a forecast model. This is a two-step process in which contributions to the forecast tendencies from individual terms of the model equations are first determined using an elaborate bookkeeping of the forecast. The second step regresses these estimates of tendencies from individual terms of the model equations against the observed total tendencies. This process is executed separately for the entire horizontal and vertical transform grid points of a global model. The summary of results based on the corrections to the physics and dynamics provided by the regression coefficients highlights the component errors of the model arising from its formulation. This study provides information on geographical and vertical distribution of forecast errors contributed by features such as nonlinear advective dynamics, the rest of the dynamics, deep cumulus convection, large-scale condensation physics, radiative processes, and the rest of physics. Several future possibilities from this work are also discussed in this paper.


2010 ◽  
Vol 67 (4) ◽  
pp. 941-962 ◽  
Author(s):  
Zhiming Kuang

Abstract An approach is presented for the construction of linear response functions of a cumulus ensemble to large-scale temperature and moisture perturbations using a cloud system–resolving model (CSRM). A set of time-invariant, horizontally homogeneous, anomalous temperature and moisture tendencies is added, one at a time, to the forcing of the CSRM. By recording the departure of the equilibrium domain-averaged temperature and moisture profiles from those of a control experiment and through a matrix inversion, a sufficiently complete and accurate set of linear response functions is constructed for use as a parameterization of the cumulus ensemble around the reference mean state represented by the control experiment. This approach is applied to two different mean state conditions in which the CSRM, when coupled with 2D gravity waves, exhibits interestingly different behaviors. With a more strongly convecting mean state forced by the large-scale vertical velocity profile taken from the Tropical Ocean and Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE), spontaneous development of convectively coupled waves requires moisture variations above the boundary layer, whereas with a mean state of radiative–convective equilibrium (RCE) not forced by large-scale vertical advection, the development of convectively coupled waves is stronger and persists even when moisture variations above the boundary layer are removed. The linear response functions were able to reproduce these behaviors of the full CSRM with some quantitative accuracy. The linear response functions show that both temperature and moisture perturbations at a range of heights can regulate convective heating. The ability for convection to remove temperature anomalies, thus maintaining convective neutrality, decreases considerably from the lower troposphere to the middle and upper troposphere. It is also found that the response of convective heating to a lower tropospheric temperature anomaly is more top-heavy in the RCE case than in the TOGA COARE case. Comparing the linear response functions with the treatment of convection in an earlier simple model by the present author indicates general consistency, lending confidence that the instability mechanisms identified in that model provide the correct explanation to the instability seen in the CSRM simulations and the instability’s dependence on the mean state.


2010 ◽  
Vol 23 (7) ◽  
pp. 1837-1853 ◽  
Author(s):  
Eric Tromeur ◽  
William B. Rossow

Abstract To better understand the interaction between tropical deep convection and the Madden–Julian oscillation (MJO), tropical cloud regimes are defined by cluster analysis of International Satellite Cloud Climatology Project (ISCCP) cloud-top pressure—optical thickness joint distributions from the D1 dataset covering 21.5 yr. An MJO index based solely on upper-level wind anomalies is used to study variations of the tropical cloud regimes. The MJO index shows that MJO events are present almost all the time; instead of the MJO event being associated with “on or off” deep convection, it is associated with weaker or stronger mesoscale organization of deep convection. Atmospheric winds and humidity from NCEP–NCAR reanalysis 1 are used to characterize the large-scale dynamics of the MJO; the results show that the large-scale motions initiate an MJO event by moistening the lower troposphere by horizontal advection. Increasingly strong convection transports moisture into the upper troposphere, suggesting a reinforcement of the convection itself. The change of convection organization shown by the cloud regimes indicates a strong interaction between the large-scale circulation and deep convection. The analysis is extended to the complete atmospheric diabatic heating by precipitation, radiation, and surface fluxes. The wave organizes stronger convective heating of the tropical atmosphere, which results in stronger winds, while there is only a passive response of the surface, directly linked to cloud radiative effects. Overall, the results suggest that an MJO event is an amplification of large-scale wave motions by stronger convective heating, which results from a dynamic reorganization of scattered deep convection into more intense mesoscale systems.


2010 ◽  
Vol 10 (13) ◽  
pp. 6371-6389 ◽  
Author(s):  
S. S. Lee ◽  
J. E. Penner

Abstract. The response of a case of thin, warm marine-boundary-layer (MBL) clouds to preindustrial (PI) and present-day (PD) conditions is simulated by a cloud-system resolving model (CSRM). Here, both the aerosol conditions and environmental conditions match those of a general circulation model (GCM). The environmental conditions are characterized by the initial condition and the large-scale forcings of humidity and temperature, as well as the surface fluxes. The response of the CSRM is compared to that simulated by the GCM. The percentage increase of liquid-water path (LWP) due to a change from the PI to PD conditions is ~3 times larger in the CSRM than that in the GCM due to the formation of cumulus clouds. The formation of cumulus clouds is controlled by a larger increase in the surface latent-heat (LH) flux in the PD environment than in the PI environment rather than by the change in aerosols. However, the aerosol increase from the PI to PD level determines the LWP response in the stratocumulus clouds, while the impacts of changes in environmental conditions are negligible for stratocumulus clouds. The conversion of cloud liquid to rain through autoconversion and accretion plays a negligible role in the CSRM in the response to aerosols, whereas it plays a role that is as important as condensation in the GCM. Also, it is notable that the explicit simulation of microphysics in the CSRM leads to a smaller LWP in the CSRM than that in the GCM using heavily parameterized microphysics for stratocumulus clouds. The smaller LWP in the CSRM is closer to an observed LWP than the LWP in the GCM for stratocumulus clouds. Supplementary simulations show that increasing aerosols increase the sensitivity of the cloud responses to the PI and PD environmental conditions. They also show that aerosol effects on clouds depend on the cloud type. The LWP of warm cumulus clouds is more sensitive to aerosols than that of stratocumulus clouds.


2006 ◽  
Vol 134 (12) ◽  
pp. 3735-3752 ◽  
Author(s):  
Robert G. Fovell ◽  
Gretchen L. Mullendore ◽  
Seung-Hee Kim

Abstract Simulations of a typical midlatitude squall line were used to investigate a mechanism for discrete propagation, defined as convective initiation ahead of an existing squall line leading to a faster propagation speed for the storm complex. Radar imagery often shows new cells appearing in advance of squall lines, suggesting a causal relationship and prompting the search for an “action-at-a-distance” mechanism to explain the phenomenon. In the simulations presented, the identified mechanism involves gravity waves of both low and high frequency generated in response to the latent heating, which subsequently propagate out ahead of the storm. The net result of the low-frequency response, combined with surface fluxes and radiative processes, was a cooler and more moist lower troposphere, establishing a shallow cloud deck extending ahead of the storm. High-frequency gravity waves, excited in response to fluctuations in convective activity in the main storm, were subsequently ducted by the storm’s own upper-tropospheric forward anvil outflow. These waves helped positively buoyant cumulus clouds to occasionally form in the deck. A fraction of these clouds persisted long enough to merge with the main line, invigorating the parent storm. Discrete propagation occurred when clouds developed into deep convection prior to merger, weakening the parent storm. The ducting conditions, as diagnosed with the Scorer parameter, are shown to be sensitive to vertical wind shear and radiation, but not to the microphysical parameterization or simulation geometry.


2007 ◽  
Vol 64 (4) ◽  
pp. 1210-1229 ◽  
Author(s):  
Stefan N. Tulich ◽  
David A. Randall ◽  
Brian E. Mapes

Abstract 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.


2012 ◽  
Vol 69 (1) ◽  
pp. 281-302 ◽  
Author(s):  
Boualem Khouider ◽  
Ying Han ◽  
Joseph A. Biello

Abstract Convective momentum transport (CMT) is the process of vertical transport of horizontal momentum by convection onto the environmental flow. The significance of CMT from mesoscale to synoptic- and planetary-scale organized cumulus convection has been established by various theoretical and observational studies. A new strategy mimicking the effect of unresolved mesoscale circulation based on the weak temperature gradient (WTG) approximation with a Gaussian profile to redistribute the heating due to parameterized cumulus convection at the subgrid scale is adopted here to construct a CMT parameterization for general circulation models (GCMs). Two main regimes of CMT are considered: an upscale squall-line regime and a downscale non-squall-line regime. An exponential probability distribution is used to select which of these two effects is active, conditional on the state of the large-scale shear. The shear itself is used as a measure of the persistence of mesoscale organized circulation due to the presence or not of tilted deep convective heating with lagged stratiform anvils. The CMT model is tested in the simple case of the multicloud model of Khouider and Majda, used here as a toy GCM. Numerical simulations are performed here for the simple case without rotation, in a parameter regime where the multicloud model exhibits packets of convectively coupled gravity waves moving in one direction, at 17 m s−1, and planetary-scale wave envelopes moving in the opposite direction, at 4–6 m s−1, reminiscent of the Madden–Julian oscillation (MJO) and the associated embedded synoptic-scale superclusters. The results herein show that the inclusion of CMT intensifies both the synoptic-scale convectively coupled waves and the manifestation of planetary-scale waves in the multicloud model. This provides evidence that the present CMT model captures the essence of the physical mechanism through which kinetic energy is transferred from the subgrid-scale mesoscale circulation to the large-scale/resolved motion. Sensitivity simulations showed that two key parameters for the CMT parameterization are the relative strength of the parameterized stratiform anvils and the dimensional threshold used in the exponential distribution for the cumulus friction and the upscale CMT forcing resulting from organized subgrid mesoscale circulation.


Author(s):  
Na Li ◽  
Baofeng Jiao ◽  
Lingkun Ran ◽  
Zongting Gao ◽  
Shouting Gao

AbstractWe investigated the influence of upstream terrain on the formation of a cold frontal snowband in Northeast China. We conducted numerical sensitivity experiments that gradually removed the upstream terrain and compared the results with a control experiment. Our results indicate a clear negative effect of upstream terrain on the formation of snowbands, especially over large-scale terrain. By thoroughly examining the ingredients necessary for snowfall (instability, lifting and moisture), we found that the release of mid-level conditional instability, followed by the release of low-level or near surface instabilities (inertial instability, conditional instability or conditional symmetrical instability), contributed to formation of the snowband in both experiments. The lifting required for the release of these instabilities was mainly a result of frontogenetic forcing and upper gravity waves. However, the snowband in the control experiment developed later and was weaker than that in the experiment without upstream terrain. Two factors contributed to this negative topographic effect: (1) the mountain gravity waves over the upstream terrain, which perturbed the frontogenetic circulation by rapidly changing the vertical motion and therefore did not favor the release of instabilities in the absence of persistent ascending motion; and (2) the decrease in the supply of moisture as a result of blocking of the upstream terrain, which changed both the moisture and instability structures leeward of the mountains. A conceptual model is presented that shows the effects of the instabilities and lifting on the development of cold frontal snowbands in downstream mountains.


2008 ◽  
Vol 8 (10) ◽  
pp. 2811-2832 ◽  
Author(s):  
K. Zhang ◽  
H. Wan ◽  
M. Zhang ◽  
B. Wang

Abstract. The radioactive species radon (222Rn) has long been used as a test tracer for the numerical simulation of large scale transport processes. In this study, radon transport experiments are carried out using an atmospheric GCM with a finite-difference dynamical core, the van Leer type FFSL advection algorithm, and two state-of-the-art cumulus convection parameterization schemes. Measurements of surface concentration and vertical distribution of radon collected from the literature are used as references in model evaluation. The simulated radon concentrations using both convection schemes turn out to be consistent with earlier studies with many other models. Comparison with measurements indicates that at the locations where significant seasonal variations are observed in reality, the model can reproduce both the monthly mean surface radon concentration and the annual cycle quite well. At those sites where the seasonal variation is not large, the model is able to give a correct magnitude of the annual mean. In East Asia, where radon simulations are rarely reported in the literature, detailed analysis shows that our results compare reasonably well with the observations. The most evident changes caused by the use of a different convection scheme are found in the vertical distribution of the tracer. The scheme associated with weaker upward transport gives higher radon concentration up to about 6 km above the surface, and lower values in higher altitudes. In the lower part of the atmosphere results from this scheme does not agree as well with the measurements as the other scheme. Differences from 6 km to the model top are even larger, although we are not yet able to tell which simulation is better due to the lack of observations at such high altitudes.


1997 ◽  
Vol 15 (8) ◽  
pp. 1048-1056 ◽  
Author(s):  
R. L. Balthazor ◽  
R. J. Moffett

Abstract. A global coupled thermosphere-ionosphere-plasmasphere model is used to simulate a family of large-scale imperfectly ducted atmospheric gravity waves (AGWs) and associated travelling ionospheric disturbances (TIDs) originating at conjugate magnetic latitudes in the north and south auroral zones and subsequently propagating meridionally to equatorial latitudes. A 'fast' dominant mode and two slower modes are identified. We find that, at the magnetic equator, all the clearly identified modes of AGW interfere constructively and pass through to the opposite hemisphere with unchanged velocity. At F-region altitudes the 'fast' AGW has the largest amplitude, and when northward propagating and southward propagating modes interfere at the equator, the TID (as parameterised by the fractional change in the electron density at the F2 peak) increases in magnitude at the equator. The amplitude of the TID at the magnetic equator is increased compared to mid-latitudes in both upper and lower F-regions with a larger increase in the upper F-region. The ionospheric disturbance at the equator persists in the upper F-region for about 1 hour and in the lower F-region for 2.5 hours after the AGWs first interfere, and it is suggested that this is due to enhancements of the TID by slower AGW modes arriving later at the magnetic equator. The complex effects of the interplays of the TIDs generated in the equatorial plasmasphere are analysed by examining neutral and ion winds predicted by the model, and are demonstrated to be consequences of the forcing of the plasmasphere along the magnetic field lines by the neutral air pressure wave.


Sign in / Sign up

Export Citation Format

Share Document