Konstantin V. Khlopenkov

SPARC: New Cloud, Snow, and Cloud Shadow Detection Scheme for Historical 1-km AVHHR Data over Canada

Abstract The identification of clear-sky and cloudy pixels is a key step in the processing of satellite observations. This is equally important for surface and cloud–atmosphere applications. In this paper, the Separation of Pixels Using Aggregated Rating over Canada (SPARC) algorithm is presented, a new method of pixel identification for image data from the Advanced Very High Resolution Radiometer (AVHRR) on board the NOAA satellites. The SPARC algorithm separates image pixels into clear-sky and cloudy categories based on a specially designed rating scheme. A mask depicting snow/ice and cloud shadows is also generated. The SPARC algorithm has been designed to work year-round (day and night) over the temperate and polar regions of North America, for current and historical AVHRR/NOAA High-Resolution Picture Transmission (HRPT) and Local Area Coverage (LAC) data with original 1-km spatial resolution. The algorithm was tested and applied to data from the AVHRR sensors flown on board NOAA-6 to NOAA-18. The met...

Generating historical AVHRR 1 km baseline satellite data records over Canada suitable for climate change studies

Generating historical AVHRR 1 km baseline satellite data records over Canada suitable for climate change studies Rasim Latifovic, Alexander P. Trishchenko, Ji Chen, William B. Park, Konstantin V. Khlopenkov, Richard Fernandes, Darren Pouliot, Calin Ungureanu, Yi Luo, Shusen Wang, Andrew Davidson, and Josef Cihlar Pages 324-346 Abstract. Satellite data are an important component of the global climate observing system (GCOS). To serve the purpose of climate change monitoring, these data should satisfy certain criteria in terms of the length of observations and the continuity and consistency between different missions and instruments. Despite the great potential and obvious advantages of satellite observations, such as frequent repeat cycles and global coverage, their use in climate studies is hindered by substantial difficulties arising from large data volumes, complicated processing, and significant computer resources required for archiving and analysis. Successful examples of satellite earth observation (EO) data in climate studies include, among others, analyses of the earths radiation budget (Earth Radiation Budget Experiment (ERBE), Scanner for Radiation Budget (ScaRaB), and Cloud and the Earths Radiant Energy System (CERES)), cloudiness (International Satellite Cloud Climatology Project (ISCCP)), vegetation research (Global Inventory Modeling and Mapping Studies (GIMMS)), and the National Oceanic and Atmospheric Administration – National Aeronautics and Space Administration (NOAA–NASA) Pathfinder Program. Despite several attempts, the great potential of the advanced very high resolution radiometer (AVHRR) 1 km satellite data for climate research remains substantially underutilized. To address this issue, the generation of a comprehensive satellite data archive of AVHRR data and products at 1 km spatial resolution over Canada for 1981–2004 (24 years) has been initiated, and a new system for processing at level 1B has been developed. This processing system was employed to generate baseline 1 day and 10 day year-round clear-sky composites for a 5700 km × 4800 km area of North America. This region is centred over Canada but also includes the northern United States, Alaska, Greenland, and surrounding ocean regions. The baseline products include top-of-atmosphere (TOA) visible and near-infrared reflectance, TOA band 4 and band 5 brightness temperature, a cloud – clear – shadow – snow and ice mask, and viewing geometry. Details of the data processing system are presented in the paper. An evaluation of the system characteristics and comparison with previous results demonstrate important improvements in the quality and efficiency of the data processing. The system can process data in a highly automated manner, both for snow-covered and snow-free scenes, and for daytime and nighttime orbits, with high georeferencing accuracy and good radiometric consistency for all sensors from AVHRR NOAA-6 to AVHRR NOAA-17. Other processing improvements include the implementation of advanced algorithms for clear sky – cloud – shadow – snow and ice scene identification, as well as atmospheric correction and compositing. At the time of writing, the assembled dataset is the most comprehensive AVHRR archive at 1 km spatial resolution over Canada that includes all available observations from AVHRR between 1981 and 2004. The archive and the processing system are valuable assets for studying different aspects of land, oceans, and atmosphere related to climate variability and climate change.

Implementation and Evaluation of Concurrent Gradient Search Method for Reprojection of MODIS Level 1B Imagery

