LSA SAF Research Articles

Evaluating a multi-step collocation approach for an ensemble climatological dataset of actual evapotranspiration over Italy

Accurate estimations of actual evapotranspiration (ET) are key in a variety of water balance applications, but divergent results can be obtained due to the large range of available methodologies. The use of an ensemble approach is a suitable alternative, as it summarizes multiple sources in an optimized strategy. In this study, an expert-based multi-step collocation (MC) approach is tested to merge six ET datasets, with the aim of reconstructing a spatiotemporally-consistent monthly dataset for Italy in the climatological period 1991–2020. The merged products are: three water balance datasets (BIG BANG, LSA SAF, and LISFLOOD), two residual surface energy balance model datasets (SSEBop, and ALEXI), and the MODIS standard ET product. The merged product is analyzed for spatio-temporal consistency and evaluated using flux observations from 11 sites. On average, the merged product has higher accuracy (mean absolute difference = 0.47 ± 0.17 mm/d, relative difference = 27.9 ± 7.5 %) than any single base dataset, and it is characterized by limited bias (mean bias error = -0.17 ± 0.26 mm/d), high correlation (r = 0.83 ± 0.10), and more uniform performance across sites. The merged ET dataset is accompanied by an estimation of the ensemble spread, which highlights large differences in ET estimates in some areas and periods characterized by severe water stress, such as in southern Italy during the summer. This large spread seems to be mostly driven by systematic differences among datasets, which affect the estimation of the reference climatology, suggesting how inter-model spread can have a defining role in further improving the merging strategies.

Evaluation of Downward Surface Longwave Flux Estimates Using Meteosat Cloud Observations

AbstractThe downward surface longwave flux (DSLF) plays a relevant role in the Earth’s surface radiative budget, which is crucial to monitor, understand and model the impact of changes at local and global scales on surface temperature and surface conditions. This study focuses on the evaluation and intercomparison of four DSLF products: (a) a recently developed all‐weather DSLF product based on the multivariate adaptive regression splines (MARS) algorithm driven by satellite cloud information from the Meteosat Second Generation (MSG) and ERA5 reanalysis screen variables; (b) the Satellite Application Facility on Land Surface Analysis (LSA SAF); (c) CERES Synoptic top‐of‐atmosphere and surface fluxes and clouds (CERES‐SYN1deg) and (d) ERA5 reanalysis. The study covers the period 2005–2021 and the MSG region focusing on monthly means. The evaluation performed against 48 ground stations from the Baseline Surface Radiation Network (BSRN) and FLUXNET2015 networks showed that the MARS product outperforms the remaining products, particularly the LSA SAF, while ERA5 and CERES show similar performance metrics. The fours products are intercompared in terms of their mean spatial variability and temporal mean annual cycles and inter‐annual variability in four selected regions, showing a high level of agreement, particularly between MARS, ERA5 and CERES. Our results highlight the clear added value of MARS with respect to LSA SAF, while providing higher spatial resolution (0.05°), constrained by satellite cloud information, when compared with ERA5 (0.25°) or CERES (1°).

Upgrade and extension of LSA-SAF land surface albedo archive from EPS Metop/AVHRR: description and quality assessment

ABSTRACT ETAL is the operational EPS Ten-Day Albedo product, produced by the EUMETSAT Satellite Application Facility for Land Surface Analysis (LSA SAF). By back-processing the full catalogue of EPS-Metop radiance data from September 2007 to June 2021, we are able to 1) extend the temporal coverage (previously the archive only went back to 2015) and 2) improve the product archive that was based on near-real time (NRT) processing; the second point is achieved by using reanalyses instead of forecasts of atmospheric conditions and by not being exposed to missing data in the NRT radiance inputs. We present this reprocessed part of the ETAL data set, called ETAL-R, and assess its quality and consistency with respect to the original archive of NRT ETAL data (for the overlapping period 2015–2021), as well as its accuracy compared to albedo from MODIS and ground stations. ETAL-R exhibits reliable long-term stability and increased homogeneity compared to the NRT archive, and the comparison against the additional reference data shows satisfactory accuracy. Overall, ETAL-R is shown to be very consistent with the ETAL NRT archive while – under certain circumstances – improving it.

