Considering the effects of canopy structure on hyperspectral radiative transfer and terrestrial photosynthesis
UNKNOWN
Frankenberg, Christian↗
Engineering topics
Wang, Yujie
Publications and source records attributed to Wang, Yujie.
Global Validation of Columnar Water Vapor Derived from EOS MODIS-MAIAC Algorithm Against the Ground-Based AERONET Observations
The water vapor is a relevant greenhouse gas in the Earth's climate system, and satellite products become one of the most effective way to characterize and monitor the columnar water vapor (CWV) content at global scale. Recently, a new product (MCD19) was released as part of MODIS (Moderate Resolution Imaging Spectroradiometer) Collection 6 (C6). This operational product from the Multi-Angle Implementation for Atmospheric Correction (MAIAC) algorithm includes a high 1-kilometer resolution CWV retrievals. This study presents the first global validation of MAIAC C6 CWV obtained from MODIS MCD19A2 product. This evaluation was performed using Aerosol Robotic Network (AERONET) observations at 265 sites (2000-2017). Overall, the results show a good agreement between MAIAC/AERONET CWV retrievals, with correlation coefficient higher than 0.95 and RMS (Root Mean Square) error lower than 0.250 centimeters. The binned error analysis revealed an underestimation (approximately 10 percent) of Aqua CWV retrievals with negative bias for CWV higher than 3.0 centimeters. In contrast, Terra CWV retrievals show a slope of regression close to unity and a low mean bias of 0.075 centimeters. While the accuracy is relatively similar between 1.0 and 5.0 centimeters for both sensor products, Terra dataset is more reliable for applications in humid tropical areas (less than 5.0 centimeters). The expected error was defined as plus or minus 15 percent, with less than 68 percent of retrievals falling within this envelope. However, the accuracy is regionally dependent, and lower error should be expected in some regions, such as South America and Oceania. Since MODIS instruments have exceeded their design lifetime, time series analysis was also presented for both sensor products. The temporal analysis revealed a systematic offset of global average between Terra and Aqua CWV records. We also found an upward trend (approximately 0.2 centimeters per decade) in Terra CWV retrievals, while Aqua CWV retrievals remain stable over time. The sensor degradation influences the ability to detect climate signals, and this study indicates the need for revisiting calibration of the MODIS bands 17-19, mainly for Terra instrument, to assure the quality of the MODIS water vapor product. Finally, this study presents a comprehensive validation analysis of MAIAC CWV over land, raising the understanding of its overall quality.
MAIAC↗
Land Surface Reflectances from Geostationary Sensors
GEONEX is a processing pipeline that produces a suite of satellite land surface products using data streams from the latest geostationary (GEO) sensors including the GOES016/ABI and the Himawari-8/AHI. The suite, created collaboratively by scientists from NASA and NOAA, includes top-of-atmosphere (TOA) reflectances, land surface reflectances (LSRs), vegetation indices, LAI/fPAR, and other downstream products. As a key component of the GEONEX product processing, we have adapted the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm to produce LSRs from the TOA data. Because the algorithm depends on building "stacks" of images, we first run internal geo-registration checks to ensure geo-spatial accuracy and consistency of the input (L1B) data before transferring them from the geostationary projection into a tile system in geographic grids. Scan-time is inferred from metadata and applied to calculate the sun-sensor angles for each grid cell. The MAIAC algorithm is run to detect clouds/shadows, estimate aerosol optical thickness (AOT), perform atmospheric corrections, and generate LSRs. We have processed 18-months (from 2016/04 onward) of AHI data over East Asia and Oceania at a 10-minute time step and 10-months (from 2018/01 onward) of ABI data over North and South Americas at a 15-minute time step. As a verification measure, we compare the GEONEX (AHI/ABI) surface reflectances with the standard MODIS products (MOD09GA) and the MODIS MAIAC products over pixels that have similar sun-view geometries. The results indicate general linear relationships between GEONEX and corresponding MODIS LSRs. In particular, the RMSEs between GEONEX and MOD09 data are comparable to those between MOD09 and MODIS MAIAC products, suggesting that the uncertainties of GEONEX LSRs fall into an acceptable range. However, direct comparisons of LSRs over pixels with different sun-view angles are not as straightforward and require more modeling efforts to correct the directional effects. Evaluation of such angular influences on the downstream products (e.g., vegetation indices) is also under investigation.
Geostationary satellite; Remote Sensing; Atmospher↗
MODIS Collection 6 MAIAC Algorithm
This paper describes the latest version of the algorithm MAIAC (Multi-Angle Implementation of Atmospheric Correction) used for processing the MODIS (Moderate-resolution Imaging Spectroradiometer) Collection6 data record. Since initial publication in 2011-2012, MAIAC has changed considerably to adapt to global processing and improve cloud/snow detection, aerosol retrievals and atmospheric correction of MODIS data. The main changes include (1) transition from a 25 to 1 km scale for retrieval of the spectral regression coefficient (SRC) which helped to remove occasional blockiness at 25 km scale in the aerosol optical depth (AOD) and in the surface reflectance, (2) continuous improvements of cloud detection, (3) introduction of smoke and dust tests to discriminate absorbing fine- and coarse mode aerosols, (4) adding over-water processing, (5) general optimization of the LUT (LookUp-Table)-based radiative transfer for the global processing, and others. MAIAC provides an interdisciplinary suite of atmospheric and land products, including cloud mask (CM), column water vapor (CWV), AOD at 0.47 and 0.55 m, aerosol type (background, smoke or dust) and fine-mode fraction over water; spectral bidirectional reflectance factors (BRF), parameters of Ross-thick Lisparse (RTLS) bidirectional reflectance distribution function (BRDF) model and instantaneous albedo. For snow-covered surfaces, we provide subpixel snow fraction and snow grain size. All products come in standard HDF4 (software library) format at 1 km resolution, except for BRF, which is also provided at 500 m resolution on a sinusoidal grid adopted by the MODIS Land team. All products are provided on per-observation basis in daily files except for the BRDF/Albedo product, which is reported every 8 days. Because MAIAC uses a time series approach, BRDF/Albedo is naturally gap-filled over land where missing values are filled-in with results from the previous retrieval. While the BRDF model is reported for MODIS Land bands 1-7 and ocean band 8, BRF is reported for both land and ocean bands 1-12. This paper focuses on MAIAC cloud detection, aerosol retrievals and atmospheric correction and describes MCD19 data products and quality assurance (QA) flags.
MAIAC Algorithm↗
Estimation of Terrestrial Global Gross Primary Production (GPP) with Satellite Data-Driven Models and Eddy Covariance Flux Data
We estimate global terrestrial gross primary production (GPP) based on models that use satellite data within a simplified light-use efficiency framework that does not rely upon other meteorological inputs. Satellite-based geometry-adjusted reflectances are from the MODerate-resolution Imaging Spectroradiometer (MODIS) and provide information about vegetation structure and chlorophyll content at both high temporal (daily to monthly) and spatial (1 km) resolution. We use satellite-derived solar-induced fluorescence (SIF) to identify regions of high productivity crops and also evaluate the use of downscaled SIF to estimate GPP. We calibrate a set of our satellite-based models with GPP estimates from a subset of distributed eddy covariance flux towers (FLUXNET 2015). The results of the trained models are evaluated using an independent subset of FLUXNET 2015 GPP data. We show that variations in light-use efficiency (LUE) with incident PAR are important and can be easily incorporated into the models. Unlike many LUE-based models, our satellite-based GPP estimates do not use an explicit parameterization of LUE that reduces its value from the potential maximum under limiting conditions such as temperature and water stress. Even without the parameterized downward regulation, our simplified models are shown to perform as well as or better than state-of-the-art satellite data-driven products that incorporate such parameterizations. A significant fraction of both spatial and temporal variability in GPP across plant functional types can be accounted for using our satellite-based models. Our results provide an annual GPP value of 140 Pg C year 1 for 2007 that is within the range of a compilation of observation-based, model, and hybrid results, but is higher than some previous satellite observation-based estimates
CO2↗
Quantitative Daily Maps of PM 2.5 Episodes for California and Other Regions: Satellite Column Water and Optical Depth as Allied Tracers of Dilution
The Western US and many regions globally present daunting difficulties in understanding PM 2.5 episodes. We evaluate extensions of a method independent of modeled source-description and transport/transformation and using several satellite remote sensing products from imaging spectrometers. The San Joaquin Valley (SJV) especially suffers few-day episodes due to shallow mixing; PM 2.5 retrieval suffers low satellite AOT (Aerosol Optical Thickness) and bright surfaces.Nevertheless, we find residual errors in our maps of of typically 5-8 micrograms per cubic meter. Episodes in the Valley reaching 60-100 micrograms per cubic meter. These maps detail pollution from Interstate 5 at the scale of a few kilometers. The maps are based on NASA's MODerate resolution Imaging Spectrometer (MODIS) data at circa 1 kilometer as processed with the Multi-Angle Implementation of Atmospheric Correction. The Bay Area Air Quality Management District has requested that we test our methods in their challenging environment characterized by multiple sub-basins defined by complex topography. Our tests suggest that nearly similar precision may be expected for wintertime conditions with high PM 2.5 . We note difficulties when measured PM 2.5 is less than 8-10 micrograms per cubic meter, but good relative precision when PM 2.5 rises above 20; i.e. in episodes of concern for morbidity and mortality. Our method stresses physically meaningful functions of MODIS-MAIAC (Multi-Angle Implementation of Atmospheric Correction)-derived AOD (Aerosol Optical Depth) and total water vapor column. A mixed-effects statistical model exploiting existing station data works powerfully to allow us daily AOT-to-PM 2.5 relationships that allow a calibration of the map. In those cases where water vapor and particles have generally similar surface sources, using the ratio of AOT / Column_water can improve the daily calibrations so as to reach our quoted precision. We briefly present some cartoon idealizations that explain this success and also the likely reasons that our mixed effects model (or "daily calibration") works; also when it should not work. The combined satellite/mixed-effects model works best for wintertime San Joaquin Valley episodes, where the meteorology of particle and H2O(v) dilution is quite appropriate. We extended and tested the methodology (a) for the Bay Area wintertime situations and (b) for smoke plume events (e.g. the October 2017 fire events of the Sonoma area). Our SJV work was evaluated using NASA's DISCOVER-AQ (Deriving Information on Surface conditions from Column and Vertically Resolved Observations Relevant to Air Quality) airborne measurements, and by season- long measurements in Fresno. If the composition and size distribution of the aerosols can be assessed for the regions we describe, retrievals should have improved accuracy.
Chatfield, Robert B.↗
GEONEX: Land Monitoring From a New Generation of Geostationary Satellite Sensors
The latest generation of geostationary satellites carry sensors such as ABI (Advanced Baseline Imager on GOES-16) and the AHI (Advanced Himawari Imager on Himawari) that closely mimic the spatial and spectral characteristics of Earth Observing System flagship MODIS for monitoring land surface conditions. More importantly they provide observations at 5-15 minute intervals. Such high frequency data offer exciting possibilities for producing robust estimates of land surface conditions by overcoming cloud cover, enabling studies of diurnally varying local-to-regional biosphere-atmosphere interactions, and operational decision-making in agriculture, forestry and disaster management. But the data come with challenges that need special attention. For instance, geostationary data feature changing sun angle at constant view for each pixel, which is reciprocal to sun-synchronous observations, and thus require careful adaptation of EOS algorithms. Our goal is to produce a set of land surface products from geostationary sensors by leveraging NASA's investments in EOS algorithms and in the data/compute facility NEX. The land surface variables of interest include atmospherically corrected surface reflectances, snow cover, vegetation indices and leaf area index (LAI)/fraction of photosynthetically absorbed radiation (FPAR), as well as land surface temperature and fires. In order to get ready to produce operational products over the US from GOES-16 starting 2018, we have utilized 18 months of data from Himawari AHI over Australia to test the production pipeline and the performance of various algorithms for our initial tests. The end-to-end processing pipeline consists of a suite of modules to (a) perform calibration and automatic georeference correction of the AHI L1b data, (b) adopt the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm to produce surface spectral reflectances along with compositing schemes and QA, and (c) modify relevant EOS retrieval algorithms (e.g., LAI and FPAR, GPP, etc.) for subsequent science product generation. Initial evaluation of Himawari AHI products against standard MODIS products indicate general agreement, suggesting that data from geostationary sensors can augment low earth orbit (LEO) satellite observations.
geostationary↗
Generating Land Surface Reflectance for the New Generation of Geostationary Satellite Sensors with the MAIAC Algorithm
The latest generation of geostationary satellite sensors, including the GOES-16/ABI and the Himawari 8/AHI, provide exciting capability to monitor land surface at very high temporal resolutions (5-15 minute intervals) and with spatial and spectral characteristics that mimic the Earth Observing System flagship MODIS. However, geostationary data feature changing sun angles at constant view geometry, which is almost reciprocal to sun-synchronous observations. Such a challenge needs to be carefully addressed before one can exploit the full potential of the new sources of data. Here we take on this challenge with Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, recently developed for accurate and globally robust applications like the MODIS Collection 6 re-processing. MAIAC first grids the top-of- atmosphere measurements to a fixed grid so that the spectral and physical signatures of each grid cell are stacked (“remembered”) over time and used to dramatically improve cloud/shadow/snow detection, which is by far the dominant error source in the remote sensing. It also exploits the changing sun-view geometry of the geostationary sensor to characterize surface BRDF with augmented angular resolution for accurate aerosol retrievals and atmospheric correction. The high temporal resolutions of the geostationary data indeed make the BRDF retrieval much simpler and more robust as compared with sun-synchronous sensors such as MODIS. As a prototype test for the geostationary-data processing pipeline on NASA Earth Exchange (GEONEX), we apply MAIAC to process 18 months of data from Himawari 8/AHI over Australia. We generate a suite of test results, including the input TOA reflectance and the output cloud mask, aerosol optical depth (AOD), and the atmospherically-corrected surface reflectance for a variety of geographic locations, terrain, and land cover types. Comparison with MODIS data indicates a general agreement between the retrieved surface reflectance products. Furthermore, the geostationary results satisfactorily capture the movement of clouds and variations in atmospheric dust/aerosol concentrations, suggesting that high quality land surface and vegetation datasets from the advanced geostationary sensors can help complement and improve the corresponding EOS products.
geostationary satellite sensors↗
The Potential Impact of Satellite-Retrieved Cloud Parameters on Ground-Level PM2.5 Mass and Composition
Satellite-retrieved aerosol optical properties have been extensively used to estimate ground-level fine particulate matter (PM2.5) concentrations in support of air pollution health effects research and air quality assessment at the urban to global scales. However, a large proportion, approximately 70%, of satellite observations of aerosols are missing as a result of cloud-cover, surface brightness, and snow-cover. The resulting PM2.5 estimates could therefore be biased due to this non-random data missingness. Cloud-cover in particular has the potential to impact ground-level PM2.5 concentrations through complex chemical and physical processes. We developed a series of statistical models using the Multi-Angle Implementation of Atmospheric Correction (MAIAC) aerosol product at 1 km resolution with information from the MODIS cloud product and meteorological information to investigate the extent to which cloud parameters and associated meteorological conditions impact ground-level aerosols at two urban sites in the US: Atlanta and San Francisco. We find that changes in temperature, wind speed, relative humidity, planetary boundary layer height, convective available potential energy, precipitation, cloud effective radius, cloud optical depth, and cloud emissivity are associated with changes in PM2.5 concentration and composition, and the changes differ by overpass time and cloud phase as well as between the San Francisco and Atlanta sites. A case-study at the San Francisco site confirmed that accounting for cloud-cover and associated meteorological conditions could substantially alter the spatial distribution of monthly ground-level PM2.5 concentrations.
MAIAC AO↗
Full-Coverage High-Resolution Daily PM(sub 2.5) Estimation using MAIAC AOD in the Yangtze River Delta of China
Satellite aerosol optical depth (AOD) has been used to assess population exposure to fine particulate matter (PM (sub 2.5)). The emerging high-resolution satellite aerosol product, Multi-Angle Implementation of Atmospheric Correction(MAIAC), provides a valuable opportunity to characterize local-scale PM(sub 2.5) at 1-km resolution. However, non-random missing AOD due to cloud snow cover or high surface reflectance makes this task challenging. Previous studies filled the data gap by spatially interpolating neighboring PM(sub 2.5) measurements or predictions. This strategy ignored the effect of cloud cover on aerosol loadings and has been shown to exhibit poor performance when monitoring stations are sparse or when there is seasonal large-scale missngness. Using the Yangtze River Delta of China as an example, we present a Multiple Imputation (MI) method that combines the MAIAC high-resolution satellite retrievals with chemical transport model (CTM) simulations to fill missing AOD. A two-stage statistical model driven by gap-filled AOD, meteorology and land use information was then fitted to estimate daily ground PM(sub 2.5) concentrations in 2013 and 2014 at 1 km resolution with complete coverage in space and time. The daily MI models have an average R(exp 2) of 0.77, with an inter-quartile range of 0.71 to 0.82 across days. The overall Ml model 10-fold cross-validation R(exp 2) (root mean square error) were 0.81 (25 gm(exp 3)) and 0.73 (18 gm(exp 3)) for year 2013 and 2014, respectively. Predictions with only observational AOD or only imputed AOD showed similar accuracy.Comparing with previous gap-filling methods, our MI method presented in this study performed bette rwith higher coverage, higher accuracy, and the ability to fill missing PM(sub 2.5) predictions without ground PM(sub 2.5) measurements. This method can provide reliable PM(sub 2.5)predictions with complete coverage that can reduce biasin exposure assessment in air pollution and health studies.
PM2.5↗
Daily Kilometer-Scale MODIS Satellite Maps of PM2.5 Describe Wintertime Episodes
The San Joaquin Valley (SJV) suffers from severe health-endangering episodes of PM2.5 aerosol loadings in wintertime; episodes last approximately 5 days and differ in geographical distribution and composition. PM2.5 stations are scattered; consequently the use of remote sensing to map variable regional patterns of these varying respirable aerosol concentrations is desirable. High-precision AOT retrievals can capture column particulate loading. However,PM2.5 mapping is challenging due to several reasons: particularly thin mixed layers (ML) and thus relatively low aerosol optical thickness (AOT) close to current measurement limits, variable and a typical composition of the aerosols, and complex surface bidirectional reflectance. However, the West does present some advantages in analysis. Air basins are isolated from long-distance transport, and experience predominant strong meteorological subsidence. Thus these Western basin regions have fewer problematic cases of overriding aerosol layers detached from the surface. To counter such local overriding, Chu et al. have described an approach for the Eastern US, and He et al have described a synoptic classification approach useful in Shanghai. The Bay Area Air Quality Management District (BAAQMD) expands our experience with the use of AOT, with lower PM2.5 and several isolated sub-basins. We have prepared daily maps of episodes in each region. We present also a sequence of increasingly detailed statistical models, AOT initially appears to contribute little information; however, inclusion of weather information reveals its utility. Lyapustin and Wang's MultiAngle Implementation of Atmospheric Correction (MAIAC) retrieval for AOT provided the most useful operational remote sensing information for these regions. It provides high (1-km) spatial resolution maps and a high percentage of availability. Empirical regression methods have found that random effects regression models (aka mixed effects models, ME) employing AOT provide good estimates of ground PM2.5 concentrations.Here, we attempt to extend these methods and evaluate the usefulness of AOT with greater physical analysis, based on DISCOVER-AQ4 experience.
Chatfield, Robert B.↗
Consistency of Vegetation Index Seasonality Across the Amazon Rainforest
Vegetation indices (VIs) calculated from remotely sensed reflectance are widely used tools for characterizing the extent and status of vegetated areas. Recently, however, their capability to monitor the Amazon forest phenology has been intensely scrutinized. In this study, we analyze the consistency of VIs seasonal patterns obtained from two MODIS products: the Collection 5 BRDF product (MCD43) and the Multi-Angle Implementation of Atmospheric Correction algorithm (MAIAC). The spatio-temporal patterns of the VIs were also compared with field measured leaf litterfall, gross ecosystem productivity and active microwave data. Our results show that significant seasonal patterns are observed in all VIs after the removal of view-illumination effects and cloud contamination. However, we demonstrate inconsistencies in the characteristics of seasonal patterns between different VIs and MODIS products. We demonstrate that differences in the original reflectance band values form a major source of discrepancy between MODIS VI products. The MAIAC atmospheric correction algorithm significantly reduces noise signals in the red and blue bands. Another important source of discrepancy is caused by differences in the availability of clear-sky data, as the MAIAC product allows increased availability of valid pixels in the equatorial Amazon. Finally, differences in VIs seasonal patterns were also caused by MODIS collection 5 calibration degradation. The correlation of remote sensing and field data also varied spatially, leading to different temporal offsets between VIs, active microwave and field measured data. We conclude that recent improvements in the MAIAC product have led to changes in the characteristics of spatio-temporal patterns of VIs seasonality across the Amazon forest, when compared to the MCD43 product. Nevertheless, despite improved quality and reduced uncertainties in the MAIAC product, a robust biophysical interpretation of VIs seasonality is still missing.
MAIAC↗
Amazon Forests Response to Droughts: A Perspective from the MAIAC Product
Amazon forests experienced two severe droughts at the beginning of the 21st century: one in 2005 and the other in 2010. How Amazon forests responded to these droughts is critical for the future of the Earth's climate system. It is only possible to assess Amazon forests' response to the droughts in large areal extent through satellite remote sensing. Here, we used the Multi-Angle Implementation of Atmospheric Correction (MAIAC) Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation index (VI) data to assess Amazon forests' response to droughts, and compared the results with those from the standard (Collection 5 and Collection 6) MODIS VI data. Overall, the MAIAC data reveal more realistic Amazon forests inter-annual greenness dynamics than the standard MODIS data. Our results from the MAIAC data suggest that: (1) the droughts decreased the greenness (i.e., photosynthetic activity) of Amazon forests; (2) the Amazon wet season precipitation reduction induced by El Niño events could also lead to reduced photosynthetic activity of Amazon forests; and (3) in the subsequent year after the water stresses, the greenness of Amazon forests recovered from the preceding decreases. However, as previous research shows droughts cause Amazon forests to reduce investment in tissue maintenance and defense, it is not clear whether the photosynthesis of Amazon forests will continue to recover after future water stresses, because of the accumulated damages caused by the droughts.
remote sensing↗
Estimation of Crop Gross Primary Production (GPP): Do Scaled (MODIS) Vegetation Indices Improve Performance? - 2
Satellite remote sensing estimates of Gross Primary Production (GPP) have routinely been made using spectral Vegetation Indices (VIs) over the past two decades. The Normalized Difference Vegetation Index (NDVI), the Enhanced Vegetation Index (EVI), the green band Wide Dynamic Range Vegetation Index (WDRVIgreen), and the green band Chlorophyll Index (CIgreen) have been employed to estimate GPP under the assumption that GPP is proportional to the product of VI and photosynthetically active radiation (PAR) (where VI is one of four VIs: NDVI, EVI, WDRVIgreen, or CIgreen). However, the empirical regressions between VI*PAR and GPP measured locally at flux towers do not pass through the origin (i.e., the zero X-Y value for regressions). Therefore they are somewhat difficult to interpret and apply. This study investigates (1) what are the scaling factors and offsets (i.e., regression slopes and intercepts) between the fraction of PAR absorbed by chlorophyll of a canopy (fAPARchl) and the VIs, and (2) whether the scaled VIs developed in (1) can eliminate the deficiency and improve the accuracy of GPP estimates. Three AmeriFlux maize and soybean fields were selected for this study, two of which are irrigated and one is rainfed. The four VIs and fAPARchl of the fields were computed with the MODerate resolution Imaging Spectroradiometer (MODIS) satellite images. The GPP estimation performance for the scaled VIs was compared to results obtained with the original VIs and evaluated with standard statistics: the coefficient of determination (R2), the root mean square error (RMSE), and the coefficient of variation (CV). Overall, the scaled EVI obtained the best performance. The performance of the scaled NDVI, EVI and WDRVIgreen was improved across sites, crop types and soil/background wetness conditions. The scaled CIgreen did not improve results, compared to the original CIgreen. The scaled green band indices (WDRVIgreen, CIgreen) did not exhibit superior performance to either the scaled EVI or NDVI in estimating crop daily GPP at these agricultural fields. The scaled VIs are more physiologically meaningful than original un-scaled VIs, but scaling factors and offsets may vary across crop types and surface conditions.
Vegetation Index↗
A New Hybrid Spatio-temporal Model for Estimating Daily Multi-year PM2.5 Concentrations Across Northeastern USA Using High Resolution Aerosol Optical Depth Data
The use of satellite-based aerosol optical depth (AOD) to estimate fine particulate matter PM(sub 2.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. 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 PM(sub 2.5) at a 1X 1 km 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 X 1 km grid predictions. We used mixed models regressing PM(sub 2.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. Our model performance was excellent (mean out-of-sample R(sup 2) = 0.88). The spatial and temporal components of the out-of-sample results also presented very good fits to the withheld data (R(sup 2) = 0.87, R(sup)2 = 0.87). In addition, our results revealed very little bias in the predicted concentrations (Slope of predictions versus withheld observations = 0.99). 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.
Air pollution↗
Estimation of Crop Gross Primary Production (GPP): I. Impact of MODIS Observation Footprint and Impact of Vegetation BRDF Characteristics
Accurate estimation of gross primary production (GPP) is essential for carbon cycle and climate change studies. Three AmeriFlux crop sites of maize and soybean were selected for this study. Two of the sites were irrigated and the other one was rainfed. The normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), the green band chlorophyll index (CIgreen), and the green band wide dynamic range vegetation index (WDRVIgreen) were computed from the moderate resolution imaging spectroradiometer (MODIS) surface reflectance data. We examined the impacts of the MODIS observation footprint and the vegetation bidirectional reflectance distribution function (BRDF) on crop daily GPP estimation with the four spectral vegetation indices (VIs - NDVI, EVI, WDRVIgreen and CIgreen) where GPP was predicted with two linear models, with and without offset: GPP = a × VI × PAR and GPP = a × VI × PAR + b. Model performance was evaluated with coefficient of determination (R2), root mean square error (RMSE), and coefficient of variation (CV). The MODIS data were filtered into four categories and four experiments were conducted to assess the impacts. The first experiment included all observations. The second experiment only included observations with view zenith angle (VZA) ≤ 35◦ to constrain growth of the footprint size,which achieved a better grid cell match with the agricultural fields. The third experiment included only forward scatter observations with VZA ≤ 35◦. The fourth experiment included only backscatter observations with VZA ≤ 35◦. Overall, the EVI yielded the most consistently strong relationships to daily GPP under all examined conditions. The model GPP = a × VI × PAR + b had better performance than the model GPP = a × VI × PAR, and the offset was significant for most cases. Better performance was obtained for the irrigated field than its counterpart rainfed field. Comparison of experiment 2 vs. experiment 1 was used to examine the observation footprint impact whereas comparison of experiment 4 vs. experiment 3 was used to examine the BRDF impact. Changes in R2, RMSE,CV and changes in model coefficients "a" and "b" (experiment 2 vs. experiment 1; and experiment 4 vs. experiment 3) were indicators of the impacts. The second experiment produced better performance than the first experiment, increasing R2 (↑0.13) and reducing RMSE (↓0.68 g C m−2 d−1) and CV (↓9%). For each VI, the slope of GPP = a × VI × PAR in the second experiment for each crop type changed little while the slope and intercept of GPP = a × VI × PAR + b varied field by field. The CIgreen was least affected by the MODIS observation footprint in estimating crop daily GPP (R2, ↑0.08; RMSE, ↓0.42 g C m−2 d−1; and CV, ↓7%). Footprint most affected the NDVI (R2, ↑0.15; CV, ↓10%) and the EVI (RMSE, ↓0.84 g C m−2 d−1). The vegetation BRDF impact also caused variation of model performance and change of model coefficients. Significantly different slopes were obtained for forward vs. backscatter observations, especially for the CIgreen and the NDVI. Both the footprint impact and the BRDF impact varied with crop types, irrigation options, model options and VI options.
Chlorophyll↗
Impacts of Light Use Efficiency and fPAR Parameterization on Gross Primary Production Modeling
This study examines the impact of parameterization of two variables, light use efficiency (LUE) and the fraction of absorbed photosynthetically active radiation (fPAR or fAPAR), on gross primary production(GPP) modeling. Carbon sequestration by terrestrial plants is a key factor to a comprehensive under-standing of the carbon budget at global scale. In this context, accurate measurements and estimates of GPP will allow us to achieve improved carbon monitoring and to quantitatively assess impacts from cli-mate changes and human activities. Spaceborne remote sensing observations can provide a variety of land surface parameterizations for modeling photosynthetic activities at various spatial and temporal scales. This study utilizes a simple GPP model based on LUE concept and different land surface parameterizations to evaluate the model and monitor GPP. Two maize-soybean rotation fields in Nebraska, USA and the Bartlett Experimental Forest in New Hampshire, USA were selected for study. Tower-based eddy-covariance carbon exchange and PAR measurements were collected from the FLUXNET Synthesis Dataset. For the model parameterization, we utilized different values of LUE and the fPAR derived from various algorithms. We adapted the approach and parameters from the MODIS MOD17 Biome Properties Look-Up Table (BPLUT) to derive LUE. We also used a site-specific analytic approach with tower-based Net Ecosystem Exchange (NEE) and PAR to estimate maximum potential LUE (LUEmax) to derive LUE. For the fPAR parameter, the MODIS MOD15A2 fPAR product was used. We also utilized fAPAR chl, a parameter accounting for the fAPAR linked to the chlorophyll-containing canopy fraction. fAPAR chl was obtained by inversion of a radiative transfer model, which used the MODIS-based reflectances in bands 1-7 produced by Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm. fAPAR chl exhibited seasonal dynamics more similar with the flux tower based GPP than MOD15A2 fPAR, especially in the spring and fall at the agricultural sites. When using the MODIS MOD17-based parameters to estimate LUE, fAPAR chl generated better agreements with GPP (r2= 0.79-0.91) than MOD15A2 fPAR (r2= 0.57-0.84).However, underestimations of GPP were also observed, especially for the crop fields. When applying the site-specific LUE max value to estimate in situ LUE, the magnitude of estimated GPP was closer to in situ GPP; this method produced a slight overestimation for the MOD15A2 fPAR at the Bartlett forest. This study highlights the importance of accurate land surface parameterizations to achieve reliable carbon monitoring capabilities from remote sensing information.
grossprimary↗
Satellite Observed Widespread Decline in Mongolian Grasslands Largely Due to Overgrazing
The Mongolian Steppe is one of the largest remaining grassland ecosystems. Recent studies have reported widespread decline of vegetation across the steppe and about 70 percent of this ecosystem is now considered degraded. Among the scientific community there has been an active debate about whether the observed degradation is related to climate, or overgrazing, or both. Here, we employ a new atmospheric correction and cloud screening algorithm (MAIAC) to investigate trends in satellite observed vegetation phenology. We relate these trends to changes in climate and domestic animal populations. A series of harmonic functions is fitted to MODIS observed phenological curves to quantify seasonal and inter-annual changes in vegetation. Our results show a widespread decline (of about 12 percent on average) in MODIS observed NDVI across the country but particularly in the transition zone between grassland and the Gobi desert, where recent decline was as much as 40 percent below the 2002 mean NDVI. While we found considerable regional differences in the causes of landscape degradation, about 80 percent of the decline in NDVI could be attributed to increase in livestock. Changes in precipitation were able to explain about 30 percent of degradation across the country as a whole but up to 50 percent in areas with denser vegetation cover (p0.05). Temperature changes, while significant, played only a minor role (r20.10, p0.05). Our results suggest that the cumulative effect of overgrazing is a primary contributor to the degradation of the Mongolian steppe and is at least partially responsible for desertification reported in previous studies.
Mongolia↗