Mellor Yamada Nakanishi Niino Research Articles

Evaluation of Nine Planetary Boundary Layer Turbulence Parameterization Schemes of the Weather Research and Forecasting Model Applied to Simulate Planetary Boundary Layer Surface Properties in the Metropolitan Region of São Paulo Megacity, Brazil

This study evaluates nine Planetary Boundary Layer (PBL) turbulence parameterization schemes from the Weather Research and Forecasting (WRF) mesoscale meteorological model, comparing hourly values of meteorological variables observed and simulated at the surface of the Metropolitan Region of São Paulo (MRSP). The numerical results were objectively compared with high-quality observations carried out on three micrometeorological platforms representing typical urban, suburban, and rural land use areas of the MRSP, during the 2013 summer and winter field campaigns as part of the MCITY BRAZIL project. The main objective is to identify which PBL scheme best represents the diurnal evolution of conventional meteorological variables (temperature, relative and specific humidity, wind speed, and direction) and unconventional (sensible and latent heat fluxes, net radiation, and incoming downward solar radiation) on the surface. During the summer field campaign and over the suburban area of the MRSP, most PBL scheme simulations exhibited a cold and dry bias and overestimated wind speed. They also overestimated sensible heat flux, with high agreement index and correlation values. In general, the PBL scheme simulations performed well for latent heat flux, displaying low mean bias error and root square mean error values. Both incoming downward solar radiation and net radiation were also accurately simulated by most of them. The comparison of the nine PBL schemes indicated the local Mellor-Yamada-Janjic (MYJ) scheme performed best during the summer period, particularly for conventional meteorological variables for the land use suburban in the MRSP. During the winter field campaign, simulation outcomes varied significantly based on the site’s land use and the meteorological variable analyzed. The MYJ, Bougeault-Lacarrère (BouLac), and Mellor-Yamada Nakanishi-Niino (MYNN) schemes effectively simulated temperature and humidity, especially in the urban land use area. The MYNN scheme also simulated net radiation accurately. There was a tendency to overestimate sensible and latent heat fluxes, except for the rural land use area where they were consistently underestimated. However, the rural area exhibited superior correlations compared to the urban area. Overall, the MYJ scheme was deemed the most suitable for representing the convectional and nonconventional meteorological variables on the surface in all urban, suburban, and rural land use areas of the MRSP.

Offshore low-level jet observations and model representation using lidar buoy data off the California coast

Abstract. Low-level jets (LLJs) occur under a variety of atmospheric conditions and influence the available wind resource for wind energy projects. In 2020, lidar-mounted buoys owned by the US Department of Energy (DOE) were deployed off the California coast in two wind energy lease areas administered by the Bureau of Ocean Energy Management: Humboldt and Morro Bay. The wind profile observations from the lidars and collocated near-surface meteorological stations (4–240 m) provide valuable year-long analyses of offshore LLJ characteristics at heights relevant to wind turbines. At Humboldt, LLJs were associated with flow reversals and north-northeasterly winds, directions that are more aligned with terrain influences than the predominant northerly flow. At Morro Bay, coastal LLJs were observed primarily during northerly flow as opposed to the predominant north-northwesterly flow. LLJs were observed more frequently in colder seasons within the lowest 250 m a.s.l. (above sea level), in contrast with the summertime occurrence of the higher-altitude California coastal jet influenced by the North Pacific High, which typically occurs at heights of 300–400 m. The lidar buoy observations also validate LLJ representation in atmospheric models that estimate potential energy yield of offshore wind farms. The European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5) was unsuccessful at identifying all observed LLJs at both buoy locations within the lowest 200 m. An extension of the National Renewable Energy Laboratory (NREL) 20-year wind resource dataset for the Outer Continental Shelf off the coast of California (CA20-Ext) yielded marginally greater captures of observed LLJs using the Mellor–Yamada–Nakanishi–Niino (MYNN) planetary boundary layer (PBL) scheme than the 2023 National Offshore Wind dataset (NOW-23), which uses the Yonsei University (YSU) scheme. However, CA20-Ext also produced the most LLJ false alarms, which are instances when a model identified an LLJ but no LLJ was observed. CA20-Ext and NOW-23 exhibited a tendency to overestimate the duration of LLJ events and underestimate LLJ core heights.

