Meteosat Second Generation Research Articles

Impact of Aerosol and Cloud on the Solar Energy Potential over the Central Gangetic Himalayan Region

We examine the impact of atmospheric aerosols and clouds on the surface solar radiation and solar energy at Nainital, a high-altitude remote location in the central Gangetic Himalayan region (CGHR). For this purpose, we exploited the synergy of remote-sensed data in terms of ground-based AERONET Sun Photometer and satellite observations from the MODerate Resolution Imaging Spectroradiometer (MODIS) and the Meteosat Second Generation (MSG), with radiative transfer model (RTM) simulations and 1 day forecasts from the Copernicus Atmosphere Monitoring Service (CAMS). Clouds and aerosols are one of the most common sources of solar irradiance attenuation and hence causing performance issues in the photovoltaic (PV) and concentrated solar power (CSP) plant installations. The outputs of RTM results presented with high accuracy under clear, cloudy sky and dust conditions for global horizontal (GHI) and beam horizontal irradiance (BHI). On an annual basis the total aerosol attenuation was found to be up to 105 kWh m−2 for the GHI and 266 kWh m−2 for BHI, respectively, while the cloud effect is much stronger with an attenuation of 245 and 271 kWh m−2 on GHI and BHI. The results of this study will support the Indian solar energy producers and electricity handling entities in order to quantify the energy and financial losses due to cloud and aerosol presence.

On the geometric accuracy and stability of MSG SEVIRI images

Monitoring the geometric quality and temporal stability of Earth Observation satellite sensors is a major task of satellite vendors. As an owner and operator of widely used meteorological satellites, EUMETSAT (European Organisation for the Exploitation of Meteorological Satellites) carries out operational activities and scientific studies on the radiometric and geometric quality assessments of geostationary Meteosat satellites among others. As part of Level 1 Product Monitoring - Evolution Studies initiated by EUMETSAT, a new methodology has been developed and implemented as a prototype software, namely the Geometric Quality Assessment (GQA) Tool; which enhances image local texture, extracts prominent features and matches them via area-based least squares, and eliminates possible outliers using a set of rules. This article presents the results of algorithm validation activities carried out by using time series images of the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) sensor aboard Meteosat Second Generation (MSG) satellites and MEdium Resolution Imaging Spectrometer (MERIS) global mosaics as reference to evaluate short- and long-term geometric temporal stability. Meteosat-11 SEVIRI time series data from a full-day and one year were employed for this purpose. In addition, the geometric quality of Meteosat-8 over Indian Ocean and Meteosat-10 in rapid scan mode was assessed for single acquisitions. The results show that the geometric accuracy of the evaluated datasets yield to the mean shifts up to one km in East-West (x) and North-South (y) directions, with larger errors in the y direction. The image inner geometric accuracy assessments revealed displacement errors with strip pattern in the three spectral bands (i.e. high-resolution visible, infra-red 10.8 and visible 0.8) evaluated within the framework.

Rain Area Detection in South-Western Kenya by Using Multispectral Satellite Data from Meteosat Second Generation.

This study presents a rain area detection scheme that uses a gradient based adaptive technique for daytime and nighttime rain area detection and correction from reflectance and infrared (IR) brightness temperatures data of the Meteosat Second Generation (MSG) satellite. First, multiple parametric rain detection models developed from MSG’s reflectance and IR data were calibrated and validated with rainfall data from a dense network of rain gauge stations and investigated to determine the best model parameters. The models were based on a conceptual assumption that clouds characterised by the top properties, e.g., high optical thickness and effective radius, have high rain probabilities and intensities. Next, a gradient based adaptive correction technique that relies on rain area-specific parameters was developed to reduce the number and sizes of the detected rain areas. The daytime detection with optical (VIS0.6) and near IR (NIR1.6) reflectance data achieved the best detection skill. For nighttime, detection with thermal IR brightness temperature differences of IR3.9-IR10.8, IR3.9-WV73 and IR108-WV62 showed the best detection skill based on general categorical statistics. Compared to the Global Precipitation Measurement (GPM) Integrated Mult-isatellitE Retrievals for GPM (IMERG) and the gauge station data from the southwest of Kenya, the model showed good agreement in the spatial dynamics of the detected rain area and rain rate.