Evapotranspiration and surface energy fluxes across Europe, Africa and Eastern South America throughout the operational life of the Meteosat second generation satellite

AbstractThe exchange of energy and water fluxes between the Earth's surface and the atmosphere is crucial to a series of processes that impact human life. Noteworthy examples are agriculture yields, water availability, intensity and extent of droughts and the ability of ecosystems to provide services to society. The relevance of these processes has motivated the Satellite Application Facility on Land Surface Analysis (LSA SAF) programme to set up an operational framework to estimate—among other variables—evapotranspiration (ET) and surface energy fluxes (SEF) on the basis of observations by the Meteosat Second Generation (MSG) satellite. The LSA SAF programme has recently launched the reprocessing of the ET and SEF datasets on the basis of the most recent version of the algorithm and homogenous forcing datasets. This article features the resulting ET/SEF dataset, a Data Record that encompasses the period from the start of the operational life of the MSG satellite (2004) till 2020 and covers the field of view of the MSG satellite (i.e. Europe, Africa and Eastern South America). Details on the algorithm and the datasets driving the ET/SEF estimates are also provided as well as a quality assessment.

Biomass burning CO, PM and fuel consumption per unit burned area estimates derived across Africa using geostationary SEVIRI fire radiative power and Sentinel-5P CO data

