Nonlinear Processes in Geophysics www.nonlin-processes-geophys.net 1023-5809 1607-7946 17 2 2010 10.5194/npg-17-187-2010 http://www.nonlin-processes-geophys.net/17/187/2010/ http://www.nonlin-processes-geophys.net/17/187/2010/npg-17-187-2010.html http://www.nonlin-processes-geophys.net/17/187/2010/npg-17-187-2010.pdf 187 200 2010-04-08 Estimating the diffusive heat flux across a stable interface forced by convective motions C. Chemel C. Staquet chantal.staquet@hmg.inpg.fr J.-P. Chollet NCAS-Weather, Centre for Atmospheric & Instrumentation Research, University of Hertfordshire, Hatfield, UK Laboratoire des Ecoulements Géophysiques et Industriels, CNRS/UJF/INPG, Grenoble, France Entrainment at the top of the convectively-driven boundary layer (CBL) is revisited using data from a high-resolution large-eddy simulation (LES). In the range of values of the bulk Richardson number RiB studied here (about 15–25), the entrainment process is mainly driven by the scouring of the interfacial layer (IL) by convective cells. We estimate the length and time scales associated with these convective cells by computing one-dimensional wavenumber and frequency kinetic energy spectra. Using a Taylor assumption, based upon transport by the convective cells, we show that the frequency and wavenumber spectra follow the Kolmogorov law in the inertial range, with the multiplicative constant being in good agreement with previous measurements in the atmosphere. We next focus on the heat flux at the top of the CBL, <IMG WIDTH="12" HEIGHT="22" ALIGN="MIDDLE" BORDER="0" src="http://www.nonlin-processes-geophys.net/17/187/2010/npg-17-187-2010-2.gif" ALT="$\mathcal{F}_i$">, which is parameterized in classical closure models for the entrainment rate we at the interface. We show that <IMG WIDTH="12" HEIGHT="22" ALIGN="MIDDLE" BORDER="0" src="http://www.nonlin-processes-geophys.net/17/187/2010/npg-17-187-2010-2.gif" ALT="$\mathcal{F}_i$"> can be computed exactly using the method proposed by Winters et al. (1995), from which the values of a turbulent diffusivity <IMG WIDTH="9" HEIGHT="9" ALIGN="BOTTOM" BORDER="0" src="http://www.nonlin-processes-geophys.net/17/187/2010/npg-17-187-2010-4.gif" ALT="$\mathcal{K}$"> across the IL can be inferred. These values are recovered by tracking particles within the IL using a Lagrangian stochastic model coupled with the LES. The relative difference between the Eulerian and Lagrangian values of <IMG WIDTH="9" HEIGHT="9" ALIGN="BOTTOM" BORDER="0" src="http://www.nonlin-processes-geophys.net/17/187/2010/npg-17-187-2010-4.gif" ALT="$\mathcal{K}$"> is found to be lower than 10%. A simple expression of we as a function of <IMG WIDTH="9" HEIGHT="9" ALIGN="BOTTOM" BORDER="0" src="http://www.nonlin-processes-geophys.net/17/187/2010/npg-17-187-2010-4.gif" ALT="$\mathcal{K}$"> is also proposed. Our results are finally used to assess the validity of the classical "first-order'' model for we. We find that, when RiB is varied, the values for we derived from the "first-order'' model with the exact computation of <IMG WIDTH="12" HEIGHT="22" ALIGN="MIDDLE" BORDER="0" src="http://www.nonlin-processes-geophys.net/17/187/2010/npg-17-187-2010-2.gif" ALT="$\mathcal{F}_i$"> agree to better than 10% with those computed directly from the LES (using its definition). The simple expression we propose appears to provide a reliable estimate of we for the largest values of RiB only. Betts, A K.: Reply to comment on the paper "Non-precipitating cumulus convection and its parameterization", Q J Roy Meteor Soc., 100, 469–471, 1974. Bretherton, C S., MacVean, M K., Bechtold, P., Chlond, A., Cotton, W R., Cuxart, J., Cuijpers, H., Khairoutdinov, M., Kosovic, B., Lewellen, D., Moeng, C.-H., Siebesma, P., Stevens, B., Stevens, D E., Sykes, I., and Wyant, M C.: An intercomparison of radiatively driven entrainment and turbulence in a smoke cloud, as simulated by different numerical models, Q J Roy Meteor Soc., 125, 391–423, 1999. Cai, X., Zhang, R., and Li, Y.: A large-eddy simulation and Lagrangian stochastic study of heavy particle dispersion in the convective boundary layer, Bound.-Lay Meteorol., 120, 413–435, 2006. Cai, X.-M. and Luhar, A K.: Fumigation of pollutants in and above the entrainment zone into a growing convective boundary layer: a large-eddy simulation, Atmos Environ., 36, 2997–3008, 2002. Champagne, F H., Friehe, C A., LaRue, J C., and Wyngaard, J C.: Flux measurements, flux estimation techniques, and fine-scale turbulence measurements in the unstable surface layer over land, J Atmos Sci., 34, 515–530, 1977. Chemel, C. and Staquet, C.: A formulation of convective entrainment in terms of mixing efficiency, J Fluid~Mech., 580, 169–178, 2007. Clarke, R H., Dyer, A J., Brook, R R., Reid, D G., and Troup, A J.: The Wangara experiment: Boundary layer data, Tech. Paper~19, CSIRO Atmospheric Research, Aspendale, Australia, 362~pp., 1971. Cohn, S A. and Angevine, W M.: Boundary layer height and entrainment zone thickness measured by lidars and wind-profiling radars, J Appl Meteorol., 39, 1233–1247, 2000. D'Asaro, E A., Winters, K B., and Lien, R.-C.: Lagrangian analysis of a convective mixed layer, J Geophys Res., 107(C5), 3040, \doi10.1029/2000JC000247, 2002. Deardorff, J W.: Stratocumulus-capped mixed layers derived from a three-dimensional model, Bound.-Lay Meteorol., 18, 495–527, 1980. Deardorff, J W., Willis, G E., and Stockton, B H.: Laboratory studies of the entrainment zone of a convectively mixed layer, J Fluid~Mech., 100, 41–64, 1980. Dosio, A., Vilá-Guerau de Arellano, J., Holtslag, A A M., and Builtjes, P J H.: Relating Eulerian and Lagrangian statistics for the turbulent dispersion in the atmospheric convective boundary layer, J Atmos Sci., 62, 1175–1191, 2005. Fedorovich, E., Conzemius, R., and Mironov, D.: Convective entrainment into a shear-free, linearly stratified atmosphere: bulk models reevaluated through large eddy simulations, J Atmos Sci., 61, 281–295, 2004. Fedorovich, E E. and Mironov, D V.: A model for a shear-free convective boundary layer with parameterized capping inversion structure, J Atmos Sci., 52, 83–95, 1995. Fernando, H J S.: Turbulent mixing in stratified fluids, Annu Rev Fluid~Mech., 23, 455–493, 1991. Gopalakrishnan, S G. and Avissar, R.: An LES study of the impacts of land surface heterogeneity on dispersion in the convective boundary layer, J Atmos Sci., 57, 352–371, 2000. Gregg, M C.: Diapycnal mixing in the thermocline: a review, J Geophys Res., 92(C5), 5249–5286, 1987. Hanna, S R.: Lagrangian and Eulerian time-scale relations in the daytime boundary layer, J Appl Meteorol., 20, 242–249, 1981. Hannoun, I A. and List, E J.: Turbulent mixing at a shear-free density interface, J Fluid~Mech., 189, 211–234, 1988. Hong, S.-Y., Noh, Y., and Dudhia, J.: A new vertical diffusion package with an explicit treatment of entrainment processes, Mon Weather~Rev., 134, 2318–2341, 2006. Hopfinger, E J.: Turbulence in stratified fluids: a review, J Geophys Res., 92(C5), 5287–5303, 1987. Hunt, J C R.: Diffusion in the stably stratified atmospheric boundary layer, J Clim Appl Meteorol., 24, 1187–1195, 1985. Hunt, J C R., Wray, A A., and Moin, P.: Eddies, stream and convergence zones in turbulent flows, Report CTR-S88, Center For Turbulence Research, Stanford, CA, USA, 193–208, 1988. Kaiser, R. and Fedorovich, E.: Turbulence spectra and dissipation rates in a wind tunnel model of the atmospheric convective boundary layer, J Atmos Sci., 55, 580–594, 1998. Kelly, M. and Wyngaard, J C.: Two-dimensional spectra in the atmospheric boundary layer, J Atmos Sci., 63, 3066–3070, 2006. Lamb, R G.: A numerical simulation of dispersion from an elevated point source in the convective planetary boundary layer, Atmos Environ., 12, 1297–1304, 1978. Lenschow, D H.: Observations of clear and cloud-capped convective boundary layers, and techniques for probing them, in: Buoyant Convection in Geophysical Flows, edited by: Plate, E J., Fedorovich, E., Viegas, D X., and Wyngaard, J C., Kluwer Academic Publishers, Dordrecht, The Netherlands, NATO ASI Series C: Mathematical and Physical Sciences, 513, 185–206, 1998. Lesieur, M. and Métais, O.: New trends in large-eddy simulations of turbulence, Annu Rev Fluid~Mech., 28, 45–82, 1996. Lewellen, D C. and Lewellen, W S.: Large-eddy boundary layer entrainment, J Atmos Sci., 55, 2645–2665, 1998. Lilly, D K.: Models of cloud-topped mixed layers under a strong inversion, Q J Roy Meteor Soc., 94, 292–309, 1968. Linden, P F.: The deepening of a mixed layer in a stratified fluid, J Fluid~Mech., 71, 385–405, 1975. Manins, P C. and Turner, J S.: The relation between the flux ratio and energy ratio in convectively mixed layers, Q J Roy Meteor Soc., 104, 39–44, 1978. Matthews, P C. and Cox, S M.: Pattern formation with a conservation law, Nonlinearity, 13, 1293–1320, 2000. Moeng, C.-H. and Wyngaard, J C.: Spectral analysis of large-eddy simulations of the convective boundary layer, J Atmos Sci., 45, 3573–3587, 1988. Noilhan, J. and Planton, S.: A simple parametrization of land surface processes for meteorological models, Mon Weather~Rev., 117, 536–549, 1989. Otte, M J. and Wyngaard, J C.: Stably stratified interfacial-layer turbulence from large-eddy simulation, J Atmos Sci., 58, 3424–3442, 2001. Peltier, L J., Wyngaard, J C., Khanna, S., and Brasseur, J G.: Spectra in the unstable surface layer, J Atmos Sci., 53, 49–61, 1996. Pleim, J E. and Xiu, A.: Development and testing of a surface flux and planetary boundary layer model for application in mesoscale models, J Appl Meteorol., 34, 16–32, 1995. Schmidt, H. and Schumann, U.: Coherent structure of the convective boundary layer derived from large-eddy simulations, J Fluid~Mech., 200, 511–562, 1989. Sorbjan, Z.: 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, 2015–2024, 1996. Stevens, B. and Lenschow, D H.: Observations, experiments, and large eddy simulation, B Am Meteorol Soc., 82, 283–294, 2001. Steyn, D G., Baldi, M., and Hoff, R M.: The detection of mixed layer depth and entrainment zone thickness from lidar backscatter profiles, J Atmos Ocean Tech., 16, 953–959, 1999. Sullivan, P P., Moeng, C.-H., Stevens, B., Lenschow, D H., and Mayor, S D.: Structure of the entrainment zone capping the convective atmospheric boundary layer, J Atmos Sci., 55, 3042–3064, 1998. Sun, J. and Wang, Y.: Effect of the entrainment flux ratio on the relationship between entrainment rate and convective Richardson number, Bound.-Lay Meteorol., 126, 237–247, 2008. Taylor, G I.: Diffusion by continuous movements, Proc London~Math Soc., 20, 196–212, 1921. Taylor, G I.: The spectrum of turbulence, P Roy Soc Lond A~Mat., 164, 476–490, 1938. Thomson, D J.: Criteria for the selection of stochastic models of particle trajectories in turbulent flows, J Fluid~Mech., 180, 529–556, 1987. Vinkovic, I., Aguirre, C., and Simoëns, S.: Large-eddy simulation and Lagrangian stochastic modeling of passive scalar dispersion in a turbulent boundary layer, J. Turbul., 7, No. 30, 2006. Weil, J C., Sullivan, P P., and Moeng, C.-H.: The use of large-eddy simulations in Lagrangian particle dispersion models, J Atmos Sci., 61, 2877–2887, 2004. Winters, K B. and D'Asaro, E A.: Diascalar flux and the rate of fluid mixing, J Fluid~Mech., 317, 179–193, 1996. Winters, K B., Lombard, P N., Riley, J J., and D'Asaro, E A.: Available potential energy and mixing in density-stratified fluids, J Fluid~Mech., 317, 115–128, 1995. Xue, M., Droegemeier, K K., and Wong, V.: The Advanced Regional Prediction System (ARPS) – A multi-scale non hydrostatic atmospheric simulation and prediction model. Part I: Model dynamics and verification, Meteorol Atmos Phys., 75, 161–193, 2000. Xue, M., Droegemeier, K K., Wong, V., Shapiro, A., Brewster, K., Carr, F., Weber, D., Liu, Y., and Wang, D.: The Advanced Regional Prediction System (ARPS) – A multi-scale non hydrostatic atmospheric simulation and prediction tool. Part II: Model physics and applications, Meteorol Atmos Phys., 76, 143–165, 2001. Yeung, P K.: Lagrangian investigations of turbulence, Annu Rev Fluid~Mech., 34, 115–142, 2002.