Next Article in Journal
Optimal Interpolation of Precipitable Water Using Low Earth Orbit and Numerical Weather Prediction Data
Next Article in Special Issue
Modeling Solar Radiation in the Forest Using Remote Sensing Data: A Review of Approaches and Opportunities
Previous Article in Journal
A Sliding Window-Based Joint Sparse Representation (SWJSR) Method for Hyperspectral Anomaly Detection
Previous Article in Special Issue
Validation of the SARAH-E Satellite-Based Surface Solar Radiation Estimates over India
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High Turbidity Solis Clear Sky Model: Development and Validation

Department F.-A. Forel for Environmental and Aquatic Sciences, Institute for Environmental Sciences, University of Geneva, 1205 Genève, Switzerland
Remote Sens. 2018, 10(3), 435; https://doi.org/10.3390/rs10030435
Submission received: 8 January 2018 / Revised: 15 February 2018 / Accepted: 8 March 2018 / Published: 10 March 2018
(This article belongs to the Special Issue Solar Radiation, Modelling and Remote Sensing)

Abstract

:
The Solis clear sky model is a spectral scheme based on radiative transfer calculations and the Lambert–Beer relation. Its broadband version is a simplified fast analytical version; it is limited to broadband aerosol optical depths lower than 0.45, which is a weakness when applied in countries with very high turbidity such as China or India. In order to extend the use of the original simplified version of the model for high turbidity values, we developed a new version of the broadband Solis model based on radiative transfer calculations, valid for turbidity values up to 7, for the three components, global, beam, and diffuse, and for the four aerosol types defined by Shettle and Fenn. A validation of low turbidity data acquired in Geneva shows slightly better results than the previous version. On data acquired at sites presenting higher turbidity data, the bias stays within ±4% for the beam and the global irradiances, and the standard deviation around 5% for clean and stable condition data and around 12% for questionable data and variable sky conditions.

Graphical Abstract

1. Introduction

Anthropogenic activities have become an important factor in climate change and continuous monitoring of the solar irradiance reaching the ground is essential to understand the impact of such changes on the environment. Unfortunately, the density of the ground measurement network is insufficient, especially on continents like Africa or Asia. To circumvent this lack of measured data, meteorological satellites are of great help and models converting the satellite images into the different radiation components are becoming increasingly useful.
The first step in evaluating the solar irradiance from satellite images is a good knowledge of the highest possible radiation transmitted by a cloudless atmosphere. This can be done with the help of clear sky model taking into account the aerosol and water vapor content of the atmosphere. Multiple clear sky models can be found in the literature, but none of them is able to evaluate solar radiation for high turbidities resulting from heating and transportation activities and encountered in countries like India or China. For example, the Gueymard CPCR2 model [1] is limited to β turbidity values lower than 0.4 (which correspond to an aerosol optical depth at 700 nm aod700 of 0.64 for rural aerosol type, Angström exponent α = 1.3), Gueymard REST2 [2] is limited to β = 1 (aod700 = 1.6, rural aerosol), Bird’s model [3] is defined for visibility values up to 23 km (which corresponds to an aod700 = 0.27), and the ESRA clear sky scheme [4,5,6] was developed for Linke turbidity values TL [7] not exceeding a value of 7 (i.e., an aerosol optical depth aod700 = 0.44 for a 2 cm water vapor column w).
The spectral Solis clear sky scheme was first developed within the Mesor European program, whose subject was the management and exploitation of solar resource knowledge. Due to the spatial and temporal range of the satellite images, the on line process to evaluate irradiance maps implies fast analytical models. In [8], Ineichen published a broadband fast analytical version of the Solis model for rural aerosol type, and in 2010 a version in the form of an Excel tool for the four types of aerosols was defined by Shettle and Fenn [9] and Shettle [10]. These versions were limited to broadband aerosol optical depth values aod700 lower than 0.45; they diverge for higher turbidity, as shown by Zhang [11]. This is illustrated in Figure 1, where the normal beam irradiance Ibn is represented versus the aerosol optical depth for the Solis 2008 model and different solar elevation angles.
This paper presents the development and validation of a new version of the analytical Solis scheme valid for the three radiation components, the global, the beam, and the diffuse, and the four aerosol types, urban, rural, maritime, and tropospheric.

2. Basis of the Solis Scheme

The original Solis clear sky model [12] was first developed within the European program Joule in 2004 in the frame of the Heliosat project, whose aim was to derive and assess the evaluation of the solar irradiance components from the geostationary meteorological satellite MSG (Meteosat Second Generation). The scheme is based on the spectral Radiative Transfer Model (RTM) LibRadTran [13,14] and the Lambert–Beer attenuation relationship:
I b n = I o · e M τ ,
where Io is the extraterrestrial irradiance, Ibn the normal beam irradiance reaching the ground, M the optical air mass, and τ the extinction coefficient. In order to adapt this relation to broadband irradiance, it has to be modified slightly. Three distinct extinction coefficients, one for each of the irradiance components, are defined, as well as three exponent coefficients in the Lambert–Beer attenuation relation. The diffuse irradiance is the result of the scattering process of the beam component during its crossing of the atmosphere. In order to minimize the number of coefficients, a common modified extraterrestrial irradiance Io′ expressed in W/m2 is defined.
The three broadband relations then take the following form:
I b n = I o · e ( τ b / sin b h )   I g h = I o · e ( τ g / sin g h ) · sin h   I d h = I o · e ( τ d / sin d h ) ,
where Ibn, Igh and Idh represent the normal beam irradiance, the horizontal global, and the horizontal diffuse components, respectively, and h is the solar elevation angle.
The model is then driven by seven parameters: the modified extraterrestrial irradiance Io′, the three extinction coefficients τb, τg and τd for the beam, global, and diffuse irradiance, respectively, and b, g, and d are the corresponding sine exponent coefficients.
The justification of these equations can be found in [1,12].