Abstract. We present the first top-down CO fire emissions inventory for Africa based on the direct relation between geostationary satellite-based fire radiative power (FRP) observations and polar-orbiting satellite observations of total column carbon monoxide (TCCO). This work significantly extends the previous Fire Radiative Energy Emissions (FREM) approach that derived total particulate matter (TPM) emission coefficients from FRP and aerosol optical depth (AOD) observations. The use of satellite-based CO observations to derive biome-specific CO emission coefficients, ECCOb, addresses key uncertainties in the use of AOD observations to estimate fire-generated CO emissions including the requirement for a smoke mass extinction coefficient in the AOD to TPM conversion and the large variation in TPM emission factors – which are used to convert TPM emissions to CO emissions. We use the FREM-derived CO emission coefficients to produce a pan-African CO fire emission inventory spanning 2004 to 2019. Regional CO emissions are in close agreement with the most recent version of GFED(v4.1s), despite the two inventories using completely different satellite datasets and methodologies. Dry matter consumed (DMC) and DMC per unit burned area are generated from our CO emission inventory – the latter using the 20 m resolution Sentinel-2 FireCCISFD burnt area (BA) product for 2019. We carry out an evaluation of our FREM-based CO emissions by using them as input in the WRF-CMAQ chemical transport model and comparing simulated TCCO fields to independent Sentinel-5P TROPOMI TCCO observations. The results of this evaluation show FREM CO emissions to generally be in good agreement with these independent measures – particularly in the case of individual fire-generated CO plumes, where modelled in-plume CO was within 5 % of satellite observations with a coefficient of determination of 0.80. Modelled and observed total CO, aggregated over the full model domain, are within 4 % of each other, though localised regions show an overestimation of modelled CO by up to 50 %. When compared to other evaluations of current state-of-the-art fire emissions inventories, the FREM CO emission inventory derived in this work shows some of the best agreement with independent observations. Updates to previously published FREM TPM emissions coefficients based on this methodology are also provided, along with a similar evaluation as conducted for CO. The methodology described in this work is forming the basis of a forthcoming near-real-time fire emissions product from Meteosat to be issued by the EUMETSAT LSA SAF (https://landsaf.ipma.pt/en/, last access: 19 December 2022).

Climatological Drought Monitoring in Switzerland Using EUMETSAT SAF Satellite Data

Climatological drought monitoring in Switzerland relies heavily on station-based precipitation and temperature data. Due to the high spatial variability and complexity of droughts, it is important to complement station-based drought indices with gridded information and to couple multiple drought indicators within the monitoring system. Here, long-term satellite-based drought parameters from the EUMETSAT SAF network are analyzed in terms of dry anomalies within their climatology’s, namely ASCAT soil water index (SWI), CM SAF land surface temperature (LST), complemented with NOAA vegetation data, and LSA SAF Meteosat evapotranspiration data. The upcoming EUMETSAT SAF climate data records on land surface temperature and evapotranspiration will cover for the first time the WMO climatological 30-year reference period. This study is the first study investigating the potential of those long-term data records for climate monitoring of droughts in Europe. The satellite datasets are compared with the standardized precipitation index (SPI), soil moisture observations from the SwissSMEX measurement network, with a modelled soil moisture index (SMI) based on observations, and with evapotranspiration measurements, focusing on the temporal dynamics of the anomalies. For vegetation and surface temperature, the dry years of 2003, 2015, and 2018 are clearly visible in the satellite data. CM SAF LSTs show strong anomalies at the beginning of the drought period. The comparison of in situ and modelled soil moisture and evapotranspiration measurements with the satellite parameters shows strong agreement in terms of anomalies. The SWI indicates high anomaly correlations of 0.56 to 0.83 with measurements and 0.63 to 0.76 with the SMI at grassland sites. The Meteosat evapotranspiration data strongly agree with the measurements, with anomaly correlations of 0.63 and 0.67 for potential and actual evapotranspiration, respectively. Due to the prevailing humid climate conditions at the considered sites, evapotranspiration anomalies during the investigated dry periods were mostly positive and thus not water limited, but were also a driver for soil moisture drought. The results indicate that EUMETSAT SAF satellite data can well complement the station-based drought monitoring in Switzerland with spatial information.

Local-scale evaluation of the simulated interactions between energy, water and vegetation in ISBA, ORCHIDEE and a diagnostic model

Abstract. The processes involved in the exchange of water, energy and carbon in terrestrial ecosystems are strongly intertwined. To accurately represent the terrestrial biosphere in land surface models (LSMs), the intrinsic coupling between these processes is required. Soil moisture and leaf area index (LAI) are two key variables at the nexus of water, energy and vegetation. Here, we evaluated two prognostic LSMs (ISBA and ORCHIDEE) and a diagnostic model (based on the LSA SAF, Satellite Application Facility for Land Surface Analysis, algorithms) in their ability to simulate the latent heat flux (LE) and gross primary production (GPP) coherently and their interactions through LAI and soil moisture. The models were validated using in situ eddy covariance observations, soil moisture measurements and remote-sensing-based LAI. It was found that the diagnostic model performed consistently well, regardless of land cover, whereas important shortcomings of the prognostic models were revealed for herbaceous and dry sites. Despite their different architecture and parametrization, ISBA and ORCHIDEE shared some key weaknesses. In both models, LE and GPP were found to be oversensitive to drought stress. Though the simulated soil water dynamics could be improved, this was not the main cause of errors in the surface fluxes. Instead, these errors were strongly correlated to errors in LAI. The simulated phenological cycle in ISBA and ORCHIDEE was delayed compared to observations and failed to capture the observed seasonal variability. The feedback mechanism between GPP and LAI (i.e. the biomass allocation scheme) was identified as a key element to improve the intricate coupling between energy, water and vegetation in LSMs.

Surface Albedo Retrieval from 40-Years of Earth Observations through the EUMETSAT/LSA SAF and EU/C3S Programmes: The Versatile Algorithm of PYALUS

Land surface albedo quantifies the fraction of the sunlight reflected by the surface of the Earth. This article presents the algorithm concepts for the remote sensing of this variable based on the heritage of several developments which were performed at Méteo France over the last decade and described in several papers by Carrer et al. The scientific algorithm comprises four steps: an atmospheric correction, a sensor harmonisation (optional), a BRDF (Bidirectional Reflectance Distribution Function) inversion, and the albedo calculation. At the time being, the method has been applied to 11 sensors in the framework of two European initiatives (Satellite Application Facility on Land Surface Analysis—LSA SAF, and Copernicus Climate Change Service—C3S): NOAA-7-9-11-14-16-17/AVHRR2-3, SPOT/VGT1-2, Metop/AVHRR-3, PROBA-V, and MSG/SEVIRI. This work leads to a consistent archive of almost 40 years of satellite-derived albedo data (available in 2020). From a single sensor, up to three different albedo products with different characteristics have been developed to address the requirements of both, near real-time (NRT) (weather prediction with a demand of timeliness of 1 h) and climate communities. The evaluation of the algorithm applied to different platforms was recently made by Lellouch et al. and Sánchez Zapero et al. in 2020 which can be considered as companion papers. After a summary of the method for the retrieval of these surface albedos, this article describes the specificities of each retrieval, lists the differences, and discusses the limitations. The plan of continuity with the next European satellite missions and perspectives of improvements are introduced. For example, Metop/AVHRR-3 albedo will soon become the medium resolution sensor product with the longest NRT data record, since MODIS is approaching the end of its life-cycle. Additionally, Metop-SG/METimage will ensure its continuity thanks to consistent production of data sets guaranteed till 2050 by the member states of the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT). In the end, the common strategy which we proposed through the different programmes may offer an unprecedented opportunity to study the temporal trends affecting surface properties and to analyse human-induced climate change. Finally, the access to the source code (called PYALUS) is provided through an open access platform in order to share with the community the expertise on the satellite retrieval of this variable.