The vortex structure and near‐surface winds of Typhoon Faxai (2019) during landfall. Part II: Evaluation of WRF simulations

AbstractThis two‐part study presents a comprehensive analysis of (1) the vortex structure, (2) the inner core planetary boundary layer (PBL) wind profile, and (3) the overland surface winds of Typhoon Faxai(2019) during landfall observationally (Part I) and in high‐resolution Weather Research and Forecasting Model (WRF) simulations (Part II). Part II presents a framework for evaluating the wind field in the WRF simulations of Typhoon Faxai (2019) by comparing them with the observations presented in Part I. First, WRF simulations with two different surface roughness tables are presented. The simulation with the default table in WRF is shown to largely overestimate the surface wind speed. The overestimation is largely mitigated by using larger based on the climatology for urban and forest areas in Japan. Next, using the modified table, WRF simulations with three different PBL parametrizations, namely, Mellor–Yamada–Janjić, Yonsei University, and Mellor–Yamada–Nakanishi–Niino schemes are presented, and their impacts on the simulated vortex structure, inner core PBL wind profile, and overland surface winds are evaluated. In particular, the simulated surface winds are in good agreement with the observations outside the inner core of Typhoon Faxai. However, during the passage of the eyewall, the maximum surface winds are underestimated in the simulations. We show that the underestimation of the inertial stability, associated with the excessively large inner core of Typhoon Faxai in the simulations, led to the overestimation of the eyewall PBL jet height. This, in turn, caused the underestimation of the peak surface winds, as based on the empirical model of the wind reduction factor proposed in Part I.

Identifying Meteorological Drivers for Errors in Modeled Winds along the Northern California Coast

Abstract An accurate wind resource dataset is required for assessing the potential energy yield of floating offshore wind farms that are expected along the California outer continental shelf. The National Renewable Energy Laboratory has developed and disseminated an updated wind resource dataset offshore of California, using the Weather Research and Forecasting Model, referred to as the CA20 dataset. As compared to buoy lidar measurements that have become available recently, the CA20 dataset showed significant positive biases for 100-m wind speeds along Northern California wind energy lease areas. To investigate the meteorological drivers for the model errors, we first consider two 1-yr simulations run with two different planetary boundary layer (PBL) parameterizations: the Mellor–Yamada–Nakanishi–Niino (MYNN) PBL scheme (the chosen configuration in the CA20 dataset) and the Yonsei University PBL scheme (which significantly reduces the bias in modeled winds). By comparing the 1-yr simulations to the concurrent lidar buoy observations, we find that errors are larger with the MYNN PBL scheme in warm seasons. We then dive deeper into the analysis by running simulations for short-term (3-day) case studies to evaluate the sensitivity of initial/boundary condition forcings on model results. By analyzing the short-term simulations, we find that during synoptic-scale northerly flows driven by the North Pacific high and inland thermal low, a coastal warm bias in the MYNN simulation is mainly responsible for the modeled wind speed bias by altering the boundary layer thermodynamics. The results of our analysis will help guide the creation of an updated version of the CA20 dataset.

Fine scale simulation of sea salt aerosol under strong wind: WRF sensitivity test on PBL schemes and LES under Typhoon Hagibis

