Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
A Microscale Linear Phased-Array Ultrasonic Transducer Based on PZT Ceramics
Next Article in Special Issue
Improved Bathymetric Mapping of Coastal and Lake Environments Using Sentinel-2 and Landsat-8 Images
Previous Article in Journal
Dissecting and Interpreting a Three-Dimensional Ground-Penetrating Radar Dataset: An Example from Northern Australia
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Remote Sensing of Daytime Water Leaving Reflectances of Oceans and Large Inland Lakes from EPIC onboard the DSCOVR Spacecraft at Lagrange-1 Point

1
Remote Sensing Division, Naval Research Laboratory, Washington, DC 20375 USA
2
NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(5), 1243; https://doi.org/10.3390/s19051243
Submission received: 13 February 2019 / Revised: 5 March 2019 / Accepted: 8 March 2019 / Published: 12 March 2019
(This article belongs to the Special Issue Advances in Quantitative Remote Sensing: Past, Present and Future)
Figure 1
<p>(<b>a</b>) The sunlit full disk color image processed from three EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 13 March 2013 at UTC 1442; (<b>b</b>) The full disk 764-nm channel image; (<b>c</b>) The portions of color image over southern part of S. America; (<b>d</b>) The portions of the 764-nm image over southern part of S. America; (<b>e</b>) Apparent reflectances as a function of wavelength for a turbid river water pixel (blue line marked with ‘+’) and a land pixel (red line marked with ‘*’). Pixels are selected in the red box marked in Panel (C) and (D).</p> ">
Figure 2
<p>(<b>a</b>) The sunlit full disk color image processed from three atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 13 March 2017 at UTC 1442; (<b>b</b>) The portions of color image over S. America; (<b>c</b>) The portions of color image of water surfaces over S. America after masking out land and cloud pixels; (<b>d</b>) Spectral reflectances as a function of wavelength for one turbid river water pixel located at the center of the small rectangle outlined in (<b>c</b>). In (<b>d</b>), the line marked with ‘+’ is the top of atmosphere apparent reflectance curve. The lines marked with ‘*’, diamond, and triangle are the retrieved water leaving reflectances for the Rayleigh only, Rayleigh plus rural aerosol model, and Rayleigh plus maritime model corrections, respectively. The inlet in (<b>c</b>) shows an enlarged portion of La Plata River.</p> ">
Figure 3
<p>(<b>a</b>) The sunlit full disk color image processed from three atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 24 September 2017 at UTC 1649; (<b>b</b>) The portion of color image over Great Lake region; (<b>c</b>) The portion of color image of water surfaces over Great Lakes after masking out land and cloud pixels; (<b>d</b>) Spectral reflectances as a function of wavelength for one chlorophyll-rich water pixel located in the left part of Lake Erie. (<b>e</b>) The 1-km resolution true color Terra MODIS image over the Great Lake area. In (<b>d</b>), the line marked with ‘+’ is the top of atmosphere apparent reflectance curve. The lines marked with ‘*’, diamond, and triangle are the retrieved water leaving reflectances for the Rayleigh only, Rayleigh plus rural aerosol model, and Rayleigh plus maritime model corrections, respectively.</p> ">
Figure 4
<p>(<b>a</b>) The sunlit full disk color image processed from three atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 26 June 2017 at UTC 1112; (<b>b</b>) The portion of image over Caspian Sea, Black Sea, Azov Sea, Persian Gulf, Red Sea, and Mediterranean Sea; (<b>c</b>) The portion of image of water surfaces over the same area as that of (<b>b</b>), but after masking out the land and cloud pixels; (<b>d</b>) Spectral reflectances as a function of wavelength for two chlorophyll-rich water pixels (Green lines) over Caspian Sea and Azov Sea and two blue water pixels (Blue lines) over eastern portions of Black Sea and Mediterranean Sea.</p> ">
Figure 5
<p>(<b>a</b>–<b>h</b>) Eight sunlit full disk color images processed from atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 26 June 2017 at UTC 0545, 0650, 0756, 0901, 1007, 1112, 1217, and 1323, respectively. The spatial locations of the Caspian Sea in these images are pointed to with red arrows.</p> ">
Figure 6
<p>(<b>a</b>–<b>h</b>) Eight color images processed from atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 26 June 2017 at UTC 0545, 0650, 0756, 0901, 1007, 1112, 1217, and 1323, respectively, and covering areas consisting of Caspian Sea, Black Sea, and Azov Sea with land and cloud features masked out.</p> ">
Review Reports Versions Notes

Abstract

