Nonlocal Planetary Boundary Layer Research Articles

Evaluation of PBL Parameterization Schemes in WRF Model Predictions during the Dry Season of the Central Amazon Basin

Planetary Boundary Layer (PBL) parameterization schemes are employed to handle subgrid-scale processes on atmospheric models, playing a key role in accurately representing the atmosphere. Recent studies have shown that PBL schemes are particularly fundamental to the depiction of PBL height (PBLH), especially over the Amazon. In the present study, we investigated the performance of PBL schemes on the representation of meteorological variables, turbulent fluxes, PBL vertical structures, and PBLH over the central Amazon basin under dry conditions, taking advantage of observations from the Observations and Modeling of the Green Ocean Amazon campaign (GoAmazon2014/5) for validation and evaluation. Numerical experiments were carried out within the WRF model using eight PBL schemes for two dry periods from 2014 (typical year) and 2015 (El-Niño year), and results from the 1-km resolution domain were directly compared to hourly in situ observations. In general, all PBL schemes present good performance to reproduce meteorological variables, with nonlocal (local) PBL schemes producing better performance in the 2014 (2015) study period. All PBL schemes in general overestimate (>100%) daytime turbulent fluxes. Thermodynamic (daytime) vertical structures are better predicted than mechanical (nocturnal) ones. The local MYNN2.5 scheme showed the overall best performance for PBLH prediction, mainly at night.

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.

Exploring the Impact of Planetary Boundary Layer Schemes on Rainfall Forecasts for Typhoon Mujigae, 2015

Sensitivity experiments were conducted on Typhoon Mujigae, which occurred in 2015, wherein the Weather Research and Forecasting Advanced Research (WRF-ARW) model was used to select two local and two nonlocal planetary boundary layer (PBL) parameterization schemes: the quasi-normal scale elimination (QNSE) and Mellor–Yamada–Janjic (MYJ) schemes, and the Yonsei University (YSU) and medium-range forecast (MRF) schemes, respectively. The differences in rainfall response in the typhoon’s inner core and outer region were evaluated by comparing the anomaly rainfall distribution, heat transmission, and mixing processes in the boundary layer among the PBL schemes. The results show that the simulated rainfall in typhoon Mujigae has large uncertainty among the PBL schemes and a significant difference between the inner and outer regions. Compared with the observation, the simulated rainfall was significantly higher in the inner core and slightly lower in the outer region. All PBL schemes accurately identified the rainfall location, although the amounts differed between the schemes. The rainfall levels in the MRF scheme were closest to the observation, followed by those in the YSU and MYJ schemes; the QNSE scheme showed the largest deviation. In general, rainfall simulation using a nonlocal boundary layer scheme such as MRF had the best results for both the inner core and the outer region.

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.

Evaluation of PBL parameterization schemes against direct observations during a land depression over Central India

Detailed observations from the Cloud Aerosol Interaction and Precipitation Enhancement Experiment (CAIPEEX) Phase III integrated observational campaign 2014 are used to validate six different planetary boundary layer (PBL) schemes in Weather Research and Forecasting (WRF) model during a land depression over central India in summer monsoon season. This work is a continuation of our previous study on a dry event associated with a heatwave over the same region (Sathyanadh et al. 2017). Four local and two nonlocal PBL schemes are evaluated with hourly observations of air temperature, eddy covariance fluxes, friction velocity, wind, and radiation components for the model performance of the complete energy budget and weather parameters. The models simulated a slightly warmer and drier surface layer with strong winds and overestimated sensible heat flux and ground heat flux during the wet event. The errors in the surface layer temperature are closely associated with errors in the sensible heat flux and incoming radiation, while the errors associated with moisture are linked to the excessive surface layer mixing. The interrelationships between various surface layer parameters and PBL heights are investigated to identify the key controlling factors of PBL evolution in the context of land depression.

Evaluating planetary boundary-layer schemes and large-eddy simulations with measurements from a 250-m meteorological mast

