Towards multidecadal consistent Meteosat surface albedo time series

Monitoring of land surface albedo dynamics is important for the understanding of observed climate trends. Recently developed multidecadal surface albedo data products, derived from a series of geostationary satellite data, provide the opportunity to study long term surface albedo dynamics at the regional to global scale. Reliable estimates of temporal trends in surface albedo require carefully calibrated and homogenized long term satellite data records and derived products. The present paper investigates the long term consistency of a new surface albedo product derived from Meteosat First Generation (MFG) geostationary satellites for the time period 1982–2006. The temporal consistency of the data set is characterized. The analysis of the long term homogeneity reveals some discrepancies in the time series related to uncertainties in the characterization of the sensor spectral response of some of the MFG satellites. A method to compensate for uncertainties in the current data product is proposed and evaluated

Surface radiation budget for climate applications

The Surface Radiation Budget (SRB) consists of the upwelling and downwelling radiation fluxes at the surface, separately determined for the broadband shortwave (SW) (0 to 5 micron) and longwave (LW) (greater than 5 microns) spectral regions plus certain key parameters that control these fluxes, specifically, SW albedo, LW emissivity, and surface temperature. The uses and requirements for SRB data, critical assessment of current capabilities for producing these data, and directions for future research are presented

Land surface albedo from MSG/SEVIRI: retrieval method, validation, and application for weather forecast

The European Meteorological Satellite Organization (EUMETSAT) maintains a number of decentralized processing centers dedicated to different scientific themes. The Portuguese Meteorological Institute hosts the Satellite Application Facility on Land Surface Analysis (LSA-SAF). The primary objective of the LSA-SAF is to provide added-value products for the meteorological and environmental science communities with main applications in the fields of climate modeling, environmental management, natural hazards management, and climate change detection. Since 2005 data from Meteosat Second Generation satellite are routinely processed in near real time by the LSA-SAF operational system in Lisbon. Presently, the delivered operational products comprise land surface albedo and temperature, shortwave and long-wave downwelling radiation fluxes, vegetation parameters and snow cover. After more than ten years (1999-2010) of research, development, and progressive operational activities, a summary of the surface albedo product characteristics and performances is presented. The relevance of LSA-SAF albedo product is analyzed through a weather forecast model (ALADIN) in order to account for the inter-annual spatial and temporal variability. Results clearly show a positive impact on the 12-hour forecast of 2m temperatures

Simulating Meteosat-7 broadband radiances using two visible channels of Meteosat-8

International audienceThis paper addresses the conversion of narrow-to-broadband radiances in the special case of the sensor SEVIRI onboard the satellite Meteosat-8 of the Meteosat Second Generation (MSG). The advent of this new program poses the problem of the transfer of operational procedures from the previous satellite Meteosat-7 to MSG. Among these operational procedures are several different implementations of the Heliosat method to convert Meteosat data into irradiance maps. An easy-to-implement fast-running method is proposed that combines the radiances of the two narrow visible bands of SEVIRI to produce broadband radiances that are almost identical to those observed in the broadband channel of Meteosat-7. A comparison between Meteosat-7 actual radiances and these simulated-radiances shows errors with a negligible bias and a relative root mean square error (RMSE) of 6%. After combination with the uncertainties found in the calibration procedures (10-13%), the total relative RMSE on the radiances amounts to 12-14%. It is concluded that the Meteosat-8 images may be used in operational procedures currently applied to Meteosat-7 images

ESTIMATING LAND SURFACE ALBEDO FROM SATELLITE DATA