:
The NASA’s Earth Polychromatic Imaging Camera (EPIC) on board the Deep Space Climate Observatory (DSCOVR) satellite has been making multiple observations of the entire sunlit Earth in a given day from the Sun-Earth Largangian L1 point since the summer of 2015. EPIC contains 10 narrow channels in the 317–780 nm solar spectral range. The data acquired with EPIC have already been used in a variety of scientific investigations, including the study of the global ozone levels, aerosol index and aerosol optical depth, UV reflectivity of clouds over land and ocean, cloud height over land and ocean, and vegetation indices. In this article, we report that EPIC data, particularly for the data measured with narrow channels centered near 443, 551, and 680 nm, can also have important applications in remote sensing of ocean color in different geographical regions. We have modified a version of a multi-channel atmospheric correction algorithm for Moderate Resolution Imaging SpectroRadiometer (MODIS) ocean color applications and adapted the algorithm for processing EPIC data. We present three case studies on water leaving reflectance retrievals from EPIC data acquired over a large turbid river, inland lakes, and oceans. We conclude that a future ocean color instrument on board a satellite at the L1 point, which provides continuous view of the full sunlit disk of the Earth, will complement and extend ocean color observations with the low Earth observing polar orbital and geostationary satellite instruments in both the spatial and time domains.

1. Introduction

The NASA’s Earth Polychromatic Imaging Camera (EPIC) [1] on board the Deep Space Climate Observatory (DSCOVR) satellite has been making multiple observations of the entire sunlit Earth in a given day since the summer of 2015. Historically, the DSCOVR mission evolved from the late 1990s NASA intended Triana Mission, which was strongly endorsed by former vice president Albert Gore [2], for climate related research. Triana was a mission designed for deployment into a stable orbit, at approximately 1.5 million kilometers from the Earth in the direction of the Sun. The orbit is known as the Lissajous orbit about the Lagrangian point 1 (L1), which is stable in the sense that the satellite remains near the Sun-Earth line and views the full sunlit disk of the Earth continuously. Triana was proposed to be an exploratory mission to investigate the scientific and technical advantages of L1 for earth observations. The continuous view of the full sunlit disk of Earth was for complementing and extending observations from low Earth orbit (LEO) or geostationary Earth orbit (GEO) satellite instruments. Later on, Triana was renamed to DSCOVR and DSCOVR was launched into space in February of 2015.
The EPIC instrument on board the DSCOVR satellite contains 10 narrow channels in the 317–780 nm solar spectral range. Descriptions of the EPIC instruments have already been given, for examples, by Herman et al. [3], Marshak et al. [1] and Yang et al. [4]. The 10 spectral channels, their full width at half maximum (FWHM), and originally anticipated primary applications are listed in Table 1. The focal plane of the EPIC system is a 2048 × 2048 pixel charge-coupled device (CCD) array. Data for the 10 channels are obtained using a 10-color filter wheel assembly and shutter [4]. The filter transmission functions are given in Herman et al. [3]. The signal to noise ratio (SNR) for the 10 channels is ~200:1 [5]. Over bright surfaces, e.g., ice and snow, or thick clouds, the maximum SNR can achieve 290:1 [3].
From Table 1, it is seen that the observation of ocean color is not listed as a primary objective. Through visual examination of EPIC images posted on public websites and detailed scientific analysis of limited number of EPIC data sets, we have found that the EPIC data can have an additional application, i.e., remote sensing of ocean color.
In this article, we report that EPIC data, particularly for the data measured with narrow channels centered near 443, 551, and 680 nm, can have important applications in remote sensing of ocean color. We have modified a version of a multi-channel atmospheric correction algorithm for Moderate Resolution Imaging SpectroRadiometer (MODIS) [6,7] ocean color applications [8] and adapted the algorithm for processing EPIC data. The data and methods are described in Section 2. Three case studies on water leaving reflectance retrievals over a large turbid river, inland lakes, and oceans are shown in Section 3. Discussions are present in Section 4. Finally, a summary is given in Section 5.

2. Data and Methods

