Global retrieval of bidirectional reflectance and albedo over land from EOS MODIS and MISR data: Theory and algorithm

This paper describes the theory and the algorithm to be used in producing a global bidirectional reflectance distribution function (BRDF) and albedo product from data to be acquired by the moderate resolution imaging spectroradiometer (MODIS) and the multiangle imaging spectroradiometer (MISR), both to be launched in 1998 on the AM-I satellite platform as part of NASA's Earth Observing System (EOS). The product will be derived using the kernel-driven semiempirical Ambrals BRDF model, utilizing five variants of kernel functions characterizing isotropic, volume and surface scattering. The BRDF and the albedo of each pixel of the land surface will be modeled at a spatial resolution of I km and once every 16 days in seven spectral bands spanning the visible and the near infrared. The BRDF parameters retrieved and recorded in the MODIS BRDF/albedo product will be intrinsic surface properties decoupled from the prevailing atmospheric state and hence suited for a wide range of applications requiring characterization of the directional anisotropy of Earth surface reflectance. A set of quality control flags accompanies the product. An initial validation of the Ambrals model is demonstrated using a variety of field-measured data sets for several different land cover types

Remote sensing of snow : Factors influencing seasonal snow mapping in boreal forest region

Monitoring of snow cover in northern hemisphere is highly important for climate research and for operational activities, such as those related to hydrology and weather forecasting. The appearance and melting of seasonal snow cover dominate the hydrological and climatic patterns in the boreal and arctic regions. Spatial variability (in particular during the spring and autumn transition months) and long-term trends in global snow cover distribution are strongly interconnected to changes in Earth System (ES). Satellite data based estimates on snow cover extent are utilized e.g. in near-real-time hydrological forecasting, water resource management and to construct long-term Climate Data Records (CDRs) essential for climate research. Information on the quantitative reliability of snow cover monitoring is urgently needed by these different applications as the usefulness of satellite data based results is strongly dependent on the quality of the interpretation. This doctoral dissertation investigates the factors affecting the reliability of snow cover monitoring using optical satellite data and focuses on boreal regions (zone characterized by seasonal snow cover). Based on the analysis of different factors relevant to snow mapping performance, the work introduces a methodology to assess the uncertainty of snow cover extent estimates, focusing on the retrieval of fractional snow cover (within a pixel) during the spring melt period. The results demonstrate that optical remote sensing is well suited for determining snow extent in the melting season and that the characterizing the uncertainty in snow estimates facilitates the improvement of the snow mapping algorithms. The overall message is that using a versatile accuracy analysis it is possible to develop uncertainty estimates for the optical remote sensing of snow cover, which is a considerable advance in remote sensing. The results of this work can also be utilized in the development of other interpretation algorithms. This thesis consists of five articles predominantly dealing with quantitative data analysis, while the summary chapter synthesizes the results mainly in the algorithm accuracy point of view. The first four articles determine the reflectance characteristics essential for the forward and inverse modeling of boreal landscapes (forward model describes the observations as a function of the investigated variable). The effects of snow, snow-free ground and boreal forest canopy on the observed satellite scene reflectance are specified. The effects of all the error components are clarified in the fifth article and a novel experimental method to analyze and quantify the amount of uncertainty is presented. The five articles employ different remote sensing and ground truth data sets measured and/or analyzed for this research, covering the region of Finland and also applied to boreal forest region in northern Europe

Multi-angle implementation of atmospheric correction for MODIS (MAIAC): 3. Atmospheric correction

This paper describes the atmospheric correction (AC) component of the Multi-Angle Implementation of Atmospheric Correction algorithm (MAIAC) which introduces a new way to compute parameters of the Ross-Thick Li-Sparse (RTLS) Bi-directional reflectance distribution function (BRDF), spectral surface albedo and bidirectional reflectance factors (BRF) from satellite measurements obtained by the Moderate Resolution Imaging Spectroradiometer (MODIS). MAIAC uses a time series and spatial analysis for cloud detection, aerosol retrievals and atmospheric correction. It implements a moving window of up to 16 days of MODIS data gridded to 1 km resolution in a selected projection. The RTLS parameters are computed directly by fitting the cloud-free MODIS top of atmosphere (TOA) reflectance data stored in the processing queue. The RTLS retrieval is applied when the land surface is stable or changes slowly. In case of rapid or large magnitude change (as for instance caused by disturbance), MAIAC follows the MODIS operational BRDF/albedo algorithm and uses a scaling approach where the BRDF shape is assumed stable but its magnitude is adjusted based on the latest single measurement. To assess the stability of the surface, MAIAC features a change detection algorithm which analyzes relative change of reflectance in the Red and NIR bands during the accumulation period. To adjust for the reflectance variability with the sun-observer geometry and allow comparison among different days (view geometries), the BRFs are normalized to the fixed view geometry using the RTLS model. An empirical analysis of MODIS data suggests that the RTLS inversion remains robust when the relative change of geometry-normalized reflectance stays below 15%. This first of two papers introduces the algorithm, a second, companion paper illustrates its potential by analyzing MODIS data over a tropical rainforest and assessing errors and uncertainties of MAIAC compared to conventional MODIS products.Keywords: MAIAC, Atmospheric correction, Time series, Multi-angle implementation of atmospheric correction, MODIS, Aerosols, Surface reflectanc

