Abstract
Both observational and numerical studies of the convective boundary layer (CBL) have demonstrated that when surface heat fluxes are small and mean wind shear is strong, convective updrafts tend to organize into horizontal rolls aligned within 10–20\(^\circ \) of the geostrophic wind direction. However, under large surface heat fluxes and weak to negligible shear, convection tends to organize into open cells, similar to turbulent Rayleigh-Bénard convection. Using a suite of 14 large-eddy simulations (LES) spanning a range of \(-z_i/L\) between zero (neutral) and 1041 (highly convective), where \(z_i\) is the CBL depth and L is the Obukhov length, the transition between roll- and cellular-type convection is investigated systematically for the first time using LES. Mean vertical profiles including velocity variances and turbulent transport efficiencies, as well the “roll factor,” which characterizes the rotational symmetry of the vertical velocity field, indicate the transition occurs gradually over a range of \(-z_i/L\); however, the most significant changes in vertical profiles and CBL organization occur from near-neutral conditions up to about \(-z_i/L \approx \) 15–20. Turbulent transport efficiencies and quadrant analysis are used to characterize the turbulent transport of momentum and heat with increasing \(-z_i/L\). It is found that turbulence transports heat efficiently from weakly to highly convective conditions; however, turbulent momentum transport becomes increasingly inefficient as \(-z_i/L\) increases.
Similar content being viewed by others
Notes
In particular, Grossman (1982) suggested that rolls and cells coexist with rolls dominating for \(5.0 \le -z_i/L \le 7.3\) and cells dominating for \(7.3 \le -z_i/L \le 21.4\).
References
Adrian R (2007) Hairpin vortex organization in wall turbulence. Phys Fluids 19(041):301
Agee E, Chen T, Dowell K (1973) A review of mesoscale cellular convection. Bull Am Meteorol Soc 54:1004–1012
Arya S, Wyngaard J (1975) Effect of baroclinicity on wind profiles and the geostrophic drag law for the convective planetary boundary layer. J Atmos Sci 32(4):767–778
Asai T (1970a) Stability of a plane parallel flow with variable vertical shear and unstable stratification. J Meteorol Soc Jpn 48:129–139
Asai T (1970b) Three-dimensional features of thermal convection in a plane Couette flow. J Meteorol Soc Jpn 48:18–29
Asai T (1972) Thermal instability of a shear flow turning the direction with height. J Meteorol Soc Jpn 50:525–532
Atkinson B, Zhang J (1996) Mesoscale shallow convection in the atmosphere. Rev Geophys 34(4):403–431
Bou-Zeid E, Meneveau C, Parlange M (2005) A scale-dependent Lagrangian dynamic model for large eddy simulation of complex turbulent flows. Phys Fluids 17(2):025–105
Brown RA (1972) On the inflection point instability of a stratified Ekman boundary layer. J Atmos Sci 29(5):850–859
Brown RA (1980) Longitudinal instabilities and secondary flows in the planetary boundary layer: a review. Rev Geophys 18(3):683–697
Brümmer B (1999) Roll and cell convection in wintertime Arctic cold-air outbreaks. J Atmos Sci 56(15):2613–2636
Businger J (1973) A note on free convection. Boundary-Layer Meteorol 4(1–4):323–326
Cancelli DM, Chamecki M, Dias NL (2014) A large-eddy simulation study of scalar dissimilarity in the convective atmospheric boundary layer. J Atmos Sci 71(1):3–15
Canuto C, Hussaini MY, Quarteroni AM, Zhang TA (2012) Spectral methods in fluid dynamics. Springer, Berlin 539 pp
Carlotti P (2002) Two-point properties of atmospheric turbulence very close to the ground: comparison of a high resolution LES with theoretical models. Boundary-Layer Meteorol 104(3):381–410
de Roode SR, Duynkerke PG, Jonker HJ (2004) Large-eddy simulation: how large is large enough? J Atmos Sci 61(4):403–421
Deardorff J (1972) Numerical investigation of neutral and unstable planetary boundary layers. J Atmos Sci 29:91–115
Dupont S, Patton E (2012) Momentum and scalar transport within a vegetation canopy following atmospheric stability and seasonal canopy changes: the CHATS experiment. Atmos Chem Phys 12(13):5913–5935
Esau I, Davy R, Outten S, Tyuryakov S, Zilitinkevich S (2013) Structuring of turbulence and its impact on basic features of Ekman boundary layers. Nonlinear Process Geophys 20(4):589–604
Faller AJ (1963) An experimental study of the instability of the laminar Ekman boundary layer. J Fluid Mech 15(163):560–576
Faller AJ (1965) Large eddies in the atmospheric boundary layer and their possible role in the formation of cloud rows. J Atmos Sci 22(2):176–184
Finnigan J (2000) Turbulence in plant canopies. Annu Rev Fluid Mech 32(1):519–571
Germano M, Piomelli U, Moin P, Cabot W (1991) A dynamic subgrid-scale eddy viscosity model. Phys Fluids A 3:1760
Gibbs JA, Fedorovich E (2014) Comparison of convective boundary layer velocity spectra retrieved from large-eddy-simulation and weather research and forecasting model data. J Appl Meteorol Clim 53(2):377–394
Grossman RL (1982) An analysis of vertical velocity spectra obtained in the BOMEX fair-weather, trade-wind boundary layer. Boundary-Layer Meteorol 23(3):323–357
Grossman RL (1984) Bivariate conditional sampling of moisture flux over a tropical ocean. J Atmos Sci 41(22):3238–3254
Hastings R (2009) Multi-platform analysis of morning transition and evolution of shallow mesoscale convection on 21 June 2002. Master’s thesis, The Pennsylvania State University
Head M, Bandyopadhyay P (1981) New aspects of turbulent boundary-layer structure. J Fluid Mech 107:297–338
Holland JZ (1973) A statistical method for analyzing wave shapes and phase relationships of fluctuating geophysical variables. J Phys Ocean 3(1):139–155
Hunt J, Carlotti P (2001) Statistical structure at the wall of the high reynolds number turbulent boundary layer. Flow Turbul Combust 66(4):453–475
Hunt JC, Morrison JF (2000) Eddy structure in turbulent boundary layers. Eur J Mech B Fluid 19(5):673–694
Hutchins N, Marusic I (2007) Large-scale influences in near-wall turbulence. Philos Trans R Soc Lond A 365(1852):647–664
Kaimal J, Wyngaard J, Haugen D, Coté O, Izumi Y, Caughey S, Readings C (1976) Turbulence structure in the convective boundary layer. J Atmos Sci 33(11):2152–2169
Kang HS, Meneveau C (2002) Universality of large eddy simulation model parameters across a turbulent wake behind a heated cylinder. J Turbul 3:26–26
Khanna S, Brasseur JG (1998) Three-dimensional buoyancy-and shear-induced local structure of the atmospheric boundary layer. J Atmos Sci 55(5):710–743
Kleissl J, Kumar V, Meneveau C, Parlange MB (2006) Numerical study of dynamic Smagorinsky models in large-eddy simulation of the atmospheric boundary layer: validation in stable and unstable conditions. Water Resour Res. doi:10.1029/2005WR004685
Kline S, Reynolds W, Schraub F, Runstadler P (1967) The structure of turbulent boundary layers. J Fluid Mech 30(04):741–773
Knaepen B, Debliquy O, Carati D (2002) Subgrid-scale energy and pseudo pressure in large-eddy simulation. Phys Fluids 14(12):4235–4241
Krishnamurti R (1970) On the transition to turbulent convection. Part 1. The transition from two-to three-dimensional flow. J Fluid Mech 42(2):295–307
Kuettner JP (1971) Cloud bands in the earth’s atmosphere: observations and theory. Tellus 23(4–5):404–426
Kumar V, Kleissl J, Meneveau C, Parlange MB (2006) Large-eddy simulation of a diurnal cycle of the atmospheric boundary layer: atmospheric stability and scaling issues. Water Resour Res 42(6):W06D09
Kuo H (1963) Perturbations of plane Couette flow in stratified fluid and origin of cloud streets. Phys Fluids 6(2):195–211
Lanotte AS, Mazzitelli IM (2013) Scalar turbulence in convective boundary layers by changing the entrainment flux. J Atmos Sci 70(1):248–265
LeMone M (1973) The structure and dynamics of horizontal roll vortices in the planetary boundary layer. J Atmos Sci 30(6):1077–1091
LeMone M (1976) Modulation of turbulence energy by longitudinal rolls in an unstable planetary boundary layer. J Atmos Sci 33:1308–1320
LeMone MA, Pennell WT (1976) The relationship of trade wind cumulus distribution to subcloud layer fluxes and structure. Mon Weather Rev 104(5):524–539
LeMone MA, Chen F, Tewari M, Dudhia J, Geerts B, Miao Q, Coulter RL, Grossman RL (2010) Simulating the IHOP_2002 fair-weather CBL with the WRF-ARW-Noah modeling system. Part ii: structures from a few kilometers to 100 km across. Mon Weather Rev 138(3):745–764
Lenschow D, Wyngaard JC, Pennell WT (1980) Mean-field and second-moment budgets in a baroclinic, convective boundary layer. J Atmos Sci 37(6):1313–1326
Li D, Bou-Zeid E (2011) Coherent structures and the dissimilarity of turbulent transport of momentum and scalars in the unstable atmospheric surface layer. Boundary-Layer Meteorol 140(2):243–262
Li D, Katul GG, Bou-Zeid E (2012) Mean velocity and temperature profiles in a sheared diabatic turbulent boundary layer. Phys Fluids 24(10):105–105
Lilly DK (1966) On the instability of Ekman boundary flow. J Atmos Sci 23(5):481–494
Mahrt L (1998) Flux sampling errors for aircraft and towers. J Atmos Ocean Technol 15(2):416–429
Mahrt L (2010) Computing turbulent fluxes near the surface: needed improvements. Agric For Meteorol 150(4):501–509
Marusic I, Mathis R, Hutchins N (2010) Predictive model for wall-bounded turbulent flow. Science 329(5988):193–196
Mathis R, Hutchins N, Marusic I (2009) Large-scale amplitude modulation of the small-scale structures in turbulent boundary layers. J Fluid Mech 628:311–337
Melfi S, Palm SP (2012) Estimating the orientation and spacing of midlatitude linear convective boundary layer features: cloud streets. J Atmos Sci 69(1):352–364
Mellado JP, van Heerwaarden CC, Garcia JR (2016) Near-surface effects of free atmosphere stratification in free convection. Boundary-Layer Meteorol 159(1):69–95
Meneveau C, Lund TS, Cabot WH (1996) A Lagrangian dynamic subgrid-scale model of turbulence. J Fluid Mech 319:353–385
Miura Y (1986) Aspect ratios of longitudinal rolls and convection cells observed during cold air outbreaks. J Atmos Sci 43(1):26–39
Moeng C, Sullivan P (1994) A comparison of shear-and buoyancy-driven planetary boundary layer flows. J Atmos Sci 51(7):999–1022
Moeng CH (1984) A large-eddy-simulation model for the study of planetary boundary-layer turbulence. J Atmos Sci 41(13):2052–2062
Nieuwstadt FTM, Mason PJ, Moeng CH, Schumann U (1993) Large-eddy simulation of the convective boundary layer: a comparison of four computer codes. Turbulent shear flows 8. Springer, Berlin, pp 343–367
Park SB, Baik JJ (2014) Large-eddy simulations of convective boundary layers over flat and urban-like surfaces. J Atmos Sci 71:18801892
Patton EG, Sullivan PP, Shaw RH, Finnigan JJ, Weil JC (2015) Atmospheric stability influences on coupled boundary-layer-canopy turbulence. J Atmos Sci. doi:10.1175/JAS-D-15-0068.1
Robinson SK (1991) Coherent motions in the turbulent boundary layer. Annu Rev Fluid Mech 23(1):601–639
Salesky ST, Chamecki M (2012) A similarity model of subfilter-scale energy for large-eddy simulations of the atmospheric boundary layer. Boundary-Layer Meteorol 145(1):69–91
Schmidt H, Schumann U (1989) Coherent structure of the convective boundary layer derived from large-eddy simulations. J Fluid Mech 200:511–562
Shah S, Bou-Zeid E (2014) Very-large-scale motions in the atmospheric boundary layer educed by snapshot proper orthogonal decomposition. Boundary-Layer Meteorol 153(3):355–387
Sharma V, Parlange M, Calaf M (2016) Perturbations to the spatial and temporal characteristics of the diurnally-varying atmospheric boundary layer due to an extensive wind farm. Boundary-Layer Meteorol. doi:10.1007/s10546-016-0195-0
Sorbjan Z (1996) Effects caused by varying the strength of the capping inversion based on a large eddy simulation model of the shear-free convective boundary layer. J Atmos Sci 53(14):2015–2024
Sullivan PP, Patton EG (2011) The effect of mesh resolution on convective boundary layer statistics and structures generated by large-eddy simulation. J Atmos Sci 68(10):2395–2415
Trier SB, Parsons DB, Clark JH (1991) Environment and evolution of a cold-frontal mesoscale convective system. Mon Weather Rev 119(10):2429–2455
Van de Boer A, Moene A, Graf A, Schüttemeyer D, Simmer C (2014) Detection of entrainment influences on surface-layer measurements and extension of Monin-Obukhov similarity theory. Boundary-Layer Meteorol 152(1):19–44
Wakimoto RM, Wilson JW (1989) Non-supercell tornadoes. Mon Weather Rev 117(6):1113–1140
Wallace JM (2016) Quadrant analysis in turbulence research: history and evolution. Annu Rev Fluid Mech 48:131–158
Wallace JM, Eckelmann H, Brodkey RS (1972) The wall region in turbulent shear flow. J Fluid Mech 54(01):39–48
Weckwerth T, Wilson J, Wakimoto R, Crook N (1997) Horizontal convective rolls: determining the environmental conditions supporting their existence and characteristics. Mon Weather Rev 125(4):505–526
Weckwerth T, Horst T, Wilson J (1999) An observational study of the evolution of horizontal convective rolls. Mon Weather Rev 127(9):2160–2179
Weckwerth TM (2000) The effect of small-scale moisture variability on thunderstorm initiation. Mon Weather Rev 128(12):4017–4030
Weckwerth TM, Wilson JW, Wakimoto RM (1996) Thermodynamic variability within the convective boundary layer due to horizontal convective rolls. Mon Weather Rev 124(5):769–784
Willmarth W, Lu S (1972) Structure of the Reynolds stress near the wall. J Fluid Mech 55(01):65–92
Wilson DK (1996) Empirical orthogonal function analysis of the weakly convective atmospheric boundary layer. part i: Eddy structures. J Atmos Sci 53(6):801–823
Wilson DK, Wyngaard JC (1996) Empirical orthogonal function analysis of the weakly convective atmospheric boundary layer. Part ii: eddy energetics. J Atmos Sci 53(6):824–841
Wilson JW, Foote GB, Cook NA, Fankhauser JC, Wade CG, Tuttle JD, Mueller CK, Krueger SK (1992) The role of boundary-layer convergence zones and horizontal rolls in the initiation of thunderstorms: a case study. Mon Weather Rev 120(9):1785–1815
Wyngaard JC (2010) Turbulence in the atmosphere. Cambridge University Press, New York 393 pp
Wyngaard JC, Moeng CH (1992) Parameterizing turbulent diffusion through the joint probability density. Boundary-Layer Meteorol 60(1–2):1–13
Xue M, Martin WJ (2006) A high-resolution modeling study of the 24 May 2002 dryline case during IHOP. Part I: numerical simulation and general evolution of the dryline and convection. Mon Weather Rev 134(1):149–171
Young G, Kristovich D, Hjelmfelt M, Foster R (2002) Rolls, streets, waves, and more: a review of quasi-two-dimensional structures in the atmospheric boundary layer. Bull Am Meteorol Soc 83(7):997–1001
Zilitinkevich SS, Hunt JCR, Esau IN, Grachev AA, Lalas DP, Akylas E, Tombrou M, Fairall CW, Fernando HJS, Baklanov AA, Joffre SM (2006) The influence of large convective eddies on the surface-layer turbulence. Q J R Meteorol Soc 132(618):1426–1456
Acknowledgements
The authors thank Margaret LeMone and three other anonymous reviewers for their insightful comments. Simulations were run on the Guillimin cluster at McGill University, managed by Calcul Quèbec and Compute Canada. The operation of this supercomputer is funded by the Canadian Foundation for Innovation (CFI), Ministère de l’Économie, de l’Innovation et des Exportations du Quèbec (MEIE), RMGA and the Fonds de Recherche du Quèbec - Nature et Technologies (FRQ-NT).
Author information
Authors and Affiliations
Corresponding author
Appendix 1: Grid convergence tests
Appendix 1: Grid convergence tests
In order to investigate the effects of grid resolution on results from our LES, we performed several additional simulations on a \(12\times 12\times 2\, \text {km}^3\) domain at six different grid resolutions ranging from \(96^3\) to \(320^3\). The simulations were configured following the case described in Sullivan and Patton (2011), where \(U_g = 1\) m s\(^{-1}\) and \(Q_0 = 0.24 \) K m s\(^{-1}\). Simulations were run for approximately 25 large-eddy turnover times \(T_{\ell } = z_i/w_*\), and averages calculated over the time interval \(10T_\ell \)–\(25T_\ell \). Parameters of the simulations for the grid convergence tests can be found in Table 2.
Mean profiles from each grid resolution can be found in Fig 13, where the mean velocity profile is displayed in Fig 13a, the mean temperature profile in Fig 13b, the total (resolved +SGS) heat flux in Fig 13c, the total vertical velocity variance \(\langle w_\mathrm{tot}^2 \rangle /w_*^2 = \langle \widetilde{w}^{\prime 2} + 2e/3 \rangle /w_*^2\) in Fig 13d (where e is the SGS energy), the total horizontal velocity variance \(\langle u_h^2 \rangle /w_*^2 = \langle \widetilde{u}^{\prime 2} + \widetilde{v}^{\prime 2} + 4e/3 \rangle /w_*^2\) in Fig 13e, and the resolved vertical velocity skewness \(\langle w^3 \rangle / \langle w^2 \rangle ^{3/2}\) in Fig 13f. The SGS energy was estimated using the similarity model of Knaepen et al. (2002), via \(e = \frac{1}{2} C_\mathrm{sim} L_{ii}\), where \(L_{ii}\) is the trace of the Leonard stress tensor \(L_{ij} = \overline{\widetilde{u}_i \widetilde{u}_j} - \overline{\widetilde{u}}_i \overline{\widetilde{u}}_j\) and \(C_\mathrm{sim} = [ (\overline{\varDelta }/\widetilde{\varDelta })^{2/3} -1 ]^{-1}\) is a similarity coefficient, where \(\overline{\varDelta }\) and \(\widetilde{\varDelta }\) denote the test and grid filters, respectively. Note that \(C_\mathrm{sim} = 1.7\) when \(\overline{\varDelta } = 2 \widetilde{\varDelta }\). The more sophisticated model of Salesky and Chamecki (2012) that accounts for the shape of the energy spectrum in the ABL yields similar results, since the assumption of a \(k^{-5/3}\) spectrum at scale \(\overline{\varDelta }\) is reasonable under highly convective conditions.
As one can see from Fig. 13a, b, mean vertical profiles of first moments (\(\langle U \rangle \) and \(\langle \varTheta \rangle \)) are in good agreement between all resolutions considered here. Profiles of the total heat flux, seen in Fig. 13c are similar for all resolutions in the lower part of the CBL; the differences become most noticeable near the inversion. Profiles of velocity variances and vertical velocity skewness are displayed in Fig. 13d–f. Here it should be emphasized that while the LES directly provides the resolved velocity and temperature fields, and the SGS heat flux (which is used in the numerical integration) through an SGS model, an a posteriori model for the SGS energy must be used to estimate the SGS contribution to \(\langle w^2\rangle _\mathrm{tot}\) and \(\langle u_h^2 \rangle _\mathrm{tot}\), since the SGS energy is included in the modified pressure, and is not necessarily modeled in LES. Thus the convergence of the total horizontal and vertical velocity variances will depend both on grid resolution and on the model for e employed. The SGS contribution to the w-skewness is unknown, due to the lack of a model to estimate it.
Regardless of the uncertainty in the true SGS contribution to the velocity variances and vertical velocity skewness, we find that, using the present approach to estimate e, all moments are nearly grid-independent for the \(256^3\) grid, or a characteristic filter width of \(\varDelta _f = 25.8\) m. Note that we find convergence at a larger value of \(\varDelta _f\) than Sullivan and Patton (2011) (they found that most mean profiles were well converged for their \(256^3\) case on a \(5120\times 5120\times 2048\) m domain, or \(\varDelta _f = 19.3\) m). We attribute the convergence at larger values of \(\varDelta _f\) in our case to the difference in the subgrid model employed in our simulations vs. Sullivan and Patton (2011). We use the Lagrangian scale-dependent dynamic model (Bou-Zeid et al. 2005), whereas Sullivan and Patton (2011) used a one-equation model for the SGS energy (e.g. Moeng 1984).
Rights and permissions
About this article
Cite this article
Salesky, S.T., Chamecki, M. & Bou-Zeid, E. On the Nature of the Transition Between Roll and Cellular Organization in the Convective Boundary Layer. Boundary-Layer Meteorol 163, 41–68 (2017). https://doi.org/10.1007/s10546-016-0220-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10546-016-0220-3