Improved estimation of global solar radiation over rugged terrains by the disaggregation of Satellite Applications Facility on Land Surface Analysis data (LSA SAF)

AbstractThis paper presents a new method to predict global solar radiation over irregular terrain, named Estimation of global solar RADiation (ERAD). The method is based on the disaggregation of Satellite Applications Facility on Land Surface Analysis (LSA SAF) data using a digital elevation model and is applied in Italy with a time step of 1 min and a spatial resolution of 200 m. A quantitative assessment of ERAD is performed in comparison with three other standard methods (Mountain Microclimate Simulation Model [MTCLIM], LSA SAF and Copernicus Atmosphere Monitoring Service [CAMS]) using measurements taken in 43 stations located in Italy or in the surrounding countries, in the years 2005–2016. Such assessment concerns the irradiance incoming on a horizontal surface, which is measured by ground radiation sensors and is summarized by means of four accuracy statistics (i.e. mean absolute error [MAE], root mean square error [RMSE], coefficient of determination [R2] and mean bias error [MBE]). Overall, the average daily global solar radiation estimates obtained by ERAD have RMSE and R2 about 25 W·m−2 and 0.943, respectively. These statistics are similar to those of LSA SAF and better than those of CAMS and, above all, MTCLIM. The bias analysis by elevation ranges shows a slight ERAD overestimation over plains and hills and a slight underestimation over mountains. An additional qualitative assessment shows how the ERAD radiation estimates are more spatially detailed than those of the other methods and are redistributed on inclined surfaces consistently with expectations.

Satellite Retrieval of Downwelling Shortwave Surface Flux and Diffuse Fraction under All Sky Conditions in the Framework of the LSA SAF Program (Part 2: Evaluation)

High frequency knowledge of the spatio-temporal distribution of the downwelling surface shortwave flux (DSSF) and its diffuse fraction (fd) at the surface is nowadays essential for understanding climate processes at the surface–atmosphere interface, plant photosynthesis and carbon cycle, and for the solar energy sector. The European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility for Land Surface Analysis operationally delivers estimation of the MDSSFTD (MSG Downwelling Surface Short-wave radiation Fluxes—Total and Diffuse fraction) product with an operational status since the year 2019. The method for retrieval was presented in a companion paper. Part 2 now focuses on the evaluation of the MDSSFTD algorithm and presents a comparison of the corresponding outputs, i.e., total DSSF and diffuse fraction (fd) components, against in situ measurements acquired at four Baseline Surface Radiation Network (BSRN) stations over a seven-month period. The validation is performed on an instantaneous basis. We show that the satellite estimates of DSSF and fd meet the target requirements defined by the user community for all-sky (clear and cloudy) conditions. For DSSF, the requirements are 20 Wm−2 for DSSF < 200 Wm−2, and 10% for DSSF ≥ 200 Wm−2. The mean bias error (MBE) and relative mean bias error (rMBE) compared to the ground measurements are 3.618 Wm−2 and 0.252%, respectively. For fd, the requirements are 0.1 for fd < 0.5, and 20% for fd ≥ 0.5. The MBE and rMBE compared to the ground measurements are −0.044% and −17.699%, respectively. The study also provides a separate analysis of the product performances for clear sky and cloudy sky conditions. The importance of representing the cloud–aerosol radiative coupling in the MDSSFTD method is discussed. Finally, it is concluded that the quality of the aerosol optical depth (AOD) forecasts currently available is accurate enough to obtain reliable diffuse solar flux estimates. This quality of AOD forecasts was still a limitation a few years ago.