Forest structure and aboveground biomass in the southwestern United States from MODIS and MISR

Red band bidirectional reflectance factor data from the NASA MODerate resolution Imaging Spectroradiometer (MODIS) acquired over the southwestern United States were interpreted through a simple geometric–optical (GO) canopy reflectance model to provide maps of fractional crown cover (dimensionless), mean canopy height (m), and aboveground woody biomass (Mg ha−1) on a 250 m grid. Model adjustment was performed after dynamic injection of a background contribution predicted via the kernel weights of a bidirectional reflectance distribution function (BRDF) model. Accuracy was assessed with respect to similar maps obtained with data from the NASA Multiangle Imaging Spectroradiometer (MISR) and to contemporaneous US Forest Service (USFS) maps based partly on Forest Inventory and Analysis (FIA) data. MODIS and MISR retrievals of forest fractional cover and mean height both showed compatibility with the USFS maps, with MODIS mean absolute errors (MAE) of 0.09 and 8.4 m respectively, compared with MISR MAE of 0.10 and 2.2 m, respectively. The respective MAE for aboveground woody biomass was ~10 Mg ha−1, the same as that from MISR, although the MODIS retrievals showed a much weaker correlation, noting that these statistics do not represent evaluation with respect to ground survey data. Good height retrieval accuracies with respect to averages from high resolution discrete return lidar data and matches between mean crown aspect ratio and mean crown radius maps and known vegetation type distributions both support the contention that the GO model results are not spurious when adjusted against MISR bidirectional reflectance factor data. These results highlight an alternative to empirical methods for the exploitation of moderate resolution remote sensing data in the mapping of woody plant canopies and assessment of woody biomass loss and recovery from disturbance in the southwestern United States and in parts of the world where similar environmental conditions prevail

Atmospheric correction of visible to middle-infrared EOS-MODIS data over land surfaces: Background, operational algorithm and validation

The NASA moderate resolution imaging spectroradiometer (MODIS) instrument will provide a global and improved source of information for the study of land surfaces with a spatial resolution of up to 250 m

A model to estimate daily albedo from remote sensing data : accuracy assessment of MODIS MCD43 product

L’albedo superficial és un paràmetre físic que afecta al clima de la Terra i, a més, suposa una de les majors incerteses radiatives en l’actual modelització climàtica. Aquest paràmetre és molt variable tant a nivell espacial com temporal degut als canvis en les propietats de les superfícies i als canvis en les condicions d’il•luminació. En conseqüència, es requereix una resolució temporal diària de l’albedo per a realitzar estudis climàtics. L’augment de la resolució espacial dels models climàtics fa necessari l’estudi de les seues característiques espacials a nivell global. La teledetecció proporciona l’única opció pràctica de proporcionar dades d’albedo a nivell global amb alta qualitat i alta resolució tant espacial com temporal. El sensor MODerate Resolution Imaging Spectroradiometer (MODIS) a bord dels satèl•lits Terra i Aqua presenta unes característiques adequades per a l’estimació d’aquest paràmetre. En el present treball realitzem diversos estudis buscant les possibles fonts d’incertesa del producte oficial d’albedo de MODIS (MCD43). A més, presentem un model que millora la resolució temporal d’aquest paràmetre.Surface albedo is a critical land physical parameter affecting the earth’s climate and is among the main radiative uncertainties in current climate modelling. This parameter is highly variable in space and time, both as a result of changes in surface properties and as a function of changes in the illumination conditions. Consequently, an albedo daily temporal resolution is required for climate studies. The increasing spatial resolution of modern climate models makes it necessary to examine its spatial features. Satellite remote sensing provides the only practical way of producing high-quality global albedo data sets with high spatial and temporal resolutions. The MODerate Resolution Imaging Spectroradiometer (MODIS) sensor on board the Terra and Aqua satellites presents the required sampling characteristics in order to derive the this parameter. In this PhD we develop several studies looking for the improvement of the official MODIS albedo product (MCD43) accuracy. Moreover, we present a model that improves the temporal resolution of this parameter