As listed in Table 1, EPIC has 10 channels in the solar spectral range between 317 and 780 nm. The ocean color information is mostly contained in images of the blue (443 nm), green (551 nm), and red (680 nm) channels. The other channels, in general, do not contain much color information about water surfaces. We use an EPIC data set acquired on 13 March 2017 at UTC 1442 to illustrate the properties of EPIC data. Figure 1a shows the sunlit portions of the full disk color image processed from the satellite Level 1B “apparent reflectance images” of three EPIC channels (Red: 680 nm; Green: 551 nm; Blue: 443 nm). At the exact nadir looking direction, the L1B data’s spatial resolution (pixel size) is approximately 8 km. For the off nadir looking pixels, such as those near the edge of the sunlit disk, the pixel sizes on the ground are coarser than 8 km.
The apparent reflectance ρ* at the satellite for a given channel is traditionally defined [9,10] as:
ρ * = π L μ 0 E 0
where L is the radiance measured by the satellite, μ0 the cosine of solar zenith angle, and E0 the extra-terrestrial solar flux. The clear water surfaces appear blue in Figure 1a. Clouds are white. However, most land surfaces appear in the distorted brownish color because the widths of the three channels are too narrow and these channels do not contain the full spectral information for the proper generation of natural-looking color images unless going through a special data processing procedure [3].
Figure 1b shows the black and white image processed from the apparent reflectances of the 764-nm channel, which is located within the strong atmospheric oxygen-A band absorption region. When generating this image, we set the pixels with ρ* (764 nm) values smaller than 0.05 to be zero. We also set the pixels with ρ* (764 nm) values greater than 1 to be 1. Clouds and land surfaces are bright in the Figure 1b image, while the blue colored water surfaces in Figure 1a are completely dark.
The shapes of the blue features in Figure 1a are approximately the same as those of black features in Figure 1b. This indicates that the lower end ρ* (764 nm) threshold value of 0.05 is appropriate for the separation of water surfaces from land surfaces and clouds. In order to further illustrate this point, we show in Figure 1c an enlarged portion of Figure 1a color image, but just covers the southern portion of S. America. The smaller area within the red-colored lines in Figure 1c contains the very turbid La Plata River, which is difficult to see. Figure 1d is the 764-nm O2-A channel image covering the same area as that of Figure 1c. The La Plata River in Figure 1d within the red-colored lines is picked up quite well. Figure 1e shows apparent reflectances as a function of wavelength for a pixel over the turbid river water (the blue line marked with ‘+’) and over a clear land surface area (the red line marked with ‘*’). For the river pixel, because of surface liquid water absorption at 764 nm and the strong atmospheric oxygen absorption at the same wavelength [11,12], the ρ* (764 nm) value is less than 0.05. Land pixels are generally brighter than water pixels, and the ρ* (764 nm) values of land pixels are generally greater than 0.05, based on our detailed study of this EPIC data set and a number of other data sets acquired over different days and seasons.
As illustrated in Figure 1e, the multi-channel apparent reflectances as a function of wavelengths for land and water pixels are affected by atmospheric molecular and aerosol scattering effects and atmospheric gas absorption effects, including ozone and oxygen. These atmospheric effects need to be removed in order to obtain land surface reflectances and water leaving reflectances from EPIC data. Previously, we developed a multi-channel atmospheric correction algorithm for MODIS ocean color applications [8]. In the MODIS algorithm, atmospheric aerosol models and optical depths are derived on the pixel by pixel basis with a spectrum-matching technique using MODIS channels centered near 865, 1240, 1640, and 2130 nm. We have modified the MODIS version of the algorithm for processing EPIC data. The modifications include the generation of atmospheric Rayleigh and aerosol scattering tables and gas absorption tables appropriate for EPIC channels. Because EPIC does not have channels centered at wavelengths longer than 800 nm, we do not retrieve aerosol information from EPIC data. Instead, we derive surface reflectances by assuming climatological aerosol models and optical depths from EPIC data acquired on very clear days.
In brief summary, we have developed an empirical but very effective technique for separating water pixels from land and cloud pixels using the 764-nm channel. Compared to using an existing surface type dataset, this technique also avoids the known geolocation accuracy issue in EPIC L1B data due to the rotation of the Earth during the imaging time frame [4]. We have also adapted a previously developed MODIS version of multi-channel atmospheric correction algorithm for processing EPIC data. Examples of water pixel masking and atmospheric corrections are present in the section below.

3. Results

We have selected three cases for the demonstration of water leaving reflectance retrievals from EPIC data. The first case is for the scene containing the very turbid La Plata River (see Figure 1) in South America. The second case is for a scene containing all the five large lakes in the Great Lake area in North America. The third case is for scenes containing highly productive Caspian Sea, Black Sea, and Azov Sea in Europe.

3.1. La Plata River, South America, March 2017

La Plata River is located in South America. Near the mouth of the river to the Atlantic Ocean, the water is very turbid and often appears in brownish color in satellite images, such as those acquired with the NASA Terra and Aqua MODIS instruments. Because the present ocean color community adopted atmospheric correction algorithms were mainly designed for processing the clear case 1 blue ocean waters, these algorithms generally do not produce reliable water leaving reflectances over turbid case 2 waters as well as over waters with chlorophyll concentration greater than about 10 mg/m3 [13]. We selected the EPIC data set acquired on March 13, 2017 at UTC 1442 in the present study. The full disk RGB image and the 764-nm O2-A channel image have already been shown in Figure 1a,b, respectively. We have made atmospheric corrections to the EPIC data set for three scenarios: Rayleigh only (no aerosol), Rayleigh plus a climatologically averaged continental aerosol model with an optical depth of 0.1 at 550 nm, and Rayleigh plus a climatologically averaged maritime aerosol model with the same optical depth of 0.1 at 550 nm.
Figure 2a is the full disk color image processed from the atmosphere-corrected data of the three EPIC channels (Red: 680 nm; Green: 551 nm; Blue: 443 nm) assuming the maritime aerosol model with an optical depth of 0.1 at 550 nm. Figure 2b is the enlarged portions of the color image for southern portions of South America. Within the red color outlined rectangle, La Plata River is still difficult to see. Figure 2c is the color image for the same portions of South America with land surfaces masked out using the technique described in Section 2. The turbid La Plata River and nearby clear ocean surfaces are shown quite well. The inlet in (c) shows an enlarged portion of La Plata River. Figure 2d shows the spectral reflectances as a function of wavelength for one turbid river water pixel located at the center of the small rectangle outlined in Figure 2c. The line marked with ‘+’ is the top of atmosphere apparent reflectance curve. The lines marked with ‘*’, diamond, and triangle are the retrieved water leaving reflectances for the Rayleigh only, Rayleigh plus rural aerosol model, and Rayleigh plus maritime model corrections, respectively. The shapes and the magnitude of the two curves for the rural aerosol model and maritime model are very similar—indicating that the exact aerosol model used in atmospheric correction for the bright turbid water pixel is not critically important. However, by careful studies of all the four curves in Figure 4d, it is seen that an accurate Rayleigh correction is very important, because the Rayleigh correction removes the largest amounts of scattering effects in the UV (388 nm), blue (443 nm), and green (551 nm) channels.

