Advanced Search
Article Contents

Shallow Convection Dataset Simulated by Three Different Large Eddy Models


doi:  10.1007/s00376-023-3106-6

  • Shallow convection plays an important role in transporting heat and moisture from the near-surface to higher altitudes, yet its parameterization in numerical models remains a great challenge, partly due to the lack of high-resolution observations. This study describes a large eddy simulation (LES) dataset for four shallow convection cases that differ primarily in inversion strength, which can be used as a surrogate for real data. To reduce the uncertainty in LES modeling, three different large eddy models were used, including SAM (System for Atmospheric Modeling), WRF (Weather Research and Forecasting model), and UCLA-LES. Results show that the different models generally exhibit similar behavior for each shallow convection case, despite some differences in the details of the convective structure. In addition to grid-averaged fields, conditionally sampled variables, such as in-cloud moisture and vertical velocity, are also provided, which are indispensable for calculation of the entrainment/detrainment rate. Considering the essentiality of the entraining/detraining process in the parameterization of cumulus convection, the dataset presented in this study is potentially useful for validation and improvement of the parameterization of shallow convection.
  • 加载中
  • Figure 1.  Location map of four shallow convection cases: RICO (red), BOMEX (orange), ATEX (green) and RICO (purple).

    Figure 2.  Large-scale forcing of (a–d) temperature, (e–h) moisture, and (i–l) vertical velocity (w). From top to bottom are for RICO, BOMEX, ATEX and ARM-SGP. In ARM-SGP, the grey, deep grey, and black lines represent large-scale forcings at the 1130, 1430 and 1730 LT, respectively.

    Figure 3.  Initial (grey) and the last 4-h mean profiles of domain-averaged potential temperature ($ \theta $, dashed line) and total water ($ {q}_{t} $, solid line) simulated by SAM (red), WRF (blue), and UCLA-LES (green), for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) the temporal mean profiles between 1130 and 1530 LT for ARM-SGP.

    Figure 4.  Initial (pink) and the last 4-h mean profiles of domain-averaged zonal ($ U $, solid line) and meridional wind ($ V $, dashed line) simulated by SAM (red), WRF (blue), and UCLA-LES (green), for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) the temporal mean profiles between 1130 and 1530 LT for ARM-SGP. The black line denotes the MMM, and the gray shading denotes the width of twice the standard deviation of the results of all LES models.

    Figure 5.  Convective core fraction for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) ARM-SGP, simulated by SAM (red), WRF (blue), and UCLA-LES (green). The black lines indicate the MMM, and the gray shading denotes the width of twice the standard deviation of the results of all LES models.

    Figure 6.  As in Fig. 5 but for convective eddy transport ($ \overline{{w}^{{'}}{q}_{t}^{{'}}} $). Solid (dashed) lines stand for full (organized) turbulent flux. See text for details.

    Figure 7.  Time–height structures of the domain-averaged cloud condensate (left) and variance of vertical velocity (right) for ARM-SGP simulated by (a, b) SAM, (c, d) WRF, and (e, f) UCLA-LES. White lines represent the top of the boundary layer.

    Figure 8.  Convective core ensemble averages of vertical velocity $ {w}_{c} $ (dashed) and total water $ {q}_{t} $ (solid) for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) ARM-SGP, simulated by SAM (red), WRF (blue), and UCLA-LES (green). The black lines indicate the MMM, and the gray shading denotes the width of twice the standard deviation of the results of all LES models.

    Figure 9.  Diagnosed fractional entrainment rate $ \epsilon $ within the convective core plume for the SAM (red), WRF (blue), and UCLA-LES (green) simulations for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) ARM-SGP. The black lines indicate the MMM, and the gray shading denotes the width of twice the standard deviation of the results of all LES models.

    Figure 10.  Convective updraft (solid lines) and cloud (dashed lines) ensemble averages of vertical velocity $ {w}_{c} $ for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) ARM-SGP, simulated by SAM (red), WRF (blue), and UCLA-LES (green).

    Figure 11.  As in Fig. 10 but for convective updraft (solid lines) and cloud (dashed lines) ensemble averages of total water $ {q}_{{t}_{c}} $

    Figure 12.  Profiles of the convective core fraction for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) ARM-SGP, simulated by SAM (red), WRF (blue), and UCLA-LES (green). Solid lines stand for simulations with the 1.5-order TKE turbulence scheme, while dashed lines stand for simulations with the Smagorinsky-type turbulence scheme in SAM and WRF.

    Figure 13.  As in Fig. 12 but for fractional entrainment rate.

    Dataset Profile
    Dataset Title Shallow Convection Dataset Simulated by Three Different Large Eddy Models
    Time Range RICO: 24 h from 0000 UTC 16 December 2004
    BOMEX: 6 h from 0000 UTC 22 June 1969
    ATEX: 8 h from 0000 UTC 7 February 1969
    ARM-SGP: 14 h from 1130 UTC 21 June 1997
    Geographical scope RICO: Northeast quadrant of a circular area ~300 km in diameter, centered around the Island of Barbuda
    BOMEX: A 500 × 500 km2 ship array in the tropical Atlantic, east of Barbados
    ATEX: Atlantic northeast trade wind region (near 12°N, 35°W)
    ARM-SGP: 140 000 km2 in Oklahoma and Kansas
    Data format netCDF4
    Data volume 1–2 MB
    Data service system https://www.zenodo.org/
    Sources of funding National Key R&D Program of China (Grant No. 2021YFC3000802) and the National Natural Science Foundation of China (Grant No. 42175165)
    Dataset composition Four netCDF files for each shallow convection case. Each file includes variables simulated by three large eddy models (SAM, WRF, and UCLA-LES): 19 2D variables, including winds, temperature, and specific humidity; and conditionally sampled variables, such as in-cloud temperature, total water and vertical velocity.
    DownLoad: CSV

    Table 1.  Surface variables and Coriolis parameter for the four convection cases.

    RICO BOMEX ATEX ARM-SGP
    1130 (LT) 1430 (LT) 1730 (LT)
    Ps (hPa) 1015.4 1015 1015 970.000
    SST (K) 299.80 300.375 298
    LHF (W m−2) 160.00 153.4 140.00 450.00 500.00 420.00
    SHF (W m−2) 7.00 9.46 8.50 140.00 140.00 100.00
    Coriolis (s−1) 0.85×10−4 0.85×10−4 0.85×10−4 0.376×10−4
    DownLoad: CSV

    Table 2.  List of the LES output variables.

    Variable Long name Unit
    U Zonal wind m s−1
    V Meridional wind m s−1
    W Vertical velocity m s−1
    W2 Variance of vertical velocity m2 s−2
    QT Total water g kg−1
    QV Water vapor g kg−1
    QC Cloud condensate g kg−1
    P Pressure hPa
    Z Height m
    A_cor Core fraction 100%
    A_upd Updraft fraction 100%
    A_cld Cloud fraction 100%
    QT_cor Mean qt in core g kg−1
    QT_upd Mean qt in updraft g kg−1
    QT_cld Mean qt in cloud g kg−1
    W_cor Mean w in core m s−1
    W_upd Mean w in updraft m s−1
    W_cld Mean w in cloud m s−1
    QTFLUX Total water flux W m−2
    Corresponding authors: Xiaocong WANG; Yimin LIU
    DownLoad: CSV
  • Bechtold, P., S. K. Krueger, W. S. Lewellen, E. van Meijgaard, C.-H. Moeng, D. A. Randall, A. van Ulden, and S. Wang, 1996: Modeling a stratocumulus-topped PBL: Intercomparison among different one-dimensional codes and with large eddy simulation. Bull. Amer. Meteor. Soc., 77, 2033−2042, https://doi.org/10.1175/1520-0477-77.9.2033.
    Bony, S., and J.-L. Dufresne, 2005: Marine boundary layer clouds at the heart of tropical cloud feedback uncertainties in climate models. Geophys. Res. Lett., 32, L20806. https://doi.org/10.1029/2005GL023851.
    Brown, A. R., and Coauthors, 2002: Large-eddy simulation of the diurnal cycle of shallow cumulus convection over land. Quart. J. Roy. Meteor. Soc., 128, 1075−1093, https://doi.org/10.1256/003590002320373210.
    De Rooy, W. C., and A. P. Siebesma, 2008: A simple parameterization for detrainment in shallow cumulus. Mon. Wea. Rev., 136, 560−576, https://doi.org/10.1175/2007MWR2201.1.
    Gu, J.-F., R. S. Plant, C. E. Holloway, and M. R. Muetzelfeldt, 2020: Pressure drag for shallow cumulus clouds: From thermals to the cloud ensemble. Geophy. Res. Lett., 47, e2020GL090460. https://doi.org/10.1029/2020GL090460.
    Ling, J., C. Y. Li, W. Zhou, X. L. Jia, and C. D. Zhang, 2013: Effect of boundary layer latent heating on MJO simulations. Adv. Atmos. Sci., 30, 101−115, https://doi.org/10.1007/s00376-012-2031-x.
    Khairoutdinov, M. F., and D. A. Randall, 2003: Cloud resolving modeling of the ARM Summer 1997 IOP: Model formulation, results, uncertainties, and sensitivities. J. Atmos. Sci., 60, 607−625, https://doi.org/10.1175/1520-0469(2003)060<0607:CRMOTA>2.0.CO;2.
    Li, Y., and M. Zhang, 2016: Cumulus over the Tibetan Plateau in the Summer Based on CloudSat-CALIPSO Data. J. Climate., 29, 1219−1230, https://doi.org/10.1175/JCLI-D-15-0492.1.
    Liu, P., J. N. Sun, and L. D. Shen, 2016: Parameterization of sheared entrainment in a well-developed CBL. Part I: Evaluation of the scheme through large-eddy simulations. Adv. Atmos. Sci., 33, 1171−1184, https://doi.org/10.1007/s00376-016-5208-x.
    Lu, C. S., Y. G. Liu, and S. J. Niu, 2011: Examination of turbulent entrainment-mixing mechanisms using a combined approach. J. Geophys. Res.: Atmos., 116, D20207. https://doi.org/10.1029/2011JD015944.
    Lu, C. S., Y. G. Liu, S. S. Yum, S. J. Niu, and S. Endo, 2012: A new approach for estimating entrainment rate in cumulus clouds. Geophys. Res. Lett., 39, L04802. https://doi.org/10.1029/2011GL050546.
    Lu, C. S., Y. G. Liu, S. J. Niu, S. Krueger, and T. Wagner, 2013: Exploring parameterization for turbulent entrainment-mixing processes in clouds. J. Geophys. Res.: Atmos., 118, 185−194, https://doi.org/10.1029/2012JD018464.
    Miao, H., X. C. Wang, Y. M. Liu, and G. X. Wu, 2021: A regime-based investigation into the errors of CMIP6 simulated cloud radiative effects using satellite observations. Geophys. Res. Lett., 48, e2021GL095399. https://doi.org/10.1029/2021GL095399.
    Miao, S. G., and W. M. Jiang, 2004: Large eddy simulation and study of the urban boundary layer. Adv. Atmos. Sci., 21, 650−661, https://doi.org/10.1007/BF02915732.
    Neggers, R. A. J., 2009: A dual mass flux framework for boundary layer convection. Part II: Clouds. J. Atmos. Sci., 66, 1489−1506, https://doi.org/10.1175/2008JAS2636.1.
    Neggers, R. A. J., J. D. Neelin, and B. Stevens, 2007: Impact mechanisms of shallow cumulus convection on tropical climate dynamics. J. Climate., 20, 2623−2642, https://doi.org/10.1175/JCLI4079.1.
    Schiro, K. A., and J. D. Neelin, 2019: Deep convective organization, moisture vertical structure, and convective transition using deep-inflow mixing. J. Atmos. Sci., 76, 965−987, https://doi.org/10.1175/JAS-D-18-0122.1.
    Schiro, K. A., J. D. Neelin, D. K. Adams, and B. R. Lintner, 2016: Deep convection and column water vapor over tropical land versus tropical ocean: A comparison between the amazon and the tropical Western Pacific. J. Atmos. Sci., 73, 4043−4063, https://doi.org/10.1175/JAS-D-16-0119.1.
    Siebesma, A. P., and J. W. M. Cuijpers, 1995: Evaluation of parametric assumptions for shallow cumulus convection. J. Atmos. Sci., 52, 650−666, https://doi.org/10.1175/1520-0469(1995)052<0650:EOPAFS>2.0.CO;2.
    Siebesma, A. P., and Coauthors, 2003: A large eddy simulation intercomparison study of shallow cumulus convection. J. Atmos. Sci., 60, 1201−1219, https://doi.org/10.1175/1520-0469(2003)60<1201:ALESIS>2.0.CO;2.
    Stevens, B., 2005: Atmospheric moist convection. Annual Review of Earth and Planetary Sciences, 33, 605−643, https://doi.org/10.1146/annurev.earth.33.092203.122658.
    Stevens, B., 2007: On the growth of layers of nonprecipitating cumulus convection. J. Atmos. Sci., 64, 2916−2931, https://doi.org/10.1175/JAS3983.1.
    Stevens, B., and Coauthors, 2001: Simulations of trade wind cumuli under a strong inversion. J. Atmos. Sci., 58, 1870−1891, https://doi.org/10.1175/1520-0469(2001)058<1870:SOTWCU>2.0.CO;2.
    Sun, J. N., 2009: On the parameterization of convective entrainment: Inherent relationships among entrainment parameters in bulk models. Adv. Atmos. Sci., 26, 1005−1014, https://doi.org/10.1007/s00376-009-7222-8.
    Tiedtke, M., 1989: A comprehensive mass flux scheme for cumulus parameterization in large-scale models. Mon. Wea. Rev., 117, 1779−1800, https://doi.org/10.1175/1520-0493(1989)117<1779:ACMFSF>2.0.CO;2.
    Vanzanten, M. C., and Coauthors, 2011: Controls on precipitation and cloudiness in simulations of trade-wind cumulus as observed during RICO. Journal of Advances in Modeling Earth Systems, 3, M06001. https://doi.org/10.1029/2011MS000056.
    Wang, S. P., and B. Stevens, 2000: Top-hat representation of turbulence statistics in cloud-topped boundary layers: A large eddy simulation study. J. Atmos. Sci., 57, 423−441, https://doi.org/10.1175/1520-0469(2000)057<0423:THROTS>2.0.CO;2.
    Wang, X. C., and M. H. Zhang, 2014: Vertical velocity in shallow convection for different plume types. Journal of Advances in Modeling Earth Systems, 12, 478−489, https://doi.org/10.1002/2014MS000318.
    Wang, X. C., and M. H. Zhang, 2015: The coupling of mixed Rossby-gravity waves with diabatic heating during the TRMM-KWAJEX field campaign. Geophys. Res. Lett., 42, 8241−8249, https://doi.org/10.1002/2015GL065813.
    Wang, Y. H., X. P. Cheng, J. F. Fei, and B. W. Zhou, 2022: Modeling the shallow cumulus-topped boundary layer at gray zone resolutions. J. Atmos. Sci., 79, 2435−2451, https://doi.org/10.1175/JAS-D-21-0339.1.
    Yamaguchi, T., and G. Feingold, 2012: Technical note: Large-eddy simulation of cloudy boundary layer with the Advanced Research WRF model. Journal of Advances in Modeling Earth Systems, 4, M09003. https://doi.org/10.1029/2012MS000164.
    Zhu, L., C. S. Lu, X. Q. Xu, X. He, J. J. Li, S. Luo, Y. Wang, and F. Wang, 2023: The probability density function related to shallow cumulus entrainment rate and its influencing factors in a large-eddy simulation. Adv. Atmos. Sci., https://doi.org/10.1007/s00376-023-2357-6.
  • [1] Lei ZHU, Chunsong LU, Xiaoqi XU, Xin HE, Junjun LI, Shi LUO, Yuan WANG, Fan WANG, 2024: The Probability Density Function Related to Shallow Cumulus Entrainment Rate and Its Influencing Factors in a Large-Eddy Simulation, ADVANCES IN ATMOSPHERIC SCIENCES, 41, 173-187.  doi: 10.1007/s00376-023-2357-6
    [2] Guojing LI, Dongxiao WANG, Changming DONG, Jiayi PAN, Yeqiang SHU, Zhenqiu ZHANG, 2024: Frontogenesis and Frontolysis of a Cold Filament Driven by the Cross-Filament Wind and Wave Fields Simulated by a Large Eddy Simulation, ADVANCES IN ATMOSPHERIC SCIENCES.  doi: 10.1007/s00376-023-3037-2
    [3] MIAO Shiguang, JIANG Weimei, 2004: Large Eddy Simulation and Study of the Urban Boundary Layer, ADVANCES IN ATMOSPHERIC SCIENCES, 21, 650-661.  doi: 10.1007/BF02915732
    [4] Hongrong SHI, Hongbin CHEN, Xiang'ao XIA, Xuehua FAN, Jinqiang ZHANG, Jun LI, Chao LING, 2017: Intensive Radiosonde Measurements of Summertime Convection over the Inner Mongolia Grassland in 2014: Difference between Shallow Cumulus and Other Conditions, ADVANCES IN ATMOSPHERIC SCIENCES, 34, 783-790.  doi: 10.1007/s00376-017-6284-2
    [5] ZHAO Qiaohua, SUN Jihua, ZHU Guangwei, 2012: Simulation and Exploration of the Mechanisms Underlying the Spatiotemporal Distribution of Surface Mixed Layer Depth in a Large Shallow Lake, ADVANCES IN ATMOSPHERIC SCIENCES, 29, 1360-1373.  doi: 10.1007/s00376-012-1262-1
    [6] ZHOU Ningfang, YU Yongqiang, QIAN Yongfu, 2009: Bimodality of the South Asia High Simulated by Coupled Models, ADVANCES IN ATMOSPHERIC SCIENCES, 26, 1226-1234.  doi: 10.1007/s00376-009-7219-3
    [7] YU Yongqiang, HE Jie, ZHENG Weipeng, and LUAN Yihua, 2013: Annual Cycle and Interannual Variability in the Tropical Pacific as Simulated by Three Versions of FGOALS, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 621-637.  doi: 10.1007/s00376-013-2184-2
    [8] Dongxiao WANG, Guojing LI, Lian SHEN, Yeqiang SHU, 2022: Influence of Coriolis Parameter Variation on Langmuir Turbulence in the Ocean Upper Mixed Layer with Large Eddy Simulation, ADVANCES IN ATMOSPHERIC SCIENCES, 39, 1487-1500.  doi: 10.1007/s00376-021-1390-6
    [9] Fang ZHOU, Hong-Li REN, 2017: Dynamical Feedback between Synoptic Eddy and Low-Frequency Flow as Simulated by BCC_CSM1.1(m), ADVANCES IN ATMOSPHERIC SCIENCES, 34, 1316-1332.  doi: 10.1007/s00376-017-6318-9
    [10] Lang ZHANG, Yaoming MA, Weiqiang MA, Binbin WANG, 2018: Comparison of Different Generation Mechanisms of Free Convection between Two Stations on the Tibetan Plateau, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 1137-1144.  doi: 10.1007/s00376-018-7195-6
    [11] Eric P. CHASSIGNET, Xiaobiao XU, 2021: On the Importance of High-Resolution in Large-Scale Ocean Models, ADVANCES IN ATMOSPHERIC SCIENCES, 38, 1621-1634.  doi: 10.1007/s00376-021-0385-7
    [12] SUN Guodong, MU Mu, 2013: Using the Lund-Potsdam-Jena Model to Understand the Different Responses of Three Woody Plants to Land Use in China, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 515-524.  doi: 10.1007/s00376-012-2011-1
    [13] LI Yan, ZHU Jiang, WANG Hui, 2013: The Impact of Different Vertical Diffusion Schemes in a Three-Dimensional Oil Spill Model in the Bohai Sea, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 1569-1586.  doi: 10.1007/s00376-012-2201-x
    [14] Peng LIU, Jianning SUN, Lidu SHEN, 2016: Parameterization of Sheared Entrainment in a Well-Developed CBL. Part I: Evaluation of the Scheme through Large-Eddy Simulations, ADVANCES IN ATMOSPHERIC SCIENCES, 33, 1171-1184.  doi: 10.1007/s00376-016-5208-x
    [15] Dezhen YIN, Fang LI, Yaqiong LU, Xiaodong ZENG, Zhongda LIN, Yanqing ZHOU, 2023: Assessment of Crop Yield in China Simulated by Thirteen Global Gridded Crop Models, ADVANCES IN ATMOSPHERIC SCIENCES.  doi: 10.1007/s00376-023-2234-3
    [16] JIANG Dabang, WANG Huijun, LANG Xianmei, 2005: Evaluation of East Asian Climatology as Simulated by Seven Coupled Models, ADVANCES IN ATMOSPHERIC SCIENCES, 22, 479-495.  doi: 10.1007/BF02918482
    [17] ZHANG Jie, ZHOU Tianjun, BAO Qing, WU Bo, 2010: The Vertical Structures of Atmospheric Temperature Anomalies Associated with El Nino Simulated by the LASG/IAP AGCM: Sensitivity to Convection Schemes, ADVANCES IN ATMOSPHERIC SCIENCES, 27, 1051-1063.  doi: 10.1007/s00376-010-9167-3
    [18] WU Chunqiang, ZHOU Tianjun, SUN De-Zheng, BAO Qing, 2011: Water Vapor and Cloud Radiative Forcings over the Pacific Ocean Simulated by the LASG/IAP AGCM: Sensitivity to Convection Schemes, ADVANCES IN ATMOSPHERIC SCIENCES, 28, 80-98.  doi: 10.1007/s00376-010-9205-1
    [19] Luo Dehai, 1999: Nonlinear Three-Wave Interaction among Barotropic Rossby Waves in a Large-scale Forced Barotropic Flow, ADVANCES IN ATMOSPHERIC SCIENCES, 16, 451-466.  doi: 10.1007/s00376-999-0023-2
    [20] Yazhou ZHANG, Zhijie LIAO, Yaocun ZHANG, Feng NIE, 2016: Characteristics of the Asian-Pacific Oscillation in Boreal Summer Simulated by BCC_CSM with Different Horizontal Resolutions, ADVANCES IN ATMOSPHERIC SCIENCES, 33, 1401-1412.  doi: 10.1007/s00376-016-5266-0
  • suppl_data.zip