3. Solis Extension to High Turbidity

To circumvent the limitation of the previous Solis scheme, we developed a new version that is applicable for aerosol optical depth values aod up to 7. This large range of turbidity values is needed to avoid unrealistic irradiance values resulting from extreme input parameters to automatic processes. However, one has to keep in mind that above aod = 2, it is not clear that the turbidity can be considered an aerosol optical depth, but a “turbidity” due to bigger particles like sand, or thin, high-altitude clouds. Nevertheless, it is important that in on line automatic production processes, the model does not diverge and still derives coherent values, even with some discrepancies from ground measurements.

3.1. Model Development

The main parameters driving the determination of the solar radiation reaching the ground for a cloudless atmosphere are the aerosol optical depth aod, the total atmospheric precipitable water vapor column w expressed in centimeters, the altitude of the considered site, and the aerosol type. Other atmospheric content such as O3, NOx, etc. have a very low influence of the atmospheric transmittivity and will be considered constant in the development of the model.
In the first step, we made spectral calculations with LibRadTran for each combination of the following parameters.
These RTM calculations permit us to generate the seven coefficients that drive the model: Io’, τb, b, τg, g, τd, and d. The next step is to develop an analytical formulation for these coefficients with the four input parameters given in Table 1.
The analysis of Io′, τb, τg dependence with the aod shows a similar behavior for these three parameters. A third-order polynomial model of the form is applicable:
I o / I o = a i · a o d 3 + b i · a o d 2 + c i · a o d + d i τ b = a b · a o d 3 + b b · a o d 2 + c b · a o d + d b τ g = a g · a o d 3 + b g · a o d 2 + c g · a o d + d g .
The behavior of Io′/Io and τg is represented in Figure 2.
We then analyzed the dependence of each of the a, b, c and d coefficients with the atmospheric water vapor content w. The behavior of the four coefficients shows a common dependence of the form:
n = n 1 · w 0.5 + n 2 · ln w + n 3 ,
where n replaces the a, b, c, and d coefficients for Io′, τb and τg.
The best fit for these coefficients is illustrated in Figure 3.
Finally, the dependence of the three above coefficients on the altitude and the normalized atmospheric pressure p/po (the pressure p at a given altitude normalized by the corresponding sea level pressure po) is analyzed. It appears that a linear regression gives good results.
The model has the following form (for ai, bi, ci, di, ab, bb, cb, db, ag, bg, cg, and dg):
a = a 1 · w 0.5 + a 2 · ln w + a 3 b = b 1 · w 0.5 + b 2 · ln w + b 3 c = c 1 · w 0.5 + c 2 · ln w + c 3 d = d 1 · w 0.5 + d 2 · ln w + d 3
and the general form ni of a, b, c and d coefficient is obtained by a linear function of p/po:
n i = n i 1 · p / p o + n i 2 .
Finally, the inputs for the model for Io′/Io, τb, and τg, consist of the aerosol optical depth at 550 nm, aod550; the water vapor content of the atmosphere, w; and the relative atmospheric pressure, p/po. The corresponding coefficients are given in a 4 × 3 × 2 matrix (see Appendix A).
For τd, the best correlation we found is a relation with τb and τg that has the form:
τ d = c t d 1 + c t d 2 · τ g + c t d 3 / τ b + c t d 4 · τ g 2 + c t d 5 / τ b 2 + c t d 6 · τ g / τ b .
The exponents of sin h in the Lambert–Beer function are best fitted as follows:
g = c a 1 + c a 2 · ln w + c a 3 · ln a o d + c a 4 ln 2 w + c a 5 · ln 2 a o d + c a 6 · ln w · ln a o d b = c b 1 + c b 2 · w + c b 3 · a o d d = c d 1 + c d 2 · ln w + c d 3 · a o d + c d 4 · a o d 2 + c d 5 · a o d 3 + c d 6 · a o d 4 + c d 7 · a o d 5 .
All coefficients for the four aerosol types are given in Appendix A, as well as in the form of an Excel tool [15].

3.2. Solis Analytical Model Versus RTM Calculated Parameter Comparison

The comparison of the seven parameters of the model is expressed as scatter plots between the Solis analytical model parameters and the corresponding LibRadTran calculated parameters. The validation points should be aligned on the 1:1 diagonal. On the graphs are also given the average parameter value, the mean bias difference mbd, the standard deviation sd and the correlation coefficient R2. An illustration is given in Figure 4 for Io′ and the g coefficient.
The average mbd and sd values cannot be taken as an absolute validation of the precision as all the values calculated from the input matrix (altitude, aod, and w) have the same weight. Nevertheless, it gives a good idea of the roughness of the parameter best fit. These values are given in Table 2.

3.3. Modeled versus RTM Calculated Irradiance Comparison