3.2. Great Lakes, North America, September 2017

Increases in the frequency and magnitude of harmful algae blooms (HABs) over the Great Lake region in North America in recent years have caused serious threats to the fresh and marine aquatic ecosystems and public health risks [14,15]. Satellite remote sensing can help to assess the water quality of the Great Lakes [16]. Because of frequent cloud contamination, it is not easy to find an EPIC scene that is clear over all the five major lakes, i.e., Lake Erie, Lake Michigan, Lake Ontario, Lake Huron, and Lake Superior. After extensive search, we have found one great EPIC scene that was clear over all the five major lakes. The EPIC data set was acquired on 24 September 2017 at UTC 1649.
Figure 3a is the full disk color image processed from the atmosphere-corrected data of the three EPIC channels (Red: 680 nm; Green: 551 nm; Blue: 443 nm) assuming the maritime aerosol model with an optical depth of 0.1 at 550 nm. Figure 3b is the enlarged portion of the color image covering the five large lakes. Figure 3c is the color image covering the same portion of the Great Lake area after masking out the land surfaces and clouds. The left portion of Lake Erie has ‘green’ color because of a major harmful algal bloom event there at the time of EPIC data acquisition [17,18]. Figure 3d shows the spectral reflectances as a function of wavelength for one chlorophyll-rich water pixel in Lake Erie and located exactly at the cross-hair point in Figure 3c. The line marked with ‘+’ is the top of atmosphere apparent reflectance curve. The lines marked with ‘*’, diamond, and triangle are the retrieved water leaving reflectances for the Rayleigh only, Rayleigh plus rural aerosol model, and Rayleigh plus maritime model corrections, respectively. The shape and the magnitude of the two curves for the rural aerosol model and maritime model are very similar, which indicates again that the exact aerosol model used in atmospheric correction for the bright and chlorophyll-rich water pixel is not very important. It is also noted that the green channel (551 nm) reflectance values after atmospheric corrections are larger than those of blue channel (443 nm) and red channel (680 nm). Therefore, we properly recovered the green color of the chlorophyll-rich pixels after atmospheric corrections.
Figure 3e is a true Terra MODIS color image covering approximately the same surface area as that of Figure 3c. The MODIS image with a spatial resolution of about 1 km was downloaded from the NASA MODIS web site. No spatial registration between the Figure 3e,c images is made. By comparing the two images, it is seen that the chlorophyll blooming feature in the left part of Lake Erie in the Figure 3e MODIS image is much sharper than that in the Figure 3c EPIC image. The lower spatial resolution EPIC image with pixel sizes of 24 km or larger smeared out the fine water features.

3.3. Caspian Sea, Black Sea, and Azov Sea, Europe, June 2017