This paper presents details regarding implementation of a novel algorithm for reprojection of Moderate Resolution Imaging Spectroradiometer (MODIS) Level 1B imagery. The method is based on a simultaneous 2-D search in latitude and longitude geolocation fields by using their local gradients. Due to the segmented structure of MODIS imagery caused by the instrument whiskbroom electrooptical design, the gradient search is realized in the following two steps: intersegment and intrasegment search. This approach resolves the discontinuity of the latitude/longitude geolocation fields caused by overlap between consecutively scanned MODIS multidetector image segments. The structure of the algorithm allows equal efficiency with nearest neighbor and bilinear interpolation. A special procedure that combines analytical and numerical schemes is designed for reprojecting imagery near the polar region, where the standard gradient search may become unstable. The performance of the method was validated by comparison of reprojected MODIS/Terra and MODIS/Aqua images with georectified Landsat-7 Enhanced Thematic Mapper Plus imagery over Canada. It was found that the proposed method preserves the absolute geolocation accuracy of MODIS pixels determined by the MODIS geolocation team. The method was implemented to reproject MODIS Level 1B imagery over Canada, North America, and Arctic circumpolar zone in the following four popular geographic projections: Plate Care (cylindrical equidistant), Lambert Conic Conformal, Universal Transverse Mercator, and Lambert Azimuthal Equal-Area. It was also found to be efficient for reprojection of Advanced Very High Resolution Radiometer and Medium Resolution Imaging Spectrometer satellite images and general-type meteorological fields, such as the North American Regional Reanalysis data sets.

A method for downscaling MODIS land channels to 250-m spatial resolution using adaptive regression and normalization

A method is proposed to derive spatially enhanced imagery for all seven Moderate Imaging Spectroradiometer (MODIS) land spectral bands at 250 m spatial resolution. Originally, only bands B1 and B2 [visible (VIS) at 0.65 μm, and near-infrared (NIR) at 0.85 μm] are available from MODIS at 250 m spatial resolution. The remaining five land channels (bands B3 to B7) are observed at 500 m resolution. The adaptive regression is constructed for each individual MODIS L1B granule of 500 m spatial resolution by splitting the area into smaller blocks and generating nonlinear regression between bands B3 to B7 and B1, B2 and NDVI. Once a set of regression coefficients is generated based on 500 m image, it is then applied to 250 m data containing only channels B1 and B2 to produce five intermediate synthetic channels (B3 to B7) at 250 m spatial resolution. The final step involves normalizing the generated 250 m images to original 500 m images to preserve radiometric consistency. It is achieved in two stages and ensures that downscaled results are unbiased relative to original observations. The developed method was applied to generate Canada-wide clear-sky composites containing all seven MODIS land spectral channels at 250 m spatial resolution over the area of North America 5700 km by 4800 km.

Achieving Subpixel Georeferencing Accuracy in the Canadian AVHRR Processing System

Precise geolocation is one of the fundamental requirements for satellite imagery to be suitable for climate applications. The Global Climate Observing System and the Committee on Earth Observing Satellites identified the requirement for the accuracy of geolocation of satellite data for climate applications as 1/3 field of view (FOV). This requirement for the series of the Advanced Very High Resolution Radiometer (AVHRR) on the National Oceanic and Atmospheric Administration platforms cannot be met without implementing the ground control point (GCP) correction, particularly for historical data, because of limited accuracy of orbit modeling and knowledge of satellite attitude angles. This paper presents a new method for precise georeferencing of the AVHRR imagery developed as part of the new Canadian AVHRR processing system (CAPS) designed for generating high-quality AVHRR satellite climate data record at 1-km spatial resolution. The method works in swath projection and uses the following: 1) the reference monthly images from Moderate Resolution Imaging Spectroradiometer at 250-m resolution; 2) orthorectification to correct for surface elevation; and 3) a novel image matching technique in swath projection to achieve the subpixel resolution. The method is designed for processing daytime data as it intensively employs observations from optical solar bands, the near-infrared channel in particular. The application of the developed processing system showed that the algorithm achieved better than 1/3 FOV geolocation accuracy for AVHRR 1-km scenes. It has very high efficiency rate (> 97%) due to the dense and uniform GCP coverage of the study area (5700 × 4800 km2 ), covering the entire Canada, the Northern U.S., Alaska, Greenland, and surrounding oceans.

Arctic circumpolar mosaic at 250 m spatial resolution for IPY by fusion of MODIS/TERRA land bands B1-B7

The first spatially enhanced Moderate Resolution Imaging Spectroradiometer (MODIS) clear‐sky mosaic for the Arctic circumpolar zone (9000 km×9000 km) is presented, as a contribution to the Canadian component of the International Polar Year (IPY) Programme. The imagery was obtained by fusion of MODIS bands B1–B2 observed at 250 m spatial resolution with bands B3–B7 observed at 500 m spatial resolution to satisfy the Global Climate Observing System (GCOS) requirement for a spatial resolution of 250 m for satellite‐based products for climate. The fusion method used adaptive regression and normalization to preserve the image radiometric properties. A new cloud and cloud shadow detection method and a clear‐sky compositing scheme were used for the 250 m multispectral data. By the end of the IPY in 2009, a decade‐long (2000–2009) time series of these data documenting the state and variability of the Arctic region at fine spatial (250 m) and temporal (10‐day) resolution will be produced if MODIS continues to operate until the end of this period. The product is generated in the Lambert Azimuthal Equal‐Area (LAEA) projection centred over the North Pole. The major intended application of the new data is mapping the surface albedo at 250 m spatial resolution. This product in turn can be used as an input for generating several other Essential Climate Variables (ECVs) as defined by the GCOS.

