This work presents a regridding procedure applied to the nitrogen dioxide (NO₂) tropospheric column data, derived from the Copernicus Sentinel 5 Precursor Tropospheric Monitoring Instrument (S5P/TROPOMI). The regridding has been performed to provide a better comparison with punctual surface observations. It will be demonstrated that TROPOMI NO₂ tropospheric column data show improved consistency with in situ surface measurements once the satellite retrievals are scaled to 1 km spatial sampling. A geostatistical technique, i.e., the ordinary kriging, has been applied to improve the spatial distribution of Level 2 TROPOMI NO₂ data, which is originally sparse and uneven because of gaps introduced by clouds, to a final spatial, regular, sampling of 1 km × 1 km. The analysis has been performed for two study areas, one in the North and the other in the South of Italy, and for May 2018-April 2020, which also covers the period January 2020-April 2020 of COVID-19 diffusion over the Po Valley. The higher spatial sampling NO₂ dataset indicated as Level 3 data, allowed us to explore spatial and seasonal data variability, obtaining better information on NO₂ sources. In this respect, it will be shown that NO₂ concentrations in March 2020 have likely decreased as a consequence of the lockdown because of COVID-19, although the far warmest winter season ever recorded over Europe in 2020 has favored a general NO₂ decrease in comparison to the 2019 winter. Moreover, the comparison between NO₂ concentrations related to weekdays and weekend days allowed us to show the strong correlation of NO₂ emissions with traffic and industrial activities. To assess the quality and capability of TROPOMI NO₂ observations, we have studied their relationship and correlation with in situ NO₂ concentrations measured at air quality monitoring stations. We have found that the correlation increases when we pass from Level 2 to Level 3 data, showing the importance of regridding the satellite data. In particular, correlation coefficients of Level 3 data, which range between 0.50–0.90 have been found with higher correlation applying to urban, polluted locations and/or cities.

This study investigates the value of satellite-based observational algorithms in supporting numerical weather prediction (NWP) for improving the alert and monitoring of extreme rainfall events. To this aim, the analysis of the very intense precipitation that affected the city of Livorno on 9 and 10 September 2017 is performed by applying three remote sensing techniques based on satellite observations at infrared/visible and microwave frequencies and by using maps of accumulated rainfall from the weather research and forecasting (WRF) model. The satellite-based observational algorithms are the precipitation evolving technique (PET), the rain class evaluation from infrared and visible observations (RainCEIV) technique and the cloud classification mask coupling of statistical and physics methods (C-MACSP). Moreover, the rain rates estimated by the Italian Weather Radar Network are also considered to get a quantitative evaluation of RainCEIV and PET performance. The statistical assessment shows good skills for both the algorithms (for PET: bias = 1.03, POD = 0.76, FAR = 0.26; for RainCEIV: bias = 1.33, POD = 0.77, FAR = 0.41). In addition, a qualitative comparison among the three technique outputs, rain rate radar maps, and WRF accumulated rainfall maps is also carried out in order to highlight the advantages of the different techniques in providing real-time monitoring, as well as quantitative characterization of rainy areas, especially when rain rate measurements from Weather Radar Network and/or from rain gauges are not available.

In this study, we compare different nowcasting techniques based upon the calculation of motion vector fields derived from spectral channels of Meteosat Second Generation—Spinning Enhanced Visible and InfraRed Imager (MSG-SEVIRI). The outputs of the nowcasting techniques are used as inputs to the Advanced Model for Estimation of Surface solar Irradiance from Satellite (AMESIS), for predicting surface solar irradiance up to 2 h in advance. In particular, the first part of the methodology consists in projecting the time evolution of each MSG-SEVIRI channel (for every pixel in the spatial domain) through extrapolation of a displacement vector field obtained by matching similar patterns within two successive MSG-SEVIRI data images. Different ways to implement the above method result in substantial differences in the predicted trajectory, leading to different performances depending on the time interval of interest. All the nowcasting techniques considered here systematically outperform the simple persistence method for all MSG-SEVIRI channels and for each case study used in this work; importantly, this occurs across the entire 2 h period of the forecast. In the second part of the algorithm, the predicted irradiance maps computed with AMESIS from the forecasted radiances, are shown to be in good agreement with irradiances derived from MSG measured radiances and improve on numerical weather model predictions, thus providing a feasible alternative for nowcasting surface solar radiation. The results show that the mean values for correlation, bias, and root mean square error vary across the time interval, ranging between 0.94, −1 W/m ${}^2$ , 61 W/m ${}^2$ after 15 min, and 0.73, −18 W/m ${}^2$ , 147 W/m ${}^2$ after 2 h, respectively.