Caspian Sea, Black Sea, and Azov Sea are located in Europe. From late June to early September in a year, the waters in the areas are usually highly productive, and chlorophyll bloom features can often be observed from satellite images. It is relatively easy to find good EPIC data sets, because the areas are often clear and free of obvious cloud contamination. In the third case study reported here, we selected the EPIC data sets acquired on 26 June 2017 with several clear overpasses over the areas during the daytime.
Figure 4a is the full disk color image processed from the atmosphere-corrected data of the three EPIC channels (Red: 680 nm; Green: 551 nm; Blue: 443 nm) assuming the maritime aerosol model with an optical depth of 0.1 at 550 nm. The EPIC data was acquired at UTC 1112. Figure 4b is the magnified portion covering the Caspian Sea, Black Sea, Azov Sea, Persian Gulf, Red Sea, and Mediterranean Sea. Figure 4c is the color image covering the same area as that of Figure 4b but with land surfaces and clouds masked out. Figure 4d shows the spectral reflectances as a function of wavelength for two chlorophyll-rich pixels—one located in the upper right portion of Caspian Sea and the other at the center of Azov Sea, one clear water pixel in eastern portion of Black Sea, and one clear water pixel in eastern part of Mediterranean Sea. It is seen that the green channel (551 nm) reflectance value for the Caspian Sea pixel is almost twice as large as that of the Azov Sea pixel. Figure 4c,d demonstrate that we are able to retrieve properly the color of water surfaces after performing the atmospheric corrections.
Figure 5a–h shows eight sunlit full disk color images processed from atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 26 June 2017 at UTC 0545, 0650, 0756, 0901, 1007, 1112, 1217, and 1323, respectively. The corrections were made by assuming a climatologically averaged maritime aerosol model with an optical depth of 0.1 at 550 nm for all the scenes. The locations of the Caspian Sea in these images are marked with red arrows. These images demonstrate that it is possible to have multi-views of the same area, such as Caspian Sea, with EPIC at the L1 point with a time interval of approximately 1 h and 5 min from the local sunrise to local sunset in a day. Figure 6a–h corresponds to Figure 5a–h, respectively, but covering much smaller areas consisting of Caspian Sea, Black Sea, and Azov Sea and with land and cloud features masked out. The spatial distributions of chlorophyll blooming features in the upper right portions of Caspian Sea and Azov Sea can vary rapidly with time. The multiple view of the areas in a day allows, in principle, the detection of rapid changes of chlorophyll bloom patterns over the areas.
Because the EPIC pixel sizes on the ground change with local view zenith angles and spatial registration errors associated with different channels in a given EPIC scene, we have not tried to make spatial registrations among the eight images in Figure 6a–h to a common latitude and longitude grid for quantitative inter-comparisons among these images. However, the eight images in Figure 6a–h have demonstrated the capability of acquiring multiple views of a given area within a day with EPIC at the L1 point. Such capability is lacking with the present low earth orbiting polar platform satellite instruments (e.g., MODIS, VIIRS).

4. Discussion

At the L1 vantage point, EPIC views the entire sunlit side of the Earth continuously and provides new opportunity for ocean color studies. As mentioned before, the EPIC pixel size at nadir is about 8 km, which is much coarser than the state-of-the-art ocean color sensors, e.g., Visible Infrared Imaging Radiometer Suite (VIIRS) [19]. Data from the L1 point with higher resolution would be much desirable. Resolution can be determined [20] by:
R ≈ 1.22 λ S/D
where R is the detector’s spatial resolution, λ the channel wavelength, S the distance between the sensor and the surface, and D the aperture diameter, which is 30 cm for EPIC [5]. At the L1 point (S = 1.5 million km), the current EPIC instrument can support a 3 km pixel size at the visible spectrum, but to achieve that, the CCD array size (currently 2048 × 2048) needs to be increased.
Another factor affecting data quality is the sensor signal to noise ratio (SNR), which is a function of many parameters, including spatial and spectral resolution, sampling time, sensor optics, detector quantum efficiency etc. The 10 EPIC channels, as listed in Table 1, are all very narrow channels with FWHMs smaller than 3 nm with designed SNR ~200 [5]. These channels were not originally designed for remote sensing of ocean color. The ocean color channels, M1–M7, implemented on the VIIRS [19] instrument currently flying on the Suomi National Polar-orbiting Partnership spacecraft, are much wider with FWHMs of approximately 20 nm. If channels similar to M1–M7 are implemented onto an EPIC-like instrument at the L1 point and keep the same scan rate as that for EPIC, with the state of the art sensor technology and an aperture size of ~100 cm, it is possible to obtain remote sensing data at a spatial resolution of ~1 km or finer with SNR similar to or better than the current EPIC instrument. Such an EPIC-like instrument, if built, would permit full sunlit earth disk observations of ocean color from sunrise to sunset with a repeating cycle of approximately 1 h.

5. Summary

The NASA’s EPIC instrument on board the DSCOVR satellite has been making multiple observations of the entire sunlit Earth in a given day from the Sun-Earth Largangian L1 point since the summer of 2015. The data acquired with EPIC have already been used for the study of the global ozone levels, aerosol index and aerosol optical depth, UV reflectivity of clouds over land and ocean, cloud height over land and ocean, and vegetation indices. Through analysis of EPIC data, we have found that the data acquired with narrow channels centered near 443, 551, and 680 nm, can also have important applications in remote sensing of ocean color. We have modified a version of a multi-channel atmospheric correction algorithm for MODIS ocean color applications and adapted the algorithm for processing EPIC data. We have shown three case studies on water leaving reflectance retrievals from EPIC data acquired over a large turbid river, inland lakes, and oceans. We conclude that a future ocean color instrument on board a satellite at the L1 point, which provides continuous view of the full sunlit disk of the earth, will complement and extend ocean color observations with the low earth observing polar orbital and geostationary satellite instruments in both the spatial and time domains.

Author Contributions

