Advanced Search
Article Contents

Numerical Simulations of a Florida Sea Breeze and Its Interactions with Associated Convection: Effects of Geophysical Representation and Model Resolution


doi: 10.1007/s00376-021-1216-6

  • The Florida peninsula in the USA has a frequent occurrence of sea breeze (SB) thunderstorms. In this study, the numerical simulation of a Florida SB and its associated convective initiation (CI) is simulated using the mesoscale community Weather Research and Forecasting (WRF) model in one-way nested domains at different horizontal resolutions. Results are compared with observations to examine the accuracy of model-simulated SB convection and factors that influence SB CI within the simulation. It is found that the WRF model can realistically reproduce the observed SB CI. Differences are found in the timing, location, and intensity of the convective cells at different domains with various spatial resolutions. With increasing spatial resolution, the simulation improvements are manifested mainly in the timing of CI and the orientation of the convection after the sea breeze front (SBF) merger into the squall line over the peninsula. Diagnoses indicate that accurate representation of geophysical variables (e.g., coastline and bay shape, small lakes measuring 10–30 km2), better resolved by the high resolution, play a significant role in improving the simulations. The geophysical variables, together with the high resolution, impact the location and timing of SB CI due to changes in low-level atmospheric convergence and surface sensible heating. More importantly, they enable Florida lakes (30 km2 and larger) to produce noticeable lake breezes (LBs) that collide with the SBFs to produce CI. Furthermore, they also help the model reproduce a stronger convective squall line caused by merging SBs, leading to more accurate locations of postfrontal convective systems.
    摘要: 美国佛罗里达半岛海风和雷暴频发。本文采用中尺度预报模式WRF,模拟研究了不同水平分辨率对佛罗里达海风及其对流触发的影响,在检验海风及其对流模拟能力的基础上,分析了影响海风及其对流触发的关键因素。结果表明,WRF模式可较好地再现出观测的海风及其对流演变。采用不同水平分辨率模拟的对流时间、位置和强度存在明显的差异。提高水平分辨率对模拟的改进主要表现在海风锋在佛罗里达半岛合并为飑线后的对流触发时间和对流移动方向上。进一步的诊断表明,提高水平分辨率对海岸线、海湾形状和尺度为10-30 km2的小尺度湖泊等地形单体的刻画更为精确,从而在模拟改进中发挥了重要的作用。更精细的地形单体可通过改变低层大气辐合和地表感热加热对海风及其对流触发的位置和时间产生重要的影响。值得注意的是,佛罗里达半岛面积为30 km2及以上的湖泊不仅会产生明显的湖风并与海风锋碰撞触发对流,还有利于模式再现出海风合并引起的更强的对流飑线,从而提高锋面后部对流系统位置的模拟。
  • 加载中
  • Figure 1.  Observed soundings for four Florida stations: (a) Tallahassee (TLH), (b) Jacksonville (JAX), (c) Tampa Bay (TPA), and (d) Miami (MIA) at 1200 UTC 6 September 2012. Images courtesy of University of Wyoming [http://weather.uwyo.edu/upperair/sounding.html]. Locations of sounding stations can be found in Fig. 3a.

    Figure 2.  Geostationary Operational Environmental Satellite-East (GOES-E) visible satellite images of convective cloud cover over the Florida peninsula at (a) 1401 UTC, (b) 1515 UTC, (c) 2001 UTC, and (d) 2031 UTC 6 September 2012. Red lines show the location of the WCSB and ECSB fronts. Yellow line in panel (d) indicates formation of the squall line due to the collision of the SBFs. Images courtesy of College of DuPage [http://www.cod.edu/] and Carl Jones [https://northflwx.wordpress.com/2012/10/01/florida-coast-to-coast-sea-breeze-convergence/].

    Figure 3.  Terrain heights and water body locations for (a) d03 and (b) d04. Red lines indicate the locations of the vertical cross sections cut west-to-east through MacDill, Air Force Base (MCF, northernmost) and Lake Okeechobee (OBE, southernmost). The colored boxes indicate locations used later in the study for a zoomed-in analysis. The locations of the 23 observing stations used to aid in the verification of the WRF simulation are also plotted. Note the increased number of lakes in (b) compared to (a). Water body names and locations are annotated in panel (c). The left corner inset of panel (c) is a zoomed-in area over Tampa Bay. (d) Locations of WRF model domains for all four-level nested domains.

    Figure 4.  WRF-simulated zonal (u-component) winds (color: red easterly, blue westerly; units: m s–1) and total wind (black arrows; unit reference vector 2.5 m s–1) at surface level for a zoomed-in region of d03 (yellow box in Fig. 3a) at: (a) 1300 UTC, (b) 1400 UTC, (c) 1500 UTC, and (d) 1600 UTC 6 September 2012.

    Figure 5.  Cross section (through northernmost red line in Fig. 3a) of d03 WRF-simulated zonal (u-component) winds (color: red easterly, blue westerly; units: m s–1) and total horizontal wind (black arrows represent horizontal winds at each height level along longitude, reference vector 10 m s–1; upward arrows mean north direction) at: (a) 1600 UTC, (b) 1700 UTC, (c) 1930 UTC, and (d) 2030 UTC 6 September 2012. Red contour lines indicate the zero-wind contour heights, which are indicative of approximate SBF heights.

    Figure 6.  18-h rainfall accumulation totals from 1200 UTC 6 September 2012 to 0600 UTC 7 September 2012 for (a) d03 WRF-simulated precipitation accumulation totals (blue contours; units: mm) and CCPA data (color; units: mm). Panel (b) is the same as (a) but for d04. WRF-simulated precipitation accumulation totals (color; units: mm) for (c) d03 and (d) d04 for the same times as (a) and (b).

    Figure 7.  Stage IV precipitation accumulation (color; units: mm) and d03 (left column) and d04 (right column) WRF-simulated precipitation accumulation (black contours; units: mm) at: (a–b) 1400 UTC, (c–d) 1800 UTC, (e–f) 2100 UTC, and (g–h) 2300 UTC 6 September 2012. The red arrows in (e) and (f) indicate squall lines.

    Figure 8.  WRF-simulated divergence (color: red convergence, blue divergence; units: s–1) and total wind (black arrows; reference vector 2.5 m s–1) at surface level for a zoomed-in region of d03 (yellow box in Fig. 3a) at: (a) 1500 UTC and (c) 1530 UTC 6 September 2012. Panels (b) and (d) are the same as a) and (c) respectively, but for the zoomed-in region of d04 (yellow box in Fig. 3b) and include max reflectivity (red contours; units: dBZ).

    Figure 9.  Temperature difference [color contours; units: oF , oF = (oC × 9/5) + 32] between the d03 and d04 WRF simulations at 0330 UTC 6 September 2012 for a zoomed-in region of d03 and d04 (yellow boxes in Figs. 3a and 3b).

    Figure 10.  Cross section (through southernmost red line in Fig. 3b) of d04 WRF-simulated zonal (u-component) winds (color: red easterly, blue westerly; units: m s–1) and total horizontal wind (black arrows represent horizontal winds at each height level along longitude, reference vector 10 m s–1; upward arrows mean north direction) at 1600 UTC 6 September 2012. Red contour line indicates the zero-wind contour heights, which indicate Lake Okeechobee’s approximate lake breeze heights. Blue arrows indicate the locations of the SBs, while the red arrows indicate the lake breezes for 1) Lake Weohyakapka, 2) Lake Kissimmee, 3) Lake Marian, and 4) Indian River.

    Figure 11.  D03 WRF-simulated divergence (color: red convergence, blue divergence; units: s–1) at surface level, max reflectivity (yellow and purple contours; units: 5 dBZ), and total wind (black arrows; reference vector 2.5 m s–1) at surface level for a zoomed-in region of d03 (purple box in Fig. 3a) at: (a) 1820 UTC, (b) 1825 UTC, (c) 1840 UTC, (d) 1845 UTC, (e) 1900 UTC, and (f) 1910 UTC 6 September 2012. Yellow reflectivity contour lines indicate the convective cell’s first appearance in the simulation and are numbered accordingly in panels (b–f); thereafter, the reflectivity contour lines are purple. Panel (a) indicates the lake effect divergence from the following water bodies: A—St John’s River, B—Crescent Lake, C—Lake George, D—Lakes Griffin, Harris, and Eustis, E—Lake Apopka. See text for more details.

    Figure 12.  WRF-simulated zonal (u-component) winds (color: red easterly, blue westerly; units: m s–1) and total wind (black arrows; reference vector 2.5 m s–1) at surface level for a zoomed-in region of d03 (yellow box in Fig. 3a) at (a) 2100 UTC, (c) 2200 UTC, and (e) 2300 UTC 6 September 2012. Panels (b), (d), and (f) are the same as (a), (c), and (e) respectively, but for the zoomed-in region of d04 (yellow box in Fig. 3b). The black circles denote the locations of squall lines.

    Table 1.  WRF model parameter settings for domains d01, d02, d03, and d04.

    Parameterd01d02d03d04
    Physical domain139 × 97 (lon. × lat.)253 × 202 (lon. × lat.)343 × 313 (lon. × lat.)502 × 661 (lon. × lat.)
    Horizontal resolution27 km9 km3 km1 km
    Vertical resolution62 levels with variable Δz (26 levels below 850 hPa): model top 50 hPa*
    Time integration120-sec time step; 36-h duration*
    Boundary conditionsDamping depth over top 500 m*
    Lateral boundary21 600-sec time interval*
    Surface layer physicsMonin-Obukhov Janjić Eta (MO-JE)*
    MicrophysicsNew Thompson*
    Longwave radiationRapid Radiative Transfer Model (RRTM)*
    Shortwave radiationMM5 (Dudhia)*
    PBL physicsMellor-Yamada-Janjić (MYJ)*
    Land-surface physicsNoah land surface model (LSM)*
    Cumulus physicsNew Kain-FritschNew Kain-FritschNoneNone
    Vertical Diffusion2nd order diffusion*
    Time-integrationRunge-Kutta 3rd order*
    Eddy coefficientHorizontal Smagorinsky 1st order closure*
    *Values are the same for all four grids.
    DownLoad: CSV

    Table 2.  Contingency table illustrating the counts used in verification statistics of dichotomous (e.g., Yes/No) forecasts and observations.

    ForecastObserved
    YesNoSum
    YesHits (H)False alarmsAll forecasted (F)
    NoMissesCorrect rejectionsNo forecasted
    (NF)
    SumAll observed (O)No observed
    (NO)
    Total (N)
    DownLoad: CSV

    Table 3.  Equitable Threat and Bias Score values for four rainfall thresholds over an 18-h forecast period (1200 UTC 6 September 2012 to 0600 UTC 7 September 2012).

    Precipitation
    Threshold (mm)
    Domain
    d03d04
    Equitable Threat Score (ETS)2.540.120.28
    6.350.10.19
    12.70.080.15
    25.400.04
    Bias Score (BS)2.540.120.26
    6.350.110.2
    12.70.120.17
    25.40.140.23
    DownLoad: CSV
  • Abbs, D. J., and W. L. Physick, 1992: Sea-breeze observations and modelling: A review. Aust. Meteor. Mag., 41, 7−19.
    Abulikemu, A., X. Xin, Y. Wang, J. F. Ding, S. S. Zhang, and W. Q. Shen, 2016: A modeling study of convection initiation prior to the merger of a sea-breeze front and a gust front. Atmospheric Research, 182, 10−19, https://doi.org/10.1016/j.atmosres.2016.07.003.
    Baker, R. D., B. H. Lynn, A. Boone, W.-K. Tao, and J. Simpson, 2001: The influence of soil moisture, coastline curvature, and land-breeze circulations on sea-breeze-initiated precipitation. Journal of Hydrometeorology, 2, 193−211, https://doi.org/10.1175/1525-7541(2001)002<0193:TIOSMC>2.0.CO;2.
    Blanchard, D. O., and R. E. López, 1985: Spatial patterns of convection in south Florida. Mon. Wea. Rev., 113, 1282−1299, https://doi.org/10.1175/1520-0493(1985)113<1282:SPOCIS>2.0.CO;2.
    Boybeyi, Z., and S. Raman, 1992: A three-dimensional numerical sensitivity study of convection over the Florida peninsula. Bound.-Layer Meteorol., 60, 325−359, https://doi.org/10.1007/BF00155201.
    Chen, F., and J. Dudhia, 2001: Coupling an advanced land surface-hydrology model with the Penn State-NCAR MM5 modeling system. Part I: Model implementation and sensitivity. Mon. Wea. Rev., 129, 569−585, https://doi.org/10.1175/1520-0493(2001)129<0569:CAALSH>2.0.CO;2.
    Clark, P. A., K. A. Browning, R. M. Forbes, C. J. Morcrette, A. M. Blyth, and H. W. Lean, 2014: The evolution of an MCS over southern England. Part 2: Model simulations and sensitivity to microphysics. Quart. J. Roy. Meteor. Soc., 140, 458−479, https://doi.org/10.1002/qj.2142.
    Crosman, E. T., and J. D. Horel, 2010: Sea and lake breezes: A review of numerical studies. Bound.-Layer Meteorol., 137, 1−29, https://doi.org/10.1007/s10546-010-9517-9.
    Dudhia, J., 1989: Numerical study of convection observed during the Winter Monsoon Experiment using a mesoscale two-dimensional model. J. Atmos. Sci., 46, 3077−3107, https://doi.org/10.1175/1520-0469(1989)046<3077:NSOCOD>2.0.CO;2.
    Etherton, B., and P. Santos, 2008: Sensitivity of WRF forecasts for south Florida to initial conditions. Wea. Forecasting, 23, 725−740, https://doi.org/10.1175/2007WAF2006115.1.
    Fankhauser, J. C., N. A. Crook, J. Tuttle, L. J. Miller, and C. G. Wade, 1995: Initiation of deep convection along boundary layer convergence lines in a semitropical environment. Mon. Wea. Rev., 123, 291−314, https://doi.org/10.1175/1520-0493(1995)123<0291:IODCAB>2.0.CO;2.
    Fovell, R. G., 2005: Convective initiation ahead of the sea-breeze front. Mon. Wea. Rev., 133, 264−278, https://doi.org/10.1175/MWR-2852.1.
    Hahmann, A. N., Y. B. Liu, and T. T. Warner, 2006: Mesoscale circulations over the Athens metropolitan area during the 2004 summer Olympic games. Proc. 6th Symposium on the Urban Environment and Forum on Managing our Physical and Natural Resources, Atlanta, GA, Amer. Meteor. Soc.
    Hou, D. C., and Coauthors, 2014: Climatology-calibrated precipitation analysis at fine scales: Statistical adjustment of stage IV toward CPC gauge-based analysis. Journal of Hydrometeorology, 15, 2542−2557, https://doi.org/10.1175/JHM-D-11-0140.1.
    Iacono, M. J., J. S. Delamere, E. J. Mlawer, M. W. Shephard, S. A. Clough, and W. D. Collins, 2008: Radiative forcing by long-lived greenhouse gases: Calculations with the AER radiative transfer models. J. Geophys. Res., 113, D13103, https://doi.org/10.1175/1520-0493(1994)122<0927:TSMECM>2.0.CO;2.
    Kain, J. S., 2004: The Kain-Fritsch convective parameterization: An update. J. Appl. Meteorol. Climatol., 43, 170−181, https://doi.org/10.1175/1520-0450(2004)043<0170:TKCPAU>2.0.CO;2.
    Knievel, J. C., D. L. Rife, J. A. Grim, A. N. Hahmann, J. P. Hacker, M. Ge, and H. H. Fisher, 2010: A simple technique for creating regional composites of sea surface temperature from MODIS for use in operational mesoscale NWP. J. Appl. Meteorol. Climatol., 49, 2267−2284, https://doi.org/10.1175/2010JAMC2430.1.
    Laird, N. F., D. A. R. Kristovich, R. M. Rauber, H. T. Ochs, and L. J. Miller, 1995: The Cape Canaveral sea and river breezes: Kinematic structure and convective initiation. Mon. Wea. Rev., 123, 2942−2956, https://doi.org/10.1175/1520-0493(1995)123<2942:TCCSAR>2.0.CO;2.
    Lin, Y., and K. E. Mitchell, 2005: The NCEP stage II/IV hourly precipitation analyses: Development and applications. Preprints, 19th Conf. on Hydrology, San Diego, CA, Amer. Meteor. Soc.
    Lock, N. A., and A. L. Houston, 2014: Empirical examination of the factors regulating thunderstorm initiation. Mon. Wea. Rev., 142, 240−258, https://doi.org/10.1175/MWR-D-13-00082.1.
    Mellor, G. L., and T. Yamada, 1974: A hierarchy of turbulence closure models for planetary boundary layers. J. Atmos. Sci., 31, 1791−1806, https://doi.org/10.1175/1520-0469(1974)031<1791:AHOTCM>2.0.CO;2.
    Mesinger, F., and Coauthors, 2006: North American regional reanalysis. Bull. Amer. Meteor. Soc., 87, 343−360, https://doi.org/10.1175/BAMS-87-3-343.
    Miller, S. T. K., B. D. Keim, R. W. Talbot, and H. Mao, 2003: Sea breeze: Structure, forecasting, and impacts. Rev. Geophys., 41, 1011, https://doi.org/10.1029/2003RG000124.
    Nicholls, M. E., R. A. Pielke, and W. R. Cotton, 1991: A two-dimensional numerical investigation of the interaction between sea breezes and deep convection over the Florida peninsula. Mon. Wea. Rev., 119, 298−323, https://doi.org/10.1175/1520-0493(1991)119<0298:ATDNIO>2.0.CO;2.
    Ogawa, S., W. M. Sha, T. Iwasaki, and Z. F. Wang, 2003: A numerical study on the interaction of a sea-breeze front with convective cells in the daytime boundary layer. J. Meteor. Soc. Japan, 81, 635−651, https://doi.org/10.2151/jmsj.81.635.
    Rao, P. A., and H. E. Fuelberg, 2000: An investigation of convection behind the Cape Canaveral sea-breeze front. Mon. Wea. Rev., 128, 3437−3458, https://doi.org/10.1175/1520-0493(2000)128<3437:AIOCBT>2.0.CO;2.
    Schwartz, C. S., and R. A. Sobash, 2019: Revisiting sensitivity to horizontal grid spacing in convection-allowing models over the central and eastern united states. Mon. Wea. Rev., 147, 4411−4435, https://doi.org/10.1175/MWR-D-19-0115.1.
    Segal, M., M. Leuthold, R. W. Arritt, C. Anderson, and J. Shen, 1997: Small lake daytime breezes: Some observational and conceptual evaluations. Bull. Amer. Meteor. Soc., 78, 1135−1148, https://doi.org/10.1175/1520-0477(1997)078<1135:SLDBSO>2.0.CO;2.
    Skamarock, W. C., and Coauthors, 2008: A description of the Advanced Research WRF version 3. NCAR Tech. Note NCAR/TN-475+STR, 113 pp.
    Thompson, G., P. R. Field, R. M. Rasmussen, W. D. Hall, 2008: Explicit Forecasts of Winter Precipitation Using an Improved Bulk Microphysics Scheme. Part II: Implementation of a New Snow Parameterization. Mon. Wea. Rev., 136, 5095−5115, https://doi.org/10.1175/2008MWR2387.1.
    Trier, S. B., 2003: Convective storms: Convective initiation. Encyclopedia of Atmospheric Sciences, J. R. Holton, Ed., Academic Press, 560−569.
    Xue, M. and W. Martin, 2006a: A high-resolution modeling study of the 24 May 2002 case during IHOP. Part I: Numerical simulation and general evolution of the dryline and convection. Mon. Wea. Rev., 134, 149−171.
    Xue, M. and W. Martin, 2006b: A high-resolution modeling study of the 24 May 2002 case during IHOP. Part II: Horizontal convective rolls and convective initiation. Mon. Wea. Rev., 134, 172−191.
  • [1] Hyo-Eun JI, Soon-Hwan LEE, Hwa-Woon LEE, 2013: Characteristics of Sea Breeze Front Development with Various Synoptic Conditions and Its Impact on Lower Troposphere Ozone Formation, ADVANCES IN ATMOSPHERIC SCIENCES, 30, 1461-1478.  doi: 10.1007/s00376-013-2256-3
    [2] HU Dingzhu, TIAN Wenshou, XIE Fei, SHU Jianchuan, and Sandip DHOMSE, , 2014: Effects of Meridional Sea Surface Temperature Changes on Stratospheric Temperature and Circulation, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 888-900.  doi: 10.1007/s00376-013-3152-6
    [3] Qiwei WANG, Ming XUE, Zhemin TAN, 2016: Convective Initiation by Topographically Induced Convergence Forcing over the Dabie Mountains on 24 June 2010, ADVANCES IN ATMOSPHERIC SCIENCES, 33, 1120-1136.  doi: 10.1007/s00376-016-6024-z
    [4] Qi-Wei WANG, Ming XUE, 2018: A High-Resolution Modeling Study of the 19 June 2002 Convective Initiation Case during IHOP_2002: Localized Forcing by Horizontal Convective Rolls, ADVANCES IN ATMOSPHERIC SCIENCES, 35, 1243-1253.  doi: 10.1007/s00376-018-7218-3
    [5] Lujun XU, Huizhi LIU, Qun DU, Yang LIU, Jihua SUN, Anlun XU, Xiaoni MENG, 2021: Characteristics of Lake Breezes and Their Impacts on Energy and Carbon Fluxes in Mountainous Areas, ADVANCES IN ATMOSPHERIC SCIENCES, 38, 603-614.  doi: 10.1007/s00376-020-0298-x
    [6] Qiwei WANG, Yi ZHANG, Kefeng ZHU, Zhemin TAN, Ming XUE, 2021: A Case Study of the Initiation of Parallel Convective Lines Back-Building from the South Side of a Mei-yu Front over Complex Terrain, ADVANCES IN ATMOSPHERIC SCIENCES, 38, 717-736.  doi: 10.1007/s00376-020-0216-2
    [7] PING Fan, GAO Shouting, WANG Huijun, 2003: A Comparative Study of the Numerical Simulation of the 1998 Summer Flood in China by Two Kinds of Cumulus Convective Parameterized Methods, ADVANCES IN ATMOSPHERIC SCIENCES, 20, 149-157.  doi: 10.1007/BF03342059
    [8] YANG Jing, BAO Qing, JI Duoying, GONG Daoyi, MAO Rui, ZHANG Ziyin, Seong-Joong KIM, 2014: Simulation and Causes of Eastern Antarctica Surface Cooling Related to Ozone Depletion during Austral Summer in FGOALS-s2, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 1147-1156.  doi: 10.1007/s00376-014-3144-1
    [9] SHENG Chunyan, Ming XUE, GAO Shouting, 2009: The Structure and Evolution of Sea Breezes During the Qingdao Olympics Sailing Test Event in 2006, ADVANCES IN ATMOSPHERIC SCIENCES, 26, 132-142.  doi: 10.1007/s00376-009-0132-y
    [10] ZHAO Haikun, WU Liguang*, and WANG Ruifang, 2014: Decadal Variations of Intense Tropical Cyclones over the Western North Pacific during 19482010, ADVANCES IN ATMOSPHERIC SCIENCES, 31, 57-65.  doi: 10.1007/s00376-013-3011-5
    [11] ZHU Guofu, CHEN Shoujun, 2003: A Numerical Case Study on a Mesoscale Convective System over the Qinghai-Xizang (Tibetan) Plateau, ADVANCES IN ATMOSPHERIC SCIENCES, 20, 385-397.  doi: 10.1007/BF02690797
    [12] Wu Zengmao, 1987: NUMERICAL STUDY OF LAKE-LAND BREEZE OVER LAKE V?TTERN SWEDEN, ADVANCES IN ATMOSPHERIC SCIENCES, 4, 198-209.  doi: 10.1007/BF02677066
    [13] Tianxue ZHENG, Yongbo TAN, Yiru WANG, 2021: Numerical Simulation to Evaluate the Effects of Upward Lightning Discharges on Thunderstorm Electrical Parameters, ADVANCES IN ATMOSPHERIC SCIENCES, 38, 446-459.  doi: 10.1007/s00376-020-0154-z
    [14] Yang Fanglin, Yuan Chongguang, 1993: Numerical Simulation of Regional Short-Range Climate Anomalies, ADVANCES IN ATMOSPHERIC SCIENCES, 10, 335-344.  doi: 10.1007/BF02658139
    [15] Xie Zhenghui, Dai Yongjiu, Zeng Qingcun, 1999: An Unsaturated Soil Water Flow Problem and Its Numerical Simulation, ADVANCES IN ATMOSPHERIC SCIENCES, 16, 183-196.  doi: 10.1007/BF02973081
    [16] DING Yihui, LI Chongyin, LIU Yanju, 2004: Overview of the South China Sea Monsoon Experiment, ADVANCES IN ATMOSPHERIC SCIENCES, 21, 343-360.  doi: 10.1007/BF02915563
    [17] Honglei Zhang, Honglei Zhang, Ming Xue, Hangfeng Shen, Xiaofan Li, Guoqing Zhai, 2023: Local Torrential Rainfall Event within a Mei-Yu Season Mesoscale Convective System: Importance of Back-Building Processes, ADVANCES IN ATMOSPHERIC SCIENCES.  doi: 10.1007/s00376-023-3033-6
    [18] LIU Yanju, DING Yihui, 2007: Sensitivity Study of the South China Sea Summer Monsoon in 1998 to Different Cumulus arameterization Schemes, ADVANCES IN ATMOSPHERIC SCIENCES, 24, 360-376.  doi: 10.1007/s00376-007-0360-y
    [19] Chen Yuejuan, Zheng Bin, Zhang Hong, 2002: The Features of Ozone Quasi-Biennial Oscillation in Tropical Stratosphere and Its Numerical Simulation, ADVANCES IN ATMOSPHERIC SCIENCES, 19, 777-793.  doi: 10.1007/s00376-002-0044-6
    [20] Jianjun LIU, Feimin ZHANG, Zhaoxia PU, 2017: Numerical Simulation of the Rapid Intensification of Hurricane Katrina (2005): Sensitivity to Boundary Layer Parameterization Schemes, ADVANCES IN ATMOSPHERIC SCIENCES, 34, 482-496.  doi: 10.1007/s00376-016-6209-5