In this study, the Meteosat Second Generation (MSG)—Spinning Enhanced Visible and Infrared Imager (SEVIRI) High Resolution Visible channel (HRV) is used in synergy with the narrow band MSG-SEVIRI channels for daytime fog detection. A new algorithm, named MSG-SEVIRI SatFog, has been designed and implemented. MSG-SEVIRI SatFog provides the indication of the presence of fog in near real time and at the high spatial resolution of the HRV channel. The HRV resolution is useful for detecting small scale daytime fog that would be missed in the MSG-SEVIRI low spatial resolution channels. By combining textural, physical and tonal tests, a distinction between fog and low stratus is performed for pixels identified as low/middle clouds or clear by the Classification-MAsk Coupling of Statistical and Physical Methods (C-MACSP) cloud detection algorithm. Suitable thresholds have been determined using a specific dataset covering different geographical areas, seasons and time of the day. MSG-SEVIRI SatFog is evaluated against METeorological Aerodrome Reports (METAR) data observations. Evaluation results in an accuracy of 69.9%, a probability of detection of 68.7%, a false alarm ratio of 31.3% and a probability of false detection of 30.0%.

Thermal infrared remote sensing measurements have greatly improved in terms of spectral, spatial, and temporal resolution. These improvements are producing a clearer picture of the land surface and Earth atmospheric composition than ever before. Nevertheless, the analysis of this big quantity of data presents important challenges due to incomplete temporal and spatial recorded information. The aim of the present paper is to discuss a methodology to retrieve missing values of some interesting geophysical variables on a spatial field retrieved from spatially scattered infrared satellite observations in order to yield level 3, regularly gridded, data. The technique is based on a 2-Dimensional (2D) Optimal Interpolation (OI) scheme and is derived from the broad class of Kalman filter or Bayesian estimation theory. The goodness of the approach has been tested on 15-min temporal resolution Spinning Enhanced Visible and Infrared Imager (SEVIRI) emissivity and surface temperature (ST) products over South Italy (land and sea), on Infrared Atmospheric Sounding Interferometer (IASI) atmospheric ammonia ( $N{H}_3$ ) concentration over North Italy and carbon monoxide ( $CO$ ), sulfur dioxide ( $S{O}_2$ ) and $N{H}_3$ concentrations over China. All these gases affect air quality. Moreover, sea surface temperature (SST) retrievals have been compared with gridded data from MODIS (Moderate-resolution Imaging Spectroradiometer) observations. For gases concentration we have considered data from 3 different emission inventories, that is, Emissions Database for Global Atmospheric Research v3.4.2 (EDGARv3.4.2), the Regional Emission inventory in ASiav3.1 (REASv3.1) and MarcoPolov0.1, plus an independent study. The results show the efficacy of the proposed strategy to better capture the daily cycle for surface parameters and to detect hotspots of severe emissions from gas sources affecting air quality such as $CO$ and $N{H}_3$ and, therefore, to yield valuable information on the variability of gas concentration to complete ground stations measurements.

A new algorithm for the estimation of atmospheric temperature (T) and water vapor (WV) vertical profiles in nonprecipitating conditions is presented. The microwave random forest temperature and water vapor (MiRTaW) profiling algorithm is based on the random forest (RF) technique and it uses microwave (MW) sounding from the Advanced Technology Microwave Sounder (ATMS) onboard the Suomi National Polar-orbiting Partnership (SNPP) satellite. Three different data sources were chosen for both training and validation purposes, namely, the ERA-Interim from the European Centre for Medium-Range Weather Forecasts (ECMWF), the Infrared Atmospheric Sounding Interferometer Atmospheric Temperature Water Vapour and Surface Skin Temperature (IASI L2 v6) from the Meteorological Operational satellites of the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT), and the radiosonde observations from the Integrated Global Radiosonde Archive (IGRA). The period from 2012 to 2016 was considered in the training dataset; particular attention was paid to the instance selection procedure, in order to reduce the full training dataset with negligible information loss. The out-of-bag (OOB) error was computed and used to select the optimal RF parameters. Different RFs were trained, one for each vertical level: 32 levels for T (within 10–1000 hPa) and 23 levels for WV (200–1000 hPa). The validation of the MiRTaW profiling algorithm was conducted on a dataset from 2017. The mean bias error (MBE) of T vertical profiles ranges within about (−0.4–0.4) K, while for the WV mixing ratio, the MBE starts at ~0.5 g/kg near the surface and decreases to ~0 g/kg at 200 hPa level, in line with the expectations.