B.-C.G. made initial observations of chlorophyll blooming features over Caspean Sea and Lake Erie from DSCOVR EPIC data sets. R.-R.L. carried out detailed data analysis. Y.Y. provided some DSCOVR EPIC data sets, suggested analysis of time series of EPIC data sets acquired over Caspean Sea, and made sensitivity studies for a possible future multi-channel instrument at the L1 point for remote sensing of ocean color.

Funding

This research is partially supported by the US Office of Naval Research and by a research grant from the Science Mission Directorate of National Aeronautics and Space Administration. The research conducted by Y. Yang is supported by the NASA DSCOVR Science Team Program.

Acknowledgments

The authors are grateful to Weidong Yang of University Space Research Association and to Pengwang Zhai of University of Maryland at Baltimore County for useful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Marshak, A.; Herman, J.; Szabo, A.; Blank, K.; Cede, A.; Carn, S.; Geogdzhaev, D.; Huang, D.; Huang, L.-K.; Knyazikhin, Y.; et al. Earth Observations from DSCOVR/EPIC Instrument. Bull. Am. Meteorol. Soc. 2018, 9, 1829–1850. [Google Scholar] [CrossRef] [PubMed]
  2. Reichhardt, T. Gore wins approval for his Triana satellite. Nature 1998, 396, 5. [Google Scholar] [CrossRef] [PubMed]
  3. Herman, J.; Huang, L.; McPeters, R.; Ziemke, J.; Cede, A.; Blank, K. Synoptic ozone, cloud reflectivity, and erythemal irradiance from sunrise to sunset for the whole earth as viewed by the DSCOVR spacecraft from the earth–sun Lagrange 1 orbit. Atmos. Meas. Technol. 2018, 11, 177–194. [Google Scholar] [CrossRef] [Green Version]
  4. Yang, Y.; Meyer, G.; Wind, Y.; Zhou, A.; Marshak, S.; Platnick, Q.; Min, A.B.; Davis, J.; Joiner, A.; Vasilkov, D.; et al. Cloud Products from Earth Polychromatic Imaging Camera (EPIC) observations: Algorithm and Initial Evaluation. Atmos. Meas. Technol. Discuss. 2018. [Google Scholar] [CrossRef]
  5. Valero, F.P.J.; Herman, J.; Minnis, P.; Collins, W.D.; Sadourny, R.; Wiscombe, W.; Lubin, D.; Ogilvie, K. Triana—A Deep Space Earth and Solar Observatory; NASA Background Report; NASA: Washington, DC, USA, 1999.
  6. Salomonson, V.V.; Barnes, W.L.; Maymon, P.W.; Montgomery, H.E.; Ostrow, H. MODIS: Advanced facility instrument for studies of the earth as a system. IEEE Trans. Geosci. Remote Sens. 1989, 27, 145–153. [Google Scholar] [CrossRef]
  7. King, M.D.; Menzel, W.P.; Kaufman, Y.J.; Tanré, D.; Gao, B.C.; Platnick, S.; Ackerman, S.A.; Remer, L.A.; Pincus, R.; Hubanks, P.A. Cloud and aerosol properties, precipitable water, and profiles of temperature and humidity from MODIS. IEEE Trans. Geosci. Remote Sens. 2003, 41, 442–458. [Google Scholar] [CrossRef]
  8. Gao, B.-C.; Montes, M.J.; Li, R.-R.; Dierssen, H.M.; Davis, C.O. An atmospheric correction algorithm for remote sensing of bright coastal waters using MODIS land and ocean channels in the solar spectral region. IEEE Trans. Geosci. Remote Sens. 2007, 45, 1835–1843. [Google Scholar] [CrossRef]
  9. Tanre, D.; Deroo, D.C.; Duhaut, P.; Herman, M.; Morcrette, J.J.; Perbos, J.; Deschamps, P.Y. Description of a computer code to simulate the satellite signal in the solar spectrum: The 5S code. Int. J. Remote Sens. 1990, 11, 659–668. [Google Scholar] [CrossRef]
  10. Gao, B.-C.; Montes, M.J.; Ahmad, Z.; Davis, C.O. Atmospheric correction algorithm for hyperspectral remote sensing of ocean color from space. Appl. Opt. 2000, 39, 887–896. [Google Scholar] [CrossRef] [PubMed]
  11. Xu, X.; Wang, J.; Wang, Y.; Zeng, J.; Torres, O.; Yang, Y.; Marshak, A.; Reid, J.; Miller, S. Passive remote sensing of altitude and optical depth of dust plumes using the oxygen A and B bands: First results from EPIC/DSCOVR at Lagrange-1 point. Geophys. Res. Lett. 2017, 44, 7544–7554. [Google Scholar] [CrossRef]
  12. Yang, Y.; Marshak, A.; Mao, J.; Lyapustin, A.; Herman, J. A Method of Retrieving Cloud Top Height and Cloud Geometrical Thickness with Oxygen A and B bands for the Deep Space Climate Observatory (DSCOVR) Mission: Radiative Transfer Simulations. J. Quant. Spectrosc. Radiat. Trans. 2013, 122, 141–149. [Google Scholar] [CrossRef]
  13. Jaelani, L.M.; Matsushita, B.; Yang, W.; Fukushima, T. Evaluation of four MERIS atmospheric correction algorithms in Lake Kasumigaura, Japan. Int. J. Remote Sens. 2013, 34, 8967–8985. [Google Scholar] [CrossRef]
  14. Pick, F.R. Blooming algae: A Canadian perspective on the rise of toxic cyanobacteria. Can. J. Fish. Aquat. Sci. 2016, 73, 1149–1158. [Google Scholar] [CrossRef]
  15. Winter, J.G.; DeSellas, A.M.; Fletcher, R.; Heintsch, L.; Morley, A.; Nakamoto, L.; Utsumi, K. Algae blooms in Ontario, Canada: Increases in reports since 1994. Lake Reserv. Manag. 2011, 27, 107–114. [Google Scholar] [CrossRef]
  16. Binding, C.E.; Greenberg, T.A.; McCullough, G.; Watson, S.B.; Page, E. An analysis of satellite-derived chlorophyll and algal bloom indices on Lake Winnipeg. J. Great Lakes Res. 2018, 44, 436–446. [Google Scholar] [CrossRef]
  17. Available online: https://earthobservatory.nasa.gov/images/91017/bloom-persists-in-lake-erie (accessed on 8 March 2019).
  18. Available online: https://www.sciencedirect.com/science/article/pii/S0380133018301874 (accessed on 8 March 2019).
  19. Xiong, X.X.; Butler, J.; Chiang, K.; Efremova, B.; Fulbright, J.; Lei, N.; McIntire, J.; Oudrari, H.; Sun, J.; Wang, Z.; et al. VIIRS on-orbit calibration methodology and performance. J. Grophys. Res. Atmos. 2013, 119, 5065–5078. [Google Scholar] [CrossRef]
  20. Kidder, S.Q.; Vonder Haar, T.H. Satellite Meteorology—An Introduction; Academic Press: San Diego, CA, USA, 1995. [Google Scholar]
