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.

The Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) provides data at 0.5° × 0.625° resolution covering a period from 1 January 1980 to the present. Natural and anthropogenic aerosols are simulated in MERRA-2, considering the Goddard chemistry, aerosol, radiation, and transport model. This model simulates the sources, sinks, and chemistry of mixed aerosol tracers: dust, sea salt, hydrophobic and hydrophilic black carbon and organic carbon, and sulfate. MERRA-2 aerosol reanalysis is a pioneering tool for investigating air quality issues, noteworthy for its global coverage and its distinction of aerosol speciation expressed in the form of aerosol optical depth (AOD). The aim of this work was to use the MERRA-2 reanalysis to study urban air pollution at a national scale by analyzing the AOD. AOD trends were evaluated for a 30-year period (1987–2017) over five Italian cities (Milan, Rome, Cagliari, Taranto, and Palermo) in order to investigate the impacts of urbanization, industrialization, air quality regulations, and regional transport on urban aerosol load. AOD evolution predicted by the MERRA-2 model in the period 2002–2017 showed a generalized decreasing trend over the selected cities. The anthropogenic signature on total AOD was between 50% and 80%, with the largest contribution deriving from sulfate.

In this study, we investigate a series of Saharan dust outbreaks toward the Mediterranean basin that occurred in late June 2021. In particular, we analyze the effect of mineral dust aerosols on radiation and cloud properties (direct, semi-direct and indirect effects), and in turn, on meteorological parameters. This is achieved by running the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) over a domain covering North Africa and the Central Mediterranean Basin. The simulations were configured using a gradual coupling strategy between the GOCART aerosol model and the Goddard radiation and microphysics schemes available in the WRF-Chem package. A preliminary evaluation of the model performances was conducted in order to verify its capability to correctly reproduce the amount of mineral dust loaded into the atmosphere within the spatial domain considered. To this purpose, we used a suite of experimental data from ground- and space-based remote sensing measurements. This comparison highlighted a model over-estimation of aerosol optical properties to the order of 20%. The evaluation of the desert dust impact on the radiation budget, achieved by comparing the uncoupled and the fully coupled (aerosol–radiation–clouds) simulation, shows that mineral dust induces a net (shortwave–longwave) cooling effect to the order of −10 W m⁻². If we consider the net dust radiative forcing, the presence of dust particles induces a small cooling effect at the top of the atmosphere (−1.2 W m⁻²) and a stronger cooling at the surface (−14.2 W m⁻²). At the same time, analysis of the perturbation on the surface energy budget yields a reduction of −7 W m⁻² when considering the FULL-coupled simulation, a positive perturbation of +3 W m⁻² when only considering microphysics coupling and −10.4 W m⁻² when only considering radiation coupling. This last result indicates a sort of “superposition” of direct, indirect and semi-direct effects of dust on the radiation budget. This study shows that the presence of dust aerosols significantly influences radiative and cloud properties and specifically the surface energy budget. This suggests (i) that dust effects should be considered in climate models in order to increase the accuracy of climate predictions over the Mediterranean region and (ii) the necessity of performing fully coupled simulations including aerosols and their effects on meteorology at a regional scale.

Carbon assimilation and wood production are influenced by environmental conditions and endogenous factors, such as species auto-ecology, age, and hierarchical position within the forest structure. Disentangling the intricate relationships between those factors is more pressing than ever due to climate change’s pressure. We employed the 3D-CMCC-FEM model to simulate undisturbed forests of different ages under four climate change (plus one no climate change) Representative Concentration Pathways (RCP) scenarios from five Earth system models. In this context, carbon stocks and increment were simulated via total carbon woody stocks and mean annual increment, which depends mainly on climate trends. We find greater differences among different age cohorts under the same scenario than among different climate scenarios under the same age class. Increasing temperature and changes in precipitation patterns led to a decline in above-ground biomass in spruce stands, especially in the older age classes. On the contrary, the results show that beech forests will maintain and even increase C-storage rates under most RCP scenarios. Scots pine forests show an intermediate behavior with a stable stock capacity over time and in different scenarios but with decreasing mean volume annual increment. These results confirm current observations worldwide that indicate a stronger climate-related decline in conifers forests than in broadleaves.

Different configurations for the Weather Research and Forecasting (WRF) model were evaluated to improve wind and temperature fields predictions in the Northern Sahara and the Mediterranean basin. Eight setups, associated with different combinations of the surface layer physical parameters, the land surface model, and the grid nudging parameters, were considered. Numerical simulations covered the entire month of November 2017. Model results were compared with surface data from meteorological stations. The introduction of the grid nudging parameters leads to a general improvement of the modeled 10 m wind speed and 2 m temperature. In particular, nudging of wind speed parameter inside the planetary boundary layer (PBL) provides the most remarkable differences. In contrast, the nudging of temperature and relative humidity parameters inside the PBL may be switched off to reduce computational time and data storage. Furthermore, it was shown that the prediction of the 10 m wind speed and 2 m temperature is quite sensitive to the choice of the surface layer scheme and the land surface model. This paper provides useful suggestions to improve the setup of the WRF model in the Northern Sahara and the Mediterranean basin. These results are also relevant for topics related with the emission of mineral dust and sea spray within the Mediterranean region.

Biogenic volatile organic compounds (BVOCs) emitted from the natural ecosystem are highly reactive and can thus impact air quality and aerosol radiative forcing. BVOC emission models (e.g., Model of Emissions of Gases and Aerosols from Nature – MEGAN) in global and regional chemical transport models still have large uncertainties in estimating biogenic trace gases because of uncertainties in emission activity factors, specification of vegetation type, and plant emission factors. This study evaluates a set of updates made to MEGAN v2.04 in the Weather Research and Forecasting model coupled with chemistry (WRF-Chem version 3.9). Our study considers four simulations for each update made to MEGAN v2.04: (i) a control run with no changes to MEGAN, (ii) a simulation with the emission activity factors modified following MEGAN v2.10, (iii) a simulation considering the changes to the plant functional type (PFT) emission factor, and (iv) a simulation with the isoprene emission factor calculated within the MEGAN module instead of being prescribed by the input database. We evaluate two regions, Europe and the southeastern United States, by comparing WRF-Chem results to ground-based monitoring observations in Europe (i.e., AirBase database) and aircraft observations obtained during the NOMADSS field campaign. We find that the updates to MEGAN v2.04 in WRF-Chem caused overpredictions in ground-based ozone concentrations in Europe and in isoprene mixing ratios compared to aircraft observations in the southeastern US. The update in emission activity factors caused the largest biases. These results suggest that further experimental and modeling studies should be conducted to address potential shortcomings in BVOC emission models.