Effect of the Aerosol Type Selection for the Retrieval of Shortwave Ground Net Radiation: Case Study Using Landsat 8 Data

This paper discusses the aerosol radiative effects involved in the accuracy of shortwave net radiation, R n.sw , with sw ∈ (400–900) nm, retrieved by the Operational Land Imager (OLI), the new generation sensor of the Landsat mission. Net radiation is a key parameter for the energy exchange between the land and atmosphere; thus, R n.sw retrieval from space is under investigation by exploiting the increased spatial resolution of the visible and near-infrared OLI data. We adopted the latest version of the Second Simulation of a Satellite Signal in the Solar Spectrum (6SV) atmospheric radiative transfer model implemented in the atmospheric correction algorithm (OLI Atmospherically-Corrected Reflectance Imagery (OLI@CRI)) developed specifically for OLI data. The values of R n.sw were obtained by varying the microphysical properties of the aerosol during the OLI@CRI retrieval of both the OLI surface reflectance, ρ oli pxl , and the incoming solar irradiance at the surface. The analysis of the aerosol effects on the R n.sw was carried out on a spectrally-homogeneous desert area located in the southwestern Nile Delta. The results reveal that the R n.sw available for energy exchange between the land and atmosphere reduces the accuracy (NRMSE 14%) when the local aerosol microphysical properties are not considered during the processing of space data. Consequently, these findings suggest that the aerosol type should be considered for variables retrieved by satellite observations concerning the energy exchange in the natural ecosystems, such as Photosynthetically-Active Radiation (PAR). This will also improve the accuracy of land monitoring and of solar energy for power generation when space data are used.