Satellite Retrieval of Downwelling Shortwave Surface Flux and Diffuse Fraction under All Sky Conditions in the Framework of the LSA SAF Program (Part 1: Methodology)

Several studies have shown that changes in incoming solar radiation and variations of the diffuse fraction can significantly modify the vegetation carbon uptake. Hence, monitoring the incoming solar radiation at large scale and with high temporal frequency is crucial for this reason along with many others. The European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility for Land Surface Analysis (LSA SAF) has operationally disseminated in near real time estimates of the downwelling shortwave radiation at the surface since 2005. This product is derived from observations provided by the SEVIRI instrument onboard the Meteosat Second Generation series of geostationary satellites, which covers Europe, Africa, the Middle East, and part of South America. However, near real time generation of the diffuse fraction at the surface level has only recently been initiated. The main difficulty towards achieving this goal was the general lack of accurate information on the aerosol particles in the atmosphere. This limitation is less important nowadays thanks to the improvements in atmospheric numerical models. This study presents an upgrade of the LSA SAF operational retrieval method, which provides the simultaneous estimation of the incoming solar radiation and its diffuse fraction from the satellite every 15 min. The upgrade includes a comprehensive representation of the influence of aerosols based on physical approximations of the radiative transfer within an atmosphere-surface associated medium. This article explains the retrieval method, discusses its limitations and differences with the previous method, and details the characteristics of the output products. A companion article will focus on the evaluation of the products against independent measurements of solar radiation. Finally, the access to the source code is provided through an open access platform in order to share the expertise on the satellite retrieval of this variable with the community.

Climate Data Records of Vegetation Variables from Geostationary SEVIRI/MSG Data: Products, Algorithms and Applications

The scientific community requires long-term data records with well-characterized uncertainty and suitable for modeling terrestrial ecosystems and energy cycles at regional and global scales. This paper presents the methodology currently developed in EUMETSAT within its Satellite Application Facility for Land Surface Analysis (LSA SAF) to generate biophysical variables from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) on board MSG 1-4 (Meteosat 8-11) geostationary satellites. Using this methodology, the LSA SAF generates and disseminates at a time a suite of vegetation products, such as the leaf area index (LAI), the fraction of the photosynthetically active radiation absorbed by vegetation (FAPAR) and the fractional vegetation cover (FVC), for the whole Meteosat disk at two temporal frequencies, daily and 10-days. The FVC algorithm relies on a novel stochastic spectral mixture model which addresses the variability of soils and vegetation types using statistical distributions whereas the LAI and FAPAR algorithms use statistical relationships general enough for global applications. An overview of the LSA SAF SEVIRI/MSG vegetation products, including expert knowledge and quality assessment of its internal consistency is provided. The climate data record (CDR) is freely available in the LSA SAF, offering more than fifteen years (2004-present) of homogeneous time series required for climate and environmental applications. The high frequency and good temporal continuity of SEVIRI products addresses the needs of near-real-time users and are also suitable for long-term monitoring of land surface variables. The study also evaluates the potential of the SEVIRI/MSG vegetation products for environmental applications, spanning from accurate monitoring of vegetation cycles to resolving long-term changes of vegetation.

Fire Activity and Fuel Consumption Dynamics in Sub-Saharan Africa