We evaluate the performance of six planetary boundary layer (PBL) schemes and a large-eddy simulation (LES) for the characterization of the neutral PBL through idealized simulations using the Weather Research and Forecasting model. The evaluation is performed by comparison with sonic anemometer measurements from a 250-m tall meteorological tower that observes close-to-homogeneous winds under the predominant westerlies. All simulations show similar behavior for the vertical temperature profile except for one that uses a non-local PBL scheme, which is known to produce excessive vertical mixing. Within the range of measurements of the tower and except for a PBL scheme that produces a too low jet, the simulations using PBL schemes, which were developed to simulate the nighttime atmosphere, generally show the largest deviations from the observed wind speeds. As expected within the surface layer, the LES shows excessive vertical shear. More importantly, we find that two PBL schemes that use their own surface-layer scheme are effectively reducing the surface roughness. When looking at the vertical profile of momentum exchange coefficient, we find very good agreement between a local and a non-local PBL scheme within the bulk of the PBL, the highest values for the simulation with the PBL scheme showing excessive vertical mixing and the lowest values for that most suitable for nighttime conditions. The comparison with the observed turbulent kinetic energy reveals a good match between the LES and a simulation using a local PBL scheme and a general underestimation (overestimation) of turbulence within the range of measurements by the PBL schemes (LES).

Addition of Multilayer Urban Canopy Models to a Nonlocal Planetary Boundary Layer Parameterization and Evaluation Using Ideal and Real Cases

AbstractThe multilayer urban canopy models (UCMs) building effect parameterization (BEP) and BEP + building energy model (BEM; a building energy model integrated in BEP) are added to the Yonsei University (YSU) planetary boundary layer (PBL) parameterization in the Weather Research and Forecasting (WRF) Model. The additions allow for the first analysis of the detailed effects of buildings on the urban boundary layer in a nonlocal closure scheme. The modified YSU PBL parameterization is compared with the other 1.5-order local PBL parameterizations that predict turbulent kinetic energy (TKE), Mellor–Yamada–Janjić and Bougeault–Lacarerre, using both ideal and real cases. The ideal-case evaluation confirms that BEP and BEP+BEM produce the expected results in the YSU PBL parameterization because the simulations are qualitatively similar to the TKE-based PBL parameterizations in which the multilayer UCMs have long existed. The modified YSU PBL parameterization is further evaluated for a real case. Similar to the ideal case, there are larger differences among the different UCMs (simple bulk scheme, BEP, and BEP+BEM) than across the PBL parameterizations when the UCM is held fixed. Based on evaluation against urban near-surface wind and temperature observations for this case, the BEP and BEP+BEM simulations are superior to the simple bulk scheme for each PBL parameterization.

Effect of Scale-Aware Nonlocal Planetary Boundary Layer Scheme on Lake-Effect Precipitation at Gray-Zone Resolutions

Abstract The effects of a nonlocal planetary boundary layer (PBL) scheme that considers scale dependency in the parameterized turbulent vertical transport are investigated for a case of wintertime lake-effect precipitation over Korea at gray-zone resolutions using a mesoscale model. An experiment using the scale-aware PBL scheme is compared with that using a conventional PBL scheme, which shows that the simulated precipitation amount at a resolution of less than 1 km is smaller with the scale-aware PBL scheme. The role of turbulent processes in simulating lake-effect precipitation is understood through interaction with microphysical processes. When the scale-aware PBL scheme is used, liquid water content is increased while ice water content is reduced. The higher cloud water content is because of enhanced condensation with stronger updrafts, attributed to the suppression of parameterized turbulent mixing. This results in higher rainwater content by enhancing autoconversion and accretion from cloud water to rainwater. The cloud ice content is reduced mainly because of the suppressed deposition and enhanced sublimation centered near the PBL top, and the snow content is reduced mainly because of the enhanced sublimation below and near the PBL top and suppressed growth of cloud ice to snow. The lower ice water content is mainly due to the drier PBL, attributed to the enhanced resolved (suppressed parameterized) turbulent moisture transport and enhanced condensation. The melting of a smaller amount of snow under dominant cold rain processes is responsible for the reduced surface precipitation with the scale-aware PBL scheme.