In the same way, the modeled irradiances are compared to the corresponding irradiances estimated from the coefficients calculated with LibRadTran. The validation is expressed as scatter plots and the usual first-order statistics. Here, again due to the unweighted input matrix, the validation values obtained do not give absolute precision. The scatter plots given in Figure 5 illustrate the behavior of the model. They represent the modeled values plotted against the corresponding values evaluated from RTM calculations, for the three components and urban type aerosols.
To obtain the diffuse component, there are two possibilities: the use of the model coefficients Io, τd and d, or to evaluate it from the global and the beam components with the use of the closure equation Idh = Igh – Ibn · sin h. The diffuse component obtained from the closure equation looks slightly better in Figure 4, but it presents negative values for very low solar elevations.
The comparison statistics are given in Table 3.
All the development and validation graphs for the rural aerosol type are given in Appendix B.

4. Model Validation against Ground Measurements

The first step in the model validation process is to apply stringent quality control to the collected ground data. Then, only data acquired under clear sky conditions should be selected; the validation is then applied on these data and the validity of the model is quantified using classical first-order statistical indicators: the mean bias difference (mbd) and the standard deviation (sd).

4.1. Quality Control of the Ground Data

The validity of the results obtained from the use of measured data is highly correlated to the quality of the data bank used as a reference. Controlling data quality is therefore the first step in the process of validating models against ground data. This should be done properly during the acquisition process and automated in order to rapidly detect significant instrumental problems like sensor failure, errors in calibration, orientation, leveling, tracking, consistency, etc. Normally, this quality control process should be done by the institution responsible for the measurements. Unfortunately, this is not always the case. Even if some quality control procedures have been implemented, they might not be sufficient to catch all errors, or the data points might not be flagged to indicate the source of the problem. A stringent control quality procedure must therefore be adopted in the present context, and its various elements are given in what follows.
If the three solar irradiance components, beam, diffuse, and global, are available, a consistency test can be applied, based on the closure equation that links them:
I g h = I d h + I b n · sin h .
A posteriori automatic quality control cannot detect all acquisition problems, however. The remaining elements to be assessed are threefold:
  • The measurement’s time stamp (needed to compute the solar geometry);
  • The sensors’ calibration coefficient used to convert the acquired data into physical values;
  • The coherence between the parameters.
It has to be noted that, for the data of Solar Village, even if the closure equation satisfies the selection conditions, as can be seen in Figure 6, a clear pattern appears due to calibration differences of up to 6–8% between the winter months and the summer months. The bias between the beam irradiance measurements and the closure equation values are negative in the summer months; it is not clear what the source of the difference is. This can be an issue of the temperature and/or solar altitude dependence of the sensors; it is difficult to find out a posteriori what component is questionable and why.
A complete description of the quality control procedure is given by Ineichen [16].

4.2. Clear Conditions Selection

To perform the selection, the criteria defined in Ineichen [17] are applied:
  • If three components are available, the closure equation must be satisfied within −50 W/m2 − 5% and +50 W/m2 + 5%; if only two components are available, all the values are kept.
  • The modified global clearness Kt′ (as defined by Perez et al. in [18]) of the measurements is higher than 0.65.
  • The stability of the global clearness index ΔKt′ is better than 0.01 (ΔKt′ is evaluated by the difference of the considered hour and the average of the considered hour, the preceding one, and the following one)
This selection minimizes the cloud contamination. It is restrictive, but for the purposes of the comparison, it ensures that only clear and stable conditions are selected. It shows its importance particularly with data from Xianghe, where often the beam component is very low or missing in the middle of the day, while the corresponding global component shows clear conditions (due to saturation?). This is illustrated in Figure 7.

4.3. Statistical Indicators and Graphical Representation

The comparison is done on an hourly basis, the model–measurements difference is computed, so that a positive value of the mean bias difference represents an overestimation of the model. The following indicators are used to quantify model performance:
  • First-order statistics for a given site: the mean bias difference (mbd) and the standard deviation (sd). In addition, qualitative visualization is made with the help of model versus measurement scatterplots.
  • The seasonal dependence of the bias and its dependence with the aerosol optical depth (aod);
  • The frequency distribution of the model–measurements differences and the corresponding cumulated frequencies.

4.4. Ground Measurements for the Validation

To assess the validity of the model for low turbidity values, we first made a validation of the same data as in [18] acquired in Geneva for the years 2004 to 2011, where the aod is pretty low.
The average aerosol optical depth in Geneva is on the order of 0.17, with a maximum of 0.5 during polluted episodes. As no specific aerosol optical depth values are measured in Geneva, we used one of the best state-of-the-art models, REST2 [5], to retrieve by retrofit the aod following the method described in [19].
The next validation step is to do a comparison with ground data acquired at a site with higher turbidity values. To conduct the validation, sites that jointly acquire the irradiance components and the aerosol optical depth are needed. They are very scarce, so that not all conditions could be analyzed. Irradiance data from three Baseline Surface Radiation Network (BSRN) sites and from one Indian National Institute of Wind Energy (NIWE) site are used for the validation. The aerosol optical depths are retrieved from the Aerosol Robotic Network (AERONET). The latitude, longitude, altitude, networks, and acquired parameters are given in Table 4.
The median values of the aod for the four sites are situated between 0.3 and 0.55, with maximum values around 1.5 to 4. The cumulated frequencies of occurrence of the aerosol optical depth and the atmospheric water vapor content are given in Figure 8.

5. Validation Results and Discussion