Seasonal Net Shortwave Radiation of Bare Arable Land in Poland and Israel According to Roughness and Atmospheric Irradiance

Tillage of arable fields, using for instance a smoothing harrow, may increase the magnitude of albedo of such soil surfaces depending on the location, the sun’s illumination and atmospheric components. As these soil surfaces absorb less shortwave radiation compared to plowed soils, the result is an atmospheric cooling and a positive effect on the Earth’s climate. This paper is the follow-on of a previous study aimed at quantifying the seasonal dynamics of net shortwave radiation reflected by bare air-dried arable land areas located in contrasting environments, i.e. Poland and Israel. Soil tillage includes a plow, a disk harrow, and a smoothing harrow. Previous work concentrated on the estimate of net shortwave radiation under clear-sky theoretical scenarios, whereas the present study deals with a realistic atmosphere throughout the year 2014. This latter is characterized by the observations of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) instrument on board the Meteosat Second Generation (MSG). The variations of the net shortwave radiation for the selected bare arable land areas were assessed in combining observations from Landsat 8 images and digital maps of land use and soil, plus model equations that calculate the diurnal variations of the broadband blue-sky albedo with roughness inclusive. The daily amount of net shortwave radiation for air-dried bare arable land in Poland and Israel for the time their spatial coverage is the largest was found to be about 40–50% and 10% lower, respectively, in cloudy-sky conditions compared to clear-sky conditions.

Validation and consistency assessment of land surface temperature from geostationary and polar orbit platforms: SEVIRI/MSG and AVHRR/Metop

Land Surface Temperature (LST) is a fundamental state variable for land surface processes, long available from satellite observations in the thermal infrared. Here we make an assessment of LST products derived from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on-board the geostationary satellite Meteosat Second Generation (MSG), and from the Advanced Very High Resolution Radiometer (AVHRR) on-board Metop polar-orbiting satellites, as provided by the Satellite Application Facility on Land Surface Analysis (LSA-SAF). The validation of LSA-SAF LST products uses ground observations taking into account the representativeness of each station at pixel scale. The in situ - satellite comparisons reveal overall accuracies of 0.13 °C for SEVIRI LST and −0.32 °C for AVHRR-based LST. Better matches are usually found for night-time observations, highlighting the influence of LST spatial and temporal variability as well as of sensor view angle on satellite daytime estimates. The averaged difference between AVHRR and SEVIRI LST products remains mostly within ±1 °C. The impact of different viewing and illumination angles is particularly pronounced for daytime matchups, leading to higher standard deviation of differences than in the case of night-time. Over heterogeneous landscapes, this primarily results from illumination-shadowing effects in daytime observations present in all LST estimates derived from infrared observations. These differences are smaller for night-time observations, or over homogeneous surfaces in general: in these cases, the dependency of LST differences on viewing geometry is smoother and likely explained by the impact of view angle on surface emissivity. Overall the LSA-SAF SEVIRI and AVHRR LST datasets are consistent and meet high accuracy standards. Their production is ensured by the LSA-SAF throughout the life-cycle of the sensors' satellite platforms, making them good candidates for long term applications and studies.

Daily grass reference evapotranspiration with Meteosat Second Generation shortwave radiation and reference ET products