Extreme rainfall event in the Northeast coast of Brazil: a numerical sensitivity study

This study investigates an extreme rainfall event which occurred in Northern Brazil (NEB) between 20 and 30th 2017 May causing several deaths and making thousands homeless. Based on a suite of microphysics and planetary boundary layer (PBL) schemes based on the WRF model, it is demonstrated that anomalous weather conditions are characterized by significant upward and eastward wind flow. Omega differences with respect to climatological conditions showed values up to − 0.04 Pa/s and wind up to 6–8 m/s in consonance with higher precipitation in May in the NEB coast. The cumulative rainfall for 11 days was higher than 500 mm in some locations, as measured by weather stations. These conditions were simulated by the WRF model under different physics parameterization schemes. In total 24 experiments with WRF were implemented. Non-local PBL demonstrated better performance than the local PBL. Moreover, the rainfall was concentrated in small portions of the region, and the local scheme limited the WRF conditions to estimate the correct maximum precipitation location. The Morrison scheme performed better compared to the other schemes. Results presented here show that the correct choices of the microphysics and PBL parameterizations are fundamental to obtain good simulation/forecast, especially for extreme rainfall events. This study demonstrates that regional modeling is crucial to provide accurate information to forecasters and decision makers to plan actions which hamper catastrophic situations such as landslides and floods in high-risk regions.

A Numerical Study of a Sea Breeze at Fuerteventura Island, Canary Islands, Spain

The genesis and evolution of a sea-breeze front associated with a shallow convection episode at Fuerteventura Island (Canaries archipelago) are investigated using the Weather Research and Forecasting (WRF) numerical model. Three local and two non-local planetary boundary layer (PBL) parametrization schemes are used, and results of the numerical simulations are compared with observations from the local meteorological surface station network. Statistical analysis shows a good agreement between simulations and measurements, in particular for the 2-m temperature. These results are used in choosing the final PBL scheme and its coupled surface-layer scheme, which is used to simulate the daytime period of the sea-breeze development and weakening. During this episode, a sea-breeze front moved inland against a prevailing offshore flow 3 h after sunrise, colliding with a downslope flow near the centre of the island. The convergence of both flows from opposing directions produced a strong updraft, forming shallow convection a few hours later.

Evaluation of WRF planetary boundary layer parameterization schemes for simulation of monsoon depressions over India

This study evaluates the fidelity of five planetary boundary layer (PBL) parameterization schemes in the advanced weather research and forecasting model for simulating monsoon depressions (MDs) over India. Five PBL schemes include; nonlocal first-order medium-range forecasting (MRF) and Yonsei University (YSU); hybrid first-order Asymmetric Convective Model version 2 (ACM2), and local one-and-a-half-order Bougeault–Lacarrere (BouLac) and Mellor–Yamada–Nakanishi–Niino (MYNN2). PBL schemes show significant impact on rainfall along with dynamical and thermodynamical parameters associated with MDs at the surface as well as at the upper levels. MRF simulates a relatively shallower, warmer and drier boundary layer compared to others. Results reveal that strong upper-level divergence and high moisture content within the lower levels are favorable for the occurrence of heavy rain associated with MDs. However, stronger wind shear within the mid-troposphere weakens the system and reduces the rain intensity. Based on the results and keeping the rainfall product in view, it is found that nonlocal PBL schemes (MRF and ACM2) have better forecast skills score than local PBL schemes (BouLac and MYNN2) over the Indian region.

Does convection-permitting simulate better rainfall distribution and extreme over Guinean coast and surroundings?