Retrieval of aerosol optical thickness over snow and ice surfaces in the Arctic using Advanced Along Track Scanning Radiometer

Aerosols in the Arctic cause radiative forcing and a variety of climatic feedbacks, which affect climate of both local and global scales. In order to assess the state of the Arctic climate, information on the aerosol type and amount is needed. Harsh conditions and remoteness of the Arctic region result in very few ground based measurements of aerosol optical thickness. Remote sensing has the potential to provide the necessary temporal and spatial coverage. A non-trivial task of aerosol retrieval over a very bright surface is being solved within the thesis; the developed retrieval consists of cloud screening over snow and two types of aerosol retrieval over snow - in the visible and infrared spectral regions. A number of validation and case studies has been performed to assess the quality of the retrieval. The developed algorithm applies to the data of Advanced Along Track Scanning Radiometer and produces maps of aerosol optical thickness over snow and ice

Monitoring Land Surface Albedo and Vegetation Dynamics Using High Spatial and Temporal Resolution Synthetic Time Series from Landsat and the MODIS BRDF/NBAR/Albedo Product

Seasonal vegetation phenology can significantly alter surface albedo which in turn affects the global energy balance and the albedo warmingcooling feedbacks that impact climate change. To monitor and quantify the surface dynamics of heterogeneous landscapes, high temporal and spatial resolution synthetic time series of albedo and the enhanced vegetation index (EVI) were generated from the 500-meter Moderate Resolution Imaging Spectroradiometer (MODIS) operational Collection V006 daily BRDF (Bidirectional Reflectance Distribution Function) / NBAR (Nadir BRDF-Adjusted Reflectance) / albedo products and 30-meter Landsat 5 albedo and near-nadir reflectance data through the use of the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM). The traditional Landsat Albedo (Shuai et al., 2011) makes use of the MODIS BRDFAlbedo products (MCD43) by assigning appropriate BRDFs from coincident MODIS products to each Landsat image to generate a 30-meter Landsat albedo product for that acquisition date. The available cloud free Landsat 5 albedos (due to clouds, generated every 16 days at best) were used in conjunction with the daily MODIS albedos to determine the appropriate 30-meter albedos for the intervening daily time steps in this study. These enhanced daily 30-meter spatial resolution synthetic time series were then used to track albedo and vegetation phenology dynamics over three Ameriflux tower sites (Harvard Forest in 2007, Santa Rita in 2011 and Walker Branch in 2005). These Ameriflux sites were chosen as they are all quite nearby new towers coming on line for the National Ecological Observatory Network (NEON), and thus represent locations which will be served by spatially paired albedo measures in the near future. The availability of data from the NEON towers will greatly expand the sources of tower albedometer data available for evaluation of satellite products. At these three Ameriflux tower sites the synthetic time series of broadband shortwave albedos were evaluated using the tower albedo measurements with a Root Mean Square Error (RMSE) less than 0.013 and a bias within the range of 0.006. These synthetic time series provide much greater spatial detail than the 500 meter gridded MODIS data, especially over more heterogeneous surfaces, which improves the efforts to characterize and monitor the spatial variation across species and communities. The mean of the difference between maximum and minimum synthetic time series of albedo within the MODIS pixels over a subset of satellite data of Harvard Forest (16 kilometers by 14 kilometers) was as high as 0.2 during the snow-covered period and reduced to around 0.1 during the snow-free period. Similarly, we have used STARFM to also couple MODIS Nadir BRDF-Adjusted Reflectances (NBAR) values with Landsat 5 reflectances to generate daily synthetic times series of NBAR and thus Enhanced Vegetation Index (NBAR-EVI) at a 30-meter resolution. While normally STARFM is used with directional reflectances, the use of the view angle corrected daily MODIS NBAR values will provide more consistent time series. These synthetic times series of EVI are shown to capture seasonal vegetation dynamics with finer spatial and temporal details, especially over heterogeneous land surfaces