Get Citation+

Export:  

Share Article

Manuscript History

Manuscript received: 07 June 2021
Manuscript revised: 10 November 2021
Manuscript accepted: 10 December 2021
通讯作者: 陈斌, bchen63@163.com
  • 1. 

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

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

Numerical Simulations of a Florida Sea Breeze and Its Interactions with Associated Convection: Effects of Geophysical Representation and Model Resolution

    Corresponding author: Zhaoxia PU, Zhaoxia.Pu@utah.edu
  • 1. Department of Atmospheric Sciences, University of Utah, Salt Lake City, Utah 84112, USA
  • 2. Air Force Institute of Technology, Wright-Patterson Air Force Base, Ohio 45433, USA

Abstract: The Florida peninsula in the USA has a frequent occurrence of sea breeze (SB) thunderstorms. In this study, the numerical simulation of a Florida SB and its associated convective initiation (CI) is simulated using the mesoscale community Weather Research and Forecasting (WRF) model in one-way nested domains at different horizontal resolutions. Results are compared with observations to examine the accuracy of model-simulated SB convection and factors that influence SB CI within the simulation. It is found that the WRF model can realistically reproduce the observed SB CI. Differences are found in the timing, location, and intensity of the convective cells at different domains with various spatial resolutions. With increasing spatial resolution, the simulation improvements are manifested mainly in the timing of CI and the orientation of the convection after the sea breeze front (SBF) merger into the squall line over the peninsula. Diagnoses indicate that accurate representation of geophysical variables (e.g., coastline and bay shape, small lakes measuring 10–30 km2), better resolved by the high resolution, play a significant role in improving the simulations. The geophysical variables, together with the high resolution, impact the location and timing of SB CI due to changes in low-level atmospheric convergence and surface sensible heating. More importantly, they enable Florida lakes (30 km2 and larger) to produce noticeable lake breezes (LBs) that collide with the SBFs to produce CI. Furthermore, they also help the model reproduce a stronger convective squall line caused by merging SBs, leading to more accurate locations of postfrontal convective systems.