This study assesses the accuracy of estimating daily grass reference evapotranspiration (PM-ETo) using daily shortwave radiation (Rs) and reference evapotranspiration (ETREF) products provided by the Meteosat Second Generation (MSG) geostationary satellite delivered by the Satellite Applications Facility on Land Surface Analysis (LSA-SAF) framework. The accuracy of using reanalysis ERA5 shortwave radiation data (Rs ERA5) provided by the European Center for Medium-Range Weather Forecasts (ECMWF) is also evaluated. The assessments were performed using observed weather variables at 37 weather stations distributed across continental Portugal, where climate conditions range from semi-arid to humid, and 12 weather stations located in Azores islands, characterized by humid, windy and often cloudy conditions. This study's use of data from a variety of climate conditions contributed to a unique and innovative assessment of the usability of LSA-SAF and ERA5 products for ETo estimation. The first assessment focused on comparing LSA-SAF estimates of Rs (Rs LSA-SAF) against ground stations (Rs ground). The results showed a good matching between the two Rs data sets for continental Portugal but a tendency for Rs LSA-SAF to under-estimate Rs ground in the cloudy islands of Azores. ETo values computed using Rs LSA-SAF data and observed temperature, humidity and wind speed (ETo LSA-SAF) were then compared with PM-ETo estimates with ground-based data, which were used as benchmark; input data of temperature and humidity needed for PM-ETo were quality checked for surface aridity effects. It was observed that ETo LSA-SAF is strongly correlated with PM-ETo (R2 > 0.97) for most locations in continental Portugal, with regression coefficient of a linear regression forced to the origin ranging between 0.95 and 1.05, mean root mean square error (RMSE) of 0.13 mm d−1, and Nash and Sutcliff efficiency of modeling (EF) above 0.95. For most Azores locations, ETo LSA-SAF over-estimated PM-ETo. This is likely a consequence of the high spatio-temporal heterogeneity of weather conditions that occur in these oceanic islands together with the different footprints of satellite (averaged over the pixel) and station observations. Reanalysis ERA5 shortwave radiation data presented similar behavior to the LSA-SAF products, however with slightly lower accuracy. The daily LSA-SAF ETREF product (ETREF LSA-SAF) was assessed and results have shown a good accuracy of this product, with acceptable RMSE and high EF values, for continental Portugal but a low accuracy for the Azores islands. A simplified bias correction approach was shown to improve both ETo derived from the LSA-SAF products, namely for Azores stations, which seem to be representative of smaller areas. The use of the FAO-PM temperature approach (PMT) was also assessed using the Rs LSA-SAF and Rs ERA5 data, which showed a superiority of the LSA-SAF product for ETo estimations (ETo PMT LSA-SAF). No significant differences (p < 0.05) were observed in terms of the median value of the RMSE when adopting ETo PMT and ETREF LSA-SAF. Differently, results showed that using the Rs LSA-SAF in the PMT approach (ETo PMT LSA-SAF) produces significantly better RMSE results than ETo PMT and ETREF LSA-SAF. Overall, the performed assessment allows concluding that the use of Rs LSA-SAF, and to a lesser extent the use of the Rs ERA5, highly improves the accuracy of computation of ETo when Rs observations are not available, including when only temperature data are accessible. The use of the ETREF LSA-SAF product is a good alternative when observed weather data are not available.

Short-Term Forecasting of Large-Scale Clouds Impact on Downwelling Surface Solar Irradiation

This study focuses on the use of cloud motion vectors (CMV) and fast radiative transfer models (FRTM) in the prospect of forecasting downwelling surface solar irradiation (DSSI). Using near-real-time cloud optical thickness (COT) data derived from multispectral images from the spinning enhanced visible and infrared imager (SEVIRI) onboard the Meteosat second generation (MSG) satellite, we introduce a novel short-term forecasting system (3 h ahead) that is capable of calculating solar energy in large-scale (1.5 million-pixel area covering Europe and North Africa) and in high spatial (5 km over nadir) and temporal resolution (15 min intervals). For the operational implementation of such a big data computing architecture (20 million simulations in less than a minute), we exploit a synergy of high-performance computing and deterministic image processing technologies (dense optical flow estimation). The impact of clouds forecasting uncertainty on DSSI was quantified in terms of cloud modification factor (CMF), for all-sky and clear sky conditions, for more generalized results. The forecast accuracy was evaluated against the real COT and CMF images under different cloud movement patterns, and the correlation was found to range from 0.9 to 0.5 for 15 min and 3 h ahead, respectively. The CMV forecast variability revealed an overall DSSI uncertainty in the range 18–34% under consecutive alternations of cloud presence, highlighting the ability of the proposed system to follow the cloud movement in opposition to the baseline persistent forecasting, which considers the effects of topocentric sun path on DSSI but keeps the clouds in “fixed” positions, and which presented an overall uncertainty of 30–43%. The proposed system aims to support the distributed solar plant energy production management, as well as electricity handling entities and smart grid operations.

An Improved Approach for Downscaling Coarse-Resolution Thermal Data by Minimizing the Spatial Averaging Biases in Random Forest