African landscape fires are widespread, recurrent and temporally dynamic. They burn large areas of the continent, modifying land surface properties and significantly affect the atmosphere. Satellite Earth Observation (EO) data play a pivotal role in capturing the spatial and temporal variability of African biomass burning, and provide the key data required to develop fire emissions inventories. Active fire observations of fire radiative power (FRP, MW) have been shown to be linearly related to rates of biomass combustion (kg s−1). The Meteosat FRP-PIXEL product, delivered in near real-time by the EUMETSAT Land Surface Analysis Satellite Applications Facility (LSA SAF), maps FRP at 3 km resolution and 15-min intervals and these data extend back to 2004. Here we use this information to assess spatio-temporal variations in fire activity across sub-Saharan Africa, and identify an overall trend of decreasing annual fire activity and fuel consumption, agreeing with the widely-used Global Fire Emissions Database (GFEDv4) based on burned area measures. We provide the first comprehensive assessment of relationships between per-fire FRE-derived fuel consumption (Tg dry matter, DM) and temporally integrated Moderate Resolution Imaging Spectroradiometer (MODIS) net photosynthesis (PSN) (Tg, which can be converted into pre-fire fuel load estimates). We find very strong linear relationships over southern hemisphere Africa (mean r = 0.96) that are partly biome dependent, though the FRE-derived fuel consumptions are far lower than those derived from the accumulated PSN, with mean fuel consumptions per unit area calculated as 0.14 kg DM m−2. In the northern hemisphere, FRE-derived fuel consumption is also far lower and characterized by a weaker linear relationship (mean r = 0.76). Differences in the parameterization of the biome look up table (BLUT) used by the MOD17 product over Northern Africa may be responsible but further research is required to reconcile these differences. The strong relationship between fire FRE and pre-fire fuel load in southern hemisphere Africa is encouraging and highlights the value of geostationary FRP retrievals in providing a metric that relates very well to fuel consumption and fire emission variations. The fact that the estimated fuel consumed is only a small fraction of the fuel available suggests underestimation of FRE by Spinning Enhanced Visible and Infrared Imager (SEVIRI) and/or that the FRE-to-fuel consumption conversion factor of 0.37 MJ kg−1 needs to be adjusted for application to SEVIRI. Future geostationary imaging sensors, such as on the forthcoming Meteosat Third Generation (MTG), will reduce the impact of this underestimation through its ability to detect even smaller and shorter-lived fires than can the current second generation Meteosat.

Land Surface Albedo Derived on a Ten Daily Basis from Meteosat Second Generation Observations: The NRT and Climate Data Record Collections from the EUMETSAT LSA SAF

Land surface albedo determines the splitting of downwelling solar radiation into components which are either reflected back to the atmosphere or absorbed by the surface. Land surface albedo is an important variable for the climate community, and therefore was defined by the Global Climate Observing System (GCOS) as an Essential Climate Variable (ECV). Within the scope of the Satellite Application Facility for Land Surface Analysis (LSA SAF) of EUMETSAT (European Organization for the Exploitation of Meteorological Satellites), a near-real time (NRT) daily albedo product was developed in the last decade from observations provided by the Spinning Enhanced Visible and Infrared Imager (SEVIRI) instrument on board the geostationary satellites of the Meteosat Second Generation (MSG) series. In this study we present a new collection of albedo satellite products based on the same satellite data. The MSG Ten-day Albedo (MTAL) product incorporates MSG observations over 31 days with a frequency of NRT production of 10 days. The MTAL collection is more dedicated to climate analysis studies compared to the daily albedo that was initially designed for the weather prediction community. For this reason, a homogeneous reprocessing of MTAL was done in 2018 to generate a climate data record (CDR). The resulting product is called MTAL-R and has been made available to the community in addition to the NRT version of the MTAL product which has been available for several years. The retrieval algorithm behind the MTAL products comprises three distinct modules: One for atmospheric correction, one for daily inversion of a semi-empirical model of the bidirectional reflectance distribution function, and one for monthly composition, that also determines surface albedo values. In this study the MTAL-R CDR is compared to ground surface measurements and concomitant albedo products collected by sensors on-board polar-orbiting satellites (SPOT-VGT and MODIS). We show that MTAL-R meets the quality requirements if MODIS or SPOT-VGT are considered as reference. This work leads to 14 years of production of geostationary land surface albedo products with a guaranteed continuity in the LSA SAF for the future years with the forthcoming third generation of European geostationary satellites.

Validation of reference evapotranspiration from Meteosat Second Generation (MSG) observations