The results of the validation against the ground data acquired in Geneva are illustrated in Figure 9 for the normal beam and the global horizontal irradiance. Done on the same dataset as in [17], it shows that the first-order statistics are slightly better than for the previous version of Solis.
The first-order statistics of the model validation against data acquired at high turbidity sites are given in Table 5 for the rural aerosol type, and in Table 6 for the urban aerosol type. The tables report the number of nb values kept for the validation, the average measured irradiance and irradiation, the relative mean bias difference, mbd, and the standard deviation, sd. Results are provided for the three irradiance components and for all sites. Due to the different climates and quality of the data, the number of selected data can be very different from one site to the other.
The main results drawn from the table are the following:
  • For all sites, choosing the rural aerosol type gives the best results in terms of mean bias and standard deviation, especially for the global component.
  • For Ilorin, the high standard deviation is due to the high variability of the sky conditions; it can also be a result of the beam retrieval from the global and the diffuse components. For Xianghe, the poor quality of the data induces this higher standard deviation.
  • The Jaipur and Solar Village datasets show the overall best results in term of mean bias and standard deviation.
  • The “daily values” represent the sum of only the hourly values kept for the validation: it is not the daily integral (the sum of all the values of the day). This is particularly the case during partially cloudy days where only a few hourly values are acquired during cloudless conditions. The consequence is that, for example, at Ilorin, where the sky conditions are highly variable and clear days from sunrise to sunset are very scarce, the average “daily value” is low.
A deeper analysis of the results leads to the following general observations:
  • Ground measurements and modeled values can be represented on the same graph in different colors. If the model faithfully reproduces the measurements, the envelopes of the data clouds should be similar. This is illustrated for the site of Solar Village in Figure 10 where the beam clearness index is plotted against the global clearness index on the left graph (a) and the diffuse fraction against the global clearness index on the right graph (b); the measurements are represented in green and the modeled values in blue.
  • The trends of the bias as a function of the aerosol optical depth show no particular pattern for all the sites, except for Xianghe, where a non-significant trend is due to the dispersion. These dependences are illustrated in Figure 11 in hourly values, for the site of Solar Village and rural aerosol type.
  • The seasonal dependence of the bias for both the global and the beam irradiance components shows no significant pattern. An illustration is given in Figure 12 for data from Solar Village in hourly values.
  • With the exception of the site of Ilorin, the distributions of the bias around the 1:1 model–measurements axis are near normal; for this reason, the first-order statistics represented by the mean bias and the standard deviation can be considered reliable. An example is given in Figure 13 for the site of Solar Village. The figure displays the hourly bias frequency distribution for both Ibn and Igh irradiance components. The cumulated frequencies (red curve) are also represented in the graphs.

6. Conclusions

When dealing with satellite images to derive the irradiance components worldwide and every 15 min, the calculation time should be as short as possible. The analytical Solis clear sky scheme was developed to fulfill this requirement but had the weakness of being limited to aerosol optical depth values lower than 0.45. The new analytical Solis scheme is valid for aod550 values up to seven, even if very high values are not realistic as an optical depth; it is probably due to bigger particles like sand, or thin, high-altitude clouds. Nevertheless, in contrary to other clear sky models, this permits us to produce coherent irradiance values, even with questionable input data. Indeed, for example, the aod, due to the difficulty of retrieving it by automatic processes, can take on out-of-range or unrealistic values.
The new Solis clear sky scheme is developed for aerosol optical depths values aod550 from 0.02 to 7, for atmospheric water vapor w content from 0.01 to 10 cm, for altitude from sea level to 7000 m, and for four aerosol types, as defined by Shettle and Fenn [2]. For the urban, rural, and tropospheric aerosol types, the validation against the original LibRadTran RTM calculations presents no bias and a standard deviation lower than 1% for the global and the beam component, and 3% for the diffuse. When dealing with maritime aerosol type, the standard deviation is 2.2%, 3.4%, and 4%, respectively, for the global, the beam, and the diffuse components.
The validation against ground measurements acquired in Geneva, Ilorin, Jaipur, Solar Village, and Xianghe gives a mean bias difference less than ±4%, a standard deviation of 3–6% for the global component, and a mean bias difference less than ±3% with a standard deviation of 5–13% for the beam component. The 12–13% standard deviations are due to either the high variability of the sky conditions or the poor quality of the ground data.
As it is very difficult to find ground measurements covering the complete range of validity of the new Solis analytical scheme, it cannot be completely validated. However, the new version faithfully reproduces the RTM LibRadTran calculations.

Nomenclature

SymbolDescriptionSymbolDescription
Ibn or Bnnormal beam irradianceMoptical air mass
Igh or Ghglobal horizontal irradiancehsolar elevation
Idh or Dhdiffuse horizontal irradiancepatmospheric pressure
Iosolar constantpoatm. pressure at sea level
Io'modified solar constant
aodaerosol optical depthnbnumber of involved values
watmospheric water vapor columnmbdmean bias difference
τextinction coefficientsdstandard deviation
RTMRadiation Transfer ModelR2correlation coefficient
BSRNBaseline Surface Radiation Network
NIWENational Institute of Wind Energy
AERONETAERosol Robotic NETwork

Appendix A. Model Coefficients

