Alexei Lyapustin

Hyperspectral remote sensing of foliar nitrogen content

A strong positive correlation between vegetation canopy bidirectional reflectance factor (BRF) in the near infrared (NIR) spectral region and foliar mass-based nitrogen concentration (%N) has been reported in some temperate and boreal forests. This relationship, if true, would indicate an additional role for nitrogen in the climate system via its influence on surface albedo and may offer a simple approach for monitoring foliar nitrogen using satellite data. We report, however, that the previously reported correlation is an artifact—it is a consequence of variations in canopy structure, rather than of %N. The data underlying this relationship were collected at sites with varying proportions of foliar nitrogen-poor needleleaf and nitrogen-rich broadleaf species, whose canopy structure differs considerably. When the BRF data are corrected for canopy-structure effects, the residual reflectance variations are negatively related to %N at all wavelengths in the interval 423–855 nm. This suggests that the observed positive correlation between BRF and %N conveys no information about %N. We find that to infer leaf biochemical constituents, e.g., N content, from remotely sensed data, BRF spectra in the interval 710–790 nm provide critical information for correction of structural influences. Our analysis also suggests that surface characteristics of leaves impact remote sensing of its internal constituents. This further decreases the ability to remotely sense canopy foliar nitrogen. Finally, the analysis presented here is generic to the problem of remote sensing of leaf-tissue constituents and is therefore not a specific critique of articles espousing remote sensing of foliar %N.

Vegetation Dynamics and Rainfall Sensitivity of the Amazon

Significance Understanding the sensitivity of tropical vegetation to changes in precipitation is of key importance for assessing the fate of the Amazon rainforest and predicting atmospheric CO2 levels. Using improved satellite observations, we reconcile observational and modeling studies by showing that tropical vegetation is highly sensitive to changes in precipitation and El Niño events. Our results show that, since the year 2000, the Amazon forest has declined across an area of 5.4 million km2 as a result of well-described reductions in rainfall. We conclude that, if drying continues across Amazonia, which is predicted by several global climate models, this drying may accelerate global climate change through associated feedbacks in carbon and hydrological cycles. We show that the vegetation canopy of the Amazon rainforest is highly sensitive to changes in precipitation patterns and that reduction in rainfall since 2000 has diminished vegetation greenness across large parts of Amazonia. Large-scale directional declines in vegetation greenness may indicate decreases in carbon uptake and substantial changes in the energy balance of the Amazon. We use improved estimates of surface reflectance from satellite data to show a close link between reductions in annual precipitation, El Niño southern oscillation events, and photosynthetic activity across tropical and subtropical Amazonia. We report that, since the year 2000, precipitation has declined across 69% of the tropical evergreen forest (5.4 million km2) and across 80% of the subtropical grasslands (3.3 million km2). These reductions, which coincided with a decline in terrestrial water storage, account for about 55% of a satellite-observed widespread decline in the normalized difference vegetation index (NDVI). During El Niño events, NDVI was reduced about 16.6% across an area of up to 1.6 million km2 compared with average conditions. Several global circulation models suggest that a rise in equatorial sea surface temperature and related displacement of the intertropical convergence zone could lead to considerable drying of tropical forests in the 21st century. Our results provide evidence that persistent drying could degrade Amazonian forest canopies, which would have cascading effects on global carbon and climate dynamics.

Global Estimates of Fine Particulate Matter Using a Combined Geophysical-Statistical Method with Information from Satellites, Models, and Monitors

We estimated global fine particulate matter (PM2.5) concentrations using information from satellite-, simulation- and monitor-based sources by applying a Geographically Weighted Regression (GWR) to global geophysically based satellite-derived PM2.5 estimates. Aerosol optical depth from multiple satellite products (MISR, MODIS Dark Target, MODIS and SeaWiFS Deep Blue, and MODIS MAIAC) was combined with simulation (GEOS-Chem) based upon their relative uncertainties as determined using ground-based sun photometer (AERONET) observations for 1998-2014. The GWR predictors included simulated aerosol composition and land use information. The resultant PM2.5 estimates were highly consistent (R(2) = 0.81) with out-of-sample cross-validated PM2.5 concentrations from monitors. The global population-weighted annual average PM2.5 concentrations were 3-fold higher than the 10 μg/m(3) WHO guideline, driven by exposures in Asian and African regions. Estimates in regions with high contributions from mineral dust were associated with higher uncertainty, resulting from both sparse ground-based monitoring, and challenging conditions for retrieval and simulation. This approach demonstrates that the addition of even sparse ground-based measurements to more globally continuous PM2.5 data sources can yield valuable improvements to PM2.5 characterization on a global scale.

A new hybrid spatio-temporal model for estimating daily multi-year PM2.5 concentrations across northeastern USA using high resolution aerosol optical depth data

