Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Classification of Personnel Targets with Baggage Using Dual-band Radar
Next Article in Special Issue
Estimation of Surface NO2 Volume Mixing Ratio in Four Metropolitan Cities in Korea Using Multiple Regression Models with OMI and AIRS Data
Previous Article in Journal
Gravitation-Based Edge Detection in Hyperspectral Images
Previous Article in Special Issue
Ground-Level NO2 Concentrations over China Inferred from the Satellite OMI and CMAQ Model Simulations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MODIS Retrieval of Aerosol Optical Depth over Turbid Coastal Water

1
Department of Chemical and Biochemical Engineering, The University of Iowa, Iowa City, IA 52242, USA
2
Center of Global and Regional Environmental Research, The University of Iowa, Iowa City, IA 52242, USA
3
Interdisciplinary Graduate Program in Informatics, The University of Iowa, Iowa City, IA 52242, USA
4
NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA
5
Marine Meteorology Division, Naval Research Laboratory, Monterey, CA 93943, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2017, 9(6), 595; https://doi.org/10.3390/rs9060595
Submission received: 23 March 2017 / Revised: 26 May 2017 / Accepted: 7 June 2017 / Published: 12 June 2017
(This article belongs to the Special Issue Remote Sensing of Atmospheric Pollution)

Abstract

:
We present a new approach to retrieve Aerosol Optical Depth (AOD) using the Moderate Resolution Imaging Spectroradiometer (MODIS) over the turbid coastal water. This approach supplements the operational Dark Target (DT) aerosol retrieval algorithm that currently does not conduct AOD retrieval in shallow waters that have visible sediments or sea-floor (i.e., Class 2 waters). Over the global coastal water regions in cloud-free conditions, coastal screening leads to ~20% unavailability of AOD retrievals. Here, we refine the MODIS DT algorithm by considering that water-leaving radiance at 2.1 μm to be negligible regardless of water turbidity, and therefore the 2.1 μm reflectance at the top of the atmosphere is sensitive to both change of fine-mode and coarse-mode AODs. By assuming that the aerosol single scattering properties over coastal turbid water are similar to those over the adjacent open-ocean pixels, the new algorithm can derive AOD over these shallow waters. The test algorithm yields ~18% more MODIS-AERONET collocated pairs for six AERONET stations in the coastal water regions. Furthermore, comparison of the new retrieval with these AERONET observations show that the new AOD retrievals have equivalent or better accuracy than those retrieved by the MODIS operational algorithm’s over coastal land and non-turbid coastal water product. Combining the new retrievals with the existing MODIS operational retrievals yields an overall improvement of AOD over those coastal water regions. Most importantly, this refinement extends the spatial and temporal coverage of MODIS AOD retrievals over the coastal regions where 60% of human population resides. This expanded coverage is crucial for better understanding of impact of anthropogenic aerosol particles on coastal air quality and climate.

Graphical Abstract

1. Introduction

Aerosols are a colloidal system of particles suspended in the atmosphere, and have significant impacts on weather, climate, and human health [1,2,3]. Because of its global observational coverage, satellite remote sensing has a critical role in quantifying these impacts. Such satellite remote sensing is being used to retrieve aerosol properties such as the aerosol optical depth (AOD), along with a well-characterized uncertainty envelope, at high spatial resolution across the globe. Indeed, since the launch of the Moderate Resolution Imaging Spectroradiometer (MODIS) on the NASA Terra satellite in 1999, AOD retrievals derived by Dark Target (DT) algorithm [4] or Deep Blue (DB) algorithm (land only) [5,6] have been used widely in the research community [7].
One persistent challenge to aerosol remote sensing is retrieval of environmental properties over coastal (or littoral) waters. Variable ocean color, along with the presence of visible sediments and sea floor along coastlines, has strong spectral and spatial variability leading to poorly constrained lower boundary conditions for aerosol retrievals. Yet, retrieving AOD in coastal waters is a much-desired part of new systems characterizing air quality and aerosol radiative effects. As ~60% of the human population lives in coastal areas [8], it is crucial to expand the satellite-remote sensing datasets to include these areas. This study aims to address this observational gap by refining the MODIS DT-ocean algorithm to retrieve AOD over turbid coastal waters.
The MODIS DT-ocean algorithm uses top of atmosphere (TOA) reflectance in six wavelength bands, ranging from 0.55 to 2.1 μm, to simultaneously retrieve AOD and Fine Mode Fraction (FMF) based on a lookup table approach. The lookup table is constructed by assuming aerosol optical properties of four fine aerosol modes and five coarse aerosol modes, coupled with the boundary conditions of molecular scattering, and a rough ocean surface (glitter, whitecaps, foam). Additionally, the lookup table assumes zero water-leaving radiance at all bands except 0.55 μm (at which a fixed water-leaving reflectance of 0.005 is used) [4,9]. The retrieval process searches the lookup table to find the best combination of fine mode and coarse mode aerosol type (out of possibly 20 combinations) such that the AOD and FMF retrievals render the best match between observed and simulated (e.g., those from the lookup table) radiances. However, while the assumption of nearly-zero water-leaving radiance may be suitable for open ocean, it is clearly not applicable to turbid coastal waters. There, the water-leaving radiances can be contributed by reflection of shallow-water sea floor (particularly in blue wavelengths) and either suspended or dissolved particulate matter in the water, especially in green to red wavelengths such as 0.55, 0.66, and 0.86 μm [10]. Hence, as part of the DT-ocean algorithm, the turbid water pixels are masked and not considered for retrieval. The method used for such masking compares the TOA reflectance at 0.55 μm with the expected counterpart from the power-law fitting using the TOA reflectance at 0.47, 1.2, 1.6, and 2.1 μm; if a significant difference (larger than 0.01) is found, the corresponding pixel is masked out for AOD retrieval [8].
We analyze the MODIS AOD retrieval unavailability over the cloud-free conditions at both global and regional scale to reveal how often the AOD is not retrieved only because of water turbidity (and not other factors such as cloud cover). As shown in Figure 1, this data availability is near total over all open ocean and decreases dramatically (by 90–100%) toward coastlines. In a global average, the unavailability of AOD is ~20% over coastal water region which is labelled as shallow ocean (within 5 km of coastline or with water depth less than 50 m) in the MODIS geolocation product [11]. In other words, in these 20% of cloud-free cases, the AOD should have been retrieved if the water was not turbid.
Based on the principle that liquid water absorption increases in shortwave infrared (SWIR), we present a new approach that uses MODIS-measured radiance at 2.1 μm to retrieve AOD over turbid coastal water no matter whether the aerosol is fine or coarse mode dominated. In the spectral range of 0.55 to 2.1 μm, the transparency of pure water to sunlight decreases rapidly with the increase of wavelength; the penetration depth (at which light attenuation is 90%) is 41 m at 0.55 μm and drops to 0.001 m at 2.1 μm [8]. Hence, unless the water is shallower than 1 mm, the water-leaving radiance contributed from the sea floor, sediments, and other contaminants in the water at 2.1 μm is nearly negligible; this vastly simplifies the lower boundary condition of the retrieval.
We present the data and the new approach in Section 2 and Section 3, respectively, followed by description of retrieval results and the validation against measurements from six Aerosol Robotic Network (AERONET) coastal sites in Section 4. Discussions and conclusions are in Section 5 and Section 6, respectively.