This article presents validation results of the grass Reference Evapotranspiration (ETo) product, which is provided operationally by the Satellite Applications Facility on Land Surface Analysis (LSA SAF). ETo is considered to be the evapotranspiration from a hypothetical extensive well-watered field covered with green grass (12 cm high and with 0.23 albedo), under the given down-welling short-wave radiation. The LSA SAF product is estimated from daily solar radiation at the surface via a methodology designed to be applicable to that well-defined reference surface. In line with its definition, the LSA SAF ETo product is based on estimates of the radiative energy available at the surface.Validation of ETo estimates is challenging, given the difficulties of finding reliable in situ measurements that meet the ETo definition. We show that for the site that matches closely the reference surface (Cabauw, The Netherlands), the LSA SAF product outperforms other commonly used methodologies, even when these use in situ observations as input. However, observations taken at other stations, located in areas that do not deviate greatly from the reference surface, put into evidence the high uncertainty in local measurements. It is shown that the LSA SAF ETo product follows well in situ values, with average differences of 0.3 mm/day or less for mid-latitude sites and of the order of 1 mm/day for the Spanish sites near Cordoba and Albacete. Local advection effects cannot be ignored in measurements performed in the latter stations, where summer conditions are mostly warm and dry. In situ ETo measurements performed in those cases with lysimeters within limited fields are higher than satellite estimates due to advection of warm dry air from the vicinity, which acts as an extra source of energy. It is shown that this effect can be parameterized as a function of near surface air temperature. However, it is argued that those local advection effects should not occur in the idealized surface referred above.In contrast to other ETo estimates, the LSA SAF product is not influenced by local aridity or advection effects, and therefore it is particularly appropriate for large scale climate assessments, including drought monitoring (e.g. by considering the ratio of actual and reference evapotranspiration). Additionally, it provides suitable estimates of irrigation requirements in support of water management.

Fire danger rating over Mediterranean Europe based on fire radiative power derived from Meteosat

Abstract. We present a procedure that allows the operational generation of daily forecasts of fire danger over Mediterranean Europe. The procedure combines historical information about radiative energy released by fire events with daily meteorological forecasts, as provided by the Satellite Application Facility for Land Surface Analysis (LSA SAF) and the European Centre for Medium-Range Weather Forecasts (ECMWF). Fire danger is estimated based on daily probabilities of exceedance of daily energy released by fires occurring at the pixel level. Daily probability considers meteorological factors by means of the Canadian Fire Weather Index (FWI) and is estimated using a daily model based on a generalized Pareto distribution. Five classes of fire danger are then associated with daily probability estimated by the daily model. The model is calibrated using 13 years of data (2004–2016) and validated against the period of January–September 2017. Results obtained show that about 72 % of events releasing daily energy above 10 000 GJ belong to the “extreme” class of fire danger, a considerably high fraction that is more than 1.5 times the values obtained when using the currently operational Fire Danger Forecast module of the European Forest Fire Information System (EFFIS) or the Fire Risk Map (FRM) product disseminated by the LSA SAF. Besides assisting in wildfire management, the procedure is expected to help in decision making on prescribed burning within the framework of agricultural and forest management practices.

Retrieval of daily gross primary production over Europe and Africa from an ensemble of SEVIRI/MSG products

The main goal of this paper is to derive a method for a daily gross primary production (GPP) product over Europe and Africa taking the full advantage of the SEVIRI/MSG satellite products from the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) sensors delivered from the Satellite Application Facility for Land Surface Analysis (LSA SAF) system. Special attention is paid to model the daily GPP response from an optimized Montheith's light use efficiency model under dry conditions by controlling water shortage limitations from the actual evapotranspiration and the potential evapotranspiration (PET). The PET was parameterized using the mean daily air temperature at 2m (Ta) from ERA-Interim data. The GPP product (MSG GPP) was produced for 2012 and assessed by direct site-level comparison with GPP from eddy covariance data (EC GPP). MSG GPP presents relative bias errors lower than 40% for the most forest vegetation types with a high agreement (r>0.7) when compared with EC GPP. For drylands, MSG GPP reproduces the seasonal variations related to water limitation in a good agreement with site level GPP estimates (RMSE=2.11gm−2day−1; MBE=−0.63gm−2day−1), especially for the dry season. A consistency analysis against other GPP satellite products (MOD17A2 and FLUXCOM) reveals a high consistency among products (RMSD<1.5gm−2day−1) over Europe, North and South Africa. The major GPP disagreement arises over moist biomes in central Africa (RMSD>3.0gm−2day−1) and over dry biomes with MSG GPP estimates lower than FLUXCOM (MBD up to −3.0gm−2day−1). This newly derived product has the potential for analysing spatial patterns and temporal dynamics of GPP at the MSG spatial resolutions on a daily basis allowing to better capture the GPP dynamics and magnitude.

Improving Land Surface Temperature Retrievals over Mountainous Regions

Algorithms for Land Surface Temperature (LST) retrieval from infrared measurements are usually sensitive to the amount of water vapor present in the atmosphere. The Satellite Application Facilities on Climate Monitoring and Land Surface Analysis (CM SAF and LSA SAF) are currently compiling a 25 year LST Climate data record (CDR), which uses water vapor information from ERA-Int reanalysis. However, its relatively coarse spatial resolution may lead to systematic errors in the humidity profiles with implications in LST, particularly over mountainous areas. The present study compares LST estimated with three different retrieval algorithms: a radiative transfer-based physical mono-window (PMW), a statistical mono-window (SMW), and a generalized split-windows (GSW). The algorithms were tested over the Alpine region using ERA-Int reanalysis data and relied on the finer spatial scale Consortium for Small-Scale Modelling (COSMO) model data as a reference. Two methods were developed to correct ERA-Int water vapor misestimation: (1) an exponential parametrization of total precipitable water (TPW) appropriate for SMW/GSW; and (2) a level reduction method to be used in PMW. When ERA-Int TPW was used, the algorithm missed the right TPW class in 87% of the cases. When the exponential parametrization was used, the missing class rate decreased to 9%, and when the level reduction method was applied, the LST corrections went up to 1.7 K over the study region. Overall, the correction for pixel orography in TPW leads to corrections in LST estimations, which are relevant to ensure that long-term LST records meet climate requirements, particularly over mountainous regions.

Assessing reference evapotranspiration at regional scale based on remote sensing, weather forecast and GIS tools

Reference evapotranspiration (ETo) is a key component in efficient water management, especially in arid and semi-arid environments. However, accurate ETo assessment at the regional scale is complicated by the limited number of weather stations and the strict requirements in terms of their location and surrounding physical conditions for the collection of valid weather data. In an attempt to overcome this limitation, new approaches based on the use of remote sensing techniques and weather forecast tools have been proposed.Use of the Land Surface Analysis Satellite Application Facility (LSA SAF) tool and Geographic Information Systems (GIS) have allowed the design and development of innovative approaches for ETo assessment, which are especially useful for areas lacking available weather data from weather stations. Thus, by identifying the best-performing interpolation approaches (such as the Thin Plate Splines, TPS) and by developing new approaches (such as the use of data from the most similar weather station, TS, or spatially distributed correction factors, CITS), errors as low as 1.1% were achieved for ETo assessment. Spatial and temporal analyses reveal that the generated errors were smaller during spring and summer as well as in homogenous topographic areas.The proposed approaches not only enabled accurate calculations of seasonal and daily ETo values, but also contributed to the development of a useful methodology for evaluating the optimum number of weather stations to be integrated into a weather station network and the appropriateness of their locations. In addition to ETo, other variables included in weather forecast datasets (such as temperature or rainfall) could be evaluated using the same innovative methodology proposed in this study.

SEVIRI PrePro: A novel software tool for the pre-processing of SEVIRI geostationary orbit EO data products

The Spinning Enhanced Visible and Infrared Imager (SEVIRI) is a geostationary orbit multispectral sensor on-board the Meteosat second Generation (MSG) platform, acquiring Earth Observation (EO) data over Earth's land surface from the optical to infrared parts of electromagnetic spectrum every 15 min. From the sensor a series of operational products are also provided to the user's community at no cost via EUMETSAT or LSA SAF portals.Herein, an open access stand-alone software product developed in Java programming language is presented for automating key pre-processing steps to all the SEVIRI operationally distributed products. The software tool, named Seviri PrePro, makes use of present day multi-core processors and is able to process very large datasets in a short time period, making it appropriate as well for use in a High Performance Computing (HPC) environment. The practical usefulness of the toolkit is also demonstrated herein using as a case study the SEVIRI evapotranspiration (ET) product.The development of SEVIRI PrePro is of significant importance to the SEVIRI users' community and is also very timely given that, to our knowledge, no similar software tool is freely distributed at present. Its use is anticipated to make a significant contribution to a large number of practical applications requiring use of SEVIRI data, including but not limited, weather forecasting and global climate monitoring at a range of geographical scales.