BACKGROUND The use of satellite-based aerosol optical depth (AOD) to estimate fine particulate matter (PM2.5) for epidemiology studies has increased substantially over the past few years. These recent studies often report moderate predictive power, which can generate downward bias in effect estimates. In addition, AOD measurements have only moderate spatial resolution, and have substantial missing data. METHODS We make use of recent advances in MODIS satellite data processing algorithms (Multi-Angle Implementation of Atmospheric Correction (MAIAC), which allow us to use 1 km (versus currently available 10 km) resolution AOD data. We developed and cross validated models to predict daily PM2.5 at a 1×1km resolution across the northeastern USA (New England, New York and New Jersey) for the years 2003-2011, allowing us to better differentiate daily and long term exposure between urban, suburban, and rural areas. Additionally, we developed an approach that allows us to generate daily high-resolution 200 m localized predictions representing deviations from the area 1×1 km grid predictions. We used mixed models regressing PM2.5 measurements against day-specific random intercepts, and fixed and random AOD and temperature slopes. We then use generalized additive mixed models with spatial smoothing to generate grid cell predictions when AOD was missing. Finally, to get 200 m localized predictions, we regressed the residuals from the final model for each monitor against the local spatial and temporal variables at each monitoring site. RESULTS Our model performance was excellent (mean out-of-sample R2=0.88). The spatial and temporal components of the out-of-sample results also presented very good fits to the withheld data (R2=0.87, R2=0.87). In addition, our results revealed very little bias in the predicted concentrations (Slope of predictions versus withheld observations = 0.99). CONCLUSION Our daily model results show high predictive accuracy at high spatial resolutions and will be useful in reconstructing exposure histories for epidemiological studies across this region.

An automatic cloud mask algorithm based on time series of MODIS measurements

[1] Quality of aerosol retrievals and atmospheric correction over land depends strongly on accuracy of the cloud mask (CM) algorithm. The heritage CM algorithms developed for AVHRR and MODIS use the latest sensor measurements of spectral reflectance and brightness temperature and perform processing at the pixel level. The algorithms are threshold-based and empirically tuned. They do not explicitly address the classical problem of cloud search, wherein the baseline clear-skies scene is defined for comparison. Here we report on a new land CM algorithm, which explicitly builds and maintains a reference clear-skies image of the surface (refcm) using a time series of MODIS measurements. The new algorithm, developed as part of the multiangle implementation of atmospheric correction (MAIAC) algorithm for MODIS, relies on the fact that clear-skies images of the same surface area have a common textural pattern, defined by the surface topography, boundaries of rivers and lakes, distribution of soils and vegetation, etc. This pattern changes slowly given the daily rate of global Earth observations, whereas clouds introduce high-frequency random disturbances. Under clear skies, consecutive gridded images of the same surface area have a high covariance, whereas in presence of clouds covariance is usually low. This idea is central to initialization of refcm, which is used to derive cloud mask in combination with spectral and brightness temperature tests. The refcm is continuously updated with the latest clear-skies MODIS measurements, thus adapting to seasonal and rapid surface changes. The algorithm is enhanced by an internal dynamic land-water-snow classification coupled with a surface change mask. An initial comparison shows that the new algorithm offers the potential to perform better than the MODIS MOD35 cloud mask in situations where the land surface is changing rapidly and over Earth regions covered by snow and ice.

Land and cryosphere products from Suomi NPP VIIRS: Overview and status

[1] The Visible Infrared Imaging Radiometer Suite (VIIRS) instrument was launched in October 2011 as part of the Suomi National Polar-Orbiting Partnership (S-NPP). The VIIRS instrument was designed to improve upon the capabilities of the operational Advanced Very High Resolution Radiometer and provide observation continuity with NASA’s Earth Observing System’s Moderate Resolution Imaging Spectroradiometer (MODIS). Since the VIIRS first-light images were received in November 2011, NASA- and NOAA-funded scientists have been working to evaluate the instrument performance and generate land and cryosphere products to meet the needs of the NOAA operational users and the NASA science community. NOAA’s focus has been on refining a suite of operational products known as Environmental Data Records (EDRs), which were developed according to project specifications under the National Polar-Orbiting Environmental Satellite System. The NASA S-NPP Science Team has focused on evaluating the EDRs for science use, developing and testing additional products to meet science data needs, and providing MODIS data product continuity. This paper presents to-date findings of the NASA Science Team’s evaluation of the VIIRS land and cryosphere EDRs, specifically Surface Reflectance, Land Surface Temperature, Surface Albedo, Vegetation Indices, Surface Type, Active Fires, Snow Cover, Ice Surface Temperature, and Sea Ice Characterization. The study concludes that, for MODIS data product continuity and earth system science, an enhanced suite of land and cryosphere products and associated data system capabilities are needed beyond the EDRs currently available from the VIIRS.

Radiative transfer codes for atmospheric correction and aerosol retrieval: intercomparison study

Results are summarized for a scientific project devoted to the comparison of four atmospheric radiative transfer codes incorporated into different satellite data processing algorithms, namely, 6SV1.1 (second simulation of a satellite signal in the solar spectrum, vector, version 1.1), RT3 (radiative transfer), MODTRAN (moderate resolution atmospheric transmittance and radiance code), and SHARM (spherical harmonics). The performance of the codes is tested against well-known benchmarks, such as Coulsons tabulated values and a Monte Carlo code. The influence of revealed differences on aerosol optical thickness and surface reflectance retrieval is estimated theoretically by using a simple mathematical approach. All information about the project can be found at http://rtcodes.ltdri.org.

Using High-Resolution Satellite Aerosol Optical Depth To Estimate Daily PM2.5 Geographical Distribution in Mexico City

Recent advances in estimating fine particle (PM2.5) ambient concentrations use daily satellite measurements of aerosol optical depth (AOD) for spatially and temporally resolved exposure estimates. Mexico City is a dense megacity that differs from other previously modeled regions in several ways: it has bright land surfaces, a distinctive climatological cycle, and an elevated semi-enclosed air basin with a unique planetary boundary layer dynamic. We extend our previous satellite methodology to the Mexico City area, a region with higher PM2.5 than most U.S. and European urban areas. Using a novel 1 km resolution AOD product from the MODIS instrument, we constructed daily predictions across the greater Mexico City area for 2004-2014. We calibrated the association of AOD to PM2.5 daily using municipal ground monitors, land use, and meteorological features. Predictions used spatial and temporal smoothing to estimate AOD when satellite data were missing. Our model performed well, resulting in an out-of-sample cross-validation R(2) of 0.724. Cross-validated root-mean-squared prediction error (RMSPE) of the model was 5.55 μg/m(3). This novel model reconstructs long- and short-term spatially resolved exposure to PM2.5 for epidemiological studies in Mexico City.

10-year spatial and temporal trends of PM 2.5 concentrations in the southeastern US estimated using high-resolution satellite data

Long-term PM2.5 exposure has been associated with various adverse health outcomes. However, most ground monitors are located in urban areas, leading to a potentially biased representation of true regional PM2.5 levels. To facilitate epidemiological studies, accurate estimates of the spatiotemporally continuous distribution of PM2.5 concentrations are important. Satellite-retrieved aerosol optical depth (AOD) has been increasingly used for PM2.5 concentration estimation due to its comprehensive spatial coverage. Nevertheless, previous studies indicated that an inherent disadvantage of many AOD products is their coarse spatial resolution. For instance, the available spatial resolutions of the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multiangle Imaging SpectroRadiometer (MISR) AOD products are 10 and 17.6 km, respectively. In this paper, a new AOD product with 1 km spatial resolution retrieved by the multi-angle implementation of atmospheric correction (MAIAC) algorithm based on MODIS measurements was used. A two-stage model was developed to account for both spatial and temporal variability in the PM2.5–AOD relationship by incorporating the MAIAC AOD, meteorological fields, and land use variables as predictors. Our study area is in the southeastern US centered at the Atlanta metro area, and data from 2001 to 2010 were collected from various sources. The model was fitted annually, and we obtained model fitting R2 ranging from 0.71 to 0.85, mean prediction error (MPE) from 1.73 to 2.50 μg m−3, and root mean squared prediction error (RMSPE) from 2.75 to 4.10 μg m−3. In addition, we found cross-validation R2 ranging from 0.62 to 0.78, MPE from 2.00 to 3.01 μgm−3, and RMSPE from 3.12 to 5.00 μgm−3, indicating a good agreement between the estimated and observed values. Spatial trends showed that high PM2.5 levels occurred in urban areas and along major highways, while low concentrations appeared in rural or mountainous areas. Our time-series analysis showed that, for the 10-year study period, the PM2.5 levels in the southeastern US have decreased by ∼20 %. The annual decrease has been relatively steady from 2001 to 2007 and from 2008 to 2010 while a significant drop occurred between 2007 and 2008. An observed increase in PM2.5 levels in year 2005 is attributed to elevated sulfate concentrations in the study area in warm months of 2005.

Sunlight mediated seasonality in canopy structure and photosynthetic activity of Amazonian rainforests

Resolving the debate surrounding the nature and controls of seasonal variation in the structure and metabolism of Amazonian rainforests is critical to understanding their response to climate change. In situ studies have observed higher photosynthetic and evapotranspiration rates, increased litterfall and leaf flushing during the Sunlight-rich dry season. Satellite data also indicated higher greenness level, a proven surrogate of photosynthetic carbon fixation, and leaf area during the dry season relative to the wet season. Some recent reports suggest that rainforests display no seasonal variations and the previous results were satellite measurement artefacts. Therefore, here we re-examine several years of data from three sensors on two satellites under a range of sun positions and satellite measurement geometries and document robust evidence for a seasonal cycle in structure and greenness of wet equatorial Amazonian rainforests. This seasonal cycle is concordant with independent observations of solar radiation. We attribute alternative conclusions to an incomplete study of the seasonal cycle, i.e. the dry season only, and to prognostications based on a biased radiative transfer model. Consequently, evidence of dry season greening in geometry corrected satellite data was ignored and the absence of evidence for seasonal variation in lidar data due to noisy and saturated signals was misinterpreted as evidence of the absence of changes during the dry season. Our results, grounded in the physics of radiative transfer, buttress previous reports of dry season increases in leaf flushing, litterfall, photosynthesis and evapotranspiration in well-hydrated Amazonian rainforests.