Solis 2017 Rural Aerosol Type
Io'/Io
111221223132
ai0.000261031−0.001810295−0.000334819−0.001533816−0.0001183430.000100560
bi0.0199082350.1958654980.0196498060.0305139080.0032462440.010001152
ci0.0357657040.7375720970.0170743080.0217126390.0048442090.005917312
di0.0728759471.0035943130.009537179−0.0001367490.0018831350.000378418
τg
111221223132
ag0.000193733−0.002501003−0.000156219−0.000277670−0.0000590400.000006851
bg−0.0011599800.0585241580.0025482290.0048297490.0007885220.000522673
cg−0.001449466−0.833692155−0.012697614−0.030033430−0.004681306−0.011682439
dg−0.141809777−0.085671381−0.015070590−0.030811277−0.007315151−0.008232404
coef a
ca1ca2ca3ca4ca5ca6
0.393927007−0.014924316−0.092174236−0.001048172−0.0091630930.006108964
τb
111221223132
ab0.000771802−0.007478508−0.000180242−0.000368091−0.0000499380.000049572
bb−0.0095651740.1487749900.0023365390.0070875300.0011026800.000284776
cb0.024325289−1.461925559−0.009964879−0.043313859−0.007992277−0.015281853
db−0.191282122−0.099096810−0.022005054−0.027680914−0.006793639−0.008818264
coef b
cb1cb2cb3
0.482261237−0.0166786720.914171831
τd
ctd1ctd2ctd3ctd4ctd5ctd6
1.6260899460.5814433000.602607926−0.0419350670.048133609−3.902240255
coef d
cd1cd2cd3cd4cd5cd6cd7
0.113144991−0.0041652290.621087873−0.3592489540.092564432−0.0112839850.000525184
Solis 2017 Urban Aerosol Type
Io'/Io
111221223132
ai0.0005180110.002795295−0.000216902−0.001116789−0.0000715070.000376597
bi0.0189992240.1806688370.0204880320.0286053980.0029689670.008819303
ci0.0350156700.5996754840.0085124090.0150891450.0043389100.004266916
di0.0720176311.0052994070.0095729120.0003813920.0019043390.000343349
τg
111221223132
ag0.000428378−0.000564884−0.000375237−0.0000727980.000128450−0.000046323
bg−0.0027330050.0362148850.0043749670.003084955−0.0010125000.000930531
cg0.003834122−0.883655173−0.015857854−0.028437392−0.001164850−0.013464145
dg−0.141741963−0.079752002−0.014920981−0.030293926−0.007709979−0.007819176
coef a
ca1ca2ca3ca4ca5ca6
0.436328716−0.015982197−0.099473876−0.001106107−0.0132733970.006225334
τb
111221223132
ab0.000565740−0.005280564−0.000003576−0.000404517−0.0000448660.000073274
bb−0.0074583420.1180926080.0006406650.0074248990.000959090−0.000036289
cb0.020375757−1.380477103−0.006939622−0.043687713−0.007132214−0.014709261
db−0.192470443−0.093585418−0.021132099−0.027442728−0.007179027−0.008451631
coef b
cb1cb2cb3
0.498959654−0.0176369160.926155270
τd
ctd1ctd2ctd3ctd4ctd5ctd6
1.8869307290.5410831410.605281729−0.0413224840.048799518−4.249335048
coef d
cd1cd2cd3cd4cd5cd6cd7
0.112629823−0.0030887670.569910110−0.3378792160.087849430−0.0107512820.000501473
Solis 2017 Tropospheric Aerosol Type
Io'/Io
111221223132
ai−0.000147911−0.008122292−0.001047420−0.001931206−0.000132858−0.000370950
bi0.0222920280.1585748110.0209479360.0256836880.0026563820.010480509
ci0.0328850910.7579046430.0143096590.0309287210.0058700700.005967166
di0.0734286210.9999001990.009230660−0.0011211320.0018549480.000311441
τg
111221223132
ag0.000217441−0.003126729−0.000132644−0.000291300−0.000002894−0.000029457
bg−0.0018092580.0680302150.0021617510.0051602920.0003105470.000821525
cg−0.004470877−0.811312676−0.010715753−0.031925980−0.004215436−0.011986451
dg−0.138496207−0.088561634−0.016862218−0.029763820−0.006341627−0.008895780
coef a
ca1ca2ca3ca4ca5ca6
0.395322938−0.015327403−0.090285920−0.001118039−0.0087275680.006173611
τb
111221223132
ab0.000953948−0.008496386−0.000088009−0.000432762−0.0000567250.000028537
bb−0.0107865340.1632752310.0015125160.0077634340.0011261510.000593655
cb0.022706050−1.435580065−0.008631610−0.044554649−0.007857069−0.016684930
db−0.192368280−0.099465438−0.020840853−0.028520351−0.007159173−0.008504057
coef b
cb1cb2cb3
0.481658656−0.0165556670.891497733
τd
ctd1ctd2ctd3ctd4ctd5ctd6
1.5804524050.5580864490.586357700−0.0456636450.045790439−3.883558603
coef d
cd1cd2cd3cd4cd5cd6cd7
0.113480292−0.0045187700.624749275−0.3585748300.092385725−0.0112770250.000525609
Solis 2017 Maritim Aerosol Type
Io'/Io
111221223132
ai0.0147871720.0812673710.0076116570.0056872110.0010572160.004652006
bi−0.0290980680.099939972−0.0044662200.0171233330.000502026−0.005634380
ci0.1350107021.0163609830.0435788390.0331172030.0094478320.024152082
di0.0620397420.9721522570.005356544−0.0011779700.001368473−0.002323682
τg
111221223132
ag−0.000258280−0.0009915980.000087922−0.000450453−0.0001032830.000036369
bg0.0038328400.038507962−0.0005399190.0066631500.0013475140.000076692
cg−0.020444433−0.882340428−0.002557914−0.034648986−0.006330049−0.009504023
dg−0.138105804−0.087198021−0.019195997−0.028342978−0.006284626−0.009009817
coef a
ca1ca2ca3ca4ca5ca6
0.370255687−0.013481692−0.096888113−0.000980585−0.0089661970.006430510
τb
111221223132
ab−0.000015662−0.0041948810.000107149−0.000489497−0.0001189760.000094402
bb−0.0038980580.107002592−0.0001869230.0075636740.001566986−0.000369907
cb0.010411340−1.597417220−0.006692200−0.040428247−0.007628685−0.012973060
db−0.198484327−0.092013892−0.019285874−0.029172125−0.007407717−0.008863425
coef b
cb1cb2cb3
0.487380474−0.0160885750.952963502
τd
ctd1ctd2ctd3ctd4ctd5ctd6
1.6330651560.6764966070.651637829−0.0339191320.053967101−3.761466993
coef d
cd1cd2cd3cd4cd5cd6cd7
0.109205787−0.0038275790.635588699−0.3726377060.095930104−0.0116452470.000539553