Figure 1. (a) The sunlit full disk color image processed from three EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 13 March 2013 at UTC 1442; (b) The full disk 764-nm channel image; (c) The portions of color image over southern part of S. America; (d) The portions of the 764-nm image over southern part of S. America; (e) Apparent reflectances as a function of wavelength for a turbid river water pixel (blue line marked with ‘+’) and a land pixel (red line marked with ‘*’). Pixels are selected in the red box marked in Panel (C) and (D).
Figure 1. (a) The sunlit full disk color image processed from three EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 13 March 2013 at UTC 1442; (b) The full disk 764-nm channel image; (c) The portions of color image over southern part of S. America; (d) The portions of the 764-nm image over southern part of S. America; (e) Apparent reflectances as a function of wavelength for a turbid river water pixel (blue line marked with ‘+’) and a land pixel (red line marked with ‘*’). Pixels are selected in the red box marked in Panel (C) and (D).
Sensors 19 01243 g001
Figure 2. (a) The sunlit full disk color image processed from three atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 13 March 2017 at UTC 1442; (b) The portions of color image over S. America; (c) The portions of color image of water surfaces over S. America after masking out land and cloud pixels; (d) Spectral reflectances as a function of wavelength for one turbid river water pixel located at the center of the small rectangle outlined in (c). In (d), the line marked with ‘+’ is the top of atmosphere apparent reflectance curve. The lines marked with ‘*’, diamond, and triangle are the retrieved water leaving reflectances for the Rayleigh only, Rayleigh plus rural aerosol model, and Rayleigh plus maritime model corrections, respectively. The inlet in (c) shows an enlarged portion of La Plata River.
Figure 2. (a) The sunlit full disk color image processed from three atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 13 March 2017 at UTC 1442; (b) The portions of color image over S. America; (c) The portions of color image of water surfaces over S. America after masking out land and cloud pixels; (d) Spectral reflectances as a function of wavelength for one turbid river water pixel located at the center of the small rectangle outlined in (c). In (d), the line marked with ‘+’ is the top of atmosphere apparent reflectance curve. The lines marked with ‘*’, diamond, and triangle are the retrieved water leaving reflectances for the Rayleigh only, Rayleigh plus rural aerosol model, and Rayleigh plus maritime model corrections, respectively. The inlet in (c) shows an enlarged portion of La Plata River.
Sensors 19 01243 g002
Figure 3. (a) The sunlit full disk color image processed from three atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 24 September 2017 at UTC 1649; (b) The portion of color image over Great Lake region; (c) The portion of color image of water surfaces over Great Lakes after masking out land and cloud pixels; (d) Spectral reflectances as a function of wavelength for one chlorophyll-rich water pixel located in the left part of Lake Erie. (e) The 1-km resolution true color Terra MODIS image over the Great Lake area. In (d), the line marked with ‘+’ is the top of atmosphere apparent reflectance curve. The lines marked with ‘*’, diamond, and triangle are the retrieved water leaving reflectances for the Rayleigh only, Rayleigh plus rural aerosol model, and Rayleigh plus maritime model corrections, respectively.
Figure 3. (a) The sunlit full disk color image processed from three atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 24 September 2017 at UTC 1649; (b) The portion of color image over Great Lake region; (c) The portion of color image of water surfaces over Great Lakes after masking out land and cloud pixels; (d) Spectral reflectances as a function of wavelength for one chlorophyll-rich water pixel located in the left part of Lake Erie. (e) The 1-km resolution true color Terra MODIS image over the Great Lake area. In (d), the line marked with ‘+’ is the top of atmosphere apparent reflectance curve. The lines marked with ‘*’, diamond, and triangle are the retrieved water leaving reflectances for the Rayleigh only, Rayleigh plus rural aerosol model, and Rayleigh plus maritime model corrections, respectively.
Sensors 19 01243 g003
Figure 4. (a) The sunlit full disk color image processed from three atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 26 June 2017 at UTC 1112; (b) The portion of image over Caspian Sea, Black Sea, Azov Sea, Persian Gulf, Red Sea, and Mediterranean Sea; (c) The portion of image of water surfaces over the same area as that of (b), but after masking out the land and cloud pixels; (d) Spectral reflectances as a function of wavelength for two chlorophyll-rich water pixels (Green lines) over Caspian Sea and Azov Sea and two blue water pixels (Blue lines) over eastern portions of Black Sea and Mediterranean Sea.
Figure 4. (a) The sunlit full disk color image processed from three atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 26 June 2017 at UTC 1112; (b) The portion of image over Caspian Sea, Black Sea, Azov Sea, Persian Gulf, Red Sea, and Mediterranean Sea; (c) The portion of image of water surfaces over the same area as that of (b), but after masking out the land and cloud pixels; (d) Spectral reflectances as a function of wavelength for two chlorophyll-rich water pixels (Green lines) over Caspian Sea and Azov Sea and two blue water pixels (Blue lines) over eastern portions of Black Sea and Mediterranean Sea.
Sensors 19 01243 g004
Figure 5. (ah) Eight sunlit full disk color images processed from atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 26 June 2017 at UTC 0545, 0650, 0756, 0901, 1007, 1112, 1217, and 1323, respectively. The spatial locations of the Caspian Sea in these images are pointed to with red arrows.
Figure 5. (ah) Eight sunlit full disk color images processed from atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 26 June 2017 at UTC 0545, 0650, 0756, 0901, 1007, 1112, 1217, and 1323, respectively. The spatial locations of the Caspian Sea in these images are pointed to with red arrows.
Sensors 19 01243 g005
Figure 6. (ah) Eight color images processed from atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 26 June 2017 at UTC 0545, 0650, 0756, 0901, 1007, 1112, 1217, and 1323, respectively, and covering areas consisting of Caspian Sea, Black Sea, and Azov Sea with land and cloud features masked out.
Figure 6. (ah) Eight color images processed from atmosphere corrected EPIC channel (Red: 680 nm; Green: 551 nm; Blue: 443 nm) images acquired on 26 June 2017 at UTC 0545, 0650, 0756, 0901, 1007, 1112, 1217, and 1323, respectively, and covering areas consisting of Caspian Sea, Black Sea, and Azov Sea with land and cloud features masked out.
Sensors 19 01243 g006
Table 1. Specifications of 10 EPIC channels and their primary applications.
Table 1. Specifications of 10 EPIC channels and their primary applications.
Spectral Channels (nm)FWHM (nm)Primary Usage
317.51Ozone, SO2
3252Ozone
3403Ozone, Aerosols
3883Aerosols, Clouds
4433Aerosols
5513Aerosols, Vegetation
6802Aerosols, Vegetation, Clouds
687.750.8Cloud Height
7641Cloud Height
779.52Clouds

Share and Cite

MDPI and ACS Style

Gao, B.-C.; Li, R.-R.; Yang, Y. Remote Sensing of Daytime Water Leaving Reflectances of Oceans and Large Inland Lakes from EPIC onboard the DSCOVR Spacecraft at Lagrange-1 Point. Sensors 2019, 19, 1243. https://doi.org/10.3390/s19051243

AMA Style

Gao B-C, Li R-R, Yang Y. Remote Sensing of Daytime Water Leaving Reflectances of Oceans and Large Inland Lakes from EPIC onboard the DSCOVR Spacecraft at Lagrange-1 Point. Sensors. 2019; 19(5):1243. https://doi.org/10.3390/s19051243

Chicago/Turabian Style

Gao, Bo-Cai, Rong-Rong Li, and Yuekui Yang. 2019. "Remote Sensing of Daytime Water Leaving Reflectances of Oceans and Large Inland Lakes from EPIC onboard the DSCOVR Spacecraft at Lagrange-1 Point" Sensors 19, no. 5: 1243. https://doi.org/10.3390/s19051243

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