A Method to Derive the Multispectral Surface Albedo Consistent with MODIS from Historical AVHRR and VGT Satellite Data

Abstract Multispectral surface albedo and bidirectional properties are required for accurate determination of the surface and atmosphere solar radiation budget. A method is developed here to obtain time series of these surface characteristics consistent with the Moderate Resolution Imaging Spectroradiometer (MODIS) using historical satellite observations with limited spectral coverage available from NOAA Advanced Very High Resolution Radiometer (AVHRR) and VEGETATION/Satellite pour l’Observation de la Terre (SPOT). A nonlinear regression model was developed that relates retrievals from four spectral channels of VEGETATION/SPOT or three spectral channels of NOAA AVHRR with retrieval from each of the seven MODIS channels designed for land applications. The model also takes into account the surface land cover type, the normalized difference vegetation index, and the seasonal cycle. It was applied to generate surface albedo and bidirectional parameters of the seven MODIS-like spectral channels at a 10-day int...

A Probabilistic Multispectral Pattern Recognition Method for Detection of Overshooting Cloud Tops Using Passive Satellite Imager Observations

AbstractDeep convective updrafts often penetrate through the surrounding cirrus anvil and into the lower stratosphere. Cross-tropopause transport of ice, water vapor, and chemicals occurs within these “overshooting tops” (OTs) along with a variety of hazardous weather conditions. OTs are readily apparent in satellite imagery, and, given the importance of OTs for weather and climate, a number of automated satellite-based detection methods have been developed. Some of these methods have proven to be relatively reliable, and their products are used in diverse Earth science applications. Nevertheless, analysis of these methods and feedback from product users indicate that use of fixed infrared temperature–based detection criteria often induces biases that can limit their utility for weather and climate analysis. This paper describes a new multispectral OT detection approach that improves upon those previously developed by minimizing use of fixed criteria and incorporating pattern recognition analyses to arriv...

Perennial snow and ice variations (2000–2008) in the Arctic circumpolar land area from satellite observations

Perennial snow and ice (PSI) extent is an important parameter of mountain environments with regard to its involvement in the hydrological cycle and the surface energy budget. We investigated interannual variations of PSI in nine mountain regions of interest (ROI) between 2000 and 2008. For that purpose, a novel MODIS data set processed at the Canada Centre for Remote Sensing at 250 m spatial resolution was utilized. The extent of PSI exhibited significant interannual variations, with coefficients of variation ranging from 5% to 81% depending on the ROI. A strong negative relationship was found between PSI and positive degree-days (threshold 0°C) during the summer months in most ROIs, with linear correlation coefficients (r) being as low as r = −0.90. In the European Alps and Scandinavia, PSI extent was significantly correlated with annual net glacier mass balances, with r = 0.91 and r = 0.85, respectively, suggesting that MODIS-derived PSI extent may be used as an indicator of net glacier mass balances. Validation of PSI extent in two land surface classifications for the years 2000 and 2005, GLC-2000 and Globcover, revealed significant discrepancies of up to 129% for both classifications. With regard to the importance of such classifications for land surface parameterizations in climate and land surface process models, this is a potential source of error to be investigated in future studies. The results presented here provide an interesting insight into variations of PSI in several ROIs and are instrumental for our understanding of sensitive mountain regions in the context of global climate change assessment.

MTSAT-1R Visible Imager Point Spread Correction Function, Part I: The Need for, Validation of, and Calibration With

The multifunctional transport satellite (MTSAT)-1R imager was launched in 2005 and is operated by the Japan Meteorological Agency (JMA). A nonlinear behavior in the MTSAT-1R visible sensor response is observed when the instrument is intercalibrated with coincident moderate resolution imaging spectroradiometer (MODIS) ray-matched radiances. Analysis reveals that the nonlinear behavior is not a result of imager navigation, sensor spectral response difference, nor scan pattern. Examination of coincident MTSAT-1R and MTSAT-2 images reveals that MTSAT-1R dark ocean radiances are affected by neighboring bright clouds, whereas large regions of dark ocean radiances are not impacted. Although the IR and visible optical paths are shared, the MTSAT-1R brightness temperatures are not affected. A dust contaminant coating the mirror, which only affects certain wavelengths, may be one explanation. To address the nonlinearity, a pixel point spread function (PSF) correction algorithm is implemented, wherein most of the radiance contribution is from the pixel field of view itself, as well as including a small contribution from all pixels within a radii of several hundred kilometers. The application of the PSF-corrected ~80% of the affected pixel radiances. After application, a near linear response is observed between the coincident MTSAT-1R and Aqua-MODIS ray-matched radiances, and the intercept is now near the predicted space count of zero. The monthly calibration gain noise is reduced by one-third when compared with the non-PSF-corrected gains. The monthly gains are the most erratic during the first two years of operation, and the MTSAT-1R visible sensor is degrading at ~1.9 % decade.