Appendix B. Validation Graphs

Remotesensing 10 00435 i001
Remotesensing 10 00435 i002
Remotesensing 10 00435 i003
Remotesensing 10 00435 i004a
Remotesensing 10 00435 i004b
Remotesensing 10 00435 i005a
Remotesensing 10 00435 i005b
Remotesensing 10 00435 i006
Remotesensing 10 00435 i007
Remotesensing 10 00435 i008
Remotesensing 10 00435 i009
Remotesensing 10 00435 i010a
Remotesensing 10 00435 i010b
Remotesensing 10 00435 i011a
Remotesensing 10 00435 i011b
Remotesensing 10 00435 i012

References

  1. Gueymard, C. A two-band model for the calculation of clear Sky Solar Irradiance, Illuminance, and Photosynthetically Active Radiation at the Earth Surface. Sol. Energy 1989, 43, 253–265. [Google Scholar] [CrossRef]
  2. Gueymard, C. Direct solar transmittance and irradiance predictions with broadband models. Part 1: Detailed theoretical performance assessment. Sol. Energy 2003, 74, 355–379. [Google Scholar] [CrossRef]
  3. Bird, R.E.; Huldstrom, R.L. Direct insolation models. Trans. ASME J. Sol. Energy Eng. 1980, 103, 182–192. [Google Scholar] [CrossRef]
  4. Wald, L. European Solar Radiation Atlas; Commission of the European Communities Presses de l’Ecole, Ecole des Mines de Paris: Paris, France, 2000. [Google Scholar]
  5. Rigollier, C.; Bauer, O.; Wald, L. On the Clear Sky Model of the ESRA—European Solar Radiation Atlas—With Respect to the Heliosat Method. Sol. Energy 2000, 68, 33–48. [Google Scholar] [CrossRef]
  6. Geiger, M.; Diabaté, L.; Ménard, L.; Wald, L. A web Service for Controlling the Quality of Measurements of Global Solar Irradiation. Sol. Energy 2002, 73, 475–480. [Google Scholar] [CrossRef]
  7. Linke, F. Transmissions-Koeffizient und Trübungsfaktor. Beiträge zur Physik der freien Atmosphäre 1992, 10, 91–103. [Google Scholar]
  8. Ineichen, P. A broadband simplified version of the Solis clear sky model. Sol. Energy 2008, 82, 758–762. [Google Scholar] [CrossRef]
  9. Shettle, E.P.; Fenn, R.W. Models Fort the Aerosol of Lower Atmosphere and the Effect of Humidity Variations on Their Optical Properties. 1979. Available online: http://www.dtic.mil/docs/citations/ADA085951 (accessed on 2 March 2018).
  10. Shettle, E.P. Models of aerosols, clouds and precipitation for atmospheric propagation studies. In Proceedings of the AGARD, Atmospheric Propagation in the UV, Visible, IR, and MM-Wave Region and Related Systems Aspects Conference, Copenhagen, Denmark, 9–13 October 1989. [Google Scholar]
  11. Zhang, T.; Stackhouse, P.W., Jr.; Chandler, W.S.; Westberg, D.J. Application of a global-to-beam irradiance model to the NASA GEWEX SRB dataset: An extension of the NASA Surface meteorology and Solar Energy datasets. Sol. Energy 2014, 110, 117–131. [Google Scholar] [CrossRef]
  12. Mueller, R.W.; Dagestad, K.F.; Ineichen, P.; Schroedter-Homscheidt, M.; Cros, S.; Dumortier, D.; Kuhlemann, R.; Olseth, J.A.; Piernavieja, G.; Reise, C.; et al. Rethinking satellite based solar irradiance modelling—The SOLIS clear-sky module. Remote Sens. Environ. 2004, 91, 160–174. [Google Scholar] [CrossRef]
  13. Mayer, B.; Kylling, A. Technical note: The libRadtran software package for radiative transfer calculations—Description and examples of use. Atmos. Phys. 2005, 5, 1855–1877. [Google Scholar] [CrossRef]
  14. Mayer, B.; Kylling, A.; Emde, C.; Buras, R.; Hamann, U. LibRadTran: Library for Radiative Transfer Calculations, Edition 1.0 for Libradtran Version 1.5-beta. 2010. Available online: https://clouds.eos.ubc.ca/~phil/courses/eosc582_2010/textfiles/libRadtran.pdf (accessed on 3 March 2018).
  15. Ineichen, P. Solis 2017 Excel Tool. Password: Solis2017. Available online: http://www.adpi.ch/Solis2017/Solis2017-tool.xlsx (accessed on 28 February 2017).
  16. Ineichen, P. Satellite Irradiance Based on MACC Aerosols: Helioclim 4 and SolarGIS, Global and Beam Components Validation. In Proceedings of the EuroSun 2014, Aix-les-Bains, France, 16–19 Septembre 2014. [Google Scholar]
  17. Ineichen, P.; Barroso, C.S.; Geiger, B.; Hollmann, R.; Marsouin, A.; Mueller, R. Satellite Application Facilities irradiance products: hourly time step comparison and validation over Europe. Int. J. Remote Sens. 2009, 30, 5549–5571. [Google Scholar] [CrossRef]
  18. Perez, R.; Ineichen, P.; Seals, R.; Zelenka, A. Making full use of the clearness index for parameterizing hourly insolation conditions. Sol. Energy 1990, 45, 111–114. [Google Scholar] [CrossRef]
  19. Ineichen, P. Validation of models that estimate the clear sky global and beam solar irradiance. Sol. Energy 2016, 132, 332–344. [Google Scholar] [CrossRef]