Identifying which factors among the planetary boundary layer (PBL) scheme, grid resolution, and source or sink terms dominantly influence the sea salt aerosol (SSA) spatial distribution is necessary. The dependence of SSA transport on PBL processes was investigated using a numerical weather prediction (NWP) model with a grid scale of ∼143 m. Size-segregated SSA transport during Typhoon Hagibis, which made landfall in Japan in 2019, was simulated using the Weather Research and Forecasting (WRF) model. The PBL schemes of Yonsei University (YSU), scale-aware YSU (Shin and Hong, 2015's scheme, SH) and Mellor-Yamada Nakanishi-Niino (MYNN), and large-eddy simulation (LES) subjected to turbulence transition enhancement via a cell perturbation method were tested to evaluate the reproducibility of observed values of wind speed, surface pressure, relative humidity, and temperature at meteorology masts. The results showed that wind speed was most influenced by the selection of the PBL scheme and LES. The distance from the coastline, as determined by the wind direction during the most intense typhoon stage, was correlated with a decrease in the time averaged SSA concentration when the storm moved inland. Despite the simulated SSA emission scheme's high sensitivity to wind speed (nearly 3rd power), the near-surface average concentration was less sensitive to the PBL schemes and LES selection than wind speed. However, the maximum arrival of the SSA differed by ∼20 % during the most vigorous stage of the typhoon. Budget analysis of the SSA transport equation highlighted that the advection term and gravity settling were most dependent on the PBL schemes and LES. The impacts of both factors on SSA transport were not limited to near-surface tendencies but were in a broader range over the atmosphere, up to 1–2 km. These results indicated that the vertical structure of the SSA field over the boundary layer depended on the PBL schemes and LES, thereby modulating the inland arrival of the SSA.

Anatomy of a Summertime Convective Event over the Arabian Region

Abstract This study investigates the structure and evolution of a summertime convective event that occurred on 14 July 2015 over the Arabian region. We use the WRF Model with 1-km horizontal grid spacing and test three PBL parameterizations: the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme; the Asymmetrical Convective Model, version 2, (ACM2) scheme; and the quasi-normal scale-elimination (QNSE) scheme. Convection initiates near the Al Hajar Mountains of northern Oman at around 1100 local time (LT; 0700 UTC) and propagates northwestward. A nonorographic convective band along the west coast of the United Arab Emirates (UAE) develops after 1500 LT as a result of the convergence of cold pools with the sea breeze from the Arabian Gulf. The model simulation employing the QNSE scheme simulates the convection initiation and propagation well. Although the MYNN and ACM2 simulations show convective initiation near the Al Hajar Mountains, they fail to simulate the development of the convective band along the UAE west coast. The MYNN run simulates colder near-surface temperatures and a weaker sea breeze, whereas the ACM2 run simulates a stronger sea breeze but a drier lower troposphere. Sensitivity simulations using horizontal grid spacings of 9 and 3 km show that lower-resolution runs develop broader convective structures and weaker cold pools and horizontal wind divergence, affecting the development of convection along the west coast of the UAE. The 1-km run using the QNSE PBL scheme realistically captures the sequence of events that leads to the moist convection over the UAE and adjacent mountains.

Determining the Impact of Boundary Layer Schemes on the Secondary Circulation of Typhoon Faxai Using Radar Observations in the Gray Zone

Abstract Using detailed radar observation data for Typhoon Faxai, which made landfall in the Tokyo metropolitan area in 2019, a sensitivity test of the boundary layer (BL) schemes for a numerical weather prediction (NWP) model was conducted for gray-zone numerical simulations with a grid spacing of 250 m. We compared the results of our simulations using an NWP model with radar observations that captured the BL and the secondary circulation structures of Faxai. We used three BL schemes based on a Reynolds-averaged model, the gray-zone model, and a large-eddy simulation (LES) model: the Mellor–Yamada–Nakanishi–Niino level 3 (MYNN3) scheme, the Anisotropic Deardorff Model (ADM) scheme, and the Deardorff (DDF) scheme, respectively. The turbulence kinetic energy was the smallest, and the inflow near Earth’s surface the strongest, in the gray-zone simulation with the DDF scheme. This simulation also produced values for BL thickness and secondary circulation that were the closest to observation and reproduced horizontal roll structures whose scale was larger than the observation. Neither experiment using the MYNN3 scheme or the ADM scheme produced rolls, but the parameterized turbulence seemed to estimate the effects of the rolls. However, their BL heights were higher than observed, suggesting that the MYNN3 and ADM schemes are not appropriate for a 250 m grid simulation of the present case. These results are also confirmed against LES with 50 m grid spacing in which the DDF scheme is used. In summary, this study provides insights into the interpretation of the properties of BL schemes in the gray zone.

Simulating Meteorological and Water Wave Characteristics of Cyclone Shaheen

The Bay of Bengal and Arabian Sea are annually exposed to severe tropical cyclones, which impose massive infrastructure damages and cause the loss of life in coastal regions. Cyclone Shaheen originally generated in the Bay of Bengal in 2021 and translated a rare east-to-west path toward the Arabian Sea. Although the cyclone’s wind field can be obtained from reanalysis datasets such as ERA5 (fifth generation European Centre for Medium-Range Weather Forecasts), the wind speed cannot be reproduced with realistic details in the regions close to the center of cyclone due to spatial resolution. In this study, to address this problem, the high-resolution advanced Weather Research and Forecasting (WRF) model is used for simulation of Shaheen’s wind field. As a critical part of the study, the sensitivity of the results to the planetary boundary layer (PBL) parameterization in terms of the track, intensity, strength and structure of the cyclone Shaheen is investigated. Five experiments are considered with five PBL schemes: Yonsei University (YSU); Mellor–Yamada–Janjić (MYJ); Mellor–Yamada–Nakanishi–Niino level 2.5 (MYNN); Asymmetric Convective Model version 2 (ACM2); Quasi-Normal Scale Elimination (QNSE). The track, intensity, and strength of the experiments are compared with the wind fields obtained from the Joint Typhoon Warning Centre (JTWC) dataset. The results imply the high dependency of the track, intensity, and strength of the cyclone to the PBL parameterization. Simulated tracks with non-local PBL schemes (YSU and ACM2) outperformed those of the local PBL schemes (MYJ, MYNN, and QNSE), especially during the rapid intensification phase of Shaheen before landfall. The YSU produced highly intensified storm, while the ACM2 results are in better agreement with the JTWC data. The most accurate track was obtained from the ERA5 data; however, this dataset overestimated the spatial size and underestimated the wind speed. The WRF model using either YSU or ACM2 overestimated the wind speed compared to that of the altimeter data. The YSU and ACM2 schemes were able to reproduce the observed increase in wind speed and pressure drop at in situ stations. The wind data from EAR5 and cyclone parametric model were applied to the SWAN model to simulate the wave regime in the Arabian Sea during the time that Shaheen was translating across the region. Janssen formulation for wind input and whitecapping dissipation source terms in combination with both ERA5 and hybrid wind were used and the minimum combined error in the prediction of significant wave height (Hs) and zero up-crossing wave period (Tz) was examined. The maximum significant wave height for hybrid wind is higher than that of ERA5, while the cyclone development was successfully inferred from the wave field of the hybrid data.

Evaluation of Model Summertime Boundary Layer Cloud Development over Complex Terrain in New York State

Abstract Weather forecasting over complex terrain with diverse land cover is challenging. Utilizing the high-resolution observations from New York State Mesonet (NYSM), we are able to evaluate the surface processes of the Weather Research Forecast (WRF) Model in a detailed, scale-dependent manner. In the study, possible impacts of land–atmosphere interaction on surface meteorology and boundary layer cloud development are investigated with different model resolutions, land surface models (LSMs), and planetary boundary layer (PBL) physical parameterizations. The High-Resolution Rapid Refresh, version 3 (HRRR), forecasting model is used as a reference for the sensitivity evaluation. Results show that over complex terrain, the high-resolution simulations (1 km × 60 vertical levels) generally perform better compared to low-resolution (3 km × 50 levels) in both surface meteorology and cloud fields. LSMs play a more important role in surface meteorology compared to PBL schemes. The NoahMP land surface model exhibits daytime warmer and drier biases compared to the Rapid Update Cycle (RUC) due to better prediction of the Bowen ratio in RUC. The PBL schemes would affect the convective strength in the boundary layer. The Shin–Hong (SH) scale-aware scheme tends to produce the strongest convective strength in the PBL, while the ACM2 PBL scheme rarely resolved convection even at 1-km resolution. By considering the radiation effect of subgrid-scale (SGS) clouds, the Mellor–Yamada–Nakanishi–Niino eddy diffusivity mass flux (MYNN-EDMF) predicted the highest cloud coverage and lowest surface solar radiation bias. The configuration of SGS clouds in MYNN-EDMF would not only significantly reduce shortwave radiation bias, but also affect the convection behaviors through land surface–cloud–radiation interaction.

How Do Planetary Boundary Layer Schemes Perform in Hurricane Conditions: A Comparison With Large‐Eddy Simulations

AbstractParameterizations of turbulent processes in planetary boundary layer (PBL) schemes impact tropical cyclone (TC) forecasts. Existing PBL schemes are mostly designed for low‐wind conditions, and assessing their uncertainties in hurricane conditions remains challenging, mostly due to very scarce observations. Using a recently developed framework based on large‐eddy simulations (LES), this study evaluates K‐profile parameterization (KPP) and high‐order PBL schemes in hurricane conditions. Among KPP PBL schemes, the Global Forecast System (GFS) scheme tends to produce excessively deep inflow layers with large values of eddy viscosity (Km). Opposite results are found for the Yonsei University (YSU) scheme. Using LES results as a benchmark, the performance of YSU and GFS schemes is improved by modifying the “shape parameter” such that Km is maximized closer to the surface, and by using a new definition of boundary layer height tailored to high‐wind conditions. The LES results also suggest an asymptotic mixing length of ∼40 m can improve the Louis‐type parameterizations of the YSU scheme that operates above the boundary layer. Among high‐order PBL schemes, the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme produces reasonably accurate vertical profiles of eddy viscosity, turbulent stress, and boundary layer winds under different high‐wind conditions. Further analysis of MYNN supports a “three‐layer” strategy for the mixing length parameterization for TCs that represents different types of turbulent regimes. In contrast, the high‐order eddy‐diffusivity mass‐flux scheme produces excessive boundary‐layer vertical mixing and a deeper inflow layer, partly attributable to a notable overestimation of the maximum allowable mixing length in the PBL code.

Sensitivity analysis of mesoscale simulations to physics parameterizations over the Belgian North Sea using Weather Research and Forecasting – Advanced Research WRF (WRF-ARW)

Abstract. The Weather Research and Forecasting (WRF) model offers a multitude of physics parameterizations to study and analyze the different atmospheric processes and dynamics that are observed in the Earth's atmosphere. However, the suitability of a WRF model setup is known to be highly sensitive to the type of weather phenomena and the type and combination of physics parameterizations. A multi-event sensitivity analysis is conducted to identify general trends and suitable WRF physics setups for three extreme weather events identified to be potentially harmful for the operation and maintenance of wind farms located in the Belgian offshore concession zone. The events considered are Storm Ciara on 10 February 2020, a low-pressure system on 24 December 2020, and a trough passage on 27 June 2020. A total of 12 WRF simulations per event are performed to study the effect of the update interval of lateral boundary conditions and different combinations of physics parameterizations (planetary boundary layer, PBL; cumulus; and microphysics). Specifically, the update interval of ERA5 lateral boundary conditions is varied between hourly and 3-hourly. Physics parameterizations are varied between three PBL schemes (Mellor–Yamada–Nakanishi–Niino, MYNN; scale-aware Shin-Hong; and scale-aware Zhang), four cumulus schemes (Kain–Fritsch, Grell–Dévényi, scale-aware Grell–Freitas, and multi-scale Kain–Fritsch), and three microphysics schemes (WRF Single-Moment five-class scheme, WSM5; Thompson; and Morrison). The simulated wind direction and wind speed are compared qualitatively and quantitatively to operational supervisory control and data acquisition (SCADA) data. Overall, a definitive best-case setup common to all three events is not identified in this study. For wind direction and wind speed, the best-case setups are identified to employ scale-aware PBL schemes. These are most often driven by hourly update intervals of lateral boundary conditions as opposed to 3-hourly update intervals, although it is only in the case of Storm Ciara that significant differences are observed. Scale-aware cumulus schemes are identified to produce better results when combined with scale-aware PBL schemes, specifically for Storm Ciara and the trough passage cases. However, for the low-pressure-system case this trend is not observed. No clear trend in utilizing higher-order microphysics parameterization considering the combinations of WRF setups in this study is found in all cases. Overall, the combination of PBL, cumulus, and microphysics schemes is found to be highly sensitive to the type of extreme weather event. Qualitatively, precipitation fields are found to be highly sensitive to model setup and the type of weather phenomena.

Performance Evaluation of Planetary Boundary Layer Schemes in Simulating Structures of Wintertime Lower Troposphere in Seoul Using One-Hour Interval Radiosonde Observation

We investigated the structures of the wintertime lower troposphere in Seoul, South Korea on 17 and 18 January 2017 by performing 1 h interval radiosonde observation and numerical simulations. In the daytime on 17 January, the height of the convective boundary layer (CBL) sharply and quickly increased when the residual layer became a part of the CBL. From the afternoon on 17 January, moist air with clouds began to substantially intrude in the lower troposphere in Seoul, and radiative heating/cooling weakened. As a result, the mixing of air in the lower troposphere was inhibited and the vertical gradients of potential temperature and water vapor mixing ratio changed little on 18 January. We evaluated the performance of four planetary boundary layer (PBL) parameterization schemes (the Yonsei University (YSU), Mellor–Yamada–Janjić (MYJ), Mellor–Yamada–Nakanishi–Niino (MYNN), and Asymmetric Convective Model version 2 (ACM2) schemes) coupled with the Weather Research and Forecasting model in simulating the structures of the lower troposphere against 1 h interval radiosonde observation. The general tendencies of the air temperature and wind speed in the lower troposphere were well-reproduced in the four simulations. However, the sharp increase in the CBL height did not appear in the four simulations, implying that the process of the residual layer becoming a part of the CBL in the daytime is not well-parameterized. Additionally, the simulated water vapor mixing ratio near the surface was smaller compared with the observation. We found that small-scale turbulence in the CBL, which mixes advected air and pre-existing air, was not reproduced well by the PBL parameterization schemes. Compared with the other simulations, the most accurate air temperature and wind speed were reproduced in the simulation with the MYJ scheme, while the CBL development and moisture advection were reproduced relatively well in the simulation with the MYNN scheme.

Sensitivity of the Tropical Cyclone Boundary Layer to Vertical Diffusion in a Turbulent Kinetic Energy‐Based Boundary Layer Parameterization Scheme at Gray‐Zone Resolution

AbstractThe impact of eddy diffusivities in the Mellor–Yamada–Nakanishi–Niino (MYNN) boundary layer parameterization on the fine‐scale turbulence and intensity of tropical cyclone (TC) simulations at gray‐zone resolution was investigated through comparison with large‐eddy simulation (LES) and observations. A series of idealized experiments were performed with varying vertical mixing lengths using the Advanced Research version of the Weather Research and Forecasting model. For enhanced eddy diffusivity cases, the storm intensity was sensitive to the eddy diffusivities and reducing the overestimation of the eddy diffusivities caused a significant increase in intensity. For the reduced eddy diffusivity cases, the TC intensity was marginally influenced by the eddy diffusivities. However, the reduction in eddy diffusivities effectively reproduced the fine‐scale turbulence. The stronger the resolved fine‐scale turbulence that could be established, the more frictional tendencies near the location of maximum wind were provided to the TC, the closer it was to the reference LES and observations. The resolved turbulence also had consequences in terms of preventing the wind gradient from increasing, which was one of the primary reasons for the discrepancies in the nonlinear intensity and fine‐scale features. The results of this study demonstrate that the modification of eddy diffusivities in the MYNN scheme can measure or guide improvements in the simulation of TC boundary layer fine‐scale features and intensity when operating at the gray‐zone resolution.

The Skill Assessment of Weather and Research Forecasting and WAVEWATCH-III Models During Recent Meteotsunami Event in the Persian Gulf

This study aims to use a fully realistic high-resolution mesoscale atmospheric and wave model to reproduce met-ocean conditions during a meteotsunami in the Persian Gulf. The atmospheric simulations were performed with the Weather and Research Forecasting (WRF) model by varying planetary boundary layer, microphysics, cumulus, and radiations parameterizations. The atmospheric results were compared to the meteorological observations (e.g., air pressure and wind speed) from the coastal and island synoptic and buoy stations of the nearest area to the meteotsunami event. The results show that using Mellor-Yamada-Nakanishi-Niino (MYNN) scheme for planetary boundary and surface layer had the best performance for stations over the water, whereas applying Mellor-Yamada-Janjic scheme for planetary boundary and Eta similarity surface layer had the best performance for stations over the land. For wave simulations, the WAVEWATCH-III model was employed with the well-known WAM-Cycle4 formulation and a more recent ST6 package. Six WRF experiments and ERA5 wind data were used to force the wave models. The new error parameter was introduced to identify the optimum wind data for wave simulation. EXP4 configuration which uses the MYNN scheme for planetary boundary and surface layer was led to minimum error, while ERA5 severely underestimated Hs and Tp parameters. For the first time, the Gaussian Quadrature Method (GQM) was implemented in the WAVEWATCH-III model and combined with a depth scale to be used in the Persian Gulf. This method is more accurate for non-linear wave-wave interaction than the default Discrete Interaction Approximation (DIA) method. Lower coefficients for dissipation term were required for GQM and the resulted bulk wave parameters were improved compared to the DIA method. The calibrated ST6 formulation with GQM resulted in a more realistic prediction of wave spectrum than the default settings of the WAVEWATCH-III.

Evaluating the Impact of Planetary Boundary Layer, Land Surface Model, and Microphysics Parameterization Schemes on Simulated GOES-16 Water Vapor Brightness Temperatures

The impact of several land surface models (LSMs) and microphysics (MP), planetary boundary layer (PBL), and surface layer schemes on the accuracy of simulated brightness temperatures (BTs) from water vapor (WV) sensitive bands was examined via comparison with observations from the GOES-16 Advanced Baseline Imager. Nine parameterization configurations were evaluated. Analysis revealed that, compared to the Thompson MP scheme, the National Severe Storms Laboratory MP scheme produced lower simulated WV BTs in the upper troposphere but higher WV BTs in the middle and lower troposphere. The configuration with the Geophysical Fluid Dynamics Laboratory MP and hybrid eddy-diffusivity mass-flux (EDMF) PBL instead of Mellor–Yamada–Nakanishi–Niino (MYNN) PBL produced higher BTs. Yet, changing the PBL from MYNN to Shin–Hong or EDMF reduced the simulated WV BTs. Changing the LSM from Noah to RUC also resulted in lower simulated WV BTs, which were further enhanced with the MYNN surface layer instead of the GFS. The location and orientation of upper-level jet streams and troughs was assessed using the location of WV gradient objects. Every configuration had an increased translation speed compared to the observations, as forecast WV gradient objects were west of the observation objects early in the forecast and then east later in the forecast.

Influence of Planetary Boundary Layer (PBL) Parameterizations in the Weather Research and Forecasting (WRF) Model on the Retrieval of Surface Meteorological Variables over the Kenyan Highlands

Regional climate models (RCMs) are crucial for climate studies and may be an alternative source of meteorological data in data-scarce regions. However, the effectiveness of the numerical weather prediction (NWP) models applied in RCMs is hampered by the parameterization of unresolved physical processes in the model. A major source of uncertainties in NWP models is the parameterization of the planetary boundary layer (PBL). This study evaluates the influence of seven PBL parameterization schemes in the Weather Research and Forecasting (WRF) model on the retrieval of four meteorological variables over the Kenyan highlands. The seven PBL schemes consist of four local schemes: the Mellor-Yamada-Janjic (MYJ), Mellor-Yamada-Nakanishi-Niino (MYNN), Bougeault-Lacarrere (BouLac), quasinormal scale elimination (QNSE), and three nonlocal schemes: asymmetrical convective model version 2 (ACM2), Shin and Hong (SHIN) and Yonsei University (YSU). The forcing data for the WRF model was obtained from the fifth generation of the European ReAnalysis (ERA5) dataset. The results were validated against observational data from the Trans-African Hydro-Meteorological Observatory (TAHMO). WRF was found to simulate surface meteorological variables with spatial details coherent with the complex topography within the Kenyan highlands, irrespective of the PBL scheme. A comparison between 2-meter temperature (T2) derived from the YSU scheme and T2 from the land component of ERA5 (ERA5-Land) indicates that surface meteorological variables derived from WRF are better suited for applications over the Kenyan highlands. The choice of the PBL scheme was found to primarily influence the simulation of the 10-meter wind speed (WS10) and rainfall as opposed to T2 and the 2-meter relative humidity (RH2). The insensitivity of the 2-meter variables to the choice of the PBL scheme is attributed to the influence of the surface layer parameterization near the surface. Results from the rainfall simulation indicate that the YSU scheme provides a more realistic depiction of PBL dynamics within the study area. Hence, the YSU scheme is best suited for simulating surface meteorological variables over the Kenyan highlands.

Stability Functions in the Stable Surface Layer Derived from the Mellor–Yamada–Nakanishi–Niino (MYNN) Scheme

Stability Functions in the Stable Surface Layer Derived from the Mellor–Yamada–Nakanishi–Niino (MYNN) Scheme

Combinational Optimization of the WRF Physical Parameterization Schemes to Improve Numerical Sea Breeze Prediction Using Micro-Genetic Algorithm

This study aims to improve the performance of the Weather Research and Forecasting (WRF) model in the sea breeze circulation using the micro-Genetic Algorithm (micro-GA). We found the optimal combination of four physical parameterization schemes related to the sea breeze system, including planetary boundary layer (PBL), land surface, shortwave radiation, and longwave radiation, in the WRF model coupled with the micro-GA (WRF-μGA system). The optimization was performed with respect to surface meteorological variables (2 m temperature, 2 m relative humidity, 10 m wind speed and direction) and a vertical wind profile (wind speed and direction), simultaneously for three sea breeze cases over the northeastern coast of South Korea. The optimized set of parameterization schemes out of the WRF-μGA system includes the Mellor–Yamada–Nakanishi–Niino level-2.5 (MYNN2) for PBL, the Noah land surface model with multiple parameterization options (Noah-MP) for land surface, and the Rapid Radiative Transfer Model for GCMs (RRTMG) for both shortwave and longwave radiation. The optimized set compared with the various other sets of parameterization schemes for the sea breeze circulations showed up to 29 % for the improvement ratio in terms of the normalized RMSE considering all meteorological variables.

The impacts of multi-physics parameterization on forecasting heavy rainfall induced by weak landfalling Typhoon Rumbia (2018)

The weak and sustained tropical cyclone (TC) Rumbia (2018) and its associated heavy precipitation after landfall are simulated by the Weather Research and Forecasting model using three microphysics parameterizations (MP) and two planetary boundary layer (PBL) schemes. Key factors regulating the heavy rainfall distribution of landed TC at three stages (landfall, inland slow-moving and recurving stage) and their sensitivities to multi-physics parameterizations are examined. Results show that heavy rainfall distribution is largely regulated by the intensity of TC itself at landfall, but more affected by environmental factors at the two inland stages, including the environmental vertical wind shear, the associated upper-level wind divergence and low-level moisture convergence. Different MP and PBL schemes mainly affect the simulation of TC thermodynamic structure and key environmental factors, leading to their different forecast skills at three stages. Specifically, the Yonsei University (YSU) PBL scheme systematically outperforms the Mellor-Yamada-Nakanishi-Niino (MYNN) scheme in simulating stronger TC and heavy rainfall intensity. With the advantageous YSU PBL scheme, the Ferrier MP scheme only shows a slight advantage in simulating the intensity of TC and heavy rainfall at landfall. However, the WRF single-moment 6-class (WSM6) and Thompson scheme produce more graupel or snow, and better simulate the key environmental factors, showing their respective advantages in simulating the heavy rainfall structure and location at the two inland stages. This implies that the use of nonlocal YSU PBL scheme and the MP schemes with sophisticated ice processes shows superiority in simulating the postlandfall heavy rainfall induced by weak TCs.

Sensitivity analysis of planetary boundary layer schemes using the WRF model in Northern Colombia during 2016 dry season

This study conducted meteorological simulations in northern Colombia by analyzing different planetary boundary layer (PBL) schemes available in the numerical Weather Research and Forecasting (WRF) model. The study area included three nested domains with horizontal resolutions of 18 km, 6 km, and 2 km, with 38 vertical levels. The evolution and structure of the PBL were analyzed during the driest months (March, April, and May 2016) and in regions with the highest particulate matter concentrations. Sensitivity analysis of the WRF model was performed with two local and two non-local PBL schemes. The results were validated using observations of the surface air temperature, relative humidity, and surface wind speed collected from three meteorological stations in the area. The PBL heights were experimentally determined using radiosonde data provided by a station located in the center of the study area. Variations in PBL heights were estimated using linear regression methods and minimization of statistical errors for the bulk Richardson number, as well as analysis of vertical temperature and wind profiles. The WRF model reliably reproduced the daily values and diurnal cycles of temperature, relative humidity, and wind speed within the PBL and accounted for the influence of topography and sea breezes. Horizontal heat advection dominates the upwelling of air masses when sea breezes are active. The onshore wind direction starts to change from east to northwest, implying a decay in the land breeze regime. All schemes overestimate the mixing height and tend to underestimate surface air temperature values at night. All show wetter conditions and underestimate wind speed. Although the non-local Yonsei University (YSU) scheme shows the best performance, it also shows the largest sources of errors when determining the behavior of the surface layer during stable conditions. Relative humidity and wind speed estimates provided by the local Mellor‐Yamada‐Nakanishi‐Niino (MYNN) scheme were closer to those recorded at the meteorological stations.