摘要: 美国佛罗里达半岛海风和雷暴频发。本文采用中尺度预报模式WRF,模拟研究了不同水平分辨率对佛罗里达海风及其对流触发的影响,在检验海风及其对流模拟能力的基础上,分析了影响海风及其对流触发的关键因素。结果表明,WRF模式可较好地再现出观测的海风及其对流演变。采用不同水平分辨率模拟的对流时间、位置和强度存在明显的差异。提高水平分辨率对模拟的改进主要表现在海风锋在佛罗里达半岛合并为飑线后的对流触发时间和对流移动方向上。进一步的诊断表明,提高水平分辨率对海岸线、海湾形状和尺度为10-30 km2的小尺度湖泊等地形单体的刻画更为精确,从而在模拟改进中发挥了重要的作用。更精细的地形单体可通过改变低层大气辐合和地表感热加热对海风及其对流触发的位置和时间产生重要的影响。值得注意的是,佛罗里达半岛面积为30 km2及以上的湖泊不仅会产生明显的湖风并与海风锋碰撞触发对流,还有利于模式再现出海风合并引起的更强的对流飑线,从而提高锋面后部对流系统位置的模拟。

    • A sea breeze (SB) is a local mesoscale circulation generated when a large temperature gradient between the warmer land surface and the colder ocean surface creates boundary-layer pressure differences, which in turn lead to a diurnal system of breezes that propagate inland (Abbs and Physick, 1992). As an SB propagates inland, its frontal boundary will force the ascent of the continental air mass above the maritime air mass through low-level convergent forcing. If the ascending air reaches its level of free convection (LFC) and maintains sufficient buoyancy through significant upward displacement, shallow cumulus (cu) or deep convective initiation (CI) may occur (Miller et al., 2003; Trier, 2003). Thus, lengthy coastal areas, such as the Florida peninsula, have an almost daily occurrence of summertime SB thunderstorms, making them ideal locations to study SBs and their associated convection. Much progress has been made in understanding CI produced from SBs; however, this and other types of warm-season convection are still among the most difficult atmospheric events to predict as far as the exact timing, location, and intensity of storm cells (Clark et al., 2014; Lock and Houston, 2014).

      The CI of Florida SBs has been investigated in numerous observational and numerical studies (e.g., Nicholls et al., 1991; Abbs and Physick, 1992; Boybeyi and Raman, 1992; Fankhauser et al., 1995; Rao and Fuelberg, 2000; Miller et al., 2003; Etherton and Santos, 2008; Crosman and Horel, 2010). It has been shown that CI is produced along an SB’s convergent boundary, which occurs at discontinuous locations along the sea breeze front (SBF), rather than in an uninterrupted line. These locations of preferred CI have often been shown to occur where the SBF collides or merges with other boundary convergence lines such as thermally driven mesoscale breezes, including sea, lake, river, or land breezes (Blanchard and López, 1985; Boybeyi and Raman, 1992; Laird et al., 1995; Baker et al., 2001), gust front (GF) or outflow boundaries (Nicholls et al., 1991; Fankhauser et al., 1995; Rao and Fuelberg, 2000), or microscale horizontal convective roll (HCR) updrafts (Fankhauser et al., 1995; Rao and Fuelberg, 2000; Ogawa et al., 2003; Fovell, 2005).

      Moreover, it has been found that CI can also occur ahead of or behind the SBF. Observations and numerical simulations by Nicholls et al. (1991) and Fankhauser et al. (1995) noted that CI could develop prior to the merger of two SBFs within an enhanced convergence zone formed as a natural consequence of declining surface heating and the inland progression of the two SBFs. In a high-resolution, three-dimensional numerical simulation conducted by Fovell (2005), an intersecting HCR’s moist updraft interacting with obstacle-generated gravity waves produced a convective roll cloud ahead of the SBF. More recently, Abulikemu et al. (2016) investigated CI prior to the merger of an SBF and gust front in North China and found that low-level convergence, conditional instability, and dynamic vertical forcing were the key ingredients for the CI produced between the two boundaries. Convection behind the stable SBF has been studied in high-resolution numerical simulations and was found to be produced from an outflow boundary intersecting a deep layer of upward motion that was enhanced and carried by microscale Kelvin–Helmholtz billows (KHBs) along the top of the SB gravity current (Rao and Fuelberg, 2000) or intersecting outflow boundaries modified and enhanced by SBF placement (Fankhauser et al., 1995).

      From these studies, it is evident that mesoscale SBs and their associated convection are highly affected by locally driven boundaries and microscale features that change on a daily and even hourly basis, making forecasts of the exact timing, location, and intensity of convective systems difficult. Operational numerical weather prediction (NWP) models have long been used to bridge the gaps in forecasters’ ability to predict warm-season CI and aid in the generation of its forecast (Miller et al., 2003; Hahmann et al., 2006; Crosman and Horel, 2010).

      Previously, numerical studies of Florida’s SBs and their associated convection have used simulations either in large mesoscale domains of several tens of kilometers or in much smaller large-eddy simulation (LES) domains of a hundred meters or less. Few studies have focused on prediction and simulation at the gray-zone scale (e.g., 1-km to 5-km horizontal resolution). Recently, numerical simulations at the gray-zone scale have been commonly used in operational mesoscale weather forecast systems. Research into the predictive performance of numerical models at different spatial resolutions at the gray-zone scale in simulating Florida SB and CI is of pressing importance in aiding future improvements to NWP forecasts in this region. In addition, increasing spatial resolution allows for the improved resolution of geographic features such as terrain height, coastline shape, lake dimensions, etc., which in turn affect geophysical properties such as instability, moisture, and lift that have been shown to be essential for CI (Lock and Houston, 2014). Although the systematic study by Schwartz and Sobash (2019) investigated 3- vs. 1-km forecast performances over the entire CONUS, comprehensive systematic examinations on sea breeze applications have not been performed. In light of the above issues, this paper examines simulation performance for a Florida mesoscale SB and associated convection with different spatial resolutions at the gray-zone scale. Special attention is paid to verifying how the model resolves the timing, location, and intensity of the SB convection as well as the effects of changing geographic features (e.g., coastline shape, water body dimensions, etc.) with different spatial resolutions at the gray-zone scale on convective ingredients and the interactions between the SB and its associated convective characteristics.

      The paper is organized as follows: section 2 briefly describes the real-data case used for the study, while section 3 briefly introduces the Weather Research and Forecasting (WRF) model configuration and setup used in the numerical simulation. Section 4 depicts the simulation results of the Florida SB CI. Section 5 discusses the response of geophysical change to different spatial resolutions at the gray-zone scale and its influence on atmospheric ingredients and CI. Finally, a summary and conclusions are given in section 6.

    2.   A brief overview of the sea breeze case
    • SBs and associated convection developed along both the east and west coasts of the Florida peninsula on 6–7 September 2012. This case was characterized by SB convergence and convective development typified by a fast-moving west coast sea breeze (WCSB) and a slow-moving east coast sea breeze (ECSB). According to Blanchard and López (1985), this is a classic “Type 3” SB pattern, characterized by widespread rainfall amounts and SB convection that is sustained throughout the late evening hours. Figure 1 shows four soundings at 1200 UTC 6 September 2012 [hereafter all times are in UTC; UTC = EDT (Eastern Daylight Time) + 4 h] from four Florida stations: Tallahassee (TLH, Fig. 1a), Jacksonville (JAX, Fig. 1b), Tampa Bay (TBW, Fig. 1c), and Miami (MFL, Fig. 1d). These soundings represent the atmospheric conditions about an hour after sunrise, several hours before the SBs and their associated convection are initiated. All four soundings show convective inhibition (CIN) values of less than –25 J kg–1 and convective available potential energy (CAPE) values ranging from almost 900 J kg–1 (TLH, Fig. 1a) to over 2500 J kg–1 (MFL, Fig. 1d). The level of free convection (LFC) of the four soundings is between 895 hPa and 850 hPa, while the lifting condensation level (LCL) of the four soundings ranges between 970 hPa and 950 hPa. Precipitable water (PW) values for the soundings are around 48 mm, except for JAX (Fig. 1b), which is slightly less than the other soundings, at 42 mm. Consequently, the soundings across the peninsula show a warm, moist, and unstable atmosphere favorable for the SB convection event, given that the weak subsidence and associated CIN over the peninsula could be overcome in the afternoon through convective heating and lifting effects.

      Figure 1.  Observed soundings for four Florida stations: (a) Tallahassee (TLH), (b) Jacksonville (JAX), (c) Tampa Bay (TPA), and (d) Miami (MIA) at 1200 UTC 6 September 2012. Images courtesy of University of Wyoming [http://weather.uwyo.edu/upperair/sounding.html]. Locations of sounding stations can be found in Fig. 3a.

      A day before this event (5 September 2012), SB convection occurred during the afternoon and extended into the evening hours (not shown). However, by 0600 UTC 6 September 2012, the sky over Florida was mostly clear and remained cloud free for the rest of the morning, allowing for plenty of solar insolation to occur over the peninsula. Both the WCSB and ECSB developed between 1320 and 1400 UTC 6 September 2012, while convection began along the WCSB at 1400 UTC (Fig. 2a) and along the ECSB at 1500 UTC (Fig. 2b) 6 September 2012, respectively. The enhanced WCSB and its associated convection traveled significantly inland, while the ECSB and convection remained adjacent to the east coast for most of this period (Fig. 2c); consequently, the two SBs collided just inland from the eastern coast of Florida at approximately 2030 UTC 6 September 2012 (Fig. 2d), producing a strong squall line of enhanced deep convection. The deep convection produced by the colliding SBs was oriented in a north-to-south line along Interstate 95 from Jacksonville (JAX) to Vero Beach (VRB), while the convection to the south was oriented in a northeast-to-southwest line across the peninsula from VRB to just inland and north of Naples (APF). Shortly following the production of this squall line at 2100 UTC 6 September 2012, there were several reports of −27 m s−1 winds along the southeast coast of Florida. The last remnants of SB convective precipitation ended at 0500 UTC 7 September 2012 in east-central portions of the Florida peninsula (not shown).

      Figure 2.  Geostationary Operational Environmental Satellite-East (GOES-E) visible satellite images of convective cloud cover over the Florida peninsula at (a) 1401 UTC, (b) 1515 UTC, (c) 2001 UTC, and (d) 2031 UTC 6 September 2012. Red lines show the location of the WCSB and ECSB fronts. Yellow line in panel (d) indicates formation of the squall line due to the collision of the SBFs. Images courtesy of College of DuPage [http://www.cod.edu/] and Carl Jones [https://northflwx.wordpress.com/2012/10/01/florida-coast-to-coast-sea-breeze-convergence/].

    3.   Simulation and data description
    • The numerical simulations of SB convection presented in this paper were performed using the Advanced Research WRF (WRF-ARW) model (version 3.7.1, Skamarock et al., 2008). A four-level, one-way nested domain configuration centered over the Florida peninsula is employed for the simulation. The outermost, two middle, and innermost horizontal domain grids consist of grid spacings of 27 km (d01), 9 km (d02), 3 km (d03), and 1 km (d04), respectively (not shown). Note that the spatial resolutions of d03 and d04 are both at the gray-zone scale. The outermost domain (d01) and d02 were chosen to cover the southeastern US and the surrounding Atlantic Ocean in order to capture dynamics that might influence the SB and convection. The innermost domain covers the Florida peninsula and its surrounding local waters (see areas covered by d03 and d04 in Figs. 3a and 3b).

      Figure 3.  Terrain heights and water body locations for (a) d03 and (b) d04. Red lines indicate the locations of the vertical cross sections cut west-to-east through MacDill, Air Force Base (MCF, northernmost) and Lake Okeechobee (OBE, southernmost). The colored boxes indicate locations used later in the study for a zoomed-in analysis. The locations of the 23 observing stations used to aid in the verification of the WRF simulation are also plotted. Note the increased number of lakes in (b) compared to (a). Water body names and locations are annotated in panel (c). The left corner inset of panel (c) is a zoomed-in area over Tampa Bay. (d) Locations of WRF model domains for all four-level nested domains.

      All four domains contain 62 vertical eta levels with 26 levels in the lower atmospheric region (below 850 hPa). The model top is set at 5 hPa. The physics options used in this study include the Unified Noah land surface model (LSM, Chen and Dudhia, 2001), Mellor-Yamada-Janjić (MYJ, Mellor and Yamada, 1974) PBL, New Thompson microphysics (Thompson et al., 2008), New Kain-Fritsch (KF, Kain, 2004) cumulus (d01 and d02 only), MM5 Dudhia shortwave radiation (Dudhia, 1989), and Rapid Radiative Transfer Model (RRTM) longwave radiation (Iacono et al., 2008). Table 1 gives a full list of the configurations used for each of the four domains. Different configurations of physical parameterization impact the onset, intensity, and merging of the SBFs. The selection of the physical parameterization options was based on sensitivity experiments with various physical schemes in the early phase of this study (not shown).

      Parameterd01d02d03d04
      Physical domain139 × 97 (lon. × lat.)253 × 202 (lon. × lat.)343 × 313 (lon. × lat.)502 × 661 (lon. × lat.)
      Horizontal resolution27 km9 km3 km1 km
      Vertical resolution62 levels with variable Δz (26 levels below 850 hPa): model top 50 hPa*
      Time integration120-sec time step; 36-h duration*
      Boundary conditionsDamping depth over top 500 m*
      Lateral boundary21 600-sec time interval*
      Surface layer physicsMonin-Obukhov Janjić Eta (MO-JE)*
      MicrophysicsNew Thompson*
      Longwave radiationRapid Radiative Transfer Model (RRTM)*
      Shortwave radiationMM5 (Dudhia)*
      PBL physicsMellor-Yamada-Janjić (MYJ)*
      Land-surface physicsNoah land surface model (LSM)*
      Cumulus physicsNew Kain-FritschNew Kain-FritschNoneNone
      Vertical Diffusion2nd order diffusion*
      Time-integrationRunge-Kutta 3rd order*
      Eddy coefficientHorizontal Smagorinsky 1st order closure*
      *Values are the same for all four grids.

      Table 1.  WRF model parameter settings for domains d01, d02, d03, and d04.

      The model initial and boundary conditions are provided by the NCEP North American Mesoscale (NAM) Forecast System analysis at 12-km grid spacing from 0000 UTC 6 September 2012 to 1200 UTC 7 September 2012 at six-hour intervals. After initialization, the model was integrated into a one-way nested mode for 36 hours until 1200 UTC 7 September 2012. The simulated SB events are verified using surface parameters obtained from local aviation routine weather reports (METARs) and aviation special weather reports (SPECIs) from 23 locations across the peninsula. A detailed comparison of the simulated results is made with NCEP’s Stage IV precipitation analysis (Lin and Mitchell, 2005) and the Climatology-Calibrated Precipitation Analysis (CCPA) (Hou et al., 2014). Furthermore, NCEP’s 32-km horizontal grid spacing/45-level vertical resolution North American Regional Reanalysis (NARR; Mesinger et al., 2006) is also used to verify the WRF simulation output.

    4.   Simulation results
    • The SB and its associated convection were examined on a five-minute basis for the WRF simulation. As both the d03 and d04 grid scales were very similar in their ability (spatially and temporally) to resolve the SB, only the results from d03 are shown in this section. It is noteworthy that the higher resolution of d04 did allow the SBF to more closely follow the coastline shape upon initiation and produced a more intense and narrower SBF, which will be discussed in later sections.

    • Figure 4 shows the initiation and evolution of the WRF-simulated ECSB and WCSB as seen by the zonal surface winds and total wind vectors. Specifically, the WRF-simulated SB is initiated between 1405 and 1435 UTC 6 September 2012 on both sides of the peninsula (Figs. 4a and 4b), as seen in the turning of the coastal winds, which is in line with the SB initiation time in the surface observations and NARR data. By 1500 UTC, both SBs are well established (Fig. 4c) and continue to strengthen and propagate through the afternoon. Due to the opposing low-level synoptic flow, the ECSB is quite discernible in the simulation, while the WCSB that is parallel to the synoptic flow becomes much harder to discern in the analysis (Fig. 4d). Simulated surface temperature also indicates that coastal areas are warmer than the nearby ocean before the SB, which leads to surface pressure decreases over land and then the land breeze turning into an SB (not shown).

      Figure 4.  WRF-simulated zonal (u-component) winds (color: red easterly, blue westerly; units: m s–1) and total wind (black arrows; unit reference vector 2.5 m s–1) at surface level for a zoomed-in region of d03 (yellow box in Fig. 3a) at: (a) 1300 UTC, (b) 1400 UTC, (c) 1500 UTC, and (d) 1600 UTC 6 September 2012.

      Figure 5 shows a temporal evolution of the zonal wind cross sections, indicating the inland progression and height of the ECSB and WCSB during the afternoon. Initially, the simulated WCSB extends higher (700 m) than the ECSB (400 m, Fig. 5a). However, by 1700 UTC 6 September 2012, the ECSB has reached 700 m at its frontal head, while the WCSB extends 900 m (Fig. 5b). The two SBFs reach a maximum height of around 1 km by 1930 UTC (Fig. 5c). The two SBs gradually approach each other (Fig. 5d) and eventually meet inland of the eastern coast of Florida at 2120 UTC 6 September 2012, which is 50 minutes after the SB frontal merger in surface observations.

      Figure 5.  Cross section (through northernmost red line in Fig. 3a) of d03 WRF-simulated zonal (u-component) winds (color: red easterly, blue westerly; units: m s–1) and total horizontal wind (black arrows represent horizontal winds at each height level along longitude, reference vector 10 m s–1; upward arrows mean north direction) at: (a) 1600 UTC, (b) 1700 UTC, (c) 1930 UTC, and (d) 2030 UTC 6 September 2012. Red contour lines indicate the zero-wind contour heights, which are indicative of approximate SBF heights.

    • Figures 6a and 6b compare the WRF-simulated accumulated precipitation to CCPA data from 1200 UTC 6 September 2012 to 0600 UTC 7 September 2012 (18-h accumulated precipitation in total) for domains d03 and d04. Both of the WRF simulation domains and the CCPA analysis observe relatively sparse areal coverage of rainfall along the west coast and in northern Florida, while distinct patterns of rainfall are observed along the east coast and in central Florida. In addition, an obvious lack of rainfall occurs northeast of Lake Okeechobee in both the CCPA analysis and the WRF simulations, although the “rain shadow” is more pronounced in d03 (Fig. 6a) than in d04. Meanwhile, light rainfall occurs across the north-central and northeast portions of the Florida peninsula in the d04 simulation (Fig. 6b), while the d03 grid scale shows almost no precipitation in those areas (Fig. 6a).

      Figure 6.  18-h rainfall accumulation totals from 1200 UTC 6 September 2012 to 0600 UTC 7 September 2012 for (a) d03 WRF-simulated precipitation accumulation totals (blue contours; units: mm) and CCPA data (color; units: mm). Panel (b) is the same as (a) but for d04. WRF-simulated precipitation accumulation totals (color; units: mm) for (c) d03 and (d) d04 for the same times as (a) and (b).

      Figure 7 further compares the hourly WRF precipitation accumulations for d03 and d04 with hourly Stage IV data. Simulated timing of the initial convection along the west coast of Florida at 1400 UTC 6 September 2012 for an inland location just northeast of Waccasassa Bay (see Fig. 3c for the location of the bay) in d03 and d04 (Figs. 7a and 7b) is concurrent with observations; however, d03 fails to resolve the convection produced east of Waccasassa Bay, contrary to the Stage IV analysis. In contrast, the higher resolution d04 is able to capture the convection east of Waccasassa Bay (Fig. 7b) but also overestimates the WCSB convection just north of Naples (APF in Figs. 3a and 3b). At this time, there are also small convective cells not associated with the SBF in the central portion of the peninsula in d04 (Fig. 7d) that are not present in d03 (Fig. 7c). By 1800 UTC 6 September 2012, d04 resolves the convection along the east coast of Florida, while d03 has only a small convective cell in the southeast corner of the peninsula associated with the ECSB (Figs. 7c and 7d). It is not until 1900 UTC 6 September 2012 that d03 begins to resolve significant convection along the ECSB.

      Figure 7.  Stage IV precipitation accumulation (color; units: mm) and d03 (left column) and d04 (right column) WRF-simulated precipitation accumulation (black contours; units: mm) at: (a–b) 1400 UTC, (c–d) 1800 UTC, (e–f) 2100 UTC, and (g–h) 2300 UTC 6 September 2012. The red arrows in (e) and (f) indicate squall lines.

      The WRF-simulated SBFs collide and merge at 2030 UTC 6 September 2012 and produce an enhanced convective squall line that is clearly visible at 2100 UTC 6 September 2012, which is markedly similar to the Stage IV analysis convective pattern (Figs. 7e and 7f). However, from 2230 UTC 6 September 2012 until 0200 UTC 7 September 2012, the orientation of the convection in the d03 WRF simulation becomes north-to-south, while the Stage IV analysis data clearly show the convective line in southern Florida maintaining a northeast-to-southwest orientation (Fig. 7g). The convection produced in d04 better captures this northeast-to-southwest convective orientation (Fig. 7h). Finally, the duration of the convection across the Florida peninsula in the d03 simulation lasts until 0530 UTC 7 September 2012, rather than 0500 UTC 7 September 2012 as seen in Stage IV and CCPA analyses, or 0400 UTC as seen in the d04 simulation (not shown). A comparison of 18-h accumulated precipitation (1200 UTC 6 September 2012 to 0600 UTC 7 September 2012) between the d03 and d04 WRF simulations (Figs. 6c and 6d) and the CCPA analysis (Figs. 6a and 6b) shows that despite good spatial coverage of the WRF-simulated convective precipitation, the maximum precipitation accumulations produced by the WRF simulation cover much smaller spatial areas than what is observed in the CCPA analysis.

    • To further examine the intensity differences between the two spatial resolutions at the gray-zone scale, a quantitative precipitation forecasting (QPF) analysis was conducted using the Equitable Threat Score (ETS), also known as the Critical Success Index (CSI), and Bias Score (BS) calculations. The values used to calculate these scores are given in a 2 × 2 contingency table (Table 2). Four specific precipitation thresholds are used in the QPF statistics calculations: 2.54 mm, 6.35 mm, 12.7 mm, and 25.4 mm. Knowing information about the forecast area (F), observed area (O), and the correctly forecasted “hits” (H), the CSI seeks to answer the question of how well the forecasted “hits” correspond to the observed “hits”. The BS seeks to answer the question of how similar the frequencies of forecasted and observed “hits” are. The BS is mathematically defined as:

      ForecastObserved
      YesNoSum
      YesHits (H)False alarmsAll forecasted (F)
      NoMissesCorrect rejectionsNo forecasted
      (NF)
      SumAll observed (O)No observed
      (NO)
      Total (N)

      Table 2.  Contingency table illustrating the counts used in verification statistics of dichotomous (e.g., Yes/No) forecasts and observations.

      The ETS score can be obtained using the following equation:

      For a given threshold, A represents the number of grid points that exceed the threshold in both the model forecast and the CCPA data; B denotes the number of grid points that exceed the threshold in the model forecast, but not in the CCPA data; and C is the number of grid points that do not reach the threshold in the model forecast, but that exceed the threshold in the CCPA analysis.

      ETS values range from 0 to 1, with 0 indicating no skill in the forecast. BS values range from 0 to infinity, with a bias of less than 1 indicating a tendency to underforecast, while a bias greater than 1 indicates a tendency to overforecast in the forecast system. The ETS and BS values for the d03 and d04 simulation results are given in Table 3. As seen by the relatively low ETS values in both domains for all four precipitation thresholds, the simulations are relatively poor at forecasting the correct intensities for the SB event, and the BS shows that the simulations underforecast the precipitation intensity events. d04 performs slightly better than d03 at nearly all precipitation thresholds.

      Precipitation
      Threshold (mm)
      Domain
      d03d04
      Equitable Threat Score (ETS)2.540.120.28
      6.350.10.19
      12.70.080.15
      25.400.04
      Bias Score (BS)2.540.120.26
      6.350.110.2
      12.70.120.17
      25.40.140.23

      Table 3.  Equitable Threat and Bias Score values for four rainfall thresholds over an 18-h forecast period (1200 UTC 6 September 2012 to 0600 UTC 7 September 2012).

      Thus, it is evident that the WRF simulation is able to capture the overall spatial locations and timing of the convective mesoscale systems (MCSs) that occurred during the 6–7 September 2012 SB event reasonably well. However, the intensities of the WRF-simulated convective precipitation are overall much weaker as compared to the CCPA and Stage IV analyses in subjective and QPF evaluations. It is also found that increasing spatial resolution at the gray-zone scale has impacts on convection activities: namely, increasing the spatial resolution at the gray-zone scale can reasonably reproduce and capture the timing of the CI, as well as the orientation of the convection after the SBF merger into the squall line over the peninsula, although it degrades the simulation of the size and organization of convective cells.

      The spatial resolution of the gray-zone scale commonly ranges from 1 km to 5 km; therefore, different spatial resolutions mean different geophysical features, which in turn affect atmospheric features. In order to further understand the geophysical features at different spatial resolutions and their impacts on the simulation of spatial, timing, and intensity of MCSs and precipitation, additional analysis of the atmospheric and convective features resolved at these two grid scales (d03 and d04 domains) is performed in the next section.

    5.   Influences of geographical features and model resolution
    • A significant difference in the ability of the two grids (d03 and d04) to resolve the observed convective precipitation occurs east of Waccasassa Bay at 1600 UTC 6 September 2012. D03 is unable to resolve this convection, while d04 is able to produce convection in this area. At 1500 UTC 6 September 2012, the WCSB is initiated in both gray-zone domains, as seen in u wind and total wind vectors. d04’s higher-resolution grid scale allows the WCSB to be initiated more closely along the coastline shape and with slightly stronger winds, leading to areas of enhanced convergence (Figs. 8b and 8d) that are absent in d03 (Figs. 8a and 8c). This narrow band of convergence produced by the SBF in d04 leads to the CI seen east of Waccasassa Bay (Fig. 8d), which is absent in d03. Simulation results also show that over the ocean surface, west of the Florida peninsula, the structures of the divergence field are similar to the less organized open boundary layer convective cells. At the same time, over the ocean surface, east of the Florida peninsula, the structures of the divergence field are similar to well-organized boundary layer horizontal convective rolls. These results agree with those discussed in Xue and Martin (2006a, b) and can be further addressed in future work.

      Figure 8.  WRF-simulated divergence (color: red convergence, blue divergence; units: s–1) and total wind (black arrows; reference vector 2.5 m s–1) at surface level for a zoomed-in region of d03 (yellow box in Fig. 3a) at: (a) 1500 UTC and (c) 1530 UTC 6 September 2012. Panels (b) and (d) are the same as a) and (c) respectively, but for the zoomed-in region of d04 (yellow box in Fig. 3b) and include max reflectivity (red contours; units: dBZ).

      In addition, the flow of the winds through the Y-shaped Tampa Bay area (see inset in Fig. 3c) to the south at this same time also shows significant enhancement in d04 due to its better resolution of the bay’s coastline. The marine air flowing through Tampa Bay is deflected to the northwest and northeast around the Interbay Peninsula into Old Tampa Bay and Hillsboro Bay, respectively, due to the SB flow. This produces strong convergent flow from the SB along the Pinellas and Interbay Peninsulas. The Tampa Bay breeze and WCSB create a narrow band of enhanced surface convergence, which aids in lifting air to its LFC and initiating convection at this time, as seen in the bright red reflectivity contours just north of the Pinellas Peninsula at 1500 UTC 6 September 2012 (Fig. 8d).

    • In addition to the bay breezes and coastlines affecting the CI in the gray-zone simulations, lake breezes (LBs) also play an important role in CI. LBs, like SBs, are frequent in the summertime when the cooler lakes create a pressure gradient effect that generates a lake-to-land flow of air. Lakes are known to significantly impact convection during summer by producing “rain shadows” on the downwind side of the lake (Segal et al., 1997). This occurs when the synoptic wind blows over the cooler lake surface and the air is cooled and moistened, generating subsidence and the suppression of the convective boundary layer, leading to more stable conditions as the wind begins to flow back over land (Segal et al., 1997). It has long been established that Lake Okeechobee, Florida’s largest lake at 1825 km2, affects convection in terms of LB and rain shadow effects, in both observations and numerical experiments (e.g., Blanchard and López, 1985; Boybeyi and Raman, 1992; Segal et al., 1997). As previously mentioned in section 3b, the Lake Okeechobee “rain shadow” is also seen in observations on 6 September 2012, as well as in the gray-zone WRF simulations (Figs. 6a and 6b), suppressing convection directly downwind of the lake. In addition, both grid spacings resolve many of the smaller Florida lakes, not merely Lake Okeechobee. Little research has been done on how the smaller Florida water bodies affect the SB, with the exception of studies done on the smaller water bodies surrounding Cape Canaveral (e.g., Laird, 1995; Rao and Fuelberg, 2000). In both gray-zone WRF simulations, midsized lakes (lakes between 30 km2 and 200 km2), the St. Johns River, and the Indian River are all resolved and effectively produce lake/river breezes. Small lakes (i.e., lakes between 10 km2 and 30 km2) are not resolved in d03, but the higher resolution d04 is able to resolve these lakes and produce associated LBs (see Fig. 9). Very small lakes (less than 10 km2) are not resolved at either gray-zone grid scale.

      Figure 9.  Temperature difference [color contours; units: oF , oF = (oC × 9/5) + 32] between the d03 and d04 WRF simulations at 0330 UTC 6 September 2012 for a zoomed-in region of d03 and d04 (yellow boxes in Figs. 3a and 3b).

      The midsized lakes (30–200 km2) in d03 and the small lakes (10–30 km2) in d04 both produce LBs of similar strength. The depths of the LBs produced by these smaller lakes are on the scale of 150–350 m. Figure 10 shows a cross section through the center of the peninsula at 1600 UTC 6 September 2012. The d04 WRF-simulated WCSB and ECSB frontal heads (indicated by blue arrows) show depths of ~600 m to 750 m, while in comparison, the LBs (indicated by red arrows) show depths from ~250 m for Lake Weohyakapka (30 km2) and Lake Marian (23 km2), to ~375 m for the southern tip of Lake Kissimmee (140 km2).

      Figure 10.  Cross section (through southernmost red line in Fig. 3b) of d04 WRF-simulated zonal (u-component) winds (color: red easterly, blue westerly; units: m s–1) and total horizontal wind (black arrows represent horizontal winds at each height level along longitude, reference vector 10 m s–1; upward arrows mean north direction) at 1600 UTC 6 September 2012. Red contour line indicates the zero-wind contour heights, which indicate Lake Okeechobee’s approximate lake breeze heights. Blue arrows indicate the locations of the SBs, while the red arrows indicate the lake breezes for 1) Lake Weohyakapka, 2) Lake Kissimmee, 3) Lake Marian, and 4) Indian River.

      In both WRF gray-zone simulations, the LBs produced by the midsized and small lakes in d03 and d04 are important areas for CI. LBs produce convergent boundary areas on the windward side of the lake that extend in an elongated semicircle around the sides of the lake (due to interactions with the opposing synoptic flow). In the WRF simulations, these convergent frontal areas produced by the LB interacting with the synoptic winds intersect with other boundary areas, such as SBs or other LBs, producing CI. A clear example of this is seen in the interactions of the ECSB, Lake George, Crescent Lake, and Lake Apopka at the d03 grid scale (Fig. 11). At 1825 UTC 6 September 2012, a strong area of convergence from Lake George and Crescent Lake perpendicularly intersects the ECSB, producing a convective cell (Cell 1 in Fig. 11). Twenty minutes later, at 1845 UTC, the convergence area between Lake George and Crescent Lake generates a convective cell (Cell 4 in Fig. 18), while just fifteen minutes later, the Crescent Lake breeze again intersects the ECSB (north of the lake), producing a convective cell at 1900 UTC (Cell 5 in Fig. 11). A final area of CI due to the LB–ECSB collision in the example occurs at 1910 UTC 6 September 2012 due to the convergent area produced by Lake Apopka (Cell 7 in Fig. 11).

      Figure 11.  D03 WRF-simulated divergence (color: red convergence, blue divergence; units: s–1) at surface level, max reflectivity (yellow and purple contours; units: 5 dBZ), and total wind (black arrows; reference vector 2.5 m s–1) at surface level for a zoomed-in region of d03 (purple box in Fig. 3a) at: (a) 1820 UTC, (b) 1825 UTC, (c) 1840 UTC, (d) 1845 UTC, (e) 1900 UTC, and (f) 1910 UTC 6 September 2012. Yellow reflectivity contour lines indicate the convective cell’s first appearance in the simulation and are numbered accordingly in panels (b–f); thereafter, the reflectivity contour lines are purple. Panel (a) indicates the lake effect divergence from the following water bodies: A—St John’s River, B—Crescent Lake, C—Lake George, D—Lakes Griffin, Harris, and Eustis, E—Lake Apopka. See text for more details.

      In addition to convection produced by the LBs intersecting with other boundary areas, the interaction of the LBs with synoptic winds also leads to CI. Figure 11f shows this interaction. As the LB produced by Lake Apopka interacts with synoptic winds, the convergence area surrounding the lake becomes more and more elongated downstream of the lake in narrow bands. As these convergent bands narrow, vertical velocities also increase, causing CI to occur on the synoptic-wind parallel sides of the lake (Cell 8 in Fig. 11f). Because the higher resolution of d04 is able to resolve more of the smaller Florida lakes, the d04 simulation shows more numerous convergent LB boundary areas, leading to more convective cells being initiated along these LB convergence lines (as seen in Fig. 7c) in central Florida, compared to the d03 simulation (Fig. 8d).

    • Besides differences in the ability of gray-zone grids to resolve the convection between the SBFs, postfrontal squall line convection is also quite different between the two simulations. In both high-resolution simulations of the innermost domains, the observed squall line produced by colliding SBs at 2130 UTC 6 September 2012 is captured reasonably well in terms of timing and spatial location (Figs. 8e and 8f). However, following the generation of this squall line from the colliding SBs (until 0100 UTC 7 September 2012), the convective system’s orientation in d03 evolves quite differently from the orientation of the observed convective system and d04 simulation results (Figs. 8g and 8h). CCPA and Stage IV analyses show that the observed post-squall line convective system tends to have a northeast–southwest orientation, while the d03 simulation creates a distinct north–south band of convective cells instead.

      The discrepancies in the post-squall line convective orientation are due to differences in the strength of the squall lines produced in the two gray-zone simulations. Figure 12 shows the zonal winds of d03 and d04, just after the SBF merger and generation of the squall line convection. It is evident that the SBF and associated convective cells produced by the d04 output are much stronger than those produced by the d03 output along the squall line. An hour later, at 2200 UTC 6 September 2012, the northern half of the d04 squall line continues to propagate eastward (Fig. 12d), while the weaker squall line in d03 remains essentially motionless in the east-central portion of the peninsula (Fig. 12c). By 2300 UTC 6 September 2012, the d03 post-squall line convection is still in the central portion of the peninsula (Fig. 12e), while the d04 post-squall line convection system has propagated eastward to Florida’s eastern coast (Fig. 12f).

      Figure 12.  WRF-simulated zonal (u-component) winds (color: red easterly, blue westerly; units: m s–1) and total wind (black arrows; reference vector 2.5 m s–1) at surface level for a zoomed-in region of d03 (yellow box in Fig. 3a) at (a) 2100 UTC, (c) 2200 UTC, and (e) 2300 UTC 6 September 2012. Panels (b), (d), and (f) are the same as (a), (c), and (e) respectively, but for the zoomed-in region of d04 (yellow box in Fig. 3b). The black circles denote the locations of squall lines.

      D04’s stronger SBF is induced in part by the higher-resolution grid’s ability to resolve a slightly higher (1°C–2°C) land-surface temperature (not shown), that is, at 1230 UTC 6 September 2012, the d04 simulation begins to show higher surface temperatures (1°C–3°C over the d03 surface temperatures) along the southern tip of Florida. By 1400 UTC 6 September 2012, d04’s higher surface temperature has extended along the western coast of Florida, and to Florida’s east coast by 1700 UTC 6 September 2012. As the atmosphere above the land and large bodies of water is sensitive to the specified temperature (Knievel et al., 2010), even this small temperature difference is enough to produce a stronger land–ocean temperature gradient in d04.

      A larger land–ocean surface temperature gradient (as seen in the d04 simulation) will lead to increased low-level flow within the SBF and increased low-level convergence at the SBF (for offshore flow). This difference in low-level convergence between the d03 and d04 resolutions can also be clearly seen in the WRF-simulated divergence fields (Fig. 8). The d03 ECSB band of frontal convergence (Figs. 8a and 8c) is much broader than the d04 ECSB band of frontal convergence (Figs. 8b and 8d), showing a much weaker d03 ECSB kinematic front (location of the SBF’s maximum near-surface wind convergence; Miller et al., 2003). Based on the wind field (Fig. 12) and divergence field (Fig. 8), it is clear that d04’s enhancement of the SBFs, as compared to d03, also leads to increased vertical velocities at the SBF, leading to stronger CI and updrafts along the front as well. When the two simulated SBFs merge into the squall line between 2030 UTC and 2130 UTC 6 September 2012, the stronger SBF and convective cells in d04 aid the propagation of the convection system eastward, rather than it stagnating the MCS in the central peninsula as seen in the d03 simulation of the post-squall line convection. It should be noted that the simulated discrepancies of the divergence fields (Fig. 8) correspond well with the differences in the resolved and unresolved lake distribution between d04 and d03, suggesting that the refined geophysical features by the grid resolution resulting from increasing spatial resolution are important for the simulation of Florida SB and CI.

    6.   Summary and concluding remarks
    • This study simulated an enhanced WCSB, or “Type 3” SB convective event (Blanchard and López, 1985), that occurred 6–7 September 2012 using the WRF model at the gray-zone grid scale with two nested inner-domain grids of 3 km and 1 km, respectively. Results indicate that the WRF model can properly simulate SBs and associated convection. Although the two spatial resolutions are both at the gray-zone scale, they have different impacts on the Florida SB and CI. In other words, the simulation improvements of d03 are manifested mainly in the organization of the convective system, while the simulation improvements of d04 are manifested mainly in the timing and location of CI along the WCSB.

      The numerical results were further investigated to see how geophysical features respond to different spatial resolutions at the gray-zone scale and how they impact convection. Results indicate that more accurately resolving geophysical features such as bays, coastlines, and small lakes (10–30 km2) by the high-resolution grids has significant impacts on the location and timing of CI along SBFs. Compared with the 3-km grid scale, the 1-km gray-zone grid scale enables better resolution (i.e., more grid points for better representation) of the Florida coastline and the smaller Florida lakes (10–30 km2), which leads to increased areas of atmospheric low-level convergence, and in turn, to stronger SBFs and LBs. The enhanced SB and LB fronts are key locations for CI; they are clearly seen in d04 but not seen in the coarser WRF grid domains (d03). In addition, the refined geophysical features in d04 also produce stronger SBFs, leading to a more intense squall line being produced following the SBF mergers. This stronger squall-line convection in d04 is key in determining more accurate postfrontal locations of the convective systems as compared to the d03 grid scale. In summary, improving the representation of geophysical features such as bays, coastlines, and small lakes (10–30 km2) leads to enhancements of low-level convergence and surface sensible heating, which are key ingredients to strengthening SBFs and aiding CI.

      The results of this study show that CI due to the Florida SB is sensitive to geophysical features that are closely related to the grid resolution of numerical models at the gray-zone scale. As this study presents one convective SB case, additional case studies are still needed, especially for simulations of Type 1 and Type 2 SB cases (Blanchard and López, 1985). In addition, an ensemble of simulations with various parameterizations or perturbed initial conditions as well as data assimilation could be useful areas for analysis and forecasts of SB convection in the future work.

      Acknowledgements. The first author gratefully acknowledges the support of the US Air Force Institute of Technology (AFIT) program during her graduate study. The views expressed in this article are those of the stated authors and do not reflect the official policy or position of the US Air Force, Department of Defense, or the US government.

      Computing resources for the simulations were provided by the Center for High-Performance Computing (CHPC) at the University of Utah. Stage IV precipitation data were provided by the NCAR/Earth Observing Laboratory (EOL) at its website http://data.eol.ucar.edu/. CCPA precipitation data were provided by the NCEP/Environmental Modeling Center (EMC) at its website http://www.emc.ncep.noaa.gov. NCEP reanalysis data were provided by the National Oceanic and Atmospheric Administration (NOAA)/Oceanic and Atmospheric Research (OAR)/Earth System Research Laboratory (ESRL) Physical Science Division (PSD) from its website, http://www.esrl.noaa.gov/psd/. ERA-Interim reanalysis data were provided by ECMWF from its website, http://apps.ecmwf.int/.

Reference

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return