Land surface albedo, defined as the ratio of the surface reflected incoming and outgoing solar radiation, is one of the key geophysical variables controlling the surface radiation budget. Surface shortwave albedo is widely used to drive climate and hydrological models. During the last several decades, remotely sensed surface albedo products have been generated through satellite-acquired data. However, some problems exist in those products due to instrument measurement inaccuracies and the failure of current retrieving procedures, which have limited their applications. More significantly, it has been reported that some albedo products from different satellite sensors do not agree with each other and some even show the opposite long term trend regionally and globally. The emergence of some advanced sensors newly launched or planned in the near future will provide better capabilities for estimating land surface albedo with fine resolution spatially and/or temporally. Traditional methods for estimating the surface shortwave albedo from satellite data include three steps: first, the satellite observations are converted to surface directional reflectance using the atmospheric correction algorithms; second, the surface bidirectional reflectance distribution function (BRDF) models are inverted through the fitting of the surface reflectance composites; finally, the shortwave albedo is calculated from the BRDF through the angular and spectral integration. However, some problems exist in these algorithms, including: 1) "dark-object" based atmospheric correction methods which make it difficult to estimate albedo accurately over non-vegetated or sparsely vegetated area; 2) the long-time composite albedo products cannot satisfy the needs of weather forecasting or land surface modeling when rapid changes such as snow fall/melt, forest fire/clear-cut and crop harvesting occur; 3) the diurnal albedo signature cannot be estimated in the current algorithms due to the Lambertian approximation in some of the atmospheric correction algorithms; 4) prior knowledge has not been effectively incorporated in the current algorithms; and 5) current observation accumulation methods make it difficult to obtain sufficient observations when persistent clouds exist within the accumulation window. To address those issues and to improve the satellite surface albedo estimations, a method using an atmospheric radiative transfer procedure with surface bidirectional reflectance modeling will be applied to simultaneously retrieve land surface albedo and instantaneous aerosol optical depth (AOD). This study consists of three major components. The first focuses on the atmospheric radiative transfer procedure with surface reflectance modeling. Instead of executing atmospheric correction first and then fitting surface reflectance in the previous satellite albedo retrieving procedure, the atmospheric properties (e.g., AOD) and surface properties (e.g., BRDF) are estimated simultaneously to reduce the uncertainties produced in separating the entire radiative transfer process. Data from the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard Terra and Aqua are used to evaluate the performance of this albedo estimation algorithm. Good agreement is reached between the albedo estimates from the proposed algorithm and other validation datasets. The second part is to assess the effectiveness of the proposed algorithm, analyze the error sources, and further apply the algorithm on geostationary satellite - the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) onboard Meteosat Second Generation (MSG). Extensive validations on surface albedo estimations from MSG/SEVIRI observations are conducted based on the comparison with ground measurements and other satellite products. Diurnal changes and day-to-day changes in surface albedo are accurately captured by the proposed algorithm. The third part of this study is to develop a spatially and temporally complete, continuous, and consistent albedo maps through a data fusion method. Since the prior information (or climatology) of albedo/BRDF plays a vital role in controlling the retrieving accuracy in the optimization method, currently available multiple land surface albedo products will be integrated using the Multi-resolution Tree (MRT) models to mitigate problems such as data gaps, systematic bias or low information-noise ratio due to instrument failure, persistent clouds from the viewing direction and algorithm limitations. The major original contributions of this study are as follows: 1) this is the first algorithm for the simultaneous estimations of surface albedo/reflectance and instantaneous AOD by using the atmospheric radiative transfer with surface BRDF modeling for both polar-orbiting and geostationary satellite data; 2) a radiative transfer with surface BRDF models is used to derive surface albedo and directional reflectance from MODIS and SEVIRI observations respectively; 3) extensive validations are made on the comparison between the albedo and AOD retrievals, and the satellite products from other sensors; 4) the slightly modified algorithm has been adopted to be the operational algorithm of Advanced Baseline Imager (ABI) in the future Geostationary Operational Environmental Satellite-R Series (GOES-R) program for estimating land surface albedo; 5) a framework of using MRT is designed to integrate multiple satellite albedo products at different spatial scales to build the spatially and temporally complete, continuous, and consistent albedo maps as the prior knowledge in the retrieving procedure

An Approach for the Long-Term 30-m Land Surface Snow-Free Albedo Retrieval from Historic Landsat Surface Reflectance and MODIS-based A Priori Anisotropy Knowledge

Land surface albedo has been recognized by the Global Terrestrial Observing System (GTOS) as an essential climate variable crucial for accurate modeling and monitoring of the Earth's radiative budget. While global climate studies can leverage albedo datasets from MODIS, VIIRS, and other coarse-resolution sensors, many applications in heterogeneous environments can benefit from higher-resolution albedo products derived from Landsat. We previously developed a "MODIS-concurrent" approach for the 30-meter albedo estimation which relied on combining post-2000 Landsat data with MODIS Bidirectional Reflectance Distribution Function (BRDF) information. Here we present a "pre-MODIS era" approach to extend 30-m surface albedo generation in time back to the 1980s, through an a priori anisotropy Look-Up Table (LUT) built up from the high quality MCD43A BRDF estimates over representative homogenous regions. Each entry in the LUT reflects a unique combination of land cover, seasonality, terrain information, disturbance age and type, and Landsat optical spectral bands. An initial conceptual LUT was created for the Pacific Northwest (PNW) of the United States and provides BRDF shapes estimated from MODIS observations for undisturbed and disturbed surface types (including recovery trajectories of burned areas and non-fire disturbances). By accepting the assumption of a generally invariant BRDF shape for similar land surface structures as a priori information, spectral white-sky and black-sky albedos are derived through albedo-to-nadir reflectance ratios as a bridge between the Landsat and MODIS scale. A further narrow-to-broadband conversion based on radiative transfer simulations is adopted to produce broadband albedos at visible, near infrared, and shortwave regimes.We evaluate the accuracy of resultant Landsat albedo using available field measurements at forested AmeriFlux stations in the PNW region, and examine the consistency of the surface albedo generated by this approach respectively with that from the "concurrent" approach and the coincident MODIS operational surface albedo products. Using the tower measurements as reference, the derived Landsat 30-m snow-free shortwave broadband albedo yields an absolute accuracy of 0.02 with a root mean square error less than 0.016 and a bias of no more than 0.007. A further cross-comparison over individual scenes shows that the retrieved white sky shortwave albedo from the "pre-MODIS era" LUT approach is highly consistent (R(exp 2) = 0.988, the scene-averaged low RMSE = 0.009 and bias = 0.005) with that generated by the earlier "concurrent" approach. The Landsat albedo also exhibits more detailed landscape texture and a wider dynamic range of albedo values than the coincident 500-m MODIS operational products (MCD43A3), especially in the heterogeneous regions. Collectively, the "pre-MODIS" LUT and "concurrent" approaches provide a practical way to retrieve long-term Landsat albedo from the historic Landsat archives as far back as the 1980s, as well as the current Landsat-8 mission, and thus support investigations into the evolution of the albedo of terrestrial biomes at fine resolution

