<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJMS</journal-id><journal-title-group><journal-title>Open Journal of Marine Science</journal-title></journal-title-group><issn pub-type="epub">2161-7384</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojms.2015.54035</article-id><article-id pub-id-type="publisher-id">OJMS-60622</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Retrieval of Aerosol Optical Depth from Oceansat-2 OCM
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>daya</surname><given-names>Lakshmi Sriperambudur</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Chandralingam</surname><given-names>Sonnati</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Nagamani</surname><given-names>Pullaiahgari Venkata</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Preethi</surname><given-names>Latha Thalathoti</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Physics, JNTUH, Hyderabad, India</addr-line></aff><aff id="aff2"><addr-line>Ocean Sciences Group, NRSC, Balanagar, Hyderabad, India</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>pvnagamani@gmail.com(NPV)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>17</day><month>09</month><year>2015</year></pub-date><volume>05</volume><issue>04</issue><fpage>443</fpage><lpage>454</lpage><history><date date-type="received"><day>1</day>	<month>September</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>24</month>	<year>October</year>	</date><date date-type="accepted"><day>27</day>	<month>October</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Aerosols are the tiny suspended particles in the atmosphere playing a major role in influencing the net radiation budget of the earth. The aerosols also affect cloud microphysics and those with land origin, possibly reduce the monsoon rainfall. Dynamic and diverse nature of the aerosols varies according to different locations. The aerosols originating due to desert dust, rural atmospheric situations, urban pollutants and marine areas have wide variations and display specific characteristics. Routine monitoring of aerosol events and their subsequent dispersal pattern are important in order to understand their role in climatic process. Hence, it is very important to study the aerosols and their retrieval strategy from the ocean colour remote sensing sensors. The satellite sensors provide platform for making observations covering large area as also their short-term and frequent repetivity. Ocean-colour sensors e.g. CZCS, SeaWiFS, MODIS, POLDAR, Oceansat-1
   
  &amp; 2 OCM have been used to study aerosols, apart from being used to study ocean-colour. Most of the ocean-colour sensors are equipped with a few additional near infrared (NIR) bands (
  λ
   
  &gt; 700 nm), which are helpful in providing vital information on atmospheric aerosols due to strong absorption by water in NIR wavelengths. The present work is an attempt to study the temporal and spatial variations of Aerosol Optical Depth (AOD) over the Bay of Bengal using Oceansat-2 Ocean Colour Monitor (OCM).
 