This research work focuses on the problem of climate simulation of rainfall over West Africa and particularly over coastal countries of the Gulf of Guinea by Regional Climate Models. The sensitivities of Weather Research and Forecasting (WRF) Model were tested for changes in horizontal resolution (convection permitting versus parameterized) on the replication of West African monsoon for the year 2014. Sensitivity test was also performed for response of rainfall to changes in microphysics (MP) and planetary boundary layer (PBL) schemes. Generally, the result shows that WRF is able to replicates rainfall distribution with an adequate representation of the dynamical features of West African monsoon system. The high-resolution (wrf-4km) shows dry bias along the coast of the Gulf of Guinea but generally outperforms wrf-24km run especially in replication of the extreme rainfall distribution. The dry bias along the coastal area is suggested to not to be only related to convection but mostly to microphysics and PBL parameterisation schemes. Differences were noticed between the dynamics of WRF and ERA-interim outputs despite the use of spectral nudging in the experiment which then suggest strong interactions between scales. These differences were observed to be restricted mainly to the low-layer monsoon flow in JJA. Both runs at 24 and 4 km hardly simulate the typical diurnal distribution of rainfall. The sensitivity of WRF to MP (only sophisticated MP were tested) and PBL reveals a stronger impact of PBL than MP on rainfall distribution and the most significant added value over the Guinean coast and surroundings area was provided by the configurations using non-local PBL scheme (as ACM2). The changes in MP and PBL schemes in general seem to have less effect on the explicit runs (wrf-4km) in the replication of the rainfall over the Gulf of Guinea and the surroundings seaboard.

WRF simulations for a heavy rainfall event in southern Italy: Verification and sensitivity tests

In this paper we study the sensitivity of the quantitative precipitation forecast (QPF) to different physical parameterization schemes and configuration settings of WRF (Weather Research and Forecasting) model for a catastrophic event occurred in Calabria, Southern Italy.The event lasted four days, from 30 October 2015 to 2 November 2015, and affected mostly the eastern side of the peninsula. Several raingauges measured >500 mm for the whole event, with the maximum recorded in Chiaravalle (739 mm for the whole event, which is about half of the yearly rainfall). Damages to railway, roads and properties were extensive.To show the sensitivity of the QPF to the model settings and physical parameterization schemes, twelve configurations of WRF are considered, in addition to the control run (CC), which is the configuration used for the operational forecast at ISAC-CNR in Lamezia Terme. In particular, we tested the sensitivity of the QPF to four different microphysical schemes, three Planetary Boundary Layer (PBL) schemes (non-local and local), two cumulus parameterization schemes, two different horizontal resolutions (1 km vs 3 km), and different initial and boundary conditions. The QPF is verified both by quantitative measures and qualitative scores, using the observations of 156 raingauges quite uniformly distributed over Calabria. The analysis refers to 12 h and 24 h rainfall, but some details are given also for shorter time scales.Results show that most of the WRF configurations are able to predict most of the features of the precipitation field. However, there are differences at the local scale (few tens of kilometers) that are important for several purposes, including civil protection. In particular, using a non-local PBL closure scheme was fundamental to correctly simulate the heavy precipitation close to the sea for this case, while other configurations predicted the rainfall more anchored to mountain peaks. The good performance of non-local PBL scheme confirms the findings of a previous study over Calabria.The results of this paper show the sensitivity of the QPF at the local scale for different WRF configurations and an ensemble approach is envisaged to improve the prediction of such catastrophic events.

Performance of the WRF Model for Surface Wind Prediction around Qatar

AbstractThe performance of the Weather Research and Forecasting (WRF) Model is examined for the region around Qatar in the context of surface winds. The wind fields around this peninsula can be complicated owing to its small size, to a complex pattern of land and sea breezes influenced by the prevailing shamal winds, and to its dry and arid nature. Modeled winds are verified with data from 19 land stations and two offshore buoys. A comparison with these data shows that nonlocal planetary boundary layer (PBL) schemes generally perform better than local schemes over land stations during the daytime, when convective conditions prevail; at nighttime, over land and over water, both schemes yield similar results. Among other parameters, modifications to standard USGS land-use descriptors were necessary to reduce model errors. The RMSE values are comparable to those reported elsewhere. Simulated winds, when used with a wave model, result in wave heights comparable to buoy measurements. Furthermore, WRF results, confirmed by data, show that at times sea breezes develop from both coasts, leading to convergence in the middle of the country; at other times, the large-scale wind impedes the formation of sea breezes on one or both coasts. Simulations also indicate greater land/sea-breeze activity in the summer than in the winter. Differences in the diurnal evolution of surface winds over land and water are found to be related to differences in the boundary layer stability. Overall, the results indicate that the WRF Model as configured here yields reliable simulations and can be used for various practical applications.