Land surface temperature (LST) plays a fundamental role in various geophysical processes at varying spatial and temporal scales. Satellite-based observations of LST provide a viable option for monitoring the spatial-temporal evolution of these processes. Downscaling is a widely adopted approach for solving the spatial-temporal trade-off associated with satellite-based observations of LST. However, despite the advances made in the field of LST downscaling, issues related to spatial averaging in the downscaling methodologies greatly hamper the utility of coarse-resolution thermal data for downscaling applications in complex environments. In this study, an improved LST downscaling approach based on random forest (RF) regression is presented. The proposed approach addresses issues related to spatial averaging biases associated with the downscaling model developed at the coarse resolution. The approach was applied to downscale the coarse-resolution Satellite Application Facility on Land Surface Analysis (LSA-SAF) LST product derived from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) sensor aboard the Meteosat Second Generation (MSG) weather satellite. The LSA-SAF product was downscaled to a spatial resolution of ~30 m, based on predictor variables derived from Sentinel 2, and the Advanced Land Observing Satellite (ALOS) digital elevation model (DEM). Quantitatively and qualitatively, better downscaling results were obtained using the proposed approach in comparison to the conventional approach of downscaling LST using RF widely adopted in LST downscaling studies. The enhanced performance indicates that the proposed approach has the ability to reduce the spatial averaging biases inherent in the LST downscaling methodology and thus is more suitable for downscaling applications in complex environments.

A Scheme to Estimate Diurnal Cycle of Evapotranspiration from Geostationary Meteorological Satellite Observations

The diurnal cycle of evapotranspiration (ET) is significant in studying the dynamics of land–atmosphere interactions. The diurnal ET cycle can be considered as an indicator of dry/wet surface conditions. However, the accuracy of current models in estimating the diurnal ET cycle is generally low. This study developed an improved scheme to estimate the diurnal cycle of ET by solving the surface energy balance equation combined with simplified parameterization, with daily ET as the constraint. Meteosat Second Generation (MSG) land surface temperature, and longwave and shortwave radiation products were the primary inputs. Daily ET was from the remote sensing-based ETMonitor model. The estimated instantaneous (30 min) ET from the improved scheme outperformed the official MSG instantaneous ET product when compared with in situ half-hourly measurements at 35 flux sites from the FLUXNET2015 dataset, and was also comparable with European Center for Medium-Range Weather Forecasts (ECMWF) ERA5 ET data, with an R2 of 0.617 and root mean square error (RMSE) of 65.8 W/m2 for the improved scheme. Results were largely improved compared with those without daily ET as the constraint. The improved method was stable for the estimation of ET’s diurnal cycle at the similar atmospheric conditions and the accuracy was comparative at different land cover surfaces. Errors in the input variables and the simplification of surface heat flux parameterization affected surface energy balance closure, which can lead to instability of the solution of constants in the simplified parameterization and further to the uncertainty of ET’s diurnal cycle estimation. Measurement errors, different source areas in measured variables, and inconsistent spatial representativeness between remote sensing and site measurements also impacted the evaluation.

Combining MWL and MSG SEVIRI Satellite Signals for Rainfall Detection and Estimation

Accurate rainfall detection and estimation are essential for many research and operational applications. Traditional rainfall detection and estimation techniques have achieved considerable success but with limitations. Thus, in this study, the relationships between the gauge (point measurement) and the microwave links (MWL) rainfall (line measurement), and the MWL to the satellite observations (area-wide measurement) are investigated for (area-wide) rainfall detection and rain rate retrieval. More precisely, we investigate if the combination of MWL with Meteosat Second Generation (MSG) satellite signals could improve rainfall detection and rainfall rate estimates. The investigated procedure includes an initial evaluation of the MWL rainfall estimates using gauge measurements, followed by a joint analysis of the rainfall estimates with the satellite signals by means of a conceptual model in which clouds with high cloud top optical thickness and large particle sizes have high rainfall probabilities and intensities. The analysis produced empirical thresholds that were used to test the capability of the MSG satellite data to detect rainfall on the MWL. The results from Kenya, during the “long rains” of 2013, 2014, and 2018 show convincing performance and reveal the potential of MWL and MSG data for area-wide rainfall detection.

Convective Initiation Proxies for Nowcasting Precipitation Severity Using the MSG-SEVIRI Rapid Scan