Figure 1. Behavior of the normal beam irradiance with the aerosol optical depth for the 2008 Solis version and different solar elevation angles.
Figure 1. Behavior of the normal beam irradiance with the aerosol optical depth for the 2008 Solis version and different solar elevation angles.
Remotesensing 10 00435 g001
Figure 2. Behavior of the enhanced solar constant Io′/Io (a) and the extinction coefficient τ (b) with the spectral aerosol optical depth at 550 nm.
Figure 2. Behavior of the enhanced solar constant Io′/Io (a) and the extinction coefficient τ (b) with the spectral aerosol optical depth at 550 nm.
Remotesensing 10 00435 g002
Figure 3. Behavior of the cubic equation coefficients a, b, c and d, respectively graph (ad), with the water vapor column w.
Figure 3. Behavior of the cubic equation coefficients a, b, c and d, respectively graph (ad), with the water vapor column w.
Remotesensing 10 00435 g003
Figure 4. Modeled Io′ coefficient (a) and modeled g coefficient (b) versus the corresponding RTM calculated value.
Figure 4. Modeled Io′ coefficient (a) and modeled g coefficient (b) versus the corresponding RTM calculated value.
Remotesensing 10 00435 g004
Figure 5. Modeled against RTM calculated irradiance comparison for the beam (a), the global (b), the diffuse (c), and the closure diffuse (d) components.
Figure 5. Modeled against RTM calculated irradiance comparison for the beam (a), the global (b), the diffuse (c), and the closure diffuse (d) components.
Remotesensing 10 00435 g005
Figure 6. Illustration of the calibration uncertainty for the site of Solar Village.
Figure 6. Illustration of the calibration uncertainty for the site of Solar Village.
Remotesensing 10 00435 g006
Figure 7. Illustration of the questionable data for the site of Xianghe.
Figure 7. Illustration of the questionable data for the site of Xianghe.
Remotesensing 10 00435 g007
Figure 8. Cumulated frequency of occurrence of the aerosol optical depth (a) and the atmospheric water vapor content (b) for the four sites.
Figure 8. Cumulated frequency of occurrence of the aerosol optical depth (a) and the atmospheric water vapor content (b) for the four sites.
Remotesensing 10 00435 g008
Figure 9. Validation against ground data acquired in Geneva for the beam irradiance (a) and the global component (b). The aod is retrofitted with REST2. For comparison, the same graphs are given for the previous version of Solis (c,d).
Figure 9. Validation against ground data acquired in Geneva for the beam irradiance (a) and the global component (b). The aod is retrofitted with REST2. For comparison, the same graphs are given for the previous version of Solis (c,d).
Remotesensing 10 00435 g009
Figure 10. Beam clearness index (a) and diffuse fraction (b) versus the global clearness index.
Figure 10. Beam clearness index (a) and diffuse fraction (b) versus the global clearness index.
Remotesensing 10 00435 g010
Figure 11. Bias versus the aerosol optical depth for normal beam irradiance (a) and global irradiance (b) for the site of Solar Village.
Figure 11. Bias versus the aerosol optical depth for normal beam irradiance (a) and global irradiance (b) for the site of Solar Village.
Remotesensing 10 00435 g011
Figure 12. Seasonal dependence of the bias for the normal beam irradiance (a) and the global irradiance (b) irradiance components for the site of Solar Village.
Figure 12. Seasonal dependence of the bias for the normal beam irradiance (a) and the global irradiance (b) irradiance components for the site of Solar Village.
Remotesensing 10 00435 g012
Figure 13. Frequency of occurrence of the bias around the 1:1 axis for the beam (a) and the global (b) irradiance components.
Figure 13. Frequency of occurrence of the bias around the 1:1 axis for the beam (a) and the global (b) irradiance components.
Remotesensing 10 00435 g013
Table 1. Aerosol types, altitudes, optical depths aod550, and water vapor columns w values for the RTM calculations with LibRadTran.
Table 1. Aerosol types, altitudes, optical depths aod550, and water vapor columns w values for the RTM calculations with LibRadTran.
Aerosol TypeAltitudeAod550W
urbansea level0.010.01
rural500 m0.030.03
maritim1000 m0.050.05
tropospheric2000 m0.11
3000 m0.150.15
4000 m0.20.2
5000 m0.40.3
6000 m0.70.5
7000 m11
1.51.5
22
33
44
56
68
710
Table 2. Average, mbd, sd and R2 for the Io′, τb, b, τg, g, τd, and d model.
Table 2. Average, mbd, sd and R2 for the Io′, τb, b, τg, g, τd, and d model.
Io'τga
aerosol typeaveragembdsdR2averagembdsdR2averagembdsdR2
rural62840930.998−1.490.000.011.0000.410.000.020.993
urban60360950.997−1.680.000.021.0000.440.000.020.991
tropospheric53800400.999−1.400.000.011.0000.410.000.020.992
maritime12,210−81340.998−1.690.000.011.0000.390.000.020.994
τbb
rural −2.210.000.011.0000.420.000.020.965
urban −2.210.000.011.0000.440.000.030.941
tropospheric −2.080.000.011.0000.400.000.020.976
maritime −2.680.000.011.0000.450.000.020.912
τdc
rural −2.750.010.040.9980.320.000.010.995
urban −3.080.010.050.9990.290.000.010.995
tropospheric −2.670.010.040.9980.330.000.010.995
maritime −2.860.010.040.9990.310.000.010.995
Table 3. Average, mbd, sd, and R2 for the irradiance components.
Table 3. Average, mbd, sd, and R2 for the irradiance components.
IghIbnIdh
aerosol typeaveragembdsdR2averagembdsdR2averagembdsdR2
rural727041.000549051.000223070.998
urban632071.000543031.000167070.997
tropospheric740051.000559071.000227070.999
maritime753−1170.997519−3180.9992662110.998
Table 4. Latitude, longitude, altitude, networks, and acquired parameters for the considered sites.
Table 4. Latitude, longitude, altitude, networks, and acquired parameters for the considered sites.
SiteNetworkLatitudeLongitudeAltitudeGlobalBeamDiffuseTaHrAod
IlorinBSRN/AERONET8.534.57350x xxxx
JaipurNIWE/AERONET26.8175.86403xxxxxx
Solar VillageBSRN/AERONET24.9146.41650xxxxxx
XiangheBSRN/AERONET39.75116.6922xxx x
Table 5. First-order statistics for the four sites and rural aerosol type. The irradiance is expressed in (W·h/m2·h) and the irradiation in (kWh/m2·day).
Table 5. First-order statistics for the four sites and rural aerosol type. The irradiance is expressed in (W·h/m2·h) and the irradiation in (kWh/m2·day).
Hourly Values
SiteNormal BeamHorizontal GlobalHorizontal Diffuse
nbAveragembdsdAveragembdsdAveragembdsd
Ilorin 1998–2005289646−2%12%8451%5%2772%22%
Jaipur 2016–20173876704%5%7153%3%1984%13%
Solar Village 1999–20023198790−2%5%784−1%4%1654%19%
Xianghe 2008–20151161762−2%13%612−5%6%132−10%33%
"daily values"
nbaveragembdsdaveragembdsdaveragembdsd
Ilorin 1998–2005392.30%8%3.01%4%0.9−1%20%
Jaipur 2016–2017693.64%4%3.83%2%1.05%10%
Solar Village 1999–20024535.5−2%3%5.5−1%3%1.24%12%
Xianghe 2008–20151784.9−2%14%3.9−5%6%0.8−10%33%
Table 6. First-order statistics for the four sites and urban aerosol type. The irradiance is expressed in (W·h/m2·h) and the irradiation in (kWh/m2·day).
Table 6. First-order statistics for the four sites and urban aerosol type. The irradiance is expressed in (W·h/m2·h) and the irradiation in (kWh/m2·day).
Hourly Values
SiteNormal BeamHorizontal GlobalHorizontal Diffuse
nbAverageMbdSdAverageMbdSdAverageMbdSd
Ilorin 1998—2005289646−3%12%845−8%7%277−20%20%
Jaipur 2016—20173876703%6%715−4%4%198−16%11%
Solar Village3198790−3%5%784−5%4%165−14%19%
Xianghe 2008—20151161762−2%13%612−8%8%132−24%30%
Daily Values
nbAverageMbdsdAverageMbdSdAverageMbdSd
Ilorin 1998—2005392.3−2%9%3.0−7%6%0.9−22%19%
Jaipur 2016—2017693.63%4%3.8−3%3%1.0−15%8%
Solar Village4535.5−3%3%5.5−5%3%1.2−13%13%
Xianghe 2008—20151784.9−2%14%3.9−8%9%0.8−24%26%

Share and Cite

MDPI and ACS Style

Ineichen, P. High Turbidity Solis Clear Sky Model: Development and Validation. Remote Sens. 2018, 10, 435. https://doi.org/10.3390/rs10030435

AMA Style

Ineichen P. High Turbidity Solis Clear Sky Model: Development and Validation. Remote Sensing. 2018; 10(3):435. https://doi.org/10.3390/rs10030435

Chicago/Turabian Style

Ineichen, Pierre. 2018. "High Turbidity Solis Clear Sky Model: Development and Validation" Remote Sensing 10, no. 3: 435. https://doi.org/10.3390/rs10030435

APA Style

Ineichen, P. (2018). High Turbidity Solis Clear Sky Model: Development and Validation. Remote Sensing, 10(3), 435. https://doi.org/10.3390/rs10030435

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

Article Metrics

Back to TopTop