Get Citation+

Export:  

Share Article

Manuscript History

Manuscript accepted: 25 September 2023
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Shallow Convection Dataset Simulated by Three Different Large Eddy Models

    Corresponding author: Xiaocong WANG, wangxc@lasg.iap.ac.cn
    Corresponding author: Yimin LIU, lym@lasg.iap.ac.cn
  • 1. College of Earth and Planetary Sciences, University of Chinese Academy of Sciences (UCAS), Beijing 100029, China
  • 2. State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics (LASG), Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences (CAS), Beijing 100029, China

Abstract: Shallow convection plays an important role in transporting heat and moisture from the near-surface to higher altitudes, yet its parameterization in numerical models remains a great challenge, partly due to the lack of high-resolution observations. This study describes a large eddy simulation (LES) dataset for four shallow convection cases that differ primarily in inversion strength, which can be used as a surrogate for real data. To reduce the uncertainty in LES modeling, three different large eddy models were used, including SAM (System for Atmospheric Modeling), WRF (Weather Research and Forecasting model), and UCLA-LES. Results show that the different models generally exhibit similar behavior for each shallow convection case, despite some differences in the details of the convective structure. In addition to grid-averaged fields, conditionally sampled variables, such as in-cloud moisture and vertical velocity, are also provided, which are indispensable for calculation of the entrainment/detrainment rate. Considering the essentiality of the entraining/detraining process in the parameterization of cumulus convection, the dataset presented in this study is potentially useful for validation and improvement of the parameterization of shallow convection.

  • Dataset Profile
    Dataset Title Shallow Convection Dataset Simulated by Three Different Large Eddy Models
    Time Range RICO: 24 h from 0000 UTC 16 December 2004
    BOMEX: 6 h from 0000 UTC 22 June 1969
    ATEX: 8 h from 0000 UTC 7 February 1969
    ARM-SGP: 14 h from 1130 UTC 21 June 1997
    Geographical scope RICO: Northeast quadrant of a circular area ~300 km in diameter, centered around the Island of Barbuda
    BOMEX: A 500 × 500 km2 ship array in the tropical Atlantic, east of Barbados
    ATEX: Atlantic northeast trade wind region (near 12°N, 35°W)
    ARM-SGP: 140 000 km2 in Oklahoma and Kansas
    Data format netCDF4
    Data volume 1–2 MB
    Data service system https://www.zenodo.org/
    Sources of funding National Key R&D Program of China (Grant No. 2021YFC3000802) and the National Natural Science Foundation of China (Grant No. 42175165)
    Dataset composition Four netCDF files for each shallow convection case. Each file includes variables simulated by three large eddy models (SAM, WRF, and UCLA-LES): 19 2D variables, including winds, temperature, and specific humidity; and conditionally sampled variables, such as in-cloud temperature, total water and vertical velocity.
    • Shallow convection frequently occurs in the trade wind regions. It redistributes energy and water by transporting heat and moisture from the near-surface to higher altitudes, thus affecting the vertical structure of thermodynamic fields. By moistening the middle and lower troposphere, shallow convection facilitates the triggering of deep convection (Neggers et al., 2007; Stevens, 2007; Schiro et al., 2016; Schiro and Neelin, 2019), which is important in the propagation of convectively coupled equatorial waves and disturbances (Ling et al., 2013; Wang and Zhang, 2014). Moreover, shallow convective clouds in the subtropics cool the Earth by increasing the reflection of the incoming solar radiation, thus greatly influencing climate change as well (Bony and Dufresne, 2005; Miao et al., 2021). Therefore, shallow convection must be adequately represented in large-scale models, yet its parameterization remains a great challenge, due to the lack of high-resolution observations (Tiedtke, 1989; Bechtold et al., 1996; Li and Zhang, 2016).

      In the last two decades, large eddy models have been widely used for verifying and improving the parameterization of shallow convection (Siebesma et al., 2003; Miao and Jiang, 2004; Wang and Zhang, 2015; Liu et al., 2016; Gu et al., 2020; Wang et al., 2022; Zhu et al., 2023). Large eddy simulation (LES), which typically runs at a horizonal resolution of 100 m, explicitly resolves cumulus scales and acts as a surrogate for real data, providing subgrid-scale (SGS) information that is unavailable from traditional observations. By utilizing LES data, Siebesma and Cuijpers (1995) evaluated assumptions used in the parameterization of shallow cumulus convection and diagnosed bulk entrainment/detrainment rates. This was later followed by De Rooy and Siebesma (2008), who improved the parameterization of entrainment/detrainment rates based on LES results. Besides, Neggers (2009) put forward a dual-mass flux framework for boundary layer convection using LES results.

      Since large eddy models are not fine enough to resolve eddies across all scales, they require parameterizations of SGS processes, such as turbulence, which inevitably introduces uncertainties in the LES modeling (Siebesma et al., 2003). In addition, shallow convection occurring over land and ocean should be distinguished, and more shallow convection cases are needed to make the parameterization more universal. In this study, four shallow convection cases with varying inversion strength were simulated by three different large eddy models.

      The remainder of the paper is structured as follows: The models and experimental design are introduced in section 2. Section 3 compares the performances of the models with respect to selected variables, including wind, the thermodynamic field, cloud, and SGS variables. Also explored in this section is the model sensitivity to the SGS turbulence scheme. Section 4 provides data records and usage notes. Conclusions are given in section 5.

    2.   Description of cases and LES models
    • Four shallow convection cases are used in this study—namely, Rain in Cumulus over the Ocean (RICO) (Vanzanten et al., 2011), the Barbados Oceanographic and Meteorological Experiment (BOMEX) (Siebesma et al., 2003), the Atlantic Trade Wind Experiment (ATEX) (Stevens et al., 2001), and the Atmospheric Radiation Measurement Program–Southern Great Plains site (ARM-SGP) (Brown et al., 2002). The first three cases are for marine equilibrium trade cumuli that differ in inversion strength, and the last case is for transient continental cumuli, thus spanning a wide range of shallow cumulus regimes. All cases come from the Global Energy and Water Cycle Experiment (GEWEX) Cloud System Studies (GCSS) Working Group on boundary layer clouds, with the locations marked in Fig. 1.

      Figure 1.  Location map of four shallow convection cases: RICO (red), BOMEX (orange), ATEX (green) and RICO (purple).

      Profiles of the large-scale forcing used to drive the LES models are shown in Fig. 2, and the surface forcings (i.e., sensible and latent heat flux) are listed in Table 1. For the three maritime cases, subsidence in the lower troposphere dries the atmosphere, which partially offsets the convective moistening. Horizontal advective cooling and drying also serve to counteract the warming and moistening caused by shallow convection, ensuring models quickly reach a steady state. For ARM-SGP, the large-scale forcing varies with time, roughly depicting the life cycle of shallow convection. Other relevant variables, such as surface pressure and temperature, are given in Table 1.

      Figure 2.  Large-scale forcing of (a–d) temperature, (e–h) moisture, and (i–l) vertical velocity (w). From top to bottom are for RICO, BOMEX, ATEX and ARM-SGP. In ARM-SGP, the grey, deep grey, and black lines represent large-scale forcings at the 1130, 1430 and 1730 LT, respectively.

      RICO BOMEX ATEX ARM-SGP
      1130 (LT) 1430 (LT) 1730 (LT)
      Ps (hPa) 1015.4 1015 1015 970.000
      SST (K) 299.80 300.375 298
      LHF (W m−2) 160.00 153.4 140.00 450.00 500.00 420.00
      SHF (W m−2) 7.00 9.46 8.50 140.00 140.00 100.00
      Coriolis (s−1) 0.85×10−4 0.85×10−4 0.85×10−4 0.376×10−4

      Table 1.  Surface variables and Coriolis parameter for the four convection cases.

    • The three large eddy models are the System for Atmospheric Modeling (SAM), developed by Khairoutdinov and Randall (2003), UCLA-LES, described in Stevens (2005), and the large-eddy model of the Weather Research and Forecasting model (WRF) (Yamaguchi and Feingold, 2012). All these models solve prognostic equations for the three components of velocity (u, v, w) and the variables specifying the thermodynamic state. The dynamical frameworks of SAM and UCLA-LES are based on the anelastic equations of motion in height coordinates, while WRF uses fully compressible equations in terrain-following hydrostatic-pressure vertical coordinates. Besides, the models differ in the way they parametrize the SGS fluxes. While SAM and WRF use 1.5-order SGS closure based on prognostic SGS turbulent kinetic energy (TKE), UCLA-LES uses Smagorinsky-type closure instead. For non-precipitating shallow convection cases in this study, the precipitation process is turned off to make sure total water is exactly conserved during model integration. Therefore, differences due to cloud microphysics are excluded.

      All LES models are configured as follows: For BOMEX and ATEX, the models use a horizontal resolution of 64 × 64, with a grid spacing of 100 m and a vertical resolution of 75 levels from the ground to 3000 m. For ARM-SGP and RICO, the horizontal domain is expanded to 128 × 128 and the vertical layers increased to 100 levels, with the horizonal and vertical resolution remaining unchanged. While SAM and UCLA-LES use a constant resolution of 40 m in altitudinal coordinates, WRF uses terrain-following hydrostatic-pressure vertical coordinates that are unequally spaced in the vertical direction. All models are integrated for 6 h for BOMEX, 8 h for ATEX, 14 h for ARM-SGP, and 24 h for RICO. Details of these experimental designs are documented in Siebesma et al. (2003) for BOMEX, Vanzanten et al. (2011) for RICO, Stevens et al. (2001) for ATEX, and Brown et al. (2002) for ARM-SGP.

    3.   Model performance
    • Figure 3 shows the initial (grey lines) and the last 4-h mean profiles of domain-averaged potential temperature $ \left(\theta \right) $ and total water ($ {q}_{t}) $. For each case, the model simulations generally bear close resemblance to each other, although differences exist in the simulated$ {q}_{t} $within the boundary layer. Compared to SAM and UCLA-LES, WRF appears to be more efficient in transporting moisture to upper levels, making the free troposphere wetter and the boundary layer drier. Such a phenomenon is, however, not observed in ARM-SGP.

      Figure 3.  Initial (grey) and the last 4-h mean profiles of domain-averaged potential temperature ($ \theta $, dashed line) and total water ($ {q}_{t} $, solid line) simulated by SAM (red), WRF (blue), and UCLA-LES (green), for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) the temporal mean profiles between 1130 and 1530 LT for ARM-SGP.

      Figure 4 shows the mean profiles of horizontal winds for the last 4 h, along with the initial one. Also superimposed is the multi-model mean (MMM; black line) and the width of twice the standard deviation of the results of all LES models (gray shading). The geostrophic wind is set to the same as the initial, which decreases with increasing height. Due to the Coriolis force, all models produce an increase in the meridional wind speed and a decrease in the zonal wind speed. The differences among the models are relatively larger than those seen in temperature and humidity.

      Figure 4.  Initial (pink) and the last 4-h mean profiles of domain-averaged zonal ($ U $, solid line) and meridional wind ($ V $, dashed line) simulated by SAM (red), WRF (blue), and UCLA-LES (green), for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) the temporal mean profiles between 1130 and 1530 LT for ARM-SGP. The black line denotes the MMM, and the gray shading denotes the width of twice the standard deviation of the results of all LES models.

    • This subsection evaluates the simulation of clouds and SGS variables. Figure 5 shows the mean profiles of the convective core fraction for the four cases, in which red, blue, green, and black lines denote results from SAM, WRF, UCLA-LES, and the MMM, respectively. The convective core here refers to the rising plumes that experience positive buoyancy and contain condensed liquid water, which is important for convective eddy transport (Siebesma and Cuijpers, 1995). The core fraction has a maximum near the cloud base, and decreases with increasing height, except for ATEX where a secondary peak forms near the cloud top due to strong inversion. All three models show similar profiles in each case. Overall, the WRF model tends to produce a higher cloud base and smaller core fraction compared with the other two models. The difference in cloud base height between WRF and the other two models is partly due to the different vertical coordinates used by the different models.

      Figure 5.  Convective core fraction for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) ARM-SGP, simulated by SAM (red), WRF (blue), and UCLA-LES (green). The black lines indicate the MMM, and the gray shading denotes the width of twice the standard deviation of the results of all LES models.

      As mentioned previously, the main role of shallow convection is to redistribute heat and moisture in the vertical direction. For any conserved variable $ \phi $, the SGS flux can be decomposed into three terms:

      where an overbar stands for the average, a prime represents deviation from the mean value within the plume, $ \alpha $ is the convective cloud fraction, and the superscripts c and e stand for the active cloudy part and the surrounding environment, respectively. A widely used assumption in cumulus convection parameterization is that in-cloud and environmental turbulence are much smaller than the organized turbulence term [third term on the right-hand side of Eq. (1)], which describes the contribution due to mean organized updrafts and compensating subsidence in the ambient environment, and thus the first two terms on the right-hand side of Eq. (1) can be dropped, also known as the top-hat approximation (Wang and Stevens, 2000). To verify the plausibility of top-hat approximation, the organized turbulent flux (dashed lines) and the full SGS flux (solid lines) of total water are shown in Fig. 6. Indeed, the organized turbulent flux is close to the full SGS flux within the convective layers. The models show similar profiles but different magnitudes for the organized turbulent fluxes. Compared to WRF and UCLA-LES, SAM produces larger fluxes because of greater simulation of the core fraction (Fig. 5).

      Figure 6.  As in Fig. 5 but for convective eddy transport ($ \overline{{w}^{{'}}{q}_{t}^{{'}}} $). Solid (dashed) lines stand for full (organized) turbulent flux. See text for details.

      The ARM-SGP case targets the evolution of the convective boundary layer (CBL). Figure 7 shows the temporal evolution of the mean profiles of cloud condensate (left) and variance of the vertical velocity (right). The CBL top is defined as the height where the buoyancy flux is minimal, marked by the white line in Fig. 7. It can be seen that the cloud base, i.e., the lifting condensation level, is always collocated with the CBL top, suggesting cumulus updrafts are firmly rooted in the CBL. After about 1100 LT, a condensation center appears above the CBL, implying cumulus updrafts reach the level of free convection. Before that, the variance of vertical velocity is much smaller above the CBL. At around 1400 LT, the variance of vertical velocity acquires a maximum above the CBL because clouds penetrate above the level of free convection and accelerate upward. Clouds then lose buoyancy quickly due to entrainment. After 1400 LT, the surface flux starts to decrease and shallow convection decays accordingly. All models produce a similar evolution of shallow convection as described above, but the detailed structure is different. For example, SAM simulated a higher cloud top than WRF and UCLA-LES. This is due to the fact that SAM produces larger variance of vertical velocity within the CBL, which helps clouds to penetrate progressively deeper before being diluted by entrainment.

      Figure 7.  Time–height structures of the domain-averaged cloud condensate (left) and variance of vertical velocity (right) for ARM-SGP simulated by (a, b) SAM, (c, d) WRF, and (e, f) UCLA-LES. White lines represent the top of the boundary layer.

    • Figure 8 shows the mean profiles of vertical velocity (dashed lines) and total water (solid lines) within the convective core plume type. In-cloud vertical velocity increases with height because of positive buoyancy by definition, whereas in-cloud total water decreases with height owing to entrainment. Note the increase of in-cloud total water near the cloud top is due to the small sample of grid points in the calculation. The results of all three model simulations are close to each other, with only considerable standard deviation near the cumulus top where the cloud fraction is small.

      Figure 8.  Convective core ensemble averages of vertical velocity $ {w}_{c} $ (dashed) and total water $ {q}_{t} $ (solid) for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) ARM-SGP, simulated by SAM (red), WRF (blue), and UCLA-LES (green). The black lines indicate the MMM, and the gray shading denotes the width of twice the standard deviation of the results of all LES models.

      Following Wang and Zhang (2014), diagnosis of the entrainment rate $ \varepsilon $ is expressed as

      where the subscripts c and e denote the active cloudy part and the surrounding environment, respectively. The diagnosed entrainment rate for each model is presented in Fig. 9, along with the MMM and the standard deviation. All models show a decrease in $ \epsilon $ with increasing height for all cases except ATEX, where a secondary entrainment peak appears near the cloud top. This is due to the upside-down convection originating from the stratocumulus top. Again, all models show similar profiles, with the difference mostly falling within the range of 0.2 km−1, in line with the in-cloud $ {q}_{t} $ behaviors shown in Fig. 8.

      Figure 9.  Diagnosed fractional entrainment rate $ \epsilon $ within the convective core plume for the SAM (red), WRF (blue), and UCLA-LES (green) simulations for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) ARM-SGP. The black lines indicate the MMM, and the gray shading denotes the width of twice the standard deviation of the results of all LES models.

      Besides convective core, two other types of plumes are provided, i.e., convective updraft, defined as plumes having liquid water and upward vertical velocity, and cloud, defined as grid points containing liquid water. Figures 10 and 11 show the mean profiles of vertical velocity and total water within convective updraft (solid lines) and cloud (dashed lines), respectively. The vertical velocity in the cloud plume either decreases or increases with height at a considerably smaller rate, which is less than half the vertical velocity in the updraft and convective core plumes. The decrease in total water with height in the cloud plume is slightly larger than that in updraft, which in turn is greater than in convective core. This is consistent with the larger entrainment rate within cloud and updraft plumes (Wang and Zhang, 2014). Comparison of Figs. 10 and 11 shows that, for different types of plumes, the difference in simulated vertical velocity is greater than the difference in total water across the models.

      Figure 10.  Convective updraft (solid lines) and cloud (dashed lines) ensemble averages of vertical velocity $ {w}_{c} $ for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) ARM-SGP, simulated by SAM (red), WRF (blue), and UCLA-LES (green).

      Figure 11.  As in Fig. 10 but for convective updraft (solid lines) and cloud (dashed lines) ensemble averages of total water $ {q}_{{t}_{c}} $

      These results suggest the need to distinguish between different types of plumes in the parameterization of shallow convection, since they have very different properties. For example, convective core may be more relevant to convective-scale transport, convective updraft may be more relevant to condensation, and the cloud plume is more relevant to radiation. These conditionally sampled variables thus shed light on the pathways toward further development of convection schemes by considering multiple plume types.

    • To see to what extent the differences in model behavior are due to the choice of SGS turbulence scheme, we carry out sensitivity experiments by using the Smagorinsky-type turbulence scheme with the Smagorinsky coefficient set to 0.15 in SAM and WRF. Figures 12 and 13 show the mean profiles of the convective core fraction and the diagnosed entrainment rate for the four cases, with the Smagorinsky-type simulations in SAM and WRF superimposed. As can be seen, the change of SGS turbulence scheme does not yield significant differences in all cases except ATEX, where the Smagorinsky-type closure tends to produce less cloud compared to the 1.5-order TKE closure (Fig. 12c). This is mainly due to the weaker mixing near the stratocumulus top. If we increase the Smagorinsky coefficient from 0.15 to 0.17, the differences between simulations using Smagorinsky-type closure and 1.5-order TKE closure are almost indistinguishable (not shown). This suggests that the difference in SGS turbulence scheme is not directly responsible for the differences in model behaviors. It is more likely that other factors, such as technological details in the framework, advection, diffusion, damping, and, more importantly, their interactions, are responsible for the differences. This also highlights the need to use different LES models to reduce model dependencies.

      Figure 12.  Profiles of the convective core fraction for the cases of (a) RICO, (b) BOMEX, (c) ATEX, and (d) ARM-SGP, simulated by SAM (red), WRF (blue), and UCLA-LES (green). Solid lines stand for simulations with the 1.5-order TKE turbulence scheme, while dashed lines stand for simulations with the Smagorinsky-type turbulence scheme in SAM and WRF.

      Figure 13.  As in Fig. 12 but for fractional entrainment rate.

    4.   Data records and usage notes
    • This dataset, simulated by the large eddy models SAM, WRF, and UCLA-LES, for four shallow convection cases, has been submitted to Zenodo and is available at https://doi.org/10.5281/zenodo.7895982. The dataset format is version 4 of the Network Common Data Form (NetCDF), which is easy to read and written by professional common software such as Climate Data Operators (CDO, https://code.mpimet.mpg.de/projects/cdo), NCAR Command Language (NCL, http://www.ncl.ucar.edu), NetCDF Operator (NCO, http://nco.sourceforge.net), and Python (https://www.python. org).

      For each convection case, there are 19 2D variables, including winds, temperature, and specific humidity. Conditionally sampled variables, such as in-cloud temperature, total water, and vertical velocity, are also provided, which are indispensable for calculating entrainment rates. Besides convective core, two other types of plumes are provided, i.e., convective updraft, defined as plumes having liquid water and upward vertical velocity, and cloud, defined as grid points containing liquid water. The dataset presented in this study can be used for validation of the parameterization of shallow convection and may provide guidance for further development of convection schemes by considering multiple plume types.

    5.   Summary
    • This paper describes a dataset simulated by the three large eddy models SAM, WRF, and UCLA-LES, for four shallow convection cases (RICO, BOMEX, ATEX, and ARM-SGP), which vary in inversion strength. The general features of the model output are presented, and the data records and usage notes are provided. Results show that the different LES models show similar performance for each case, albeit with noticeable differences manifested in the horizontal winds. The MMM and the standard deviation of the results of all LES models are also given. If all models produce similar features, this means the standard deviation between the models is small, so their average can be taken as the “ground truth”. If the models show greater diversity in some features, then their mean plus/minus the standard deviation yields an acceptable range for that field. Sensitivity experiments show that the differences of the SGS turbulence scheme are not directly responsible for the differences in model behavior. Other factors, such as technological details in the model framework, are more likely to matter.

      Compared to SAM and UCLA-LES, the WRF model appears to be more efficient in transporting moisture to upper levels, leading to a wetter free atmosphere and drier boundary layer. SAM tends to simulate a higher cloud top than WRF and UCLA-LES, due to the fact that SAM produces a larger variance of vertical velocity within the CBL, which helps clouds to penetrate progressively deeper. Conditionally sampled variables, such as in-cloud moisture and vertical velocity, are also analyzed, which are indispensable for calculating entrainment/detrainment rates. All models show that the entrainment rate decreases with height in all cases except ATEX, where a secondary peak appears near the cloud top. Given that the entraining/detraining process is central to the parameterization of cumulus convection (Sun, 2009; Lu et al., 2011, 2012, 2013), the dataset presented in this study is potentially useful for validating and improving such parameterization schemes.

      Acknowledgements. The authors thank the two anonymous reviewers for their constructive comments, which certainly helped to clarify and improve the original paper. This work was jointly supported by the National Key R&D Program of China (Grant No. 2021YFC3000802) the National Natural Science Foundation of China (Grant No. 42175165), and the National Key Scientific and Technological Infrastructure project “Earth System Numerical Simulation Facility” (EarthLab).

      Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

      Variable Long name Unit
      U Zonal wind m s−1
      V Meridional wind m s−1
      W Vertical velocity m s−1
      W2 Variance of vertical velocity m2 s−2
      QT Total water g kg−1
      QV Water vapor g kg−1
      QC Cloud condensate g kg−1
      P Pressure hPa
      Z Height m
      A_cor Core fraction 100%
      A_upd Updraft fraction 100%
      A_cld Cloud fraction 100%
      QT_cor Mean qt in core g kg−1
      QT_upd Mean qt in updraft g kg−1
      QT_cld Mean qt in cloud g kg−1
      W_cor Mean w in core m s−1
      W_upd Mean w in updraft m s−1
      W_cld Mean w in cloud m s−1
      QTFLUX Total water flux W m−2
      Corresponding authors: Xiaocong WANG; Yimin LIU

      Table 2.  List of the LES output variables.

Reference

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return