The Performance of a Scale-Aware Nonlocal PBL Scheme for the Subkilometer Simulation of a Deep CBL over the Taklimakan Desert

Although realistic representation of the convective boundary layer (CBL) in the desert region in Northwest China is important for weather forecasts and climate simulations, evaluations of the performance of various planetary boundary layer (PBL) schemes in simulating the CBL in the region are rare. In this study, the performance of a scale-aware PBL scheme newly implemented into the Weather Research and Forecasting (WRF) model in simulating the CBL in the Taklimakan desert is evaluated based on a comparison with both the WRF-LES simulations and observations, with the focus on scale dependencies of the simulations compared to the conventional PBL scheme. A series of simulations are performed with a scale-aware PBL scheme (Shin-Hong) and the conventional PBL scheme (YSU) for a deep CBL observed at Tazhong station in the central Taklimakan on 1 July 2016. The CBL was over 5000 m deep with wider and deeper rolls than in a shallow boundary layer. The results showed that the vertical structure simulated with the Shin-Hong scheme was closer to that in both the WRF-LES (large-eddy-simulation) and observations than that simulated with the YSU. The simulation with the scale-aware scheme reproduced cellular rolls similar to those in the WRF-LES, while the conventional PBL scheme struggled to trigger intense convective cells rather than cellular rolls. The results strongly suggest that the scale-aware nonlocal PBL scheme can be used to adequately reproduce the scale and evolution of the observed rolls in the deep CBL in Taklimakan desert at subkilometer resolutions.

PBL parameterization schemes sensitivity analysis on WRF modeling of a tornadic event environment in Skala Lakonia in September 2015

The aim of the present study is the evaluation of the thermodynamic characteristics of a high impact tornadic event in Skala Lakonia, southern Peloponnese, on September 21, 2015. The study was carried out using the Weather Research and Forecasting (AR-WRF) mesoscale model allowing high spatial resolution for the analysis. For the particular model setup used in this study, different planetary boundary layer (PBL) parameterization schemes were tested. The current study downscales ECMWF reanalysis initial and lateral boundary conditions to 1km spatial resolution to assess the mixing layer convective available potential energy (MLCAPE – MCAPE hereafter), the reflectivity perceived by radar (dBZ), the storm relative helicity (SRH) and the energy helicity index (EHI) for the different PBL schemes. The validation was conducted using in situ and remote sensing observations of the surrounding area. The results illustrate deep convective activity with MCAPE values indicating a strong instability in the region of interest using both local and non-local closure PBL schemes. The study evaluates the affectability of the PBL physical processes with regard to the thermodynamic variables on the model performance. Finally, Mellor Yamada Janjic local (MYJ) closure PBL scheme is the best PBL scheme for studying this tornado since its simulation results display deep convection activity along with the highest validation scores.

Evaluation of WRF PBL parameterization schemes against direct observations during a dry event over the Ganges valley