In this study, we investigate the ability of several convective initiation predictors based on satellite infrared observations to distinguish convective weak precipitation events from those leading to intense rainfall. The two types of precipitation are identified according to hourly rainfall, respectively less than 10 mm and greater than 30 mm. The analysis is conducted on a representative dataset containing 92 severe and weak precipitation events collected over the Italian peninsula in the period 2016–2019 over June-September. The events are selected to be short-lived (i.e., less than 12 h) and localized (i.e., less than 50×50km2). Italian National Radar Network products, namely the Vertical Maximum Intensity (VMI) and the Surface Rain Total (SRT) variables (from Dewetra Platform by CIMA, Italian Civil Protection Department), are used as indicators of convection (i.e., VMI greater than 35 dBZ echo intensity) and cumulated rainfall, respectively. The considered predictors are linear combinations of spectral infrared channels measured with the Rapid Scan Service (RSS) Spinning Enhanced Visible and InfraRed Imager (SEVIRI) aboard Meteosat Second Generation (MSG) geostationary satellites. We select a 5×5 SEVIRI pixel-box centered on the storm core and perform a statistical analysis of the predictors up to 2.5 h around the event occurrence. We demonstrate that some of the proxies—describing growth and glaciation storm properties—show few degrees contrast between severe and nonsevere precipitation cases, hence carrying significant information to help discriminate the two types. We design a threshold scheme based on the three most informative predictors to distinguish weak and strong precipitation events. This analysis yields accuracy higher than 0.6 and the probability of false detection lower than 0.26; in terms of reducing false alarms, this method shows slight better performances compared to related works, at the expense of a lower probability of detection. The overall results, however, show limited capability for these infrared proxies as stand-alone predictors to distinguish severe from nonsevere precipitation events. Nonetheless, these may serve as additional tools to reduce the false alarm ratio in nowcasting algorithms for convective orographic storms. This study also provides further insight into the correlation between early infrared fields signatures prior to convection and subsequent evolution of the storms, extending previous works in this field.

Reconstruction of daytime land surface temperatures under cloud-covered conditions using integrated MODIS/Terra land products and MSG geostationary satellite data

There is considerable demand for satellite observations that can support spatiotemporally continuous mapping of land surface temperature (LST) because of its strong relationships with many surface processes. However, the frequent occurrence of cloud cover induces a large blank area in current thermal infrared-based LST products. To effectively fill this blank area, a new method for reconstructing the cloud-covered LSTs of Terra Moderate Resolution Imaging Spectroradiometer (MODIS) daytime observations is described using random forest (RF) regression approach. The high temporal resolution of the Meteosat Second Generation (MSG) LST product assisted in identifying the temporal variations in cloud cover. The cumulative downward shortwave radiation flux (DSSF) was estimated as the solar radiation factor for each MODIS pixel based on the MSG DSSF product to represent the impact from cloud cover on incident solar radiation. The RF approach was used to fit an LST linking model based on the datasets collected from clear-sky pixels that depicted the complicated relationship between LST and the predictor variables, including the surface vegetation index (the normalized difference vegetation index and the enhanced vegetation index), normalized difference water index, solar radiation factor, surface albedo, surface elevation, surface slope, and latitude. The fitted model was then used to reconstruct the LSTs of cloud-covered pixels. The proposed method was applied to the Terra/MODIS daytime LST product for four days in 2015, spanning different seasons in southwestern Europe. A visual inspection indicated that the reconstructed LSTs thoroughly captured the distribution of surface temperature associated with surface vegetation cover, solar radiation, and topography. The reconstructed LSTs showed similar spatial pattern according to the comparison with clear-sky LSTs from temporally adjacent days. In addition, evaluations against Global Land Data Assimilation System (GLDAS) NOAH 0.25° 3-h LST data and reference LST data derived based on in-situ air temperature measurements showed that the reconstructed LSTs presented a stable and reliable performance. The coefficients of determination derived with the GLDAS LST data were all above 0.59 on the four examined days. These results indicate that the proposed method has a strong potential for reconstructing LSTs under cloud-covered conditions and can also accurately depict the spatial patterns of LST.

Investigating ground-based radar and spaceborne infrared radiometer synergy for lightning areal prediction in complex orography