Surface solar irradiance from SCIAMACHY measurements: algorithm and validation

Broadband surface solar irradiances (SSI) are, for the first time, derived from SCIAMACHY (SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY) satellite measurements. The retrieval algorithm, called FRESCO (Fast REtrieval Scheme for Clouds from the Oxygen A band) SSI, is similar to the Heliosat method. In contrast to the standard Heliosat method, the cloud index is replaced by the effective cloud fraction derived from the FRESCO cloud algorithm. The MAGIC (Mesoscale Atmospheric Global Irradiance Code) algorithm is used to calculate clear-sky SSI. The SCIAMACHY SSI product is validated against globally distributed BSRN (Baseline Surface Radiation Network) measurements and compared with ISCCP-FD (International Satellite Cloud Climatology Project Flux Dataset) surface shortwave downwelling fluxes (SDF). For one year of data in 2008, the mean difference between the instantaneous SCIAMACHY SSI and the hourly mean BSRN global irradiances is −4 W m<sup>−2</sup> (−1 %) with a standard deviation of 101 W m<sup>−2</sup> (20 %). The mean difference between the globally monthly mean SCIAMACHY SSI and ISCCP-FD SDF is less than −12 W m<sup>−2</sup> (−2 %) for every month in 2006 and the standard deviation is 62 W m<sup>−2</sup> (12 %). The correlation coefficient is 0.93 between SCIAMACHY SSI and BSRN global irradiances and is greater than 0.96 between SCIAMACHY SSI and ISCCP-FD SDF. The evaluation results suggest that the SCIAMACHY SSI product achieves similar mean bias error and root mean square error as the surface solar irradiances derived from polar orbiting satellites with higher spatial resolution

Estimation of the aerosol radiative forcing at ground level, over land, and in cloudless atmosphere, from METEOSAT-7 observation: method and first results

International audienceA new method is proposed to estimate the spatial and temporal variability of the solar radiative flux reaching the surface (DSSF) over land, as well as the Aerosol Radiative Forcing (ARF), in cloud-free atmosphere. The objective of global applications of the method is fulfilled by using the visible broadband of METEOSAT-7 satellite which scans Europe and Africa on a half-hourly basis. The method relies on a selection of best correspondence between METEOSAT-7 radiance and DSSF computed with a radiative transfer code. The validation of DSSF is performed comparing retrievals with ground-based measurements acquired in two contrasted environments, i.e. an urban site near Paris and a continental background site in South East of France. The study is concentrated on aerosol episodes occurring around the 2003 summer heat wave, providing 42 cases of comparison for variable solar zenith angle (from 59° to 69°), variable aerosol type (biomass burning emissions and urban pollution), and variable aerosol optical thickness (a factor 6). The method reproduces measurements of DSSF within an accuracy assessment of 20 Wm-2 (5% in relative) in 70% of the cases, and within 40 Wm-2 in 90% of the cases. Considering aerosol is the main contributor in changing the measured radiance at the top of the atmosphere, DSSF temporal variability is assumed to be caused only by aerosols, and consequently the ARF at ground level and over land is also retrieved: ARF is computed as the difference between DSSF and a parameterised aerosol-free reference level. Retrievals are linearly correlated with the ground-based measurements of the aerosol optical thickness (AOT): sensitivity is included between 120 and 160 Wm-2 per unity of AOT at 440 nm. AOT being an instantaneous measure indicative of the aerosol columnar amount, we therefore prove the feasibility to infer instantaneous aerosol radiative impact at the ground level over land with METEOSAT-7 visible channel