</p></abstract><kwd-group><kwd>AOD</kwd><kwd> OCM-2</kwd><kwd> Atmospheric Correction</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Atmospheric aerosol plays a significant role in the Earth’s radiation budget through radiative forcing and chemical perturbations. The net effect of aerosols is to cool the climate system by reflecting sunlight [<xref ref-type="bibr" rid="scirp.60622-ref1">1</xref>] . Quantifying the net effect requires accurate information on the distribution of aerosol properties that have to be estimated from satellite observations. Estimating aerosol properties is also one of the first and important steps in generating high-level ocean colour products from satellite observations. The major problem in their characterization is on account of their short lifetime because of which they have high spatial and temporal variability [<xref ref-type="bibr" rid="scirp.60622-ref2">2</xref>] . The transport processes may bring in aerosols from other locations and affect the local climate there. These factors reinforce the necessity of aerosol monitoring on a larger spatial scale than can be provided by the ground based measurements (IPCC, 2001). Satellite based observations can provide detailed knowledge in this regard on a long timescale covering a large spatial area [<xref ref-type="bibr" rid="scirp.60622-ref3">3</xref>] . Aerosol monitoring from space based instruments consists in extracting the atmospheric contribution from the total signal measured by the satellite sensor. Aerosol monitoring from previous sensors was limited to studies over oceans which have a distinct advantage in that the total measured signal is not much affected by reflectance from ocean surface away from sun-glint area [<xref ref-type="bibr" rid="scirp.60622-ref4">4</xref>] . Aerosol retrieval over oceans is thus more accurate and reliable.</p><p>Aerosol particles are important to scientists because they represent an area of great uncertainty in their efforts to understand the Earth’s climate system. Depending upon their size, type, and location, aerosols can either cool the surface, or warm it. They can help clouds to form, or they can inhibit cloud formation. And if inhaled, aerosols can be harmful to people’s health. There are many applications for aerosol optical thickness data: 1) atmospheric correction of remotely sensed surface features; 2) monitoring of sources and sinks of aerosols; 3) monitoring of volcanic eruptions and forest fire; 4) radiative transfer model; 5) air quality; 6) health and environment; 7) earth radiation budget; 8) climate change [<xref ref-type="bibr" rid="scirp.60622-ref5">5</xref>] .</p></sec><sec id="s2"><title>2. OCM-2 Instrument</title><p>OCM-2 sensor is the follow on sensor of OCM-1 and is designed to measure the spectral variability of the water leaving radiance that is related to the concentration of phytoplankton pigments, suspended matter and coloured dissolved organic matter in coastal and ocean waters, and to characterize the atmospheric aerosols. Oceansat-2, OCM (OCM-2) is identical to OCM-1, except minor spectral shift in band 6 and 7. In OCM-1, band-6 centered at 670 nm is shifted to 620-nm in OCM-2 for better quantification of suspended sediments and band 7 centered at 765-nm is shifted to 740-nm to avoid oxygen absorption as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. The technical details and algorithms developed for retrieving the geophysical data products from OCM-2 can be obtained from [<xref ref-type="bibr" rid="scirp.60622-ref6">6</xref>] . The coverage for LAC and GAC is shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>. The technical specifications of OCM-1 and OCM-2 are presented in <xref ref-type="table" rid="table1">Table 1</xref>.</p></sec><sec id="s3"><title>3. Geophysical Parameters from OCM-2 Data</title><p>The geophysical products from Oceansat-II OCM data at basic spatial resolution of 360 meters along with their specifications are shown in the <xref ref-type="table" rid="table2">Table 2</xref>. The normalized water-leaving radiance in 412, 443, 490, 510 and</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Spectral changes in Oceansat-2 ocean colour monitor (OCM)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1470216x5.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Pictorial representation of (a) local area overage (b) global area coverage (GAC) for Oceansat-2 ocean colour monitor</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1470216x6.png"/></fig><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Specifications of the two ocean colour sensors OCM-1 and OCM-2</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Parameters</th><th align="center" valign="middle" >OCM-2 Specifications</th><th align="center" valign="middle" >OCM-2 Specifications</th></tr></thead><tr><td align="center" valign="middle" >IGFOV at nominal altitude (m)</td><td align="center" valign="middle" >360 &#215; 250</td><td align="center" valign="middle" >360 &#215; 250</td></tr><tr><td align="center" valign="middle" >Swath (km)</td><td align="center" valign="middle" >1420</td><td align="center" valign="middle" >1420</td></tr><tr><td align="center" valign="middle" >Acquisition Mode</td><td align="center" valign="middle" >Push Broom</td><td align="center" valign="middle" >Push Broom</td></tr><tr><td align="center" valign="middle" >3. No. of spectral bands</td><td align="center" valign="middle" >8</td><td align="center" valign="middle" >8</td></tr><tr><td align="center" valign="middle" >4. Spectral range (nm)</td><td align="center" valign="middle" >402 - 885</td><td align="center" valign="middle" >402 - 885</td></tr><tr><td align="center" valign="middle" >5. Spectral bands</td><td align="center" valign="middle" >B1: 402 - 422 nm B2: 433 - 453 nm B3: 480 - 500 nm B4: 500 - 520 nm B5: 545 - 565 nm B6: 660 - 680 nm B7: 745 - 785 nm B8: 845 - 885 nm</td><td align="center" valign="middle" >B1: 404 - 424 nm B2: 431 - 451 nm B3: 476 - 496 nm B4: 500 - 520 nm B5: 546 - 566 nm B6: 610 - 630 nm B7: 725 - 755 nm B8: 845 - 885 nm</td></tr><tr><td align="center" valign="middle" >6. Quantization Bits</td><td align="center" valign="middle" >12</td><td align="center" valign="middle" >12</td></tr><tr><td align="center" valign="middle" >7. Along track steering</td><td align="center" valign="middle" >&#177;200</td><td align="center" valign="middle" >&#177;200</td></tr><tr><td align="center" valign="middle" >8. Data acquisition modes</td><td align="center" valign="middle" >Local Area Coverage (LAC) ? 360 m &#215; 236 m Resolution</td><td align="center" valign="middle" >Local Area Coverage (LAC) ? 360 m &#215; 236 m Resolution Global Area Coverage (GAC) ? 1 km &#215; 1 km Resolution</td></tr><tr><td align="center" valign="middle" >Data Formats</td><td align="center" valign="middle" >Super Structure</td><td align="center" valign="middle" >HDF</td></tr><tr><td align="center" valign="middle" >Data Processing</td><td align="center" valign="middle" >In-house</td><td align="center" valign="middle" >SeaDAS</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Geophysical products from Oceansat-II OCM data</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >S. No.</th><th align="center" valign="middle" >Parameter</th><th align="center" valign="middle" >Variable Range</th><th align="center" valign="middle" >Targeted error budget</th></tr></thead><tr><td align="center" valign="middle" >1.</td><td align="center" valign="middle" >Normalized water leaving radiance (nL<sub>w</sub>) in 412, 443, 490, 510, 555 and 620 nm</td><td align="center" valign="middle" >0.0 - 5.0 W∙cm<sup>−2</sup>∙nm<sup>−1</sup>∙sr<sup>−1 </sup> (Variable for different bands)</td><td align="center" valign="middle" >&lt;5%</td></tr><tr><td align="center" valign="middle" >2.</td><td align="center" valign="middle" >Chlorophyll-a concentration</td><td align="center" valign="middle" >0.05 - 30.0 mg∙m<sup>−3 </sup> (For case―1 waters)</td><td align="center" valign="middle" >&lt;30%</td></tr><tr><td align="center" valign="middle" >3.</td><td align="center" valign="middle" >Vertical Diffuse attenuation coefficient (K<sub>d</sub>) at 490-nm</td><td align="center" valign="middle" >0.01 - 0.50 m<sup>−1</sup></td><td align="center" valign="middle" >&lt; 15%</td></tr><tr><td align="center" valign="middle" >4.</td><td align="center" valign="middle" >Total Suspended matter (TSM) concentration</td><td align="center" valign="middle" >0.0 - 200 mg L<sup>−1</sup></td><td align="center" valign="middle" >&lt;20 %</td></tr><tr><td align="center" valign="middle" >5.</td><td align="center" valign="middle" >Aerosol Optical Depth (AOD) over oceans at 865-nm</td><td align="center" valign="middle" >0.0 - 1.0</td><td align="center" valign="middle" >&lt;20%</td></tr></tbody></table></table-wrap><p>555-nm bands will be the output of the atmospheric correction procedure and will be subsequently used in the estimation of bio-geo-physical variables. The chlorophyll-α product can be used as a proxy to the presence of phytoplankton in the seawater as chlorophyll-α is the dominant pigment found in most of the phytoplankton. The diffuse attenuation coefficient (K<sub>d</sub>) product is an apparent optical property, which defines the rate of decrease of downwelling irradiance falling with depth in the water column. The total suspended matter (TSM) product will provide quantitative measure of the inorganic particulate matter present in the suspended form mainly in the coastal waters. The aerosol optical depth (AOD) is a measure of the atmospheric turbidity, which will be characterized at 865-nm, as this wavelength is almost insensitive to reflectance from ocean waters.</p></sec><sec id="s4"><title>4. Retrieval Strategy</title><p>The geophysical retrieval procedures will be applied only to cloud free water pixels for daytime conditions over open-ocean and coastal waters that are not affected by sun glint. If a cloud mask, land mask or sun glint mask is not set for a given pixel geophysical retrieval for OCM-2 top-of-atmosphere (TOA) radiances is performed. To perform atmospheric correction, non-OCM-2 data sets such as total ozone amount and sea surface wind will also be required. The output of the atmospheric correction algorithm will be water-leaving reflectance or remote sensing reflectance (R<sub>rs</sub>), which will form input to the in-water bio-optical algorithms. The atmospheric correction algorithm over shallow water or turbid water will be performed by “borrowing” the aerosol properties over the nearest deep water.</p></sec><sec id="s5"><title>5. Cloud Masking</title><p>Seawater, and its constituents alter spectral composition of visible (400 - 700 nm) radiation reaching the satellite sensor. However, near infrared (NIR) radiation is strongly absorbed by relatively clear water, resulting in a uniformly low albedo. Clouds, in contrast, have a wide range of reflectivity at the wavelength measured by OCM sensor. An albedo-based cloud masking approach has been adopted for OCM processing by making uses of NIR band centred at 865 nm. The albedo at 865-nm, α<sub>865</sub> is calculated using the following equations:</p><disp-formula id="scirp.60622-formula1513"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x7.png"  xlink:type="simple"/></disp-formula><p>where, L<sub>t</sub>(865) is the total radiance at the sensor at 865-nm, L<sub>i</sub>(865) is the incident light at 865-nm, t(865, λ) is the diffuse transmittance between the surface and sensor at 865-nm, and θ is the satellite zenith angle [<xref ref-type="bibr" rid="scirp.60622-ref7">7</xref>] . The diffuse transmittance between the surface and sensor at 865 nm is calculated as</p><disp-formula id="scirp.60622-formula1514"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x8.png"  xlink:type="simple"/></disp-formula><p>where τ<sub>r</sub>, τ<sub>oz</sub> are the Rayleigh and ozone optical thicknesses respectively, at 865-nm. The incident light at 865-nm is computed as</p><disp-formula id="scirp.60622-formula1515"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x9.png"  xlink:type="simple"/></disp-formula><p>where F<sub>0</sub> is the extra-terrestrial solar irradiance at 865-nm and θ<sub>0</sub> is the solar zenith angle. The diffuse transmittance between the sun and the ocean surface at 865-nm follows from Equation (2), is calculated as;</p><disp-formula id="scirp.60622-formula1516"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x10.png"  xlink:type="simple"/></disp-formula><p>It was found that OCM image pixels of 865-nm band having albedo greater than 1.1% would be masked successfully for clouds/thin clouds and haze. This threshold value is also found to be effective for masking land and sun glint affected pixels as well.</p></sec><sec id="s6"><title>6. Sun Glint Mask and Correction</title><p>Sun glint refers to the phenomenon of incoming solar radiation directly reflected from the ocean surface to the sensor. For an absolutely flat ocean surface, the sun glint occurs at one point where the zenith angles of sun and sensor are identical and their azimuth angles are opposite. However, the ocean surface is never absolutely flat. The wind-derived surface roughness enlarges the sun glint area. For remote sensing of the ocean and atmosphere optical properties, the measurement of radiances affected by sun glint has to be avoided and masked out. There are usually no meaningful retrievals in regions significantly contaminated by sun glint. The Oceansat-2 OCM is capable of operationally tilting the sensor &#177; 20˚ away from the nadir to minimize sun glint contamination. An operational scheme is developed and optimal tilt angles have been identified to minimize the sun-glint in different seasons for the Arabian Sea and Bay of Bengal for Oceansat-1 OCM data [<xref ref-type="bibr" rid="scirp.60622-ref8">8</xref>] . Apart from implementing this technique some OCM data showed sun-glint in the months of April and August. A method proposed by [<xref ref-type="bibr" rid="scirp.60622-ref9">9</xref>] to determine when sun glint needs to be masked out. The probability of a pixel being contaminated by sun glitter for the incoming sun direction (θ<sub>0</sub>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x11.png" xlink:type="simple"/></inline-formula>) and the observing direction (θ,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x12.png" xlink:type="simple"/></inline-formula>) can be written as</p><disp-formula id="scirp.60622-formula1517"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x13.png"  xlink:type="simple"/></disp-formula><p>where, σ<sup>2</sup> is the mean square surface slope [<xref ref-type="bibr" rid="scirp.60622-ref10">10</xref>] which is a function of wind speed, W<sub>s </sub></p><disp-formula id="scirp.60622-formula1518"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x14.png"  xlink:type="simple"/></disp-formula><p>A threshold value of 1.5 percent for the probability to mask out the sun glint area will be used.</p></sec><sec id="s7"><title>7. Atmospheric Correction Algorithm Description</title><sec id="s7_1"><title>7.1. Processing Outline</title><p>The Ocean colour data processing for the retrieval of the geophysical parameters is a two-step procedure. The first procedure involves the atmospheric correction of the top of the atmosphere (TOA) sensor detected radiances to estimate the sea surface reflectivity. The output of this atmospheric correction procedure will be subsequently used for the estimation of in-water and atmospheric variables.</p><p><xref ref-type="fig" rid="fig3">Figure 3</xref> presents a schematic flowchart for the Gordon-Wang atmospheric correction algorithm, adapted from [<xref ref-type="bibr" rid="scirp.60622-ref11">11</xref>] and also subsequently used for OCM data processing by [<xref ref-type="bibr" rid="scirp.60622-ref12">12</xref>] . The total radiance L<sub>t</sub>, measured at the top of the atmosphere in each of the visible to near-infrared bands is divided by the extraterrestrial solar irradiance F<sub>0</sub> to obtain the measured reflectance, ρ<sub>t</sub>. The extraterrestrial solar irradiance (F<sub>0</sub>) values will be computed for the OCM spectral response functions (SRF). The reflectance contributed by whitecaps is estimated from the surface wind speed W and subtracted from ρ<sub>t</sub>. Corrections for ozone absorption and Rayleigh reflectance ρ<sub>r</sub> are applied to obtain the ρ<sub>t</sub> − ρ<sub>r</sub> term, which is considered equivalent to aerosol reflectance in the two NIR bands of OCM-2 data i.e. 740-nm and 865-nm. The algorithm then calculates the spectral dependency of the candidate aerosol type based on the value of epsilon parameter, which is ratio of aerosol reflectance in the 740 and 865-nm bands, assuming that the water-leaving reflectance in each of these bands is zero due to strong absorption of light by water beyond 700-nm. Extrapolation of the aerosol reflectance in shorter wavelength for each pixel is done using an exponential model for aerosol spectral behavior with epsilon parameter as an exponent. After subtraction of the aerosol contribution, the water-leaving reflectance is obtained in each of the visible bands by dividing by the diffuse atmospheric transmittance.</p></sec><sec id="s7_2"><title>7.2. Theoretical Description of the Atmospheric Correction Algorithm</title><p>The reflectance backscattered from the atmosphere and/or sea surface is typically at least an order of magnitude</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Flow chart for the algorithm for the atmospheric correction of OCM data</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1470216x15.png"/></fig><p>larger than the desired radiance scattered out of the water. The contributions of the water-leaving reflectance to the TOA reflectance decreases with the increase of the viewing angle because of the reduction of the diffuse transmittance. Therefore, it is mandatory to correct for atmospheric effects, to retrieve any quantitative parameter from space. This section describes the basis of the algorithm for removing the atmospheric effects from the OCM-2 imagery over the ocean to derive normalised water-leaving radiance in the visible range of electromagnetic spectrum. The normalised water-leaving radiance, [L<sub>w</sub>]<sub>N</sub>, has been defined by [<xref ref-type="bibr" rid="scirp.60622-ref13">13</xref>] in a following manner,</p><disp-formula id="scirp.60622-formula1519"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x16.png"  xlink:type="simple"/></disp-formula><p>where, L<sub>w</sub>(λ) is the radiance backscattered out of the water at a wavelength λ, τ<sub>r</sub>(λ) and τ<sub>oz</sub>(λ) are the optical thickness of the atmosphere associated with molecular (Rayleigh) scattering and Ozone absorption, respectively and θ<sub>0</sub> is the solar zenith angle. The normalised water leaving radiance is approximately the radiance that would exit the ocean in the absence of the atmosphere with the sun at the zenith.</p><p>The radiance received by a ocean colour sensor at the top of the atmosphere (TOA) in a spectral band centred at wavelength λ, L<sub>t</sub>(λ), can be divided into the following components: L<sub>path</sub>(λ) the radiance generated along the optical path by scattering in the atmosphere and by specular reflection of the atmospherically scattered light from the sea surface; L<sub>g</sub>(λ) the contribution arising from specular reflection from the sea surface (sun glitter), and, L<sub>w</sub>(λ) the desired water-leaving radiance; i.e.,</p><disp-formula id="scirp.60622-formula1520"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x17.png"  xlink:type="simple"/></disp-formula><p>Converting the Equation (8) in terms of reflectance, where, the reflectance ρ is defined as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x18.png" xlink:type="simple"/></inline-formula>, where L is the radiance in a given solar and viewing geometry, F<sub>0</sub> is extra-terrestrial solar irradiance, and μ<sub>0</sub> is the cosine of the solar zenith angle, the Equation (8) can be rewritten as</p><disp-formula id="scirp.60622-formula1521"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x19.png"  xlink:type="simple"/></disp-formula><p>where, ρ<sub>a</sub> and ρ<sub>r</sub> are the reflectance generated along the optical path by scattering in the atmosphere due to aerosol and Rayleigh scattering, ρ<sub>ra</sub>(λ) is a multiple interaction term between molecules and aerosols, ρ<sub>g</sub> is the specular reflection or sun glitter component and ρ<sub>w</sub> is the desired water leaving reflectance. In Equation (9), T and t are the direct and diffuse transmittance terms, respectively. The direct component of the transmittance T is associated with the sun glitter term, which is highly directional. The direct transmittance is given by</p><disp-formula id="scirp.60622-formula1522"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x20.png"  xlink:type="simple"/></disp-formula><p>where μ<sub>v</sub> is the satellite viewing angle and τ<sub>r</sub>, τ<sub>a</sub>, τ<sub>oz</sub> are, the Rayleigh, aerosol, and Ozone optical thicknesses. The diffuse transmittance t, associated with water leaving reflectance term in Equation (9) is given by [<xref ref-type="bibr" rid="scirp.60622-ref14">14</xref>] as</p><disp-formula id="scirp.60622-formula1523"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x21.png"  xlink:type="simple"/></disp-formula><p>The Rayleigh reflectance path radiance, ρ<sub>r</sub>, is computed by the following equation,</p><disp-formula id="scirp.60622-formula1524"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x22.png"  xlink:type="simple"/></disp-formula><p>where, ω<sub>0r</sub> is single scattering albedo for molecular scattering, τ<sub>r</sub> is the spectral Rayleigh optical depth, F<sub>0</sub> is extra-terrestrial solar flux, P<sub>r</sub> is the molecular scattering phase function and m<sub>v</sub> = cosq<sub>v</sub>. The Rayleigh scattering phase function is given as [<xref ref-type="bibr" rid="scirp.60622-ref15">15</xref>]</p><disp-formula id="scirp.60622-formula1525"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x23.png"  xlink:type="simple"/></disp-formula><p>The forward/backward scattering angle g<sup>&#177;</sup> in Equation (13) is computed as</p><disp-formula id="scirp.60622-formula1526"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x24.png"  xlink:type="simple"/></disp-formula><p>where f is the relative azimuth (sensor azimuth―sun’s azimuth), θ<sub>v</sub> is the sensor viewing angle and θ<sub>0</sub> is the sun zenith angle (illumination directions). The variability of the mean extraterrestrial solar flux, F<sub>0</sub>, for seasonal variation of the Earth-Sun distance is calculated by following equation</p><disp-formula id="scirp.60622-formula1527"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x25.png"  xlink:type="simple"/></disp-formula><p>where, d<sup>2</sup> is the factor accounting for the actual Sun-Earth distance as a function of the Julian day of the year and μ<sub>0</sub> is the cosine of the solar zenith angle. The optical depth τ<sub>oz</sub>(λ) can be assumed constant over time and space around Indian region and nominal values of τ<sub>oz</sub> for OCM bands 1 to 8 are 0, 0.00163, 0.0090, 0.0193, 0.0364, 0.0405, 0.0040, and 0, respectively.</p>The Single Scattering Approximation<p>In Equation (9), ρ<sub>ra</sub> is the interaction term between molecular and aerosol scattering [<xref ref-type="bibr" rid="scirp.60622-ref16">16</xref>] . The term ρ<sub>ra</sub> accounts for the interaction between Rayleigh and aerosol scattering, e.g., photons first scattered by the air then scattered by aerosols, or photons first scattered by aerosols then air, etc. This term is considered as zero in the single scattering case, in which photons are only scattered once, and it can be ignored as long as the amount of multiple scattering is small, i.e., at small Rayleigh and aerosol optical thicknesses and [<xref ref-type="bibr" rid="scirp.60622-ref17">17</xref>] has shown that it is useful to consider path radiance term in the limit that the optical thickness of the atmosphere is =1. This is referred as the single-scattering limit. Formulas for the reflectance in this limit are referred as the single-scattering approximation. For a single scattering approximation Equation (9) can be rewritten as</p><disp-formula id="scirp.60622-formula1528"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x26.png"  xlink:type="simple"/></disp-formula><p>where, aerosol contribution ρ<sub>as</sub> is provided by:</p><disp-formula id="scirp.60622-formula1529"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x27.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x28.png" xlink:type="simple"/></inline-formula> is the aerosol scattering phase function for a scattering angle γ, ω<sub>a</sub> is the aerosol single scattering albedo, and r is Fresnel reflectance of the interface for a given incident angle.</p><p>Following the work of [<xref ref-type="bibr" rid="scirp.60622-ref13">13</xref>] , who have shown that for near-infrared channels the water leaving radiance coming out of clear open ocean waters can be assumed to be near zero because of strong water absorption, and if the effect of sun glitter is also assumed to be negligible, as in the case of OCM sensor which has provision for tilting mechanism to avoid sun glint [<xref ref-type="bibr" rid="scirp.60622-ref8">8</xref>] , then Equation (16) can be written as:</p><disp-formula id="scirp.60622-formula1530"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x29.png"  xlink:type="simple"/></disp-formula><p>Therefore, the effects of aerosols<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x30.png" xlink:type="simple"/></inline-formula>, in the imagery can be estimated at the two NIR bands of 740 and 856-nm from the sensor-measured radiances and the computed Rayleigh scattering reflectance. This quantity is then extrapolated and removed in the visible bands of OCM data. Let us assume that the path reflectance at two bands in the NIR at <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x31.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x32.png" xlink:type="simple"/></inline-formula>, where, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x33.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x34.png" xlink:type="simple"/></inline-formula>, can be estimated from the Equation</p><p>(18). Given that <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x35.png" xlink:type="simple"/></inline-formula> can be precisely calculated given the surface atmospheric pressure and therefore <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x36.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x37.png" xlink:type="simple"/></inline-formula> can be determined from Equation (18) at <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x38.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x39.png" xlink:type="simple"/></inline-formula>. This allows estimation of the parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x40.png" xlink:type="simple"/></inline-formula> as:</p><disp-formula id="scirp.60622-formula1531"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x41.png"  xlink:type="simple"/></disp-formula><p>This <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x42.png" xlink:type="simple"/></inline-formula> parameter is subsequently used to compute spectral variation, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x43.png" xlink:type="simple"/></inline-formula>, of aerosol properties in the shorter wavelengths. The value of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x44.png" xlink:type="simple"/></inline-formula> for the OCM band <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x45.png" xlink:type="simple"/></inline-formula> is calculated from<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x46.png" xlink:type="simple"/></inline-formula>, this yields<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x47.png" xlink:type="simple"/></inline-formula>, which when combined with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x48.png" xlink:type="simple"/></inline-formula>, provides the desired path radiance in λ<sub>i</sub> wavelength. Clearly, the key to this procedure is the estimation of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x49.png" xlink:type="simple"/></inline-formula> from<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x50.png" xlink:type="simple"/></inline-formula>. The aerosol reflectance, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x45.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x46.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x51.png" xlink:type="simple"/></inline-formula>, in the i<sup>th</sup> wavelength is calculated by using the following equation:</p><disp-formula id="scirp.60622-formula1532"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x52.png"  xlink:type="simple"/></disp-formula><p>Extrapolation of aerosol properties from <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x53.png" xlink:type="simple"/></inline-formula> to<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x54.png" xlink:type="simple"/></inline-formula>, which involves more than a factor of two in wavelengths has been studied by [<xref ref-type="bibr" rid="scirp.60622-ref11">11</xref>] . This study has provided detailed insight into the possible spectral behaviour of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x55.png" xlink:type="simple"/></inline-formula>. They have computed <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x56.png" xlink:type="simple"/></inline-formula> for several aerosol models. Three basic aerosol models namely; the tropospheric model with no oceanic contribution; the maritime model for which 99% of the particles have the tropospheric characteristics and 1% the Oceanic; and the Coastal model for which 99.5% of the particles have the tropospheric characteristics and 0.5% the oceanic. This model was introduced to represent the aerosol over the oceans nearer the coast (less Oceanic contribution). The properties of all the three aerosol models depend on the wavelength and relative humidity. Sample results from this study for<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x57.png" xlink:type="simple"/></inline-formula>, where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x58.png" xlink:type="simple"/></inline-formula> is taken to be 865 nm (OCM band 8), are redrawn from [<xref ref-type="bibr" rid="scirp.60622-ref11">11</xref>] in <xref ref-type="fig" rid="fig4">Figure 4</xref>. These results suggest that there is a strong variation of ε with aerosol model and relative humidity. The increase in the particle size (due to swelling) with increasing relative humidity clearly reduces the spectral variation of ε. They have shown that over the range 412 - 865 nm <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x59.png" xlink:type="simple"/></inline-formula> can be considered to be an exponential function of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x60.png" xlink:type="simple"/></inline-formula>, for the [<xref ref-type="bibr" rid="scirp.60622-ref18">18</xref>] aerosol models and [<xref ref-type="bibr" rid="scirp.60622-ref11">11</xref>] have used these results for the extension of CZCS atmospheric correction algorithm for use with SeaWiFS and MODIS. The best-fit equation can be written as:</p><disp-formula id="scirp.60622-formula1533"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x61.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.60622-formula1534"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x62.png"  xlink:type="simple"/></disp-formula><p>where, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x63.png" xlink:type="simple"/></inline-formula>is the Rayleigh and ozone corrected aerosol radiances computed from the satellite image data. This algorithm is mainly valid for open ocean (case 1 water) where the assumption of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x64.png" xlink:type="simple"/></inline-formula> for NIR bands holds well. However, this is not true in the case of turbid coastal waters in which the water leaving reflectance at NIR bands are often not negligible. Further the estimated water-leaving reflectance are converted to normalised</p><p>water-leaving reflectance, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x65.png" xlink:type="simple"/></inline-formula>, which is earlier defined as</p><disp-formula id="scirp.60622-formula1535"><label>(23)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x66.png"  xlink:type="simple"/></disp-formula><p>where, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x67.png" xlink:type="simple"/></inline-formula>is the atmospheric diffuse transmittance in the solar direction with the solar zenith angle of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x68.png" xlink:type="simple"/></inline-formula>.</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> ε(λ, 865) for nadir viewing with θ<sub>0</sub> = 60˚ for the Maritime, Coastal, and Tropospheric aerosol models. For each model, the relative humidity values are 50%, 80%, and 98% from the upper to the lower curves (after Gordon and Wang, 1994)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1470216x69.png"/></fig><p>The estimated value of normalized water-leaving reflectance values in different spectral bands can then be used to generate ocean-colour products such as chlorophyll concentration.</p></sec><sec id="s7_3"><title>7.3. Atmospheric Correction Output</title><p>Output of the atmospheric correction algorithm will be used for bio-optical variable generation and output will be the normalized water-leaving radiance/remote-sensing reflectance, R<sub>rs</sub>, in the OCM-2 spectral bands of 412, 443, 490, 510, 555, and 620-nm. Quality indices or flags, such as for negative water-leaving/remote sensing reflectance, clouds, sun-glint mask and turbid water, will also be produced.</p></sec></sec><sec id="s8"><title>8. Aerosol Optical Depth (AOD) over Oceans at 865-nm</title><p>Aerosols are responsible for a number of physical effects in the atmosphere, the most important being modifying the atmospheric radiation balance [<xref ref-type="bibr" rid="scirp.60622-ref19">19</xref>] by reflecting away the incoming solar radiation and retaining the outgoing terrestrial radiation (greenhouse effect). Aerosols also reduce visibility, cause air pollution and influence the cloud formation microphysics by acting as condensation nuclei. Satellites, with their capability for large area coverage and short-term repeativity, are the most ideal means for acquiring global information on aerosols [<xref ref-type="bibr" rid="scirp.60622-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.60622-ref3">3</xref>] . The currently orbiting ocean colour sensors like Oceansat-1 OCM, SeaWiFS, MODIS and MERIS though primarily meant for ocean colour remote sensing, have the additional capability for monitoring of global distribution of marine aerosols.</p><p>Space borne ocean-colour remote sensor detected radiance is heavily contaminated by solar radiation backscattered by the atmosphere air molecules &amp; aerosols. This radiance is called the atmospheric radiance. So for the detection of the oceanic constituents first step is therefore the removal of the atmospheric contribution from the sensor radiances. For this NIR channel (l &gt; 700 nm) is used as ocean surface acts as a dark background due to the high absorption by water. Therefore the sensor-detected radiances can be considered as the atmospheric radiance, which can be treated just the sum of the Rayleigh &amp; Aerosol path radiance produced by the scattering of light by air molecules &amp; aerosols. So according to [<xref ref-type="bibr" rid="scirp.60622-ref15">15</xref>] radiances detected by a space borne sensor at top of atmosphere (TOA) at wavelength l &gt; 700 nm can be split into</p><disp-formula id="scirp.60622-formula1536"><label>(24)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x70.png"  xlink:type="simple"/></disp-formula><p>where,</p><p>L<sub>t</sub> = Sensor detected radiance;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x71.png" xlink:type="simple"/></inline-formula>is Aerosol path radiance.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/6-1470216x72.png" xlink:type="simple"/></inline-formula>is Rayleigh path radiance.</p><p>where, F<sub>o</sub> is Extra terrestrial solar flux.; J<sub>v</sub> is Satellite viewing angle; w<sub>oa</sub> is Aerosol single scattering albedo; w<sub>or</sub> is Rayleigh single scattering albedo (~1.0); t<sub>a</sub> is Aerosol optical depth; t<sub>r</sub> is Rayleigh optical depth; P<sub>a</sub>(J) function related to Aerosol scattering Phase function; P<sub>r</sub>(J) function related to Rayleigh scattering Phase function. This P<sub>a</sub>(J)/P<sub>r</sub>(J) related to the Aerosol/Rayleigh scattering Phase function according to [<xref ref-type="bibr" rid="scirp.60622-ref15">15</xref>] is given as</p><disp-formula id="scirp.60622-formula1537"><label>(25)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x73.png"  xlink:type="simple"/></disp-formula><p>where,</p><p>R(J<sub>v</sub>) Frensel reflectance of the water surface along (J<sub>v</sub>) satellite viewing angle.</p><p>R(J<sub>s</sub>) Frensel reflectance of the water surface along (J<sub>s</sub>) solar zenith angle.</p><p>J<sup>&#177;</sup> represents the Forward (+)/Backward (−) scattering angle.</p><p>The scattering angles in direction to the sensor and in direction to the sensor via the air-sea interface is given as</p><disp-formula id="scirp.60622-formula1538"><label>(26)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x74.png"  xlink:type="simple"/></disp-formula><p>where, f is the azimuth difference between the sensor viewing and solar illumination direction.</p><p>The Rayleigh Phase function is computed by following equation</p><disp-formula id="scirp.60622-formula1539"><label>(27)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x75.png"  xlink:type="simple"/></disp-formula><p>The aerosol phase function can be approximated by the two-term Henyey-Greenstein phase function:</p><disp-formula id="scirp.60622-formula1540"><label>(28)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x76.png"  xlink:type="simple"/></disp-formula><p>with α = 0.985, g<sub>1</sub> = 0.8, g<sub>2</sub> = 0.5 for marine aerosols [<xref ref-type="bibr" rid="scirp.60622-ref15">15</xref>] .</p><p>For calculating the specular reflectance at the air-sea interface R(J<sub>v</sub>) and R(J<sub>s</sub>), we need the law of refraction (Snell’s law):</p><disp-formula id="scirp.60622-formula1541"><label>(29)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x77.png"  xlink:type="simple"/></disp-formula><p>where n is the refractive index of water, i the incident zenith angle J<sub>v</sub> or J<sub>s</sub> and j the zenith angle of the refracted beam under water and n the refractive index, which is for water 4/3, and the Fresnel law for unpolarized light which is:</p><disp-formula id="scirp.60622-formula1542"><label>(30)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x78.png"  xlink:type="simple"/></disp-formula><p>Assuming w<sub>oa</sub> &#187; 1.0 for marine aerosols, we can determine the aerosol optical depth (AOD) from the sensor detected radiance as</p><disp-formula id="scirp.60622-formula1543"><label>(31)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/6-1470216x79.png"  xlink:type="simple"/></disp-formula><p>This algorithm was also successfully used by [<xref ref-type="bibr" rid="scirp.60622-ref20">20</xref>] with Oceansat-1 OCM data for the estimation of AOD over the oceanic areas and now the same algorithm is used with Oceansat-2 OCM data. <xref ref-type="fig" rid="fig5">Figure 5</xref> shows the flow chart for the computation of the bio-optical and atmospheric parameters from the Oceansat-2 OCM data.</p><p>As OCM-2 has two day repetivity two day composite of the successive two day passes are sufficient to cover the entire north Indian Ocean (NIO). However, it needs one week passes to generate a single global coverage AOD product to cover the entire global oceans. As a example of weekly (8-day) composite and monthly composites covering the north Indian Ocean and global oceans are shown in <xref ref-type="fig" rid="fig6">Figure 6</xref>(a) and <xref ref-type="fig" rid="fig6">Figure 6</xref>(b) respectively.</p></sec><sec id="s9"><title>9. Conclusion</title><p>In this present study the methodology and the processing scheme for estimating the aerosol optical depth (AOD)</p><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Flow chart for the computation of the bio-geo-physical and atmospheric parameters from Oceansat-2 OCM data</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1470216x80.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Variability of the aerosol optical depth (AOD) (a) over the north Indian Ocean (weekly/8-day composite) and (b) over the global oceans (monthly composite)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/6-1470216x81.png"/></fig><p>at 865 nm from Oceansat-2 Ocean Colour Monitor (OCM) has been explained. From this study it is clear that OCM-2 derived AOD products are highly useful for studying the atmospheric turbidity and plays an important role in generating the accurate geophysical products. The LAC and GAC AOD products of OCM-2 show a very good spatio-temporal variability over the north Indian Ocean and over the global oceans respectively. Further these AOD products over the north Indian Ocean needs to be intercompared and validated with the existing contemporary ocean colour sensors.</p></sec><sec id="s10"><title>Acknowledgements</title><p>The authors would like to thank the ISRO’s Earth Observations (EO) missions for providing the consistent ocean colour data and SeaDAS development group at NASA GSFC for providing the standard SeaDASv6.4 and the recent version 7 for processing OCM-2 data. We Thank Dr. Prakash Chauhan, Space Applications Centre, Ahmedabad for providing valuable information on Indian Ocean Colour sensors and the retrieval strategy.</p></sec><sec id="s11"><title>Cite this paper</title><p>Udaya LakshmiSriperambudur,ChandralingamSonnati,Nagamani PullaiahgariVenkata,Preethi LathaThalathoti, (2015) Retrieval of Aerosol Optical Depth from Oceansat-2 OCM. Open Journal of Marine Science,05,443-454. doi: 10.4236/ojms.2015.54035</p></sec></body><back><ref-list><title>References</title><ref id="scirp.60622-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Kaufman, Y.J., Tanré, D. and Boucher, O. (2002) A Satellite View of Aerosols in the Climate System. Nature, 419, 215-223. http://dx.doi.org/10.1038/nature01091</mixed-citation></ref><ref id="scirp.60622-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Seinfeld, J.H. and Pandis, S.N. (2012) Atmospheric Chemistry and Physics: From Air Pollution to Climate Change. John Wiley &amp; Sons.</mixed-citation></ref><ref id="scirp.60622-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Kaufman, Y.J., Tanre, D., Gordon, H.R., Nakajima, T., Lenoble, J., Frouin, R., Grassl, H., Herman, B.M., King, M.D. and Teillet, P.M. (1997) Passive Remote Sensing of Tropospheric Aerosol and Atmospheric Correction for the Aerosol Effect. Journal of Geophysical Research, 102, 16815-16830. http://dx.doi.org/10.1029/97JD01496</mixed-citation></ref><ref id="scirp.60622-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">King, M.D., Kaufman, Y.J., Tanré, D. and Nakajima, T. (1999) Remote Sensing of Tropospheric Aerosols from Space: Past, Present, and Future.</mixed-citation></ref><ref id="scirp.60622-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Bhaskaran, S., Phillip, N., Rahman, A. and Mallick, J. (2011) Applications of Satellite Data for Aerosol Optical Depth (AOD) Retrievals and Validation with AERONET Data. Atmospheric and Climate Sciences, 1, 61-67.   
http://dx.doi.org/10.4236/acs.2011.12007</mixed-citation></ref><ref id="scirp.60622-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Prakash, C., and Nagamani, P.V. (2007) Algorithm Theoretical Basis Document for Geophysical Parameter Retrieval using Oceansat-2 OCM, OCEANSAT-II-UP-3/SAC-2007.</mixed-citation></ref><ref id="scirp.60622-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Eckstien, B.A. and Simpson, J.J. (1981) Cloud Screening Coastal Zone Colour Scanner Images Using Channel 5. International Journal of Remote Sensing, 12, 2359-2377.</mixed-citation></ref><ref id="scirp.60622-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Mohan, M. and Chauhan, P. (2001) Simulations for Optimal Payload Tilt to Avoid Sunglint in IRS-P4 Ocean Colour Monitor (OCM) Data around the Indian Subcontinent. International Journal of Remote Sensing, 22, 185-190.  
http://dx.doi.org/10.1080/014311601750038929</mixed-citation></ref><ref id="scirp.60622-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">McClain, C.R. and Yeh, E.-N. (1994) Sun Glint Flag Sensitivity Study. SeaWiFS Technical Report Series, Vol. 13, NASA Technical Memorandum 104566.</mixed-citation></ref><ref id="scirp.60622-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Coakley, J.A., Cess, R.D. and Yurevich, F.B. (1983) The Effect of Tropospheric Aerosols on the Earth’s Radiation Budget: A Parameterisation for Climate Models. Journal of Atmospheric Science, 40, 116-138.  
http://dx.doi.org/10.1175/1520-0469(1983)040&lt;0116:TEOTAO&gt;2.0.CO;2</mixed-citation></ref><ref id="scirp.60622-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Cox, C. and Munk, W. (1954) Measurement of the Roughness of the Sea Surface from Photographs of the Sun Glitter. Journal of Optical Society of America, 44, 838-850. http://dx.doi.org/10.1364/JOSA.44.000838</mixed-citation></ref><ref id="scirp.60622-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Gordon, H.R. and Wang, M. (1994) Retrieval of Water Leaving Radiance and Aerosol Optical Thickness over the Ocean with SeaWiFS: A Preliminary Algorithms. Applied Optics, 33, 443-452.  
http://dx.doi.org/10.1364/AO.33.000443</mixed-citation></ref><ref id="scirp.60622-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Chauhan, P., Mohan, M., Sarngi, R.K., Kumari, B., Nayak, S. and Matondkar, S.G.P. (2002) Surface Chlorophyll A Estimation Using IRS-P4 OCM Data in the Arabian Sea. International Journal of Remote Sensing, 23, 1663-1676. 
http://dx.doi.org/10.1080/01431160110075866</mixed-citation></ref><ref id="scirp.60622-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Gordon, H.R. and Clarke, G.L. (1981) Clear Water Radiance for Atmospheric Correction of CZCS Imagery. Applied Optics, 20, 4157-4180. http://dx.doi.org/10.1364/AO.20.004175</mixed-citation></ref><ref id="scirp.60622-ref15"><label>15</label><mixed-citation publication-type="book" xlink:type="simple">Gordon, H.R. and Morel, A. (1983) Remote Assessment of Ocean Color for Interpretation of Satellite Visible Imagery: A Review. In: Barber, R.T., Mooers, N.K., Bowman, M.J. and Zeitzschel, B., Eds., Lecture Notes on Coastal and Estuarine Studies, Springer-Verlag, New York, 114. http://dx.doi.org/10.1029/ln004</mixed-citation></ref><ref id="scirp.60622-ref16"><label>16</label><mixed-citation publication-type="book" xlink:type="simple">Doerffer, R. (1992) Imaging Spectroscopy for Detection of Chlorophyll and Suspended Matter. In: Toselli, F. and Bodechtel, J., Eds., Imaging Spectroscopy: Fundamentals and Prospective Applications, Kluwer Academic Publishers, Dordrecht, Boston and London, 215-257.</mixed-citation></ref><ref id="scirp.60622-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Deschamps, P.Y., Herman, M. and Tanre, D. (1983) Modeling of the Atmospheric Effects and Its Applications to the Remote Sensing of Ocean Color Remote Sensing. Applied Optics, 34, 2068-2080.  
http://dx.doi.org/10.1364/ao.22.003751</mixed-citation></ref><ref id="scirp.60622-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Gordon, H.R. (1997) Atmospheric Correction of Ocean-Colour Imagery in the Earth Observing System Era. Journal of Geophysical Research, 102, 17081-17106. http://dx.doi.org/10.1029/96JD02443</mixed-citation></ref><ref id="scirp.60622-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Shettle, E.P. and Fenn, R.W. (1979) Models for the Aerosols of the Lower Atmosphere and the Effects of Humidity Variations on Their Optical Properties. AFGL-TR-79-0214, 675, 94 p.</mixed-citation></ref><ref id="scirp.60622-ref20"><label>20</label><mixed-citation publication-type="book" xlink:type="simple">Joseph, J.H. (1984) Gaps in the Knowledge of Aerosols—Analysis and Remedies. In: Gerber, H.E. and Deepak, A., Eds., Aerosols and Their Climatic Effects, A. Deepak Publishing, Hampton, 275-278.</mixed-citation></ref></ref-list></back></article>