A new multi-sensor approach, named PoLCast (Probability of Lightning foreCast), for predicting the lightning activity in a complex orography geographical area is proposed and discussed. The PoLCast input information are the ground-based weather radar horizontally polarized reflectivity factor and the atmospheric instability indexes, derived from the Spin Enhanced Visible Infrared Imager (SEVIRI) radiometer onboard the Meteosat Second Generation (MSG) satellite. The weather radar data are used to calculate the probability of lightning following a direct relation between the maximum values of the reflectivity factor and the lightning occurrence, whereas the atmospheric instability indexes from SEVIRI are used to constrain the probability of lighting, derived from weather radar data, and enhance such probability in cases of more unstable troposphere. To test the PoLCast methodology, the output of the Blitzortung and the “Sistema Italiano Rilevamento Fulmini” (SIRF) ground-based lightning sensor networks are used together with the C-band Mt. Midia weather radar within its 180 km diameter coverage over the Central Italy area. Both satellite and radar data are pre-processed into PoLCast to obtain a single time series in terms of the areal probability of lightning (PoL). PoLCast performances are evaluated in terms of statistical scores using 12 heterogeneous case studies over Central Italy. Even though the number of available cases is relatively limited, quantitative results show high areal PoL (from 0 to 100%) with a case-by-case variability of false alarm rate from 3 to 72%. The advantage of a multi-sensor technique, such as PoLCast, with respect to an approach using weather radar data only, becomes more evident when lightning activity is not present and the leading time of lightning forecast exceeds 2.5 h.

Near-Real-Time Application of SEVIRI Aerosol Optical Depth Algorithm

Within the framework of the Satellite-based Monitoring Initiative for Regional Air quality (SAMIRA) project, the near-real-time (NRT) operation has been documented for an in-house developed algorithm used for the retrieval of aerosol optical depth (AOD) maps from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) sensor onboard the Meteosat Second Generation (MSG). With the frequency of 15 min at a spatial resolution of roughly 5.5 × 5.5 km the AOD maps are provided for the country domains of Poland, the Czech Republic, Romania, and Southern Norway. A significant improvement has been reported in terms of modification of the existing prototype algorithm that it suits the operational NRT AOD retrieval for an extended area. This is mainly due to the application of the optimal interpolation method for the AOD estimation on reference days with the use of ground-based measurements of the Aerosol Robotic Network (AERONET) and the Aerosol Research Network (PolandAOD-NET) as well as simulations of the Copernicus Atmosphere Monitoring Service (CAMS). The main issues that have been addressed regarding surface reflectance estimation, cloud screening and uncertainty calculation. Exemplary maps of the NRT retrieval have been presented.

Near Real-Time Monitoring of the Christmas 2018 Etna Eruption Using SEVIRI and Products Validation

On the morning of 24 December 2018, an eruptive event occurred at Etna, which was followed the next day by a strong sequence of shallow earthquakes. The eruptive episode lasted until 30 December, ranging from moderate strombolian to lava fountain activity coupled with vigorous ash/gas emissions and a lava flow effusion toward the eastern volcano flank of Valle del Bove. In this work, the data collected from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) instruments on board the Meteosat Second Generation (MSG) geostationary satellite are used to characterize the Etna activity by estimating the proximal and distal eruption parameters in near real time. The inversion of data indicates the onset of eruption on 24 December at 11:15 UTC, a maximum Time Average Discharge Rate (TADR) of 8.3 m3/s, a cumulative lava volume emitted of 0.5 Mm3, and a Volcanic Plume Top Height (VPTH) that reached a maximum altitude of 8 km above sea level (asl). The volcanic cloud ash and SO2 result totally collocated, with an ash amount generally lower than SO2 except on 24 December during the climax phase. A total amount of about 100 and 35 kt of SO2 and ash respectively was emitted during the entire eruptive period, while the SO2 fluxes reached peaks of more than 600 kg/s, with a mean value of about 185 kg/s. The SEVIRI VPTH, ash/SO2 masses, and flux time series have been compared with the results obtained from the ground-based visible (VIS) cameras and FLux Automatic MEasurements (FLAME) networks, and the satellite images collected by the MODerate resolution Imaging Spectroradiometer (MODIS) instruments on board the Terra and Aqua- polar satellites. The analysis indicates good agreement between SEVIRI, VIS camera, and MODIS retrievals with VPTH, ash, and SO2 estimations all within measurement errors. The SEVIRI and FLAME SO2 flux retrievals show significant discrepancies due to the presence of volcanic ash and a gap of data on the FLAME network. The results obtained in this study show the ability of geostationary satellite systems to characterize eruptive events from the source to the atmosphere in near real time during the day and night, thus offering a powerful tool to mitigate volcanic risk on both local population and airspace and to give insight on volcanic processes.