Accurate representations of the planetary boundary layer (PBL) are important in all weather forecast systems, especially in simulations of turbulence, wind and air quality in the lower atmosphere. In the present study, detailed observations from the Cloud Aerosol Interaction and Precipitation Enhancement Experiment - Integrated Ground based Observational Campaign (CAIPEEX-IGOC) 2014 comprising of the complete surface energy budget and detailed boundary layer observations are used to validate Advanced Research Weather Research and Forecasting (WRF) model simulations over a diverse terrain over the Ganges valley region, Uttar Pradesh, India. A drying event in June 2014 associated with a heat wave is selected for validation.Six local and nonlocal PBL schemes from WRF at 1km resolution are compared with hourly observations during the diurnal cycle. Near-surface observations of weather parameters, radiation components and eddy covariance fluxes from micrometeorological tower, and profiles of variables from microwave radiometer, and radiosonde observations are used for model evaluations. Models produce a warmer, drier surface layer with higher wind speed, sensible heat flux and temperature than observations. Layered boundary layer dynamics, including the residual layer structure as illustrated in the observations over the Ganges valley are missed in the model, which lead to deeper mixed layers and excessive drying.Although it is difficult to identify any single scheme as the best, the qualitative and quantitative analyses for the entire study period and overall reproducibility of the observations indicate that the MYNN2 simulations describe lower errors and more realistic simulation of spatio-temporal variations in the boundary layer height.

Constraining the sulfur dioxide degassing flux from Turrialba volcano, Costa Rica using unmanned aerial system measurements

Observed sulfur dioxide (SO2) mixing ratios onboard unmanned aerial systems (UAS) during March 11–13, 2013 are used to constrain the three-day averaged SO2 degassing flux from Turrialba volcano within a Bayesian inverse modeling framework. A mesoscale model coupled with Lagrangian stochastic particle backward trajectories is used to quantify the source-receptor relationships at very high spatial resolutions (i.e., <1km). The model shows better performance in reproducing the near-surface meteorological properties and observed SO2 variations when using a first-order closure non-local planetary boundary layer (PBL) scheme. The optimized SO2 degassing fluxes vary from 0.59±0.37 to 0.83±0.33ktd−1 depending on the PBL scheme used. These fluxes are in good agreement with ground-based gas flux measurements, and correspond to corrective scale factors of 8–12 to the posteruptive SO2 degassing rate in the AeroCom emission inventory. The maximum a posteriori solution for the SO2 flux is highly sensitive to the specification of prior and observational errors, and relatively insensitive to the SO2 loss term and temporal averaging of observations. Our results indicate relatively low degassing activity but sustained sulfur emissions from Turrialba volcano to the troposphere during March 2013. This study demonstrates the utility of low-cost small UAS platforms for volcanic gas composition and flux analysis.

Implementation of non-local boundary layer schemes in the Regional Atmospheric Modeling System and its impact on simulated mesoscale circulations

This paper proposes the implementation of different non-local Planetary Boundary Layer schemes within the Regional Atmospheric Modeling System (RAMS) model. The two selected PBL parameterizations are the Medium-Range Forecast (MRF) PBL and its updated version, known as the Yonsei University (YSU) PBL. YSU is a first-order scheme that uses non-local eddy diffusivity coefficients to compute turbulent fluxes. It is based on the MRF, and improves it with an explicit treatment of the entrainment. With the aim of evaluating the RAMS results for these PBL parameterizations, a series of numerical simulations have been performed and contrasted with the results obtained using the Mellor and Yamada (MY) scheme, also widely used, and the standard PBL scheme in the RAMS model. The numerical study carried out here is focused on mesoscale circulation events during the summer, as these meteorological situations dominate this season of the year in the Western Mediterranean coast. In addition, the sensitivity of these PBL parameterizations to the initial soil moisture content is also evaluated. The results show a warmer and moister PBL for the YSU scheme compared to both MRF and MY. The model presents as well a tendency to overestimate the observed temperature and to underestimate the observed humidity, considering all PBL schemes and a low initial soil moisture content. In addition, the bias between the model and the observations is significantly reduced moistening the initial soil moisture of the corresponding run. Thus, varying this parameter has a positive effect and improves the simulated results in relation to the observations. However, there is still a significant overestimation of the wind speed over flatter terrain, independently of the PBL scheme and the initial soil moisture used, even though a different degree of accuracy is reproduced by RAMS taking into account the different sensitivity tests.