Introduction
Natural and anthropogenic aerosols directly affect the energy exchange in the Earth-atmosphere coupled system by scattering and absorbing the solar radiation and indirectly by changing the lifetime and the microphysical properties of the clouds [1][2][3].There are still uncertainties regarding the radiative effects of the atmospheric aerosol, and many efforts are underway to overcome the resulting errors on the energy exchange in natural ecosystems [4].Concerning the direct effect, aerosols determine the direct and diffuse components of the radiation available at the surface for land-atmosphere energy exchange.A key parameter for energy exchange studies is the net radiation, R n , which is defined by the difference between the incoming and the outgoing irradiance, which in turn depend on surface parameters.The R n is the total of shortwave (R n.sw ) and longwave (R n.lw ) radiation.The R n.sw depends on the capability of the surface to reflect solar radiation (surface albedo), while the R n.lw depends on the temperature and emissivity of the surface [5].
Remote sensing provides the spatial distribution of surface net radiation for both land use management and climate modeling.
To date, remote sensing data have restricted net radiation analysis to low spatial resolutions, as reported in [6].On a global scale, in the 1985, the Earth Radiation Budget Experiment (ERBE), instruments on the Earth Radiation Budget Satellite (ERBS) [7] and then the Geostationary Earth Radiation Budget (GERB) [8] were designed with low spatial resolution specifically for studies on broadband net planetary energy.Later, the Clouds and the Earth's Radiant Energy System (CERES) was designed to follow the highly successful ERBE by increasing the spatial resolution to 20 km.The CERES surface radiation budget is retrieved by using the cloud and aerosol properties derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument [9,10], which is on-board the Earth Observing System (EOS) platforms.The POLarizationand Directionality of the Earth's Reflectance (POLDER) is another useful sensor for this kind of study at a low spatial resolution (about 6 km) [11].
The most recent generation of Earth Observation sensors with improvements in terms of spectral and spatial resolution provides a new opportunity for the analysis of energy exchange in the Earth-atmosphere coupled system.Thus, the temporal and spatial combination of the new generation of the European Space Agency (ESA) [12] and the National Aeronautics and Space Administration (NASA) [13] sensors introduces a unique Earth observation system for systematic coverage at a medium resolution with a high revisit frequency of up to five days.Masek et al. [14] reported that the latest Landsat mission of NASA meets the requirements defined by the Global Climate Observing System (GCOS) for studying and monitoring land and coastal ecosystems.Indeed, the Operational Land Imager (OLI) sensor on-board Landsat 8 with a spatial resolution of 30 m in the shortwave spectral domain (visible to middle-infrared) is able to detect modifications of ecosystems by anthropic or natural sources [13].
Consequently, concerning the monitoring of ecosystems or the study of the energy exchange at the surface in the shortwave domain, remote sensing applications require the accurate simulation of the radiative effects of the aerosol on the shortwave net radiation, R n.sw , available at the ground level for energy exchange between the surface and atmosphere [15].Mira et al. [5] highlight the role of the surface albedo for R n.sw retrieval using remote sensing data.Thus, the Bidirectional Reflectance Distribution Function (BRDF) of the surface is required as primary input for characterizing the anisotropy of the surface albedo depending on the zenith and the azimuth angles along illumination and viewing directions.The BRDF can be determined from data acquired by multi-view sensors, such as POLDER, to reconstruct the reflectance function directly from the measurements [16] or by coupling a measurement with kernel-driven models in the case of a sensor with a single view acquisition, such as MODIS [10].The operational BRDF and albedo nadir reflectance products from MODIS data can be downloaded at no cost from the EOS Data Gateway (http://eos.nasa.gov/imswelcome/)and are described in Schaaf et al. [10] and Wanner et al. [17].In all cases, R n.sw is strictly connected to the surface and the aerosol physical properties, as explained in Mielonen et al. [18].Both the surface albedo and aerosol influence the energy exchange between the biosphere and the atmosphere [19].
In order to retrieve the shortwave net radiation from space, the sensor radiance firstly needs to be converted to directional surface reflectance by the atmospheric correction processing [20,21].This processing is mandatory for the assessment of spectral indices retrieved by remote sensing, such as the Photosynthetically-Active Radiation (PAR) and the Leaf Area Index (LAI) driven by the energy exchange and aimed at ecosystem productivity [22,23].The accurate directional surface reflectance from remote sensing data is successfully achieved by physically-based atmospheric correction algorithms developed using radiative transfer models, such as the open-source 6SV1.1 [24,25], the latest version of the Second Simulation of a Satellite Signal in the Solar Spectrum (6S [26]) and MODTRAN (Moderate Resolution Atmospheric Transmittance and Radiance) [27].These models provide an accurate simulation of the solar radiation through the Earth-atmosphere coupled system by solving the radiative transfer equation [26,28].If a ground station of the the AERONET worldwide network [29] is within the remote sensing image, the atmospheric correction can exploit the inverse products of the station by representing the microphysical properties of the aerosol [30,31], which improve the accuracy of the simulated radiative effects due to the aerosol [15,32,33].The surface reflectance can thus be completely and accurately derived.
The accuracy of the surface reflectance retrieved from space very much depends on the correct modeling of the atmospheric scattering and absorption due to the aerosol [21].The principal aerosol parameter of the radiative transfer model is the optical thickness of the aerosol at 550 nm, τ 550 , with robust and well-known retrieval algorithms for all of the above-mentioned sensors [24,[34][35][36][37].Other crucial parameters to simulate the radiative effects of the aerosol are the microphysical properties of the aerosol: the size distribution, dV(r)/dln(r), and the real, Re λ and the imaginary, Ri λ , parts of the refractive index, integrated along the atmospheric column [38,39].Unlike τ 550 , there is still no robust algorithm for the retrieval of the microphysical properties of the aerosol from space.Levy et al. [40] presents the global aerosol optical properties applied to MODIS for aerosol retrieval.Specifically, for the reflectance retrieval of dark surfaces, new aerosol models had to be introduced to improve the accuracy of the surface reflectance retrieved by atmospheric correction processing, as reported in Ahmad et al. [41].Furthermore, in the case of strongly-absorbing aerosols with high aerosol optical thickness over a dark surface, also the aerosol profile needs to be considered for the high accuracy of the atmospheric correction products [42] by using data available from the several infrastructures [43][44][45][46][47].
We investigated the influence of the aerosol types with the default exponential vertical aerosol profile of the 6SVmodel, on the shortwave net radiation from OLI data by considering the OLI image atmospherically corrected with aerosol dust-like component as the benchmark, after validation with the reflectance MODIS product, ρ modis .The results highlight that an accurate study and monitoring of the ecosystem is provided by new generation Earth Observation data if aerosol properties are taken into account.Thus, the aerosol radiative effects in terms of light absorption and scattering [28] could also be considered to analyze the uncertainty of the satellite-derived surface albedo products [19].Our case study shows how a refined aerosol characterization reduces the uncertainties of the surface albedo and the shortwave net radiation from the new generation of the satellite remote sensing data.

Study Area
The study area was a very high reflective desert located in the southwestern Nile Delta, 70 km from Cairo (Figure 1), where an AERONET station was in operation.This desert area has a rectangular shape of 7.85 km × 5.5 km, with the centroid in 30°40 6.6 E and 29°53 57 N with an average height of 137 m.The target is of particular interest due to the Calibration and Validation (CAL/VAL) of the remote sensing products, such as surface reflectance.The geographical position also ensure the maximum illumination of the site, which is useful for the analysis of energy exchange.The closeness of the study area to an urban environment could be advantageous for the specific analysis of the radiative effects of anthropic aerosol on energy exchange.

OLI and MODIS Data
The remote sensing data employed were available on EarthExplorer, the long-term archive of USGS (http://earthexplorer.usgs.gov/).Two sun-synchronous sensors was considered (Table 1): the Moderate Resolution Imaging Spectroradiometer (MODIS) on board EOS AM-1 (Terra) and the Optical Land Imager (OLI) on board the Landsat 8 satellite.These sensors differ in terms of the revisit time (1 day for MODIS vs. 16 days for OLI) and the spatial resolution of 250 m (Bands 1-2), 500 m (Bands 3-7) for MODIS and 30 m for the OLI (Band 1-7), respectively.6 December 2015 was selected to ensure the concurrent availability of OLI at 8:30 GMT and MODIS products derived from the image acquired at 9:35 GMT.The meteorological conditions when Landsat passed over the study area showed low cloud cover with a temperature of 13°, wind speed of 4.6 m/s and direction ENE.The OLI data were collected in Path 177 Row 039 at the level of processing L1T (terrain corrected) and were in the reference system WGS84 projection UTM 36.The MOD09GA and MCD43A1 products of the MODIS sensor were downloaded considering the tile h20 and v06.MOD09GA is the surface reflectance in a daily gridded L2G product with spatial resolution of 500 m and sinusoidal projection (LPDAAC).The MCD43A1 product contains the BRDF parameters to derive the albedo by a three-dimensional system composed of isotropic, volumetric and surface scattering components (iso, vol, geo) for each MODIS band.The BRDF parameters are derived over a 16-day period describing the surface anisotropy of each pixel with a spatial resolution of 500 m and sinusoidal projection.

AERONET Data
The AERONET network [29] provides direct and inverse products (http://aeronet.gsfc.nasa.gov/)derived from the radiative measurements performed by the automatic Cimel Sun sky-radiometer [30].These products are the integrated columnar values of atmospheric parameters, which are adopted to simulate the atmospheric radiative effects on the solar radiation propagating through the atmosphere by radiative transfer models, such as 6SV [25].In this study, the direct and inverse products from the AERONET stations in Cairo and El Farafra 400 km from the test site were used as the closest atmospheric state measurements to the considered OLI image in space and time.The two AERONET stations are taken into consideration for aerosol characterization for the distance (El Cairo) and the environment (El Farafra).El Cairo is the closest station to the site and located in a urban environment, whereas El Farafra station is far from the study area, but in the same environment (i.e., desert).
Table 2 reports the direct products retrieved from the Sun-radiance measurements (at Level 1.5).The aerosol optical thickness at 550 nm, τ 550 , is the principal aerosol parameter investigated from space [24,[34][35][36]48,49]. τ 550 was derived by applying exponential interpolation between the two nearest aerosol optical thicknesses τ 440 at 440 nm and τ 870 at 870 nm available from the two AERONET stations at the OLI overpass time.Table 2 shows low values of τ 550 in the two AERONET stations, revealing that aerosol events did not occur during the Landsat overpass.wv and o3 represent the columnar content of water vapor and ozone, respectively.Both atmospheric parameters principally influence the radiative transfer within the absorption bands of the water vapor and ozone, which are outside the OLI channels considered in this study.
The plots in Figure 2 show the inverse products of the two AERONET stations [32] (at Level 1.5) and the microphysical properties of the basic components included in the 6SV model (dust-like, water-soluble, oceanic and soot) [25].On the left, the plot shows the volumetric size distribution of the aerosol types considered in this study, dV(r)/dln(r).As expected for a desert site, a coarse mode-dominated distribution was measured in El Farafra.The trimodal distribution in Cairo is combined by the nuclei and the accumulation modes (fine particles) and the coarse mode, the usual pattern of size distribution in urban sites, as explained on in Kokhanovsky et al. [49].
The centered plot in Figure 2 shows the real part of the refractive index, Re, for each aerosol type considered.The real parts for the dust-like and water-soluble components are similar as in the case of the inverse products of the Cairo and El Farafra stations (Re = 1.6).The Re values of soot and oceanic are furthest away from the Cairo and El Farafra products and probably less representative of the local aerosol.
Concerning the imaginary parts of the refractive index, Ri, reported on the right in Figure 2, all of the aerosol types presented the Ri as being low sensitive to wavelength as in the case of the Re, with the exception of the Ri of the oceanic component.The oceanic component has low Ri, leading to negligible absorption of the solar radiation at the visible spectral domain.2) and the microphysical properties of the aerosol basic components reported in [25].On the left, the volumetric size distribution; in the center, the real part of the refractive index; on the right, the imaginary part of the refractive index.

Methodology
The radiative effects of the aerosol on the R n.sw retrieved from space were investigated as outlined in Figure 3.The core component is the surface albedo achieved pixel-by-pixel from the OLI image.The surface albedo, α oli pxl , was assessed following Shuai et al. [50] combining Landsat surface reflectance, ρ oli , with MODIS surface anisotropy products (BRDF) and black and white albedo (α moli bs , α moli ws ) under the assumption of a spectral homogeneous desert area.The surface reflectance was retrieved by a physically-based algorithm developed for the atmospheric correction of OLI images, OLI Atmospherically-Corrected Reflectance Imagery (OLI@CRI) using the 6SV model.The algorithm was implemented to work with the microphysical properties of the aerosol available as AERONET inverse products or with the microphysical properties of four aerosol basic component included in the 6SV model (dust-like, water-soluble, oceanic, soot) [25] considering an exponential vertical distribution of the aerosol with a scale height of 2 (default aerosol profile of 6SV model).
Using this algorithm, the OLI image was atmospherically corrected with different aerosol types.The OLI corrected images were compared to the concurrent MODIS surface reflectance product in order to validate the OLI reflectance, as suggested in [51].The albedo was achieved by combining OLI surface reflectance with MODIS surface anisotropy information [50].Finally, the R n.sw was obtained by the equation expressed in Mira et al. [5] for all of the aerosol types within the OLI channels in the spectral domain 400-900 nm.

OLI@CRI Atmospheric Correction Algorithm
The OLI@CRI algorithm was developed for the atmospheric correction processing of the OLI images specifically over land following the procedure described in [38,48] adapting the vector version of the Second Simulation of a Satellite Signal in the Solar Spectrum (6SV) radiative transfer model.This model is the improved version of the 6S used for the new generation of atmospheric correction algorithms, which was successfully adopted by [52].The algorithm was implemented with the open-source GNU data language (GDL) [53] and works with the OLI image available in L1T format.As with the previous algorithms [38,39,48], OLI@CRI can be applied by considering the standard aerosol basic components (dust-like, water-soluble, oceanic and soot [25]) provided by default by the 6SV model or the microphysical properties of the aerosol provided by an AERONET station as inverse products.The OLI@CRI algorithm overcomes the limits of software such as Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS) [14], where a continental aerosol model is assumed for atmospheric correction processing and an over-correction of the surface reflectance is expected.LEDAPS therefore restricts the potential of the remote sensing data by excluding the microphysical properties of the aerosol.Indeed, the aerosol optical properties required for accurate radiative modeling during the atmospheric correction are derived from the observed microphysical properties by the Mie theory solved in the 6SV model [25] (p.108).
As usual, the ground reflectance at the λ wavelength with solar illumination (θ s , φ s ) and sensor viewing (θ v , φ v ) directions is expressed by [21,24,26]: where ρ TOA = πL v /µ s E s is the sensor reflectance expressed by the TOA (Top Of Atmosphere) solar irradiance, E s .The µ v,s = cos(θ v,s ) are the cosine of the view (v) and Sun (s) zenith angles to define the illumination and viewing directions.The θ v,s is the view (v) and Sun (s) azimuth angles.t g is the gas transmittance; ρ atm is the atmospheric reflectance; S is the atmospheric spherical albedo [28]; t S = t ↑ + t ↓ represents the total transmittance obtained from the downward and upward directions; and finally, t ↑,↓ = e −τ(λ)/µ v,s + t ↑,↓ di f is the transmittance of up and down directions composed of direct and diffuse component depending on the total optical thickness due to the aerosol and molecules (τ).The 6SV model [24] is implemented in OLI@CRI to simulate all of the radiative quantities of Equation ( 1) in the 400-900-nm spectral domain with a sampling of 2.5 nm.The OLI@CRI algorithm computes the convolution of the radiative quantities on the spectral response of the OLI sensor reported in Barsi et al. [54], following the procedure described in [48].
The surface reflectance, ρ oli pxl (λ, θ s , φ s , θ v , φ v ), is retrieved by removing the effect of the neighboring pixels into the pixel viewed from the sensor (adjacency effect) by the empirical formula [21,24,26]: where the < ρ gnd > is the mean of the pixel adjacent to the viewing pixel covering the entire swath of the OLI sensor.

Albedo
The albedo, α oli (λ, θ s ), was obtained following Shuai et al. [50], where the surface anisotropy is solved at OLI spatial resolution starting from the MODIS spatial resolution.The albedo was considered in the OLI channels in the visible spectral domain, which is more sensitive to the radiative impact of the aerosol.The homogeneous area was chosen at MODIS and OLI spatial resolutions to meet the hypothesis in Shuai et al. [50].The generation of the Landsat albedo [50] is based on the kernel-driven reciprocal RossThick-LiSparseBRDF model used in the MODIS algorithm for albedo retrieval [10,17,25].
where K vol,λ (θ v , θ s , φ) is the volume-scattering Ross-Thick kernel and K geo,λ (θ v , θ s , φ) is the surface-scattering LiSparse kernel.f iso,λ ; f vol,λ and f geo,λ are the parameters of the three-dimensional MCD43A1 product available for each MODIS band.These parameters were reprojected in the same reference system (WGS84 projection UTM fuse 36) of the OLI image.Based on Equation (3), both the albedo due to the direct solar beam without diffuse component, black-sky albedo α modis bs (θ, λ) and the albedo due to the diffuse component of the solar beam without the direct component, white-sky albedo α modis ws (λ), are generated by the following expression [10]: where f iso,λ ; f vol,λ and f geo,λ are explained in Equation (3) and θ s is the solar zenith angle at Landsat acquisition time as explained in [50].The g ik terms in Equation ( 4) are reported in Table 3. Shuai et al. [50] introduced the ratio between the albedo and reflectance, with the geometrical acquisition of Landsat, resampled to MODIS spatial resolution for albedo retrieval at Landsat resolution: where ρ moli (λ, θ s , θ v , φ) is the OLI surface reflectance obtained from Equation (2) resampled to MODIS spatial resolution.Under the hypothesis of a homogeneous area, verified by the desert area chosen for this study, Shuai, Y. et al. [50] affirmed that the ratio at MODIS spatial resolution is similar to the ratio at Landsat spatial resolution, a moli bs (λ, θ s , θ v , φ) ∼ a bs (λ, θ s , θ v , φ) and a moli ws (λ, θ s , θ v , φ) ∼ a ws (λ, θ s , θ v , φ).Thus, white-sky and black-sky albedo at Landsat resolution are: Finally, the albedo at Landsat spatial resolution for each pixel, α oli pxl (θ s , λ), is retrieved following [15]: where α bs (θ, λ) is the black-sky albedo and α ws (λ) is the white-sky albedo.The Snet 6SV (θ s , τ 550 , aer.phys λ ) is the ratio of the diffuse albedo on the total component of solar radiation simulated by the 6SV model at Landsat acquisition time (solar zenith angle θ s ), with the aerosol optical thickness at 550 nm of the Cairo AERONET station and considering the microphysical properties of the aerosol, aer.phys λ , as downloaded from AERONET stations in Cairo and El Farafra and referring to the aerosol basic components (dust-like, water-soluble, oceanic and soot [25]).

Shortwave Net Radiation, R n.sw
The instantaneous shortwave net radiation at ground level in the spectral domain 400-900 nm, R n.sw , is achieved by the equation introduced in Mira et al. [5] limited to the visible and near-infrared components where the radiative effects of the aerosol are stronger than any of those in other spectral domains: R n.sw,λ = (1 where R ↓ sw is the at ground total component of the downward solar irradiance, including direct and diffuse components; α oli (θ s , λ) is the surface albedo at λ wavelength and illuminated at θ s angle.The R ↓ sw can be measured by using specific instruments or simulated by using an atmospheric radiative transfer model, such as 6SV.As reported in Shuai et al. [50], the aerosol is a source of errors for R ↓ sw assessment.The α oli (θ s , λ) can be very accurately and reliably derived from the surface reflectance if the microphysical properties of the aerosol are taken into account during the atmospheric correction of the remote sensing data [39,55].

OLI Surface Reflectance, ρ oli pxl
The OLI@CRI algorithm was applied to the OLI image acquired over the desert area under six atmospheric conditions defined by the AERONET inverse products available from the Cairo and El Farafra stations and by the size distribution and refractive index of the four aerosol basic components included in 6SV (dust-like, water-soluble, oceanic and soot [25]).
The τ 550 , the water vapor and the ozone columnar contents were fixed to the direct products of the AERONET station in Cairo since this was the closest to the study area.These atmospheric parameters were changed to the direct products of the AERONET station only when the OLI@CRI algorithm was applied to the OLI image with El Farafra atmospheric conditions.

Validation with MODIS Reflectance Product
The six atmospherically-corrected OLI images were resampled at MODIS spatial resolution to validate the OLI surface reflectance with the MODIS surface reflectance product (MOD09GA).The OLI resampling was performed in the selected area (Desert B type in Figure 4), which verified the spectral homogeneity conditions for the generation of Landsat albedo [50] and the condition of lower cloud cover with maximum insolation to evaluate the radiative impact of the aerosol on R n.sw .The homogeneous area was selected by the K-means unsupervised classification, which identified a spectral cluster of pixels depending on the land cover classes.The Desert B type (Figure 4) was selected as the highest reflective class in the image.
Figure 5 shows the six resampled OLI surface reflectances, ρ moli (λ, θ s , θ v , φ), and the MODIS reflectance products of the study area.The MODIS reflectance was considered as the benchmark to identify the most suitable aerosol type to validate the OLI reflectance.The nearest OLI reflectance to the MODIS product is obtained in the case of dust-like microphysical properties, while the case of Cairo inverse products is the most distant from the MODIS reflectance.This effect is similar to the aerosol condition described in Bas et al. [39], where the role of a tri-modal size distribution on the accuracy of retrieved surface reflectance from space needs to be investigated.Reflectance Imagery (OLI@CRI) algorithm considering the aerosol microphysical properties as the inverse products of Cairo and El Farafra AERONET stations and as the four basic component included in the 6SVmodel (dust-like, water-soluble, oceanic and soot [25]).
Table 4 shows the Normalized Root Mean Square Error (NRMSE) of the OLI surface reflectance of Figure 5 considering the MODIS surface reflectance of the study area as the benchmark.The NRMSE confirmed that the OLI surface reflectance nearest to the MODIS one is obtained with the dust-like In fact, the NRMSE for dust-like aerosol (NRMSE-dl) is lower than the others.The NRMSE for oceanic (NRMSE-oc) close to the NRMSE-dl confirms the crucial role of the size distribution of the radiative impact over a bright surface.The NRMSE for water-soluble (NRMSE-ws) is less than for soot (NRMSE-so).Table 4. Normalized Root Mean Square Error (NRMSE) of the six OLI surface reflectance obtained by OLI@CRI with the available microphysical properties of the aerosol: Cairo, El Farafra and the four basic components (dl for dust-like; ws for water-soluble; oc for oceanic; so for soot) with respect to the MOD09GA MODIS product on the study area.

OLI Albedo
The surface albedo, α oli pxl (θ s , λ), was retrieved using the ρ oli (λ, θ s , θ v , φ) and the MODIS products, which are the BRDF parameters, for each one of the six atmospherically-corrected OLI images.
Table 5 reports the surface albedo averaged by the pixels over the study area, for the OLI channels overlapping the MODIS ones.Table 5 showed that the second, third and fourth channels (from the first to third column) in the El Cairo station, water-soluble and oceanic are equal to the lowest values, while soot has the highest values.El Farafra has similar values to soot in the forth channel.The fifth band has the highest variability of values; in this channel, El Farafra has the highest albedo.These results highlight that dust-like is different in the spectral range 482-654 nm, while at 864 nm, it has values comparable to soot.Moreover, the results obtained by AERONET microphysical properties supply different results, confirming that the different microphysical conditions at a long distance provide a difference in albedo values.The microphysical properties of the soot clearly led to a different capability of the surface to reflect the solar radiation with regard to the capability obtained with the other aerosol types.
Table 5. Surface albedo averaged by the pixels over the study area at Landsat resolution for the second to the fifth OLI channels.The values refer to the surface albedo obtained from the OLI image atmospherically corrected with AERONET inverse products and the microphysical properties of the four basic components of the aerosol included in 6SV (dust-like, water-soluble, oceanic, soot [25]).R n.sw was solved pixel-by-pixel by Equation (10) after the retrieval of the Snet 6SV (θ s , τ 550 , aer.phys λ ).Table 6 reports the Snet 6SV (θ s , τ 550 , aer.phys λ ) over the study area for the second to the fifth OLI channels, simulated by the 6SV model for the aerosol microphysical properties as the inverse products of the Cairo and El Farafra AERONET stations and as the basic components of the aerosol included in 6SV (dust-like, water-soluble, oceanic, soot [25]).The lowest Snet 6SV (θ s , τ 550 , aer.phys λ ) is for the microphysical properties of the soot, whereas the highest values are for the El Farafra inverse products.

Aerosol
Table 7 shows R n.sw averaged by the pixels over the study area at Landsat resolution for the OLI channel overlapping the MODIS channels (as in Table 6).R n.sw was achieved for the six images from the atmospheric correction of the OLI image using the six aerosol types defined by the microphysical properties considered in this study: the AERONET inverse products available from the Cairo and El Farafra stations and the four basic components of the aerosol included in 6SV (dust-like, water-soluble, oceanic, soot [25]).Table 6.Snet 6SV (θ s , τ 550 , aer.phys λ ) over the study area for the second to the fifth OLI channel, simulated by the 6SV model for the aerosol microphysical properties as the inverse products of the Cairo and El Farafra AERONET stations and as the basic components of the aerosol included in 6SV (dust-like, water-soluble, oceanic, soot [25]).

Aerosol Types
Snet The accuracy of the OLI shortwave net radiation, R n.sw , and surface reflectance, ρ oli pxl , products was evaluated by the NRMSE of the products obtained for each aerosol type condition with respect to the dust-like basic component values validated with the MODIS reflectance product.The ρ oli was retrieved from the ρ oli pxl averaged by the pixels over the study area.Table 8 reports both NRMSEs achieved using all of the OLI channels (the second to the fifth) overlapping the MODIS channels.

Discussion and Conclusions
This paper discusses the aerosol impacts on the net radiation in the shortwave, R n.sw , from space applying the procedure for R sw.n retrieval by Equation (10) [5] over the study area, which satisfied the condition of maximum illumination and minimum cloud cover on a bright surface, as well as the requirement of spectral homogeneity (Figure 4) for Landsat albedo generation [50].The procedure for the albedo retrieval [50] (Section 3.2) is based on Equation ( 9) and involved concurrently MODIS products (BRDF, α ws and α bs in Section 3.2) and OLI surface reflectance (ρ oli pxl in Section 3.1).The OLI surface reflectance was retrieved by applying the OLI@CRI algorithm described Section 4.1.The atmospheric correction of the OLI image was performed using the inverse products of the two AERONET stations and the microphysical properties of the four aerosol basic components included in 6SV (dust-like, water-soluble, oceanic and soot) [25].The NRMSE of the six OLI surface reflectances obtained with the six aerosol conditions was achieved with respect to the MODIS image for the overlapping MODIS and OLI channels as reported in Table 4.The surface reflectance retrieved by considering the microphysical properties of the dust-like case produced lower NRMSE with respect to the MODIS surface reflectance (Table 4).Consequently, the influence of the aerosol conditions on the accuracy of ρ oli pxl in the study area was performed considering ρ oli pxl in the dust-like case as the validated surface reflectance.
The results highlight that the radiative impact of the aerosol plays a crucial role in the accuracy of the surface reflectance, ρ oli pxl , and the shortwave net radiation, R n.sw , achieved through OLI data.The use of these remote sensing products needs to take into account the aerosol type considered during the data processing.Indeed, the microphysical properties are the main issue for aerosol characterizations with respect to their radiative impact on remote sensing data, specifically in the atmospheric window.Table 2 shows low values of τ 550 in the two AERONET stations, revealing that aerosol events did not occur during the Landsat overpass.Hence, a low variability in the radiative impact due to different aerosol types is expected, whereas a great influence of wv and o3 is not expected, because their respective absorption bands are outside the OLI channels considered in this study.
The NRMSE of the OLI surface reflectance over the study area (NRSME.ρoli in Table 8 in Section 4.3) shows that the accuracy of the bright surface reflectance depends principally on the volumetric size distribution, dV(r)/dln(r).The lower NRSME.ρoli in Table 8 is in the case of the coarse components, oceanic and El Farafra aerosol in Figure 2. Oceanic is the aerosol basic component composed of the coarse mode as the dust-like and the bi-modal size distribution of the aerosol measured in El Farafra consist mainly of the coarse mode.The other NRSME.ρoli in Table 8 increases in the case of fine aerosol, water-soluble and soot components.The maximum NRSME.ρoli is achieved when a tri-modal size distribution characterizes the aerosol Cairo in Figure 2).As expected, the NRMSE of the OLI reflectance obtained with the AERONET products cannot be similar, since the stations are located in different environments (urban and desert) with their respective microphysical properties of the aerosol (cf. Figure 2) Concerning the absorption of solar radiation, the imaginary part of the refractive index led to a separation between the solar radiation absorption by the aerosol type with low NRSME.ρ oli (oceanic component and El Farafra).The oceanic component is non-absorbing aerosol, as confirmed by the very low imaginary part of the refractive index (right plot in Figure 2).On the other hand, the dust-like component and El Farafra aerosol are absorbing aerosol with a spectrally flat imaginary part of the refractive index in the visible spectral domain (Ri = 0.008 in Vermote et al. [25,28]).Thus, the refractive index seems to be the most negligible aerosol property for ρ oli retrieval in the case of a bright surface.
Nevertheless, the absorption of solar radiation plays a significant role in the retrieval of remotely-sensed products connected with land-atmosphere energy exchange, such as shortwave net radiation.The higher NRMSE of R n.sw (NRSME.R n.sw in Table 8 in Section 4.3) is for the soot and oceanic components.Their imaginary part of the refractive index differs significantly from the R i of the dust-like, the basic component selected as a reference from the reflectance validation; whereas for radiative remote sensing products, dV(r)/dln(r) seems to be negligible.The lower NRSME.R n.sw in Table 8 for R n.sw retrieved by the water-soluble component shows that a fine-mode aerosol (water-soluble) and a coarse-mode aerosol (dust-like) with similar imaginary parts of the refractive index (right plot of Figure 2) have the same radiative impact on the satellite-derived products.In the case study, a rough selection of the aerosol type can produce a maximum percent error of 14% (Table 8) in the retrieval of the R sw.n for low values of τ 550 products (Table 2).
The radiative impact of the aerosol on the whole procedure is well described by the size distribution, and the real and imaginary parts of the refractive index become a requirement for the retrieval of land-atmosphere energy exchange, especially if the visible spectral domain is involved.Thus, the aerosol microphysical properties improve the accuracy of the spectral indices, such as Photosynthetically-Active Radiation (PAR) and the Leaf Area Index (LAI), driven by the energy exchange and aimed at the productivity of land and coastal ecosystems.Also worth mentioning is the increasing in the accuracy of the total net radiation, R n , which is the key driving force for evapotranspiration coupled with the soil heat flux [56], if the longwave net radiation, R n.lw , is combined with the R n.sw achieved by considering the aerosol types.
We plan to complete the analysis of the radiative impact of aerosols on satellite-derived products by applying this procedure in other aerosol loading conditions taking into account the vertical distribution of the aerosol.This procedure is also the first step in involving the new generation of Earth Observation sensors for the retrieval of shortwave net radiation at a higher spatial resolution, although further research is needed to overcome the spectral homogeneity conditions.

Figure 1 .
Figure 1.The rectangle is the study area located in the southwestern Nile Delta.The geographical coordinates in the reference system WGS84 are shown along the border of the map.The AERONETstations are shown by red triangles.

Figure 2 .
Figure 2. The microphysical properties of the aerosol available from AERONET stations in Cairo and El Farafra measured close in time to the Landsat overpass (Table2) and the microphysical properties of the aerosol basic components reported in[25].On the left, the volumetric size distribution; in the center, the real part of the refractive index; on the right, the imaginary part of the refractive index.

Figure 3 .
Figure 3. Processing chain for retrieval of the R n.sw at OLI spatial resolution.

Figure 4 .
Figure 4. Unsupervised classification of the study area with main land cover clusters.

Figure 5 .
Figure 5.The OLI surface reflectance retrieved by applying of the OLI Atmospherically-Corrected Reflectance Imagery (OLI@CRI) algorithm considering the aerosol microphysical properties as the inverse products of Cairo and El Farafra AERONET stations and as the four basic component included in the 6SVmodel (dust-like, water-soluble, oceanic and soot[25]).

Table 1 .
Spectral and spatial characteristics of MODIS and OLI sensors.The MODIS Bands 1 and 2 have a

Table 2 .
Direct products of the AERONET stations in Cairo and El Farafra closest in time to the OLI image considered in this study(6 December 2015).

Table 7 .
R n.sw averaged by pixels at Landsat spatial resolution over the study area.The R n.sw is reported for the OLI channels (the second to the fifth) overlapping the MODIS channels.

Table 8 .
NRMSE of the shortwave net radiation, R n.sw , and the surface reflectance, ρ oli , obtained for each aerosol condition with respect to the values in the dust-like reference case over the study area.The NRMSEs were achieved considering the second to the fifth OLI channels.