A Combined IR-GPS Satellite Analysis for Potential Applications in Detecting and Predicting Lightning Activity

Continuous estimates of the vertical integrated precipitable water vapor content from the tropospheric delay of the signal received by the antennas of the global positioning system (GPS) are used in this paper, in conjunction with the measurements of the Meteosat Second Generation (MSG) spinning enhanced visible and infrared imager (SEVIRI) radiometer and with the lightning activity, collected here by the ground-based lightning detection network (LINET), in order to identify links and recurrent patterns useful for improving nowcasting applications. The analysis of a couple of events is shown here as an example of more general behavior. Clear signs appear before the peak of lightning activity on a timescale from 2 to 3 h. In particular, the lightning activity is generally preceded by a period in which the difference between SEVIRI brightness temperature (TB) at channel 5 and channel 6 (i.e., ∆TB) presents quite constant values around 0 K. This trend is accompanied by an increase in precipitable water vapor (PWV) values, reaching a maximum in conjunction with the major flash activity. The results shown in this paper evidence good potentials of using radiometer and GPS measurements together for predicting the abrupt intensification of lightning activity in nowcasting systems.

3D-VAR Data Assimilation of SEVIRI Radiances for the Prediction of Solar Irradiance in Italy Using WRF Solar Mesoscale Model—Preliminary Results

Solar power generation is highly fluctuating due to its dependence on atmospheric conditions. The integration of this variable resource into the energy supply system requires reliable predictions of the expected power production as a basis for management and operation strategies. This is one of the goals of the Solar Cloud project, funded by the Italian Ministry of Economic Development (MISE)—to provide detailed forecasts of solar irradiance variables to operators and organizations operating in the solar energy industry. The Institute of Methodologies for Environmental Analysis of the National Research Council (IMAA-CNR), participating to the project, implemented an operational chain that provides forecasts of all the solar irradiance variables at high temporal and horizontal resolution using the numerical weather prediction Advanced Research Weather Research and Forecasting (WRF-ARW) Solar version 3.8.1 released by the National Center for Atmospheric Research (NCAR) in August 2016. With the aim of improving the forecast of solar irradiance, the three-dimensional (3D-Var) data assimilation was tested to assimilate radiances from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) aboard the Meteosat Second Generation (MSG) geostationary satellite into WRF Solar. To quantify the impact, the model output is compared against observational data. Hourly Global Horizontal Irradiance (GHI) is compared with ground-based observations from Regional Agency for the Protection of the Environment (ARPA) and with MSG Shortwave Solar Irradiance estimations, while WRF Solar cloud coverage is compared with Cloud Mask by MSG. A preliminary test has been performed in clear sky conditions to assess the capability of the model to reproduce the diurnal cycle of the solar irradiance. The statistical scores for clear sky conditions show a positive performance of the model with values comparable to the instrument uncertainty and a correlation of 0.995. For cloudy sky, the solar irradiance and the cloud cover are better simulated when the SEVIRI radiances are assimilated, especially in the short range of the simulation. For the cloud cover, the Mean Bias Error one hour after the assimilation time is reduced from 41.62 to 20.29 W/m2 when the assimilation is activated. Although only two case studies are considered here, the results indicate that the assimilation of SEVIRI radiance improves the performance of WRF Solar especially in the first 3 hour forecast.

Near real-time estimation of Sea and Land surface temperature for MSG SEVIRI sensors

