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.