In this study, the Weather Research and Forecasting model with online coupled chemistry (WRF-Chem) is applied to simulate an intense Saharan dust outbreak event that took place over the Mediterranean in May 2014. Comparison of a simulation using a physics-based desert dust emission scheme with a numerical experiment using a simplified (minimal) emission scheme is included to highlight the advantages of the former. The model was found to reproduce well the synoptic meteorological conditions driving the dust outbreak: an omega-like pressure configuration associated with a cyclogenesis in the Atlantic coasts of Spain. The model performances in reproducing the atmospheric desert dust load were evaluated using a multi-platform observational dataset of aerosol and desert dust properties, including optical properties from satellite and ground-based sun photometers and lidars, plus in situ particulate matter mass concentration (PM) data. This comparison allowed us to investigate the model ability in reproducing both the horizontal and the vertical displacement of the dust plume, as well as its evolution in time. The comparison with satellite (MODIS-Terra) and sun photometers (AERONET) showed that the model is able to reproduce well the horizontal field of the aerosol optical depth (AOD) and its evolution in time (temporal correlation coefficient with AERONET of 0.85). On the vertical scale, the comparison with lidar data at a single site (Rome, Italy) confirms that the desert dust advection occurs in several, superimposed "pulses" as simulated by the model. Cross-analysis of the modeled AOD and desert dust emission fluxes further allowed for the source regions of the observed plumes to be inferred. The vertical displacement of the modeled dust plume was in rather good agreement with the lidar soundings, with correlation coefficients among aerosol extinction profiles up to 1 and mean discrepancy of about 50 %. The model–measurement comparison for PM₁₀ and PM_2.5 showed a good temporal matching, although it revealed a marked overestimation of PM₁₀ and PM_2.5 (of the order of 70 % during the dust peak). For PM₁₀, it was also possible to investigate the accordance between the model- and the measurement-based dust PM₁₀, this confirming the model PM₁₀ overestimation to be related to over-predicted dust mass up to a factor of 140 %. In all the model–measurement comparisons performed, the enhanced capabilities of the physics-based emission scheme with respect to its simplified, minimal version were evident and are documented.

Between 19 and 22 January 2014, a baroclinic wave moving eastward from the Atlantic Ocean generated a cut-off low over the Strait of Gibraltar and was responsible for the subsequent intensification of an extra-tropical cyclone. This system exhibited tropical-like features in the following stages of its life cycle and remained active for approximately 80 h, moving along the Mediterranean Sea from west to east, eventually reaching the Adriatic Sea. Two different modeling approaches, which are comparable in terms of computational cost, are analyzed here to represent the cyclone evolution. First, a multi-physics ensemble using different microphysics and turbulence parameterization schemes available in the WRF (weather research and forecasting) model is employed. Second, the COAWST (coupled ocean–atmosphere wave sediment transport modeling system) suite, including WRF as an atmospheric model, ROMS (regional ocean modeling system) as an ocean model, and SWAN (simulating waves in nearshore) as a wave model, is used. The advantage of using a coupled modeling system is evaluated taking into account air–sea interaction processes at growing levels of complexity. First, a high-resolution sea surface temperature (SST) field, updated every 6 h, is used to force a WRF model stand-alone atmospheric simulation. Later, a two-way atmosphere–ocean coupled configuration is employed using COAWST, where SST is updated using consistent sea surface fluxes in the atmospheric and ocean models. Results show that a 1D ocean model is able to reproduce the evolution of the cyclone rather well, given a high-resolution initial SST field produced by ROMS after a long spin-up time. Additionally, coupled simulations reproduce more accurate (less intense) sea surface heat fluxes and a cyclone track and intensity, compared with a multi-physics ensemble of standalone atmospheric simulations.

Occasionally, storms that share many features with tropical cyclones, including the presence of a quasi-circular “eye” a warm core and strong winds, are observed in the Mediterranean. Generally, they are known as Medicanes, or tropical-like cyclones (TLC). Due to the intense wind forcings and the consequent development of high wind waves, a large number of sea spray droplets—both from bubble bursting and spume tearing processes—are likely to be produced at the sea surface. In order to take into account this process, we implemented an additional Sea Spray Source Function (SSSF) in WRF-Chem, model version 3.6.1, using the GOCART (Goddard Chemistry Aerosol Radiation and Transport) aerosol sectional module. Traditionally, air-sea momentum fluxes are computed through the classical Charnock relation that does not consider the wave-state and sea spray effects on the sea surface roughness explicitly. In order to take into account these forcing, we implemented a more recent parameterization of the sea surface aerodynamic roughness within the WRF surface layer model, which may be applicable to both moderate and high wind conditions. The implemented SSSF and sea surface roughness parameterization have been tested using an operative model sequence based on COAWST (Coupled Ocean Atmosphere Wave Sediment Transport) and WRF-Chem. The third-generation wave model SWAN (Simulating Waves Nearshore), two-way coupled with the WRF atmospheric model in the COAWST framework, provided wave field parameters. Numerical simulations have been integrated with the WRF-Chem chemistry package, with the aim of calculating the sea spray generated by the waves and to include its effect in the Charnock roughness parametrization together with the sea state effect. A single case study is performed, considering the Medicane that affected south-eastern Italy on 26 September 2006. Since this Medicane is one of the most deeply analysed in literature, its investigation can easily shed some light on the feedbacks between sea spray and drag coefficients.