Land and Sea Surface Temperatures (LST and SST) are both recognized as Essential Climate Variables, and are routinely retrieved by a wealth of satellites. However, for validated approaches, the latest data are usually not available to the general public. We offer to bridge this gap, by using Meteosat Second Generation (MSG) Spinning Enhanced Visible and InfraRed Imager (SEVIRI), with its 15 min temporal resolution. Here, we present generic algorithms for the retrieval of both LST and SST, valid for the SEVIRI instrument onboard MSG platforms 8–11, which we validate using hourly data of 4 ground stations and 11 buoys in Spain over the years 2015 to 2018. These validations show that in the best conditions of surface homogeneity (cloud-free summer nights), errors in our LST estimation are below 1.5 K for stations with good thermal homogeneity. Comparison with LSA-SAF (Land Surface Analysis - Satellite Application Facility) LST shows differences below 2 K for most of SEVIRI disk, with higher differences in arid areas and during daytime. As for SST retrieval, the average error amount to 0.67 K for cloud-free buoy data. These algorithms have been implemented in a near-real time processing chain, which provide actualized LST and SST maps every 15 min within 5 min of image reception. These maps, along with other products, can be freely consulted from a dedicated webpage (https://www.uv.es/iplsat).

A Sea Surface Temperature data record (2004–2012) from Meteosat Second Generation satellites

The Ocean and Sea-Ice Satellite Application Facility (OSI-SAF) of the EUropean Organisation for the Exploitation of METeorological Satellites (EUMETSAT) has performed a reprocessing of Sea Surface Temperature (SST) from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard Meteosat Second Generation (MSG) archive (2004–2012).The retrieval method consists of a non-linear split-window algorithm and an algorithm correction relying on simulations of infrared brightness temperatures performed using atmospheric profiles of water vapour and temperature from a Numerical Weather Prediction model, and a radiative transfer model. The cloud mask used is the Climate SAF reprocessing of the MSG/SEVIRI archive which is consistent over the period considered. Atmospheric Saharan dust has a strong impact on the retrieved SST in the Atlantic and Mediterranean regions, they are taken into consideration through the computation of the Saharan Dust Index (SDI) which is then used to determine an empirical correction applied to SST.The reprocessing has benefited from the experience of the OSI SAF team in operational near real time processing of MSG/SEVIRI data, and the methods have been improved to provide a higher quality SST. The MSG/SEVIRI SST reprocessing dataset consists of hourly level 3 composites of sub-skin temperature projected onto a regular 0.05° grid over the region delimited by 60N,60S and 60W,60E. It has been thoroughly validated against drifting buoys and moored buoys measurements.Results of this validation have shown that the reprocessed data record is of significantly better quality than the OSI SAF operational processing (for instance the day-time robust standard deviation is 0.45 K for the operational processing and 0.35 K for the reprocessed dataset). The data record has been used to characterize the diurnal variability of SST over large temporal and spatial scales.

Cloud Detection Based on High Resolution Stereo Pairs of the Geostationary Meteosat Images

Due to the considerable impact of clouds on the energy balance in the atmosphere and on the earth surface, they are of great importance for various applications in meteorology or remote sensing. An important aspect of the cloud research studies is the detection of cloudy pixels from the processing of satellite images. In this research, we investigated a stereographic method on a new set of Meteosat images, namely the combination of the high resolution visible (HRV) channel of the Meteosat-8 Indian Ocean Data Coverage (IODC) as a stereo pair with the HRV channel of the Meteosat Second Generation (MSG) Meteosat-10 image at 0° E. In addition, an approach based on the outputs from stereo analysis was proposed to detect cloudy pixels. This approach is introduced with a 2D-scatterplot based on the parallax value and the minimum intersection distance. The mentioned scatterplot was applied to determine/detect cloudy pixels in various image subsets with different amounts of cloud cover. Apart from the general advantage of the applied stereography method, which only depends on geometric relationships, the cloud detection results are also improved because: (1) The stereo pair is the HRV bands of the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) sensor, with the highest spatial resolution available from the Meteosat geostationary platform; and (2) the time difference between the image pairs is nearly 5 s, which improves the matching results and also decreases the effect of cloud movements. In order to prove this improvement, the results of this stereo-based approach were compared with three different reflectance-based target detection techniques, including the adaptive coherent estimator (ACE), constrained energy minimization (CEM), and matched filter (MF). The comparison of the receiver operating characteristics (ROC) detection curves and the area under these curves (AUC) showed better detection results with the proposed method. The AUC value was 0.79, 0.90, 0.90, and 0.93 respectively for ACE, CEM, MF, and the proposed stereo-based detection approach. The results of this research shall enable a more realistic modelling of down-welling solar irradiance in the future.