Toward More Accurate Modeling of Canopy Radiative Transfer and Leaf Electron Transport in Land Surface Modeling
Explore the source record for details and available documents.
Yujie Wang↗
Engineering topics
Yujie Wang
Publications and source records attributed to Yujie Wang.
Calibration of the SNPP and NOAA 20 VIIRS Sensors for Continuity of the MODIS Climate Data Records
Accurate long-term sensor calibration and periodic re-processing to ensure consistency and continuity of atmospheric, land and ocean geophysical retrievals from space within the mission period and across different missions is a major requirement of climate data records. In this work, we applied the Multi-Angle Implementation of Atmospheric Correction (MAIAC)-based vicarious calibration technique over Libya-4 desert site to perform calibration analysis of Visible Infrared Imaging Radiometer Suite (VIIRS) on Suomi National Polar-orbiting Partnership (SNPP) and NOAA-20 satellites. For both VIIRS sensors we characterized residual linear calibration trends and cross-calibrated both sensors to MODerate resolution Imaging Spectroradiometer (MODIS) Aqua regarded as a calibration standard. The relative spectral response (RSR) differences were accounted for using the German Aerospace Center (DLR) Earth Sensing Imaging Spectrometer (DESIS) hyperspectral surface reflectance data. Our results agree with independent vicarious calibration results of both the MODIS/VIIRS Characterization Support Team as well as the CERES Imager and Geostationary Calibration Group within estimated uncertainty of 1–2%. Analysis of MAIAC geophysical products with the new calibration shows a high level of agreement of MAIAC aerosol, surface reflectance and NDVI records between MODIS and VIIRS. Excluding high aerosol optical depth (AOD), all three sensors agree in AOD with mean difference (MD) less than 0.01 and residual mean squared difference rmsd ∼ 0.04. Spectral geometrically normalized surface reflectance agrees within rmsd of 0.003–0.005 in the visible and 0.01–0.012 at longer wavelengths. The residual surface reflectance differences are fully explained by differences in spectral filter functions. Finally, difference in NDVI is characterized by rmsd ∼ 0.02 and MD less than 0.003 for NDVI based on VIIRS imagery bands I1/I2 and less than 0.01 for NDVI based on VIIRS radiometric bands M5/M7. In practical sense, these numbers indicate consistency and continuity in MAIAC records ensuring the smooth transition from MODIS to VIIRS.
MAIAC↗
A Novel Atmospheric Correction Algorithm to Exploit the Diurnal Variability in Hypertemporal Geostationary Observations
This study developed a new atmospheric correction algorithm, GeoNEX-AC, that is independent from the traditional use of spectral band ratios but dedicated to exploiting information from the diurnal variability in the hypertemporal geostationary observations. The algorithm starts by evaluating smooth segments of the diurnal time series of the top-of-atmosphere (TOA) reflectance to identify clear-sky and snow-free observations. It then attempts to retrieve the Ross-Thick–Li-Sparse (RTLS) surface bi-directional reflectance distribution function (BRDF) parameters and the daily mean atmospheric optical depth (AOD) with an atmospheric radiative transfer model (RTM) to optimally simulate the observed diurnal variability in the clear-sky TOA reflectance. Once the initial RTLS parameters are retrieved after the algorithm’s burn-in period, they serve as the prior information to estimate the AOD levels for the following days and update the surface BRDF information with the new clear-sky observations. This process is iterated through the full time span of the observations, skipping only totally cloudy days or when surface snow is detected. We tested the algorithm over various Aerosol Robotic Network (AERONET) sites and the retrieved results well agree with the ground-based measurements. This study demonstrates that the high-frequency diurnal geostationary observations contain unique information that can help to address the atmospheric correction problem from new directions.
atmospheric correction↗
A Comparison of Multi-Angle Implementation of Atmospheric Correction and MOD09 Daily Surface Reflectance Products from MODIS
This study presents the first systematic comparison of MAIAC Collection 6 MCD19A1 daily surface reflectance (SR) product with standard MODIS SR (MOD/MYD09). The study was limited to four tiles located in mid-Atlantic United States (H11V05), Canada (H12V03),central Amazon (H11V09), and North-Eastern China (H27V05) and used over 5000 MODIS granules in 2018. Overall, there is a remarkable agreement between the best quality pixels of the two products, in particular in the Red and NIR bands. Over selected tiles, the evaluation found that MAIAC provides from 4 to 25% more high-quality retrievals thanMOD09 annually, with the largest difference in tropical regions, confirming results of the previous studies. The comparison of spectral characteristics showed a systematic MAIAC-MOD09 difference increasing from NIR to Blue, typical of biases of a Lambertian assumption in MOD09 algorithm. Over the North-Eastern China, MCD19A1 SR is found more stable at wide range of aerosol optical depth (AOD) variations, whereasMOD09 SR shows a consistent positive bias increasing with AOD and at shorter wavelengths. The observed SR differences can be attributed to differences in cloud detection, aerosol retrieval and in atmospheric correction which is performed using an accurate BRDF-coupled radiative transfer model in MAIAC and a Lambertian surface model in MOD09. While this study is not representative of the global performance because of its limited geographical coverage, it should help the land community to better understand the differences between the two products.
MAIAC↗
Evaluation and Intercomparison of Wildfire Smoke Forecasts from Multiple Modeling Systems for the 2019 Williams Flats Fire
Wildfire smoke is one of the most significant concerns of human and environmental health, associated with its substantial impacts on air quality, weather, and climate. However, biomass burning emissions and smoke remain among the largest sources of uncertainties in air quality forecasts. In this study, we evaluate the smoke emissions and plume forecasts from 12 state-of-the-art air quality forecasting systems during the Williams Flats fire in Washington State, US, August 2019, which was intensively observed during the Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) field campaign. Model forecasts with lead times within 1 d are intercompared under the same framework based on observations from multiple platforms to reveal their performance regarding fire emissions, aerosol optical depth (AOD), surface PM2.5, plume injection, and surface PM2.5 to AOD ratio. The comparison of smoke organic carbon (OC) emissions suggests a large range of daily totals among the models, with a factor of 20 to 50. Limited representations of the diurnal patterns and day-to-day variations of emissions highlight the need to incorporate new methodologies to predict the temporal evolution and reduce uncertainty of smoke emission estimates. The evaluation of smoke AOD (sAOD) forecasts suggests overall underpredictions in both the magnitude and smoke plume area for nearly all models, although the high-resolution models have a better representation of the fine-scale structures of smoke plumes. The models driven by fire radiative power (FRP)-based fire emissions or assimilating satellite AOD data generally outperform the others. Additionally, limitations of the persistence assumption used when predicting smoke emissions are revealed by substantial underpredictions of sAOD on 8 August 2019, mainly over the transported smoke plumes, owing to the underestimated emissions on 7 August. In contrast, the surface smoke PM2.5 (sPM2.5) forecasts show both positive and negative overall biases for these models, with most members presenting more considerable diurnal variations of sPM2.5. Overpredictions of sPM2.5 are found for the models driven by FRP-based emissions during nighttime, suggesting the necessity to improve vertical emission allocation within and above the planetary boundary layer (PBL). Smoke injection heights are further evaluated using the NASA Langley Research Center's Differential Absorption High Spectral Resolution Lidar (DIAL-HSRL) data collected during the flight observations. As the fire became stronger over 3–8 August, the plume height became deeper, with a day-to-day range of about 2–9 km a.g.l. However, narrower ranges are found for all models, with a tendency of overpredicting the plume heights for the shallower injection transects and underpredicting for the days showing deeper injections. The misrepresented plume injection heights lead to inaccurate vertical plume allocations along the transects corresponding to transported smoke that is 1 d old. Discrepancies in model performance for surface PM2.5 and AOD are further suggested by the evaluation of their ratio, which cannot be compensated for by solely adjusting the smoke emissions but are more attributable to model representations of plume injections, besides other possible factors including the evolution of PBL depths and aerosol optical property assumptions. By consolidating multiple forecast systems, these results provide strategic insight on pathways to improve smoke forecasts.
AOD↗
Retrievals of Aerosol Optical Depth and Spectral Absorption from DSCOVR EPIC
A new algorithm is described for joint retrievals of the aerosol optical depth and spectral absorption from EPIC observations in the UV—Vis spectral range. The retrievals are illustrated on examples of the wildfire smoke events over North America, and dust storms over greater Sahara region in 2018. An initial evaluation of single scattering albedo (SSA) at 443 nm over these regions shows a good agreement with AERONET data, generally within the uncertainty of AERONET SSA of ± 0.03. A particularly good agreement is achieved for dust with R∼0.62, rmse∼0.02, negligible bias, and 85% points within the expected error. This new capability is part of version 2 MAIAC EPIC algorithm. The v2 algorithm has recently completed reprocessing of the EPIC record covering the period of 2015–2020.
EPIC↗
The AOD Sensitivity Comparison between MODIS Multi-Angle Implementation of Atmospheric Correction (MAIAC) and Standard MODIS Surface Reflectance
This study compares the Multi-Angle Implementation of Atmospheric Correction (MAIAC) and standard MODIS surface reflectance (SR) products (MOD09)for various aerosol optical depth (AOD) levels. Data from one MODIS tile in northern China for the whole year 2018 is compared. The results reveal that MAIAC SR has good stability for the full range of AOD for which SR is reported, whereas MOD09 SR shows increasing high bias with AOD increases. The cross-comparison of spectral characteristics between MAIAC SR and MOD09 SR shows a systematic MAIAC-MOD09 difference increasing from NIR to Blue. This pattern is consistent with the bias caused by Lambertian assumption used in MOD09 SR algorithm.
MAIAC↗
An AeroCom–AeroSat study: intercomparison of satellite AOD datasets for aerosol model evaluation
To better understand and characterize current uncertainties in the important observational constraint of climate models of aerosol optical depth (AOD), we evaluate and intercompare 14 satellite products, representing nine different retrieval algorithm families using observations from five different sensors on six different platforms. The satellite products (super-observations consisting of 1°×1° daily aggregated retrievals drawn from the years 2006, 2008 and 2010) are evaluated with AErosol RObotic NETwork (AERONET) and Maritime Aerosol Network (MAN) data. Results show that different products exhibit different regionally varying biases (both under- and overestimates) that may reach ±50 %, although a typical bias would be 15 %–25 % (depending on the product). In addition to these biases, the products exhibit random errors that can be 1.6 to 3 times as large. Most products show similar performance, although there are a few exceptions with either larger biases or larger random errors. The intercomparison of satellite products extends this analysis and provides spatial context to it. In particular, we show that aggregated satellite AOD agrees much better than the spatial coverage (often driven by cloud masks) within the 1°×1° grid cells. Up to ∼50 % of the difference between satellite AOD is attributed to cloud contamination. The diversity in AOD products shows clear spatial patterns and varies from 10 % (parts of the ocean) to 100 % (central Asia and Australia). More importantly, we show that the diversity may be used as an indication of AOD uncertainty, at least for the better performing products. This provides modellers with a global map of expected AOD uncertainty in satellite products, allows assessment of products away from AERONET sites, can provide guidance for future AERONET locations and offers suggestions for product improvements. We account for statistical and sampling noise in our analyses. Sampling noise, variations due to the evaluation of different subsets of the data, causes important changes in error metrics. The consequences of this noise term for product evaluation are discussed.
Aerosol↗
Evaluating impacts of snow, surface water, soil and vegetation on empirical vegetation and snow indices for the Utqiaġvik tundra ecosystem in Alaska with the LVS3 model
Satellite observations for the Arctic and boreal region may contain information of vegetation, soil, snow, snowmelt, and/or other surface water bodies. We investigated the impacts of vegetation, soil, snow and surface water on empirical vegetation/snow indices on a tundra ecosystem area located around Utqiaġvik (formerly Barrow) of Alaska with the Moderate Resolution Imaging Spectrometer (MODIS) images in 2001–2014. Empirical vegetation indices, such as normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), the index of near infrared of vegetation (NIRv), and modified EVI (EVI2), have been used to monitor vegetation. Normalized difference snow index (NDSI) has been widely applied to monitor snow. The vegetation cover fraction (VGCF), the soil cover fraction (SOILCF), the snow cover fraction (SNOWCF), the surface water body cover fraction (WaterBodyCF), the fractional absorption of photosynthetically active radiation (PAR) by vegetation chlorophyll (fAPARchl), the fractional absorption of PAR by non-chlorophyll components of the vegetation (fAPARnon-chl), and the fractional absorption of PAR by the entire canopy (fAPARcanopy) are retrieved with the MODIS images and a coupled Leaf-Vegetation-Soil-Snow-Surface water body radiative transfer model, LVS3. The vegetation indices (NDVI, EVI, EVI2 and NIRv) differ from VGCF, fAPARchl, fAPARnon-chl, and fAPARcanopy. In addition to vegetation, we find that soil, snow and surface water also have impacts on vegetation indices NDVI, EVI (EVI2), and NIRv. Presence of snow makes lower the observed values of NDVI, EVI2 and NIRv. After snowmelt is gone, the vegetation indices (NDVI, EVI, EVI2 and NIRv) linearly decrease with SOILCF and WaterBodyCF, and WaterBodyCF has stronger impacts on these vegetation indices than SOILCF. The relationship between EVI and snow is complicated. NDSI non-linearly increases with SNOWCF, but linearly increases with sum of SNOWCF and WaterBodyCF (sum = 0.5893 × NDSI +0.4342, R2 = 0.976). NDSI linearly decreases with VGCF, and the relationship between NDSI and SOILCF is complex. Retrievals of VGCF, fAPARchl, fAPARnon-chl and fAPARcanopy with the LVS3 model provide alternatives for vegetation monitoring and ecological modeling.
Qingyuan Zhang↗
First Provisional Land Surface Reflectance Product from Geostationary Satellite Himawari-8 AHI
A provisional surface reflectance (SR) product from the Advanced Himawari Imager (AHI) on-board the new generation geostationary satellite (Himawari-8) covering the period between July 2015 and December 2018 is made available to the scientific community. The Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm is used in conjunction with time series Himawari-8 AHI observations to generate 1-km gridded and tiled land SR every 10 minutes during day time. This Himawari-8 AHI SR product includes retrieved atmospheric properties (e.g., aerosol optical depth at 0.47μm and 0.51μm), spectral surface reflectance (AHI bands 1–6), parameters of the RTLS BRDF model, and quality assurance flags. Product evaluation shows that Himawari-8 AHI data on average yielded 35% more cloud-free, valid pixels in a single day when compared to available data from the low earth orbit (LEO) satellites Terra/Aqua with MODIS sensor. Comparisons of Himawari-8 AHI SR against corresponding MODIS SR products (MCD19A1) over a variety of land cover types with the similar viewing geometry show high consistency between them, with correlation coefficients (r) being 0.94 and 0.99 for red and NIR bands, respectively. The high-frequency geostationary data are expected to facilitate studies of ecosystems on daily to diurnal time scales, complementing observations from networks such as the FLUXNET.
Himawari-8 AHI↗
MAIAC Thermal Technique for Smoke Injection Height from MODIS
We present a new algorithm to derive smoke plume height (Ha) using thermal contrast from the rising mixture of aerosol and emitted gases in the MODIS 11m channel. Validation shows good agreement with wind-corrected MISR-MINX values, with about 60% of the MODIS Terra thermal retrievals within 500m of MISR Ha, and 450m lower on average. The bias is expected because the thermal technique represents an effective rather than a top plume height from MISR MINX. Comparison of MODIS Aqua retrievals with CALIOP CALIPSO shows similar statistics, with standard deviation of 458m for the mean plume height and 216m lower on average. Ha is part of the MAIAC MODIS Collection 6 suite of products (MCD19), accessible via the Land Product Distributed Active Archive Center (LP DAAC). Aerosol injection height is reported in the daily MAIAC atmospheric product MCD19A2 along with the cloud mask, column water vapor, aerosol optical depth, AOD uncertainty and aerosol type, at 1km resolution on global Sinusoidal grid. Despite some limitations, the vastly increased coverage from MODIS observations makes it a valuable dataset complementing the established MISR and CALIOP products.
atmosphere↗
Impacts of Snow and Cloud Covers on Satellite-Derived PM 2.5 Levels
Satellite aerosol optical depth (AOD) has been widely employed to evaluate ground fine particle (PM 2.5 ) levels, whereas snow/cloud covers often lead to a large proportion of non-random missing AOD. As a result, the fully covered and unbiased PM 2.5 estimates will be hard to generate. Among the current approaches to deal with the data gap issue, few have considered the cloud-AOD relationship and none of them have considered the snow-AOD relationship. This study examined the impacts of snow and cloud covers on AOD and PM 2.5 and made full-coverage PM 2.5 predictions with the consideration of these impacts. To estimate the missing AOD, daily gap-filling models with snow/cloud fractions and meteorological covariates were developed using the random forest algorithm. By using these models in New York State, a daily AOD data set with a 1-km resolution was generated with a complete coverage. The“out-of-bag” R 2 of the gap-filling models averaged 0.93 with an interquartile range from 0.90 to 0.95. Subsequently, a random forest-based PM 2.5 prediction model with the gap-filled AOD and covariates was built to predict fully covered PM 2.5 estimates. A ten-fold cross-validation for the prediction model showed a good performance with an R 2 of 0.82. In the gap-filling models, the snow fraction was of higher significance in the snow season compared with the rest of the year. The prediction models fitted with/without the snow fraction also suggested the discernible changes in PM 2.5 patterns, further confirming the significance of this parameter. Compared with the methods without considering snow and cloud covers, our PM 2.5 prediction surfaces showed more spatial details and reflected small-scale terrain-driven PM 2.5 patterns. The proposed methods can be generalized to the areas with extensive snow/cloud covers and large proportions of missing satellite AOD for predicting PM 2.5 levels with high resolutions and complete coverage.
AOD↗
An Example of Aerosol Pattern Variability Over Bright Surface Using High Resolution MODIS MAIAC: The Eastern and Western Areas of the Dead Sea and Environs
The extreme rate of evaporation of the Dead Sea (DS) has serious implications for the surrounding area, including atmospheric conditions. This study analyzes the aerosol properties over the western and eastern parts of the DS during the year 2013, using MAIAC (Multi-Angle Implementation of Atmospheric Correction) for MODIS, which retrieves aerosol optical depth (AOD) data at a resolution of 1km. The main goal of the study is to evaluate MAIAC over the study area and determine, for the first time, theprevailing aerosol spatial patterns. First, the MAIAC-derived AOD data was compared with data from three nearby AERONET sites (Nes Ziona - an urban site, and Sede Boker and Masada - two arid sites), and with the conventional Dark Target (DT) and Deep Blue (DB) retrievals for the same days and locations, on a monthly basis throughout 2013. For the urban site, the correlation coefficient (r) for DTDB products showed better performance than MAIAC (r = 0.80, 0.75, and 0.64 respectively) year-round. However, in the arid zones, MAIAC showed better correspondence to AERONET sites than the conventional retrievals (r = 0.58 - 0.60 and 0.48 - 0.50 respectively). We investigated the difference in AOD levels, and its variability, between the Dead Sea coasts on a seasonal basis and calculated monthlyseasonal AOD averages for presenting AOD patterns over arid zones. Thus, we demonstrated that aerosol concentrations show a strong preference for the western coast, particularly during the summer season. This preference, is most likely a result of local anthropogenic emissions combined with the typical seasonal synoptic conditions, the Mediterranean Sea breeze, and the region complex topography. Our results also indicate that a large industrial zone showed higher AOD levels compared to an adjacent reference-site, i.e., 13 percent during the winter season.
AOD↗
Seasonal Monitoring and Estimation of Regional Aerosol Distribution over Po Valley, Northern Italy, Using a High-Resolution MAIAC Product
In this work, the new 1-km-resolved Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm is employed to characterize seasonal AOD-PM10 correlations over northern Italy. The accuracy of the new dataset is assessed versus the widely used Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 5.1 Aerosol Optical Depth (AOD) data, retrieved at 0.55 microns with spatial resolution of 10 km (MYD04). We focused on evaluating the ability of these two products to characterize both temporal and spatial distributions of aerosols within urban and suburban areas. Ground PM10 measurements were obtained from 73 of the Italian Regional Agency for Environmental Protection (ARPA) monitoring stations, spread across northern Italy, for a three-year period from 2010 to 2012. The Po Valley area (northern Italy) was chosen as the study domain because of severe urban air pollution, resulting from the highest population and industrial manufacturing density in the country, being located in a valley where two surrounding mountain chains favor the stagnation of pollutants. We found that the global correlations between PM10 and AOD are R(sup 2) = 0.83 and R(sup 2) = 0.44 for MYD04_L2 and for MAIAC, respectively, suggesting for a greater sensitiveness of the high-resolution product to small-scale deviations. However, the introduction of Relative Humidity (RH) and Planetary Boundary Layer (PBL) depth corrections gave a significant improvement to the PM AOD correlation, which led to similar performance: R(sup 2) = 0.96 for MODIS and R(sup 2) = 0.95 for MAIAC. Furthermore, the introduction of the PBL information in the corrected AOD values was found to be crucial in order to capture the clear seasonal cycle shown by measured PM10 values. The study allowed us to define four seasonal linear correlations that estimate PM10 concentrations satisfactorily from the remotely sensed MAIAC AOD retrieval. Overall, the results show that the high resolution provided by MAIAC retrieval data is much more relevant than 10km MODIS data to characterize PM10 in this region of Italy which has a pretty limited geographical domain, but a broad variety of land usages and consequent particulate concentrations.
MODIS↗
Spatiotemporal Prediction of Fine Particulate Matter Using High-Resolution Satellite Images in the Southeastern US 2003-2011
Numerous studies have demonstrated that fine particulate matter (PM(sub 2.5), particles smaller than 2.5 micrometers in aerodynamic diameter) is associated with adverse health outcomes. The use of ground monitoring stations of PM(sub 2.5) to assess personal exposure, however, induces measurement error. Land-use regression provides spatially resolved predictions but land-use terms do not vary temporally. Meanwhile, the advent of satellite-retrieved aerosol optical depth (AOD) products have made possible to predict the spatial and temporal patterns of PM(sub 2.5) exposures. In this paper, we used AOD data with other PM(sub 2.5) variables, such as meteorological variables, land-use regression, and spatial smoothing to predict daily concentrations of PM(sub 2.5) at a 1 sq km resolution of the Southeastern United States including the seven states of Georgia, North Carolina, South Carolina, Alabama, Tennessee, Mississippi, and Florida for the years from 2003 to 2011. We divided the study area into three regions and applied separate mixed-effect models to calibrate AOD using ground PM(sub 2.5) measurements and other spatiotemporal predictors. Using 10-fold cross-validation, we obtained out of sample R2 values of 0.77, 0.81, and 0.70 with the square root of the mean squared prediction errors of 2.89, 2.51, and 2.82 cu micrograms for regions 1, 2, and 3, respectively. The slopes of the relationships between predicted PM2.5 and held out measurements were approximately 1 indicating no bias between the observed and modeled PM(sub 2.5) concentrations. Predictions can be used in epidemiological studies investigating the effects of both acute and chronic exposures to PM(sub 2.5). Our model results will also extend the existing studies on PM(sub 2.5) which have mostly focused on urban areas because of the paucity of monitors in rural areas.
PM2.5↗
Estimation of Crop Gross Primary Production (GPP): Fapar(sub Chl) Versus MOD15A2 FPAR
Photosynthesis (PSN) is a pigment level process in which antenna pigments (predominately chlorophylls) in chloroplasts absorb photosynthetically active radiation (PAR) for the photochemical process. PAR absorbed by foliar non-photosynthetic components is not used for PSN. The fraction of PAR absorbed (fAPAR) by a canopy/vegetation (i.e., fAPAR(sub canopy)) derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) images, referred to as MOD15A2 FPAR, has been used to compute absorbed PAR (APAR) for PSN (APARPSN) which is utilized to produce the standard MODIS gross primary production (GPP) product, referred to as MOD17A2 GPP. In this study, the fraction of PAR absorbed by chlorophyll throughout the canopy (fAPAR(sub chl)) was retrieved from MODIS images for three AmeriFlux crop fields in Nebraska. There are few studies in the literature that compare the performance of MOD15A2 FPAR versus fAPAR(sub chl) in GPP estimation. In our study MOD15A2 FPAR and the retrieved fAPAR(sub chl) were compared with field fAPAR(sub canopy) and the fraction of PAR absorbed by green leaves of the vegetation (fAPAR(sub green)). MOD15A2 FPAR overestimated field fAPAR(sub canopy) in spring and in fall, and underestimated field fAPAR(sub canopy) in midsummer whereas fAPAR(sub chl) correctly captured the seasonal phenology. The retrieved fAPAR(sub chl) agreed well with field fAPAR(sub green) at early crop growth stage in June, and was less than field fAPAR(sub green) in late July, August and September. GPP estimates with fAPAR(sub chl) and with MOD15A2 FPAR were compared to tower flux GPP. GPP simulated with fAPAR(sub chl) was corroborated with tower flux GPP. Improvements in crop GPP estimation were achieved by replacing MOD15A2 FPAR with fAPAR(sub chl) which also reduced uncertainties of crop GPP estimates by 1.12-2.37 g C m(sup −2) d(sup −1).
GPM↗
Uncertainty Analysis of the GeoNEX Top-of-Atmospheric Reflectance Products Generated from the Third-Generation Geostationary Satellite Sensors
The GeoNEX (Geostationary-NASA Earth eXchange) Level-1G products consist of top-of-atmosphere (TOA) bi-directional reflectance factor (BRF) and brightness temperature generated with data streams from the latest geostationary (GEO) sensors including GOES-16/17 ABI, Himawari-8/9 AHI, and GK-2A AMI on a global tiled common grid (60oN-60o and 180oW-180oE) in geographic coordinates. With their 16 spectral bands, 0.01o/0.02o nadir spatial resolution, and 10-minute temporal resolutions, these products provide exciting opportunity to monitor Earth surface processes. However, the unique Sun-Target-Satellite geometry of geostationary sensors demands special attention in analyzing/interpreting these datasets. In this study we present a systematic analysis on the relationship between the radiometric uncertainties of the GeoNEX TOA reflectance and the corresponding solar/satellite zenith angles. We show that the signal-to-noise ratio (SNR) of the BRF are positively proportional to the square roots of the cosine of solar illuminating zenith angles. That is, the BRF data are noisier earlier in the morning or later in the afternoon than in the mid of the day. The cosine of satellite viewing zenith angles do not directly influence the SNR of the TOA BRF. However, they positively regulate the relative importance of the surface component in the TOA BRF. This means that variations in surface reflectance are more difficult to detect for pixels with larger view zenith angles, even when the SNR of the TOA BRF is the same. We are developing metrics to specify such illumination-view geometry related uncertainties in the GeoNEX L1G TOA BRF products so that this key information can be easily accessed by the user community.
Geostationary satellite↗
Generation of Land Surface Reflectance with Combined Geo-KOMPSAT-2A AMI and Himawari 8 AHI Observations
The latest generation of geostationary satellites has opened a new era of Earth observations with unprecedented spatiotemporal resolution and spectral range. Together with GOES 16/17 ABI, FY4-A AGRI, and Himawari-8 AHI, a new Korean geostationary satellite (Geo-KOMPSAT-2A AMI) has operationally collected a full-disk image in 16 channels every ten minutes since July 2019, allowing diurnal land surface monitoring over a large proportion of Asia and all of Oceania. Retrieving accurate surface reflectance (SR) over land from GK-2A/AMI is a challenging but high priority objective. One of the challenges is the absence of a spectral band in the 2.2 m SWIR range from AMI, which is required by many atmospheric correction algorithms to retrieve atmospheric aerosol properties. To remedy this issue, we adopt a strategy that combines concurrent GK-2A/AMI and Himawari 8/AHI observations in order to derive AMI SR. We have adapted the NASA Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm to process the data stream from Himawari 8/AHI. The advantages of the MAIAC algorithm is its capability to exploit the high temporal frequency and varying illumination geometry of the geostationary data for advanced cloud/snow detection, aerosol retrieval, and characterization of surface bidirectional reflectance properties. Leveraging the similarities of spectral bands and the sun-target-sensor geometry between AMI and AHI, we are able to create denser time series of observations and enhanced BRDF samples over most of the spatial coverage of AMI (and AHI). The combined stereo-type observations not only help derive SR for AMI but also enhance retrievals of the corresponding AHI surface products. We evaluate the resulting AMI SR using ground (AERONET) observations and corresponding MODIS products. Further, we discuss potential challenges in utilizing the geostationary satellite data for land surface monitoring.
geostationary satellite↗