2. Data

2.1. MODIS Data

MODIS is an earth-viewing sensor on board the Terra and Aqua polar-orbiting satellites, launched in December 1999 and May 2002, respectively. MODIS has 36 channels spanning the spectral range from 0.41 to 15 μm with spatial resolutions of 250 m (2 channels), 500 m (5 channels), and 1 km (29 channels). Its 2330-km swath width enables it to provide near-global converge daily. Terra and Aqua cross the equator from north to south (descending node) at approximately 10:30 a.m. local time and from south to north (ascending node) at approximately 1:30 p.m. local time, respectively.
Here, all MODIS data products are labeled as MxDNN, where x is substituted by O for Terra and Y for Aqua, respectively, and NN is the serial number of a specific product. King et al. [12] presented a general description of MODIS atmosphere data processing architecture and products. In this study, the MODIS-calibrated TOA reflectance product (MxD02) (http://mcst.gsfc.nasa.gov/content/l1b-documents), atmospheric profile product (MxD07), geolocation product (MxD03), and aerosol product (MxD04) [4] are used for retrieving aerosol optical depth over turbid coastal water. TOA reflectance of 2.1 μm with spatial resolution of 500 m from MxD02 is used to retrieve AOD through a lookup table (LUT) method. Gas absorption by water vapor and ozone column are corrected from the TOA reflectance with spatial resolution of 5 km from MxD07 product [4]. The MxD03 product has a spatial resolution of 1 km and is used to distinguish between water and land pixels. All other ancillary information from MxD04 aerosol product for each valid AOD retrieval (at spatial resolution of 10 km at nadir) are also used as input to the retrieval algorithm in this study; the information includes aerosol mode index, reflectance weighting parameter, and National Centers for Environmental Prediction (NCEP) analysis of wind speed (2 m above the surface). In addition, a cloud mask with a spatial resolution of 500 m from MxD04 aerosol product is also used to ensure only AOD during cloud-free conditions over the turbid coastal region is retrieved by this study.

2.2. AERONET Data

AOD measurements from the ground-based AERONET sun photometers are commonly used for validating MODIS retrievals. Here we use data from six coastal water sites for this study (Table 1), and their locations are marked with green dots in global map in Figure 1f. All AERONET sun photometers (SP) measure direct solar radiation at 0.44 μm, 0.67 μm, 0.87 μm, 0.94 μm and 1.02 μm, and these measurements are used to infer AOD through Beer-Lambert-Bouguer law with quality at Level 1.0 (unscreened), Level 1.5 (cloud screened), and Level 2.0 (cloud-screened and quality-assured) [13,14]. We evaluate MODIS AOD (both MxD04 product and retrieval of this study) at 0.55 μm against an AERONET counterpart that is derived through linearly interpolating AERONET AOD at 0.44 and 0.67 μm in the logarithm domain.
Dalma and MVCO are two sites over the ocean, with distance to coastline being 48 km and 29 km, respectively, while the rest of the sites are over land within 6 km from the coastline. According to our analysis, most AERONET sites are more than 10 km away from the coast lines, and lack dedicated long-term continuous measurements of AOD over the turbid coastal water. Here, we use the SP AOD data from August 2015, December 2015, May 2016, August 2004, March 2014, and September 2004 at MVCO Bhola, Anmyon, Dalma, Karachi, and MAARCO sites, respectively, because these time periods have the most measurements available at their corresponding sites. AERONET Level 2.0 data are used for most sites except Bhola and Anmyon where only Level 1.5 data are available and utilized here. The monthly mean AERONET 0.44–0.87 μm Ångström exponent ranges from 0.597 to 1.842 over these six sites. Hence, these sites represent a wide range of atmospheric conditions, ranging from coarse mode dominated to fine mode dominated cases [15].

2.3. Data Extraction Procedure

MxD04 [4] AOD and the AOD retrieval from the new algorithm are evaluated against AERONET measurements. The spatio-temporal matching approach by Ichoku et al. [16] is applied to collocate AERONET AOD measurements and MODIS retrievals (MxD04 and/or new algorithm) for comparison. AERONET measurements within ±30 min of the MODIS overpass time are averaged and compared against MODIS retrievals averaged within a 50-km diameter circular region centered over the AERONET sites. MxD04 products include DT, DB, and DT/DB merged retrievals, corresponding to 0.55 μm AOD from Dark Target algorithms (DT-land and DT-ocean algorithms), Deep Blue algorithm, and combination of the two retrievals, respectively [4]. DT retrievals cover both vegetated land and ocean while DB is limited to land. As the DT-land algorithm is not designed to retrieve AOD over bright desert, there are significant retrieval gaps over land. The DB algorithm was originally designed to complement the gaps, based on the principal that desert pixels are relatively darker in deep-blue bands [5]. However, now DB has been extended to vegetated land surface as well. Thus, for this study, we use the DB/DT merged product (Quality flag = 1, 2, 3 over ocean, and Quality flag = 3 over land) from MxD04. Specifically, DB/DT merged product over ocean is indeed DT.

3. Retrieval Algorithm

3.1. Retrieval Principal and Sensitivity Analysis

Like the existing DT algorithm, our new algorithm is based on a lookup table (LUT) approach, meaning that it is attempting to match simulated TOA reflectance to the observed TOA reflectance. The best-match solution represents the AOD and other properties of the aerosol. Over the open ocean, the atmospheric signal tends to dominate that of clear water. Near coastlines, however, sediments and turbid waters can dominate the signal, making aerosol retrieval impossible. However, in longer wavelengths, such as 2.1 μm (at which the imaginary part of refractive index for liquid water is several of magnitude larger than that in the visible), the penetration depth is so small, that contribution from the water should be nearly negligible.
To test this assumption, we compare the sensor sensitivity at 2.1 μm at which water-leaving radiance is negligible regardless of water turbidity with that at shorter wavelength which has large water-leaving radiance over turbid coastal water (taking 0.65 μm as an example). Figure 2 shows 2.1 μm TOA reflectance has better sensitivity to aerosol than 0.65 μm in both fine and coarse mode aerosol situations when 0.65 μm surface reflectance is large (turbid coastal water). When AOD is small (less than 0.15 at 0.55 μm, Figure 2b), the TOA reflectance (Figure 2c,d) is nearly the surface reflectance. Surface reflectance of turbid coastal water is up to 0.4 at 0.65 μm (Figure 2c) while it is very small (less than 0.0035) at 2.1 μm (Figure 2d). Figure 2e–h present a simulation of TOA reflectance through UNL-VRTM model [17] with various surface reflectances in fine aerosol (average TOA reflectance of the four fine aerosol modes defined in the DT-ocean LUT) and coarse aerosol (average TOA reflectance of the coarse aerosol modes defined in the DT-ocean LUT) situations. In a fine aerosol situation, 0.65 μm TOA reflectance is an increasing function of AOD when surface reflectance is small (0.03, clear water) (Figure 2e). However, it first decreases and then increases slowly when surface reflectance is large (0.3 or 0.4, turbid coastal water) (Figure 2e). The gradient of 2.1 μm MODIS digital signal (defined as change of MODIS digital count, d n * * , with respect to 0.55 μm AOD, or ( d n * * ) 2.1 / ( A O D ) 0.55 ) is ~25 regardless of AOD values (Figure 2j) while ( d n * * ) 0.65 / ( A O D ) 0.55 is smaller than 25 when AOD is less than 0.5 (1.2) and surface reflectance is 0.3 (0.4). This means that it is better to use TOA reflectance at 2.1 μm rather than 0.65 μm to retrieve fine AOD in turbid coastal water situation. Additionally, while the sensitivity of digital signal at 2.1 μm to the AOD change is a factor of 3–8 smaller than the counterpart at 0.65 μm (e.g., contrasts of blue lines between Figure 2i,j) in the clear water situation, 2.1 μm still has significant sensitivity to the change of fine-mode AOD, and its detection limit for fine-mode AOD is ~0.04 (i.e., the inverse of ( d n * * ) 2.1 / ( A O D ) 0.55 ). The reason that 2.1 μm still has reasonably good sensitivity for fine-mode aerosols is because the ocean surface is nearly black at 2.1 μm, albeit fine aerosol extinction decrease significantly.
In a coarse aerosol situation, ( d n * * ) 0.65 / ( A O D ) 0.55 is similar to that in a fine aerosol situation and is a little smaller than ( d n * * ) 2.1 / ( A O D ) 0.55 when 0.65 μm surface reflectance is small (Figure 2k,l). However, ( d n * * ) 0.65 / ( A O D ) 0.55 is a factor of 6–10 smaller than ( d n * * ) 2.1 / ( A O D ) 0.55 when water is turbid (e.g., contrast of red/green curves between Figure 2k,l). Overall, over turbid water, it is better to use TOA reflectance at 2.1 μm rather than 0.65 μm to retrieve AOD, regardless of whether fine or coarse mode dominates. In other words, 2.1 μm should be used to retrieve coarse-mode AOD regardless of water turbidity.

3.2. Algorithm Implementation and Steps

The retrieval algorithm is designed to supplement the existing DT algorithm. Hence, the lookup table (LUT) from the MODIS DT-ocean AOD retrieval algorithm is also used here [4,9]. That LUT is created by using Ahmad and Fraser’s [18] radiative transfer code, assuming aerosol optical properties from four “fine” (effective radius < 1 μm) and five “coarse” (effective radius > 1 μm) lognormal mode. The fine modes include two water-insoluble and two soluble modes, whereas the coarse modes are separated into three soluble sea-salt like, and two insoluble dust-like modes (http://darktarget.gsfc.nasa.gov/algorithm/ocean/aerosol-models). This LUT is defined at various 0.55 μm AOD values in the range of 0 to 3 and at different Sun-Earth-satellite geometries [4,9]. Furthermore, as water surface reflectance depends on surface wind speed, the LUT is constructed at four wind speeds, specifically 2 m s−1, 6 m s−1, 10 m s−1, and 14 m s−1.
Like the DT algorithm, combination of a fine mode and a coarse mode is required for AOD retrieval. In the retrieval procedure, MxD02 TOA reflectance at 2.1 μm is used to fit the corresponding LUT value ρ 2.1 LUT ( τ 0.55 t o t ) which is a weighted sum of pure fine mode LUT value ρ 2.1 f ( τ 0.55 t o t ) and pure coarse mode LUT ρ 2.1 c ( τ 0.55 t o t ) at a given 0.55 μm AOD value τ 0.55 t o t .
ρ 2.1 LUT ( τ 0.55 t o t ) = η ρ 2.1 f ( τ 0.55 t o t ) + ( 1 η ) ρ 2.1 c ( τ 0.55 t o t )
The value of reflectance weighting parameter η equals the fraction of total AOD at 0.55 μm contributed by the fine mode [9]. In addition, the LUT is calculated under the assumption of “gas free” (gas absorption is not included), thus correction [4] of MxD02 TOA reflectance is required to match the LUT before retrieving AOD.
Figure 3 shows the steps of our algorithm to retrieve AOD over coastal water, and hereafter, we call this algorithm the Coastal Water (CW) algorithm. Each step of CW is described below with a note that these steps essentially follow DT-Ocean algorithm except CW uses the 2.1 μm to retrieve AOD for those cloud-free turbid coastal water scenes. Hence, for each valid cloud-free turbid coastal water scene (i.e., a box of 20 × 20 pixels at 500 m resolution or at 10 km resolution at nadir, as available in the standard MODIS aerosol product), the following steps are implemented.
  • Collect and organize 20 × 20 pixels at 500 m resolution, remove pixels that are defined by land/sea mask as “land”, designated by ice/snow mask to be “ice”, designated by the cloud mask to be “cloud”, or removed by other tests.
  • Discard the brightest 25% and darkest 25% pixels defined with 0.86 μm reflectance.
  • Conduct gas (H2O, CO2, and O3) correction [4] for the remaining pixels.
  • Calculate the mean 2.1 μm TOA reflectance if there are still no less than 10 pixels. Otherwise retrieval is not conducted. Calculate the sun glint angle [9]. If the sun glint angle is less than 40°, the retrieval is not conducted.
  • Prescribe single scattering properties of the aerosol. By only using 2.1 μm reflectance to retrieve AOD, there is no sensitivity to aerosol optical properties. Figure 4 shows that AOD can differ up to 0.2 in 100 km from the coast, but FMF differs only by 0.08 in 100 km from the coast. Thus, we assume the single scattering properties (including FMF) and surface wind speed for a turbid coastal water pixel is the same as those used for the AOD retrieval by the standard MODIS algorithm over its closest open-ocean pixel (within 100 km radius). The assumption that aerosol type does not change over moderate spatial scale is reasonable and was used in the atmospheric correction of SeaWiFS imagery over turbid coastal waters [19].
  • Use the mean 2.1 μm TOA reflectance and lookup table determined by Equation (1) to retrieve AOD over the turbid coastal water where MxD04 product is unavailable. In application of Equation (1), all ancillary information (aerosol mode selection and FMF) is obtained from step 5.

4. Results

Figure 5 shows the example of CW applied to an Aqua-observed scene (10 December 2015) over the Bay of Bengal. According to the standard DT retrieval (MYD04), AOD reaches 1.6 in the center of the granule and decreases gradually southward. MYD04 does not provide retrievals over the coastal water at north part of the Bay of Bengal (Figure 5c). This non-retrieval region is in cloud-free conditions (Figure 5a), so AOD was not retrieved due to turbid water or underlying sea-floor. Figure 5b shows non-retrievals are in 36.4% of this cloud-free, coastal area (within 5 km of coastline or with water depth less than 50 m).
The CW algorithm enhances DT-ocean, by retrieving AOD over those turbid water conditions (Figure 5d). These “new” retrievals are consistent, in that there appears to be smooth transition in retrieved AOD (e.g., land → CW → ocean). In addition, the new CW AOD retrievals compare well with corresponding values observed by AERONET (overlaid as filled circles in Figure 5c,d).
To validate, we follow the standard protocol, in that AERONET measurements within ±30 min of the MODIS overpass time are averaged and compared against MODIS retrievals averaged within a 25-km radius circular region centered over the AERONET site. This means that there are three situations: MxD04 retrievals only, CW retrievals only, and cases where both are available within the 25-km of the AERONET site.
AOD (at 0.55 μm) retrievals from MxD04 product and CW algorithm are evaluated against AERONET data (as shown in Figure 6). We use a Venn diagram in Figure 6 to represent how MODIS AOD retrievals (MxD04 and/or CW) and AERONET observations are collocated. Figure 6a is a scatter plot of MxD04 retrievals versus AERONET measurements in the cases that have retrievals only from MxD04, and these cases can be divided into two categories: (a) all cloud-free pixels within 25-km radius proximity of AERONET site are retrieved through MODIS operational algorithm or (b) only part of cloud-free pixels are retrieved through MODIS DT algorithm, but the AOD of the rest of the cloud-free pixels cannot be retrieved with the CW algorithm because these pixels are over turbid water and no open-ocean pixel is close enough (i.e., within 100 km) to provide aerosol single scattering property for CW algorithm. Figure 6b is a scatter plot of CW retrievals versus AERONET measurements in the cases where all cloud-free pixels within 25-km radius proximity of AERONET site cannot be retrieved through MODIS DT algorithm but some or all of these cloud-free pixels can be retrieved through CW algorithm. Comparison of Figure 6a,b shows that normalized mean bias (NMB, N M B = i = 1 N ( ( τ i M O D I S τ i A E R O E N T ) / τ i A E R O E N T ) / N ) and root mean square error (RMSE) of CW algorithm (12.0% and 0.141) are smaller than that of MxD04 (15.9% and 0.213). The percentage of collocated pairs within expected error envelope (+(0.04 + 10%), −(0.02 + 10%), asymmetric) and correlation coefficient increase from 35.5% and 0.72 in MxD04 only situation (Figure 6a) to 39.3% and 0.94 in CW only situation (Figure 6b), respectively.
In addition to the situations where only MxD04 or only CW retrievals are available within the 25-km radius proximity of AERONET site, there is a third one where some of the cloud-free pixels are retrieved through the operational MODIS algorithm while some are retrieved through the CW algorithm. Figure 6c,d are similar to Figure 6a,b, respectively, but both include the AOD retrievals from their respective counterparts of this third situation. Overall, for all possible retrievals by each algorithm, CW retrievals are comparable to MxD04 in quality, although the MODIS operational algorithm has slightly more samples.
The basis of MODIS-AERONET collocation is that air masses are always in motion and the average of MODIS AOD retrievals in a certain area which encompass an AERONET site should be comparable to the temporal statistics of the AERONET measurements [16]. In the situation that both MxD04 and CW retrievals are available within the 25-km radius proximity of AERONET site, if only MxD04 or CW retrievals are collocated with AERONET observation, it will be biased comparison. Therefore, we further show the inter-comparison between either MxD04 or CW AOD with AERONET AOD at this situation (Figure 6f–h). The percentage of collocated pairs in the EE for (MxD04 + CW) combined vs. AERONET (53.1%, Figure 6h) is larger than that of MxD04-AERONET (49.2%, Figure 6f) and CW-AERONET (40.6%, Figure 6g). The RMSE of (MxD04 + CW)-AERONET (0.110) is smaller than that of MxD04-AERONET (0.114) and CW-AERONET (0.148).
Inter-comparison of MxD04 and CW merged AOD with AERONET AOD is shown in Figure 6e; such merged AOD include all data points in Figure 6a,b,h. The number of collocated pairs increases from 190 (62 in MxD04 only situation plus 128 in the situation that both MxD04 and CW retrievals are available) to 218 (28 in CW only situation are added). In addition, the inter-comparison statistics (Figure 6e) is still comparable to or better than those retrieved from MxD04 (Figure 6a,c) in quality. Overall, CW retrievals supplement MODIS DT, and improve the AOD retrievals both spatially and temporally without degrading (and sometimes increasing) the DT-Ocean AOD retrieval quality.
In addition to evaluation of MODIS AOD diagnostic error in Figure 6, prognostic error is presented in Figure 7. Two thirds of AERONET-MODIS collocations are within the expected error envelope (y = ±(0.05 + 0.15x)) (Figure 7c), thus 0.05 + 0.15 (MODIS AOD) can be considered as its prognostic error. Figure 7d–f show RMSE as a function of MODIS; RMSE goes up as MODIS AOD increases and combined retrievals (Figure 7f) are comparable to MxD04 (Figure 7d).

5. Discussion

The CW algorithm uses 2.1 μm to retrieve AOD over turbid coastal water as water-leaving radiance is negligible at that band. As only one band is used in the algorithm, aerosol single scattering properties need to be prescribed. We assume the single scattering properties for a turbid coastal water pixel are the same as those used for the AOD retrieval by the standard MODIS algorithm over its closest open-ocean pixel (within 100 km radius). When no pixel is within 100 km, we do not conduct retrieval. Further evaluation can be targeted at the assumptions regarding the use of aerosol single scattering properties from adjacent non-turbid water pixels within certain threshold distance.
The 6 AERONET sites used for this analysis were located in polluted regions. We expect to evaluate how the algorithm performs in an unpolluted zone if there is a new AERONET site or other field campaign that satisfies all the conditions: (1) it is close to coast line; (2) coastal water is turbid; and (3) the region is unpolluted.
The CW algorithm potentially can also be used for the other aerosol retrieval algorithm over the ocean such as the SeaWiFS Ocean Aerosol Retrieval (SOAR) algorithm developed by Sayer et al. [20]. In SOAR, TOA reflectance measured by the SeaWiFS bands centered near 0.51 μm, 0.67 μm, and 0.86 μm are used to retrieve AOD regardless of water turbidity. Hence, as pointed out by Sayer et al. [20], the performance of SOAR retrieval is expected to be poorer in turbidity cases and a double of TOA reflectance as a result of water turbidity causes a positive error in retrieved 0.55 μm AOD of 0.25 [20]. Indeed, in the evaluation of SOAR algorithm, AOD retrievals that are proximity to coast are not included [20].
As MODIS sensors aboard on Terra and Aqua have been providing data for more than one decade and will be decommissioned by the early 2020s, its aerosol record is expected to be continued by the Visible Infrared Imaging Radiometer Suite (VIIRS) on board Suomi-NPP (S-NPP) which was launched in late 2011. Thus, future investigation will also include the application of the approach in this study to VIIRS. VIIRS has 2.26 μm band at which water leaving radiance is also negligible like 2.1 μm band on MODIS and thus the method can be likewise applied.

6. Conclusions

The MODIS Dark Target (DT) algorithm has been applied to retrieve AOD over land and ocean since early 2000. Although there have been significant improvements to the algorithm in the past decade, it has not yet been able to retrieve AOD over turbid coastal water due to its high water-leaving radiance from the ocean bottom or water color. We designed the Coastal Water AOD retrieval algorithm (CW algorithm) for these regions to supplement the current DT algorithm. The AOD retrieval algorithm for turbid coastal water takes advantage of the fact that water-leaving radiance is negligible at 2.1 μm, hence this band is only used to retrieve while other auxiliary information such as aerosol single scattering properties are obtained from DT retrievals over the nearby non-turbid water surfaces. In other words, the aerosol modes and reflectance weighting parameter of the pixel to be retrieved is substituted by the closest counterpart from the MxD04 ocean AOD product.
The CW algorithm not only fills the gaps of MxD04 AOD retrievals over turbid coastal water but also improves their comparison with AERONET measurements. CW AOD retrievals are validated against measurements of six AERONET sites that are located at coastal regions. The new algorithm yields ~18% more of MODIS-AERONET collocated pairs, and CW AOD retrievals are comparable or better than MxD04 in quality; the RMSE of MxD04 and CW product is 0.154 and 0.147, respectively. In the situation that both AOD retrievals from MxD04 and CW are available within a 50 km diameter circular regions centered at AERONET sites, the merged AOD retrievals show better agreement with AERONET AOD either of them alone, in terms of the percentage of collocated pairs in the EE (49.2%, 40.6%, and 53.1% for MxD04, CW, and merged product, respectively) and RMSE (0.114, 0.148, 0.110 for MxD04, CW, and merged product, respectively). In addition, 0.05 + 0.15(MODIS AOD) can be considered as prognostic error of the merged product.

Acknowledgments

Funding for this study was provided by the NASA Earth Science Division as part of NASA's Terra/Aqua and VIIRS science team program (grant numbers: NNX16AT63G and NNX17AC94A) and NASA’s GEO-CAPE and KORUS-AQ program (grant number: NNX16AT82G), and was also supported by Office of Naval Research (ONR's) Multidisciplinary University Research Initiatives (MURI) Program under the award N00014-16-1-2040.

Author Contributions

Yi Wang, Jun Wang and Robert C. Levy conceived and designed the experiments; Yi Wang and Xiaoguang Xu performed the experiments and analyzed the data; Yi Wang and Jun Wang wrote the paper; Jeffrey S. Reid provided guidance to the project, reviewed and edited the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Haywood, J.; Boucher, O. Estimates of the direct and indirect radiative forcing due to tropospheric aerosols: A review. Rev. Geophys. 2000, 38, 513–543. [Google Scholar] [CrossRef]
  2. Rosenfeld, D.; Lohmann, U.; Raga, G.B.; Dowd, C.D.; Kulmala, M.; Fuzzi, S.; Reissell, A.; Andreae, M.O. Flood or drought: How do aerosols affect precipitation? Science 2008, 321, 1309. [Google Scholar] [CrossRef] [PubMed]
  3. Lim, S.S.; Vos, T.; Flaxman, A.D.; Danaei, G.; Shibuya, K.; Adair-Rohani, H.; AlMazroa, M.A.; Amann, M.; Anderson, H.R.; Andrews, K.G.; et al. A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990–2010: A systematic analysis for the global burden of disease study 2010. Lancet 2012, 380, 2224–2260. [Google Scholar] [CrossRef]
  4. Levy, R.C.; Mattoo, S.; Munchak, L.A.; Remer, L.A.; Sayer, A.M.; Patadia, F.; Hsu, N.C. The collection 6 modis aerosol products over land and ocean. Atmos. Meas. Tech. 2013, 6, 2989–3034. [Google Scholar] [CrossRef]
  5. Hsu, N.C.; Jeong, M.J.; Bettenhausen, C.; Sayer, A.M.; Hansell, R.; Seftor, C.S.; Huang, J.; Tsay, S.C. Enhanced deep blue aerosol retrieval algorithm: The second generation. J. Geophys. Res. 2013, 118, 9296–9315. [Google Scholar] [CrossRef]
  6. Sayer, A.M.; Hsu, N.C.; Bettenhausen, C.; Jeong, M.J. Validation and uncertainty estimates for modis collection 6 “deep blue” aerosol data. J. Geophys. Res. 2013, 118, 7864–7872. [Google Scholar] [CrossRef]
  7. IPCC. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change. In Climate Change 2013: The Physical Science Basis; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013; p. 1535. [Google Scholar]
  8. Li, R.-R.; Kaufman, Y.J.; Gao, B.-C.; Davis, C.O. Remote sensing of suspended sediments and shallow coastal waters. IEEE Trans. Geosci. Remote Sens. 2003, 41, 559–566. [Google Scholar]
  9. Remer, L.A.; Kaufman, Y.J.; Tanré, D.; Mattoo, S.; Chu, D.A.; Martins, J.V.; Li, R.R.; Ichoku, C.; Levy, R.C.; Kleidman, R.G.; et al. The modis aerosol algorithm, products, and validation. J. Atmos. Sci. 2005, 62, 947–973. [Google Scholar] [CrossRef]
  10. Anderson, J.C.; Wang, J.; Zeng, J.; Leptoukh, G.; Petrenko, M.; Ichoku, C.; Hu, C. Long-term statistical assessment of aqua-modis aerosol optical depth over coastal regions: Bias characteristics and uncertainty sources. Tellus B 2013, 65. [Google Scholar] [CrossRef]
  11. Wolfe, R.E.; Nishihama, M.; Fleig, A.J.; Kuyper, J.A.; Roy, D.P.; Storey, J.C.; Patt, F.S. Achieving sub-pixel geolocation accuracy in support of modis land science. Remote Sens. Environ. 2002, 83, 31–49. [Google Scholar] [CrossRef]
  12. King, M.D.; Menzel, W.P.; Kaufman, Y.J.; Tanre, D.; Bo-Cai, G.; Platnick, S.; Ackerman, S.A.; Remer, L.A.; Pincus, R.; Hubanks, P.A. Cloud and aerosol properties, precipitable water, and profiles of temperature and water vapor from modis. IEEE Trans. Geosci. Remote Sens. 2003, 41, 442–458. [Google Scholar] [CrossRef]
  13. Smirnov, A.; Holben, B.N.; Eck, T.F.; Dubovik, O.; Slutsker, I. Cloud-screening and quality control algorithms for the aeronet database. Remote Sens. Environ. 2000, 73, 337–349. [Google Scholar] [CrossRef]
  14. Holben, B.N.; Eck, T.F.; Slutsker, I.; Tanré, D.; Buis, J.P.; Setzer, A.; Vermote, E.; Reagan, J.A.; Kaufman, Y.J.; Nakajima, T.; et al. Aeronet—A federated instrument network and data archive for aerosol characterization. Remote Sens. Environ. 1998, 66, 1–16. [Google Scholar] [CrossRef]
  15. O’Neill, N.T.; Eck, T.F.; Smirnov, A.; Holben, B.N.; Thulasiraman, S. Spectral discrimination of coarse and fine mode optical depth. J. Geophys. Res. 2003, 108, 4559. [Google Scholar] [CrossRef]
  16. Ichoku, C.; Chu, D.A.; Mattoo, S.; Kaufman, Y.J.; Remer, L.A.; Tanré, D.; Slutsker, I.; Holben, B.N. A spatio-temporal approach for global validation and analysis of modis aerosol products. Geophys. Res. Lett. 2002, 29. [Google Scholar] [CrossRef]
  17. Wang, J.; Xu, X.; Ding, S.; Zeng, J.; Spurr, R.; Liu, X.; Chance, K.; Mishchenko, M. A numerical testbed for remote sensing of aerosols, and its demonstration for evaluating retrieval synergy from a geostationary satellite constellation of geo-cape and goes-r. J. Quant. Spectrosc. Radiat. Transf. 2014, 146, 510–528. [Google Scholar] [CrossRef]
  18. Ahmad, Z.; Fraser, R.S. An iterative radiative transfer code for ocean-atmosphere systems. J. Atmos. Sci. 1982, 39, 656–665. [Google Scholar] [CrossRef]
  19. Hu, C.; Carder, K.L.; Muller-Karger, F.E. Atmospheric correction of seawifs imagery over turbid coastal waters: A practical method. Remote Sens. Environ. 2000, 74, 195–206. [Google Scholar] [CrossRef]
  20. Sayer, A.M.; Hsu, N.C.; Bettenhausen, C.; Ahmad, Z.; Holben, B.N.; Smirnov, A.; Thomas, G.E.; Zhang, J. Seawifs ocean aerosol retrieval (soar): Algorithm, validation, and comparison with other data sets. J. Geophys. Res. Atmos. 2012, 117. [Google Scholar] [CrossRef]
Figure 1. Regional (ae and gi) and global (f) distribution of AOD unavailable ratio which is defined as the ratio between number of non-retrieval pixels and all pixels in cloud-free conditions. MYD04 product in 2016 are used to conduct the statistics. Zoom in on the red-box regions of (f) are shown in (ae and gi). Lime points in (f) are AERONET sites used for validation and more information is in Table 1. Numbers in the pink boxes are AOD unavailable ratio over coastal water in clear sky condition.
Figure 1. Regional (ae and gi) and global (f) distribution of AOD unavailable ratio which is defined as the ratio between number of non-retrieval pixels and all pixels in cloud-free conditions. MYD04 product in 2016 are used to conduct the statistics. Zoom in on the red-box regions of (f) are shown in (ae and gi). Lime points in (f) are AERONET sites used for validation and more information is in Table 1. Numbers in the pink boxes are AOD unavailable ratio over coastal water in clear sky condition.
Remotesensing 09 00595 g001
Figure 2. (ad) are MODIS Aqua true color image, 0.55 μm AOD, 0.65 μm TOA reflectance, and 2.1 μm TOA reflectance at Bohai Sea, China on 23 February 2017, respectively. Red box is the region where water is turbid and atmosphere is clean. (ef) are simulation (solar zenith angle is 24°, view zenith angle 54°, and relative azimuth angle 60°) of 0.65 μm and 2.1 μm TOA reflectance with fine aerosol model, respectively. (gh) are similar to (ef) but for coarse aerosol model. (ij) are gradient of 0.65 μm and 2.1 μm MODIS digital signal (dn**) with respect to 0.55 μm AOD with fine aerosol model. (kl) are similar to (ij) but for coarse aerosol model. Values in the legends are surface reflectance used to simulate TOA reflectance at corresponding wavelength.
Figure 2. (ad) are MODIS Aqua true color image, 0.55 μm AOD, 0.65 μm TOA reflectance, and 2.1 μm TOA reflectance at Bohai Sea, China on 23 February 2017, respectively. Red box is the region where water is turbid and atmosphere is clean. (ef) are simulation (solar zenith angle is 24°, view zenith angle 54°, and relative azimuth angle 60°) of 0.65 μm and 2.1 μm TOA reflectance with fine aerosol model, respectively. (gh) are similar to (ef) but for coarse aerosol model. (ij) are gradient of 0.65 μm and 2.1 μm MODIS digital signal (dn**) with respect to 0.55 μm AOD with fine aerosol model. (kl) are similar to (ij) but for coarse aerosol model. Values in the legends are surface reflectance used to simulate TOA reflectance at corresponding wavelength.
Remotesensing 09 00595 g002
Figure 3. Flowchart of retrieving AOD over coastal water.
Figure 3. Flowchart of retrieving AOD over coastal water.
Remotesensing 09 00595 g003
Figure 4. (a) are monthly mean 0.55 μm AOD and fine mode fraction (FMF) from MOD04 (Terra) in December, 2015. (b) is average absolute AOD difference as a function a distance with respect to reference points (black solid circles in (a)). (c) are monthly mean 0.55 μm AOD and fine mode fraction (FMF) from MOD04 (Terra) in December, 2015. (d) is similar to (b), but for FMF.
Figure 4. (a) are monthly mean 0.55 μm AOD and fine mode fraction (FMF) from MOD04 (Terra) in December, 2015. (b) is average absolute AOD difference as a function a distance with respect to reference points (black solid circles in (a)). (c) are monthly mean 0.55 μm AOD and fine mode fraction (FMF) from MOD04 (Terra) in December, 2015. (d) is similar to (b), but for FMF.
Remotesensing 09 00595 g004
Figure 5. (a) Aqua MODIS true color image on 10 December 2015. (b) Blue and red represents coastal water where MYD AOD retrievals are available (blue) and unavailable (red), respectively in clear sky condition. (c) MYD04 0.55 μm DT/DB merged AOD retrievals. (d) is similar to (c), but the coastal water gaps where MYD04 AOD retrievals are not available in clear sky condition (red box) are filled by AOD retrievals from CW algorithm. AERONET AOD measurement overlaps on (c) and (d).
Figure 5. (a) Aqua MODIS true color image on 10 December 2015. (b) Blue and red represents coastal water where MYD AOD retrievals are available (blue) and unavailable (red), respectively in clear sky condition. (c) MYD04 0.55 μm DT/DB merged AOD retrievals. (d) is similar to (c), but the coastal water gaps where MYD04 AOD retrievals are not available in clear sky condition (red box) are filled by AOD retrievals from CW algorithm. AERONET AOD measurement overlaps on (c) and (d).
Remotesensing 09 00595 g005
Figure 6. Scatter plots of 0.55 μm AOD retrievals from MODIS (MxD04 and/or CW) versus AERONET observations of 6 sites in one month. Venn diagram on the upper left represents number of collocated data sets for MxD04-AERONET (blue), CW-AERONET (red), (MxD04 + CW)-AERONET (purple). Color part of Venn diagram on the top of each scatter plot (ah) represents which collocated data sets are plotted; see details in the text for further description of the Venn diagram. (a) is MxD04-AERONET in the situation that only MxD04 retrievals are available within the 25-km of AERONET site. (b) is CW-AERONET in the situation that only CW retrievals are available. (cd) are similar to (ab), respectively, but both additionally include AOD retrievals from their respective counterparts of the situation that both MxD04 and CW retrievals are available. (e) is (MxD04 + CW)-AERONET in all the three situations. (fh) are MxD04-AERONET, CW-AERONET, (MxD04 + CW)-AERONET, respectively, in the situation that both MxD04 and CW retrievals are available. 1:1 lines and expected error (EE) envelopes (+(0.04 + 10%), −(0.02 + 10%), asymmetric) are plotted as dot and dashed lines. The number of collocated pairs (N), root mean square error (RMSE), normalized mean bias (NMB), and linear correlation coefficients are also shown.
Figure 6. Scatter plots of 0.55 μm AOD retrievals from MODIS (MxD04 and/or CW) versus AERONET observations of 6 sites in one month. Venn diagram on the upper left represents number of collocated data sets for MxD04-AERONET (blue), CW-AERONET (red), (MxD04 + CW)-AERONET (purple). Color part of Venn diagram on the top of each scatter plot (ah) represents which collocated data sets are plotted; see details in the text for further description of the Venn diagram. (a) is MxD04-AERONET in the situation that only MxD04 retrievals are available within the 25-km of AERONET site. (b) is CW-AERONET in the situation that only CW retrievals are available. (cd) are similar to (ab), respectively, but both additionally include AOD retrievals from their respective counterparts of the situation that both MxD04 and CW retrievals are available. (e) is (MxD04 + CW)-AERONET in all the three situations. (fh) are MxD04-AERONET, CW-AERONET, (MxD04 + CW)-AERONET, respectively, in the situation that both MxD04 and CW retrievals are available. 1:1 lines and expected error (EE) envelopes (+(0.04 + 10%), −(0.02 + 10%), asymmetric) are plotted as dot and dashed lines. The number of collocated pairs (N), root mean square error (RMSE), normalized mean bias (NMB), and linear correlation coefficients are also shown.
Remotesensing 09 00595 g006
Figure 7. (ac) are scatter plots of AERONET 0.55 μm AOD observations versus retrievals from MODIS (MxD04 and/or CW) for the 6 monthly cases. Sampling method is same as in Figure 6. 1:1 lines and expected error (EE) envelopes (+(0.05 + 15%), −(0.05 + 15%), symmetric) are plotted as dotted and dashed lines. (df) are root mean square error (RMSE) as a function of MODIS AOD in the 0.2 interval. Numbers above the averaged dot show how many AERONET-MODIS collocations are available to calculate RMSE.
Figure 7. (ac) are scatter plots of AERONET 0.55 μm AOD observations versus retrievals from MODIS (MxD04 and/or CW) for the 6 monthly cases. Sampling method is same as in Figure 6. 1:1 lines and expected error (EE) envelopes (+(0.05 + 15%), −(0.05 + 15%), symmetric) are plotted as dotted and dashed lines. (df) are root mean square error (RMSE) as a function of MODIS AOD in the 0.2 interval. Numbers above the averaged dot show how many AERONET-MODIS collocations are available to calculate RMSE.
Remotesensing 09 00595 g007
Table 1. Information of AERONET sites used for validation.
Table 1. Information of AERONET sites used for validation.
SiteLocation *PeriodData LevelMonthly Mean 0.44–0.87 μm Ångström Exponent
MVCO, New England41.3°N 70.6°WAugust 20152.01.842
Bhola, Bangladesh22.2°N 90.8°EDecember 20151.51.206
Anmyon, S. Korea36.5°N 126.3°EMay 20161.51.076
Dalma, UAE24.5°N 52.3°EAugust 20042.00.711
Karachi, Pakistan24.9°N 67.0°EMarch 20142.00.701
MAARCO, UAE24.7°N 54.7°ESeptember 20042.00.597
* All the locations are shown as lime points in Figure 1f.

Share and Cite

MDPI and ACS Style

Wang, Y.; Wang, J.; Levy, R.C.; Xu, X.; Reid, J.S. MODIS Retrieval of Aerosol Optical Depth over Turbid Coastal Water. Remote Sens. 2017, 9, 595. https://doi.org/10.3390/rs9060595

AMA Style

Wang Y, Wang J, Levy RC, Xu X, Reid JS. MODIS Retrieval of Aerosol Optical Depth over Turbid Coastal Water. Remote Sensing. 2017; 9(6):595. https://doi.org/10.3390/rs9060595

Chicago/Turabian Style

Wang, Yi, Jun Wang, Robert C. Levy, Xiaoguang Xu, and Jeffrey S. Reid. 2017. "MODIS Retrieval of Aerosol Optical Depth over Turbid Coastal Water" Remote Sensing 9, no. 6: 595. https://doi.org/10.3390/rs9060595

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop