Derivatives of Spectral Aerosol Optical Depth for Partitioning Type and Loading

Abstract

A spectral method is provided for partitioning type and loading with aerosol optical depth. Based on multi-spectral optical aerosol depth, particle-size distribution and refractive index are derived by normalizing first- and second-order derivatives for processing quantitative calibration of main components. According to the optical feature parameters of various aerosol types, a radiation theory is applied to simulate multi-spectral optical depth for each density, including those of mixed types. The intrinsic parameters of aerosol types are figured out by constructing normalized derivative aerosol indices (NDAI). The clear characteristic differences between aerosol types are used to figure out main components of aerosols and their mixing ratios. The simulation result of the normalized index of various aerosol type is in good agreement with the ground observation data of Aerosol Robotic Network. It shows that NDAI is quite practicable in quantitative calibration of main components of atmospheric aerosol.

Claims

1. A method of spectral derivatives of aerosol optical depth (AOD) for partitioning type and loading, comprising steps of: (a) first step: based on optical feature parameters of various aerosol types, obtaining a model of Second Simulation of a Satellite Signal in the Solar Spectrum (6S model) to calculate spectral AODs of said various aerosol types, wherein said various aerosol types comprises dust (DS), biomass burning (BB), anthropogenic pollutants (AP), and various mixtures of DS, BB, and AP; and (b) second step: based on said spectral AODs of said various aerosol types, processing calculation with normalized derivative aerosol indices (NDAI) to obtain particle-size distributions and complex refractive indices derived from normalized first- and second-order derivatives of said spectral AODs of said various aerosol types to obtain intrinsic parameters of said various aerosol types to calculate main components of aerosols and mixing ratios thereof to identify each single type of aerosol and quantitatively distinguish main components of mixed aerosol.

2. The method according to claim 1, wherein a main component of said BB is black carbon and main components of said AP are sulfate and nitrate.

3. The method according to claim 1, wherein said optical feature parameters of said various aerosol types are based on said particle-size distributions and said complex refractive indices provided by the World Meteorological Organization (WMO).

4. The method according to claim 1, wherein, with the normalization of said first-order derivatives of said spectral AODs of said various aerosol types, said particle-size distributions of DS, BB, and AP are clearly distinguished by an equation as follows: NDAI.sub.(λ.sub.1.sub.,λ.sub.2.sub.)≡∇τ.sub.(λ.sub.1.sub.,λ.sub.2.sub.)/τ.sub.λ.sub.ref, wherein λ is a wavelength (micrometer, μm); ∇τ.sub.(λ.sub.1.sub.,λ.sub.2.sub.) is a first-order derivative of spectral AOD of gap between λ.sub.1 and λ.sub.2; τ.sub.λ.sub.ref is a normalization reference of loading of spectral AOD; and NDAI.sub.(λ.sub.1.sub.,λ.sub.2.sub.) is a spectral derivative of gap between λ.sub.1 and λ.sub.2 as being normalized with τ.sub.λ.sub.ref.

5. The method according to claim 1, wherein, with the normalization of said second-order derivatives of said spectral AODs of said various aerosol types, features of and differences between said various aerosol types on scattering and absorption are obtained to distinguish and identify said various aerosol types by an equation as follows: $\frac{{\partial }^2\tau }{\partial {\lambda }^2}/{\tau }_{{\lambda }_{\mathrm{ref}}}\cong {\nabla }^2{\tau }_{\left({\lambda }_1,{\lambda }_2,{\lambda }_3\right)}/{\tau }_{{\lambda }_{\mathrm{ref}}}=\frac{\left({\tau }_{{\lambda }_1}-2{\tau }_{{\lambda }_2}+{\tau }_{{\lambda }_3}\right)/{\tau }_{{\lambda }_{\mathrm{ref}}}}{\left({\lambda }_1-{\lambda }_2\right)\left({\lambda }_2-{\lambda }_3\right)},$ wherein τ.sub.λ.sub.ref is a normalization reference of loading of spectral AOD and is selected from τ.sub.λ.sub.1, τ.sub.λ.sub.2, and τ.sub.λ.sub.3 to bind a dynamic range to spectral derivative of said various aerosol type; and λ is a wavelength.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0010] The present invention will be better understood from the following detailed description of the preferred embodiment according to the present invention, taken in conjunction with the accompanying drawings, in which

[0011] FIG. 1 is the view showing the preferred embodiment according to the present invention;

[0012] FIG. 2 is the view showing the spectral distributions of the spectral aerosol optical depths (AOD);

[0013] FIG. 3A is the view showing the first-order derivatives of the unnormalized AODs;

[0014] FIG. 3B is the view showing the first-order derivatives of the normalized AODs;

[0015] FIG. 3C is the view showing the second-order derivatives of the unnormalized AODs;

[0016] FIG. 3D is the view showing the second-order derivatives of the normalized AODs;

[0017] FIG. 4A is the view showing the unnormalized first-order parameters;

[0018] FIG. 4B is the view showing the normalized first-order parameters;

[0019] FIG. 4C is the view showing the unnormalized second-order parameters;

[0020] FIG. 4D is the view showing the normalized second-order parameters;

[0021] FIG. 4E is the view showing the unnormalized second-order derivatives of the AOD intrinsic features;

[0022] FIG. 4F is the view showing the normalized second-order derivatives of the AOD intrinsic features;

[0023] FIG. 5 is the view showing the result integrating theoretical simulation and actual observation;

[0024] FIG. 6 is the view showing the result integrating theoretical simulation and actual observation by using the normalized aerosol indices;

[0025] FIG. 7 is the view showing the normalized first- and second-order derivatives of the data of ground observation and theoretical simulation; and

[0026] FIG. 8A and FIG. 8B are the views showing the component proportions of the three representative aerosols.

DESCRIPTION OF THE PREFERRED EMBODIMENT

[0027] The following description of the preferred embodiment is provided to understand the features and the structures of the present invention.

[0028] Please refer to FIG. 1 to FIG. 8B, which are a view showing a preferred embodiment according to the present invention; a view showing spectral distributions of AODs; a view showing first-order derivatives of unnormalized AODs; a view showing first-order derivatives of normalized AODs; a view showing second-order derivatives of unnormalized AODs; a view showing second-order derivative of normalized AODs; a view showing unnormalized first-order parameters; a view showing normalized first-order parameters; a view showing unnormalized second-order parameters; a view showing normalized second-order parameters; a view showing unnormalized second-order derivatives of AOD intrinsic features; a view showing normalized second-order derivatives of AOD intrinsic features; a view showing a result integrating theoretical simulation and actual observation; a view showing a result integrating theoretical simulation and actual observation by normalized aerosol indices; a view showing normalized first- and second-order derivatives of data of ground observation and theoretical simulation; and a view showing component proportions of three representative aerosols. As shown in the figures, the present invention is a method of spectral AOD derivatives for partitioning type and loading, where a normalized derivative aerosol index (NDAI) is used to integrate data of theoretical simulation and actual observation for examining various aerosol types with the relationships of particle-size distributions and complex refractive indices together with first- and second-order derivatives of spectral AODs and constructing optical intrinsic parameters of dust (DS), biomass burning (BB), and anthropogenic pollutants (AP); and, thus, each single type of aerosol is identified and main components of each mixed aerosol are quantitatively distinguished. The present invention comprises the following steps:

[0029] (a) Processing theoretical simulation s1: Regarding a theoretical simulation, based on optical feature parameters of various aerosol types, a model of Second Simulation of a Satellite Signal in the Solar Spectrum (6S model) is used to calculate spectral aerosol optical depths (AOD) of the various aerosol types. The various aerosol types comprises DS, BB, AP, and various mixtures of DS, BB, and AP, where the main component of BB is black carbon and the main components of AP are sulfate and nitrate. Therein, the optical feature parameters of the various aerosol types are based on particle-size distributions and complex refractive indices of aerosols provided by the World Meteorological Organization (WMO). As listed in Table 1, n.sub.r and n.sub.i are the real number part and the imaginary number part of the complex refractive index, respectively; R.sub.mean is a geometric mean radius; and R.sub.std is a geometric standard deviation.

TABLE-US-00001 TABLE 1 λ DS AP BB (micrometer, μm) n.sub.r n.sub.i n.sub.r n.sub.i n.sub.r n.sub.i 0.400 1.53 8.00E−03 1.53 5.00E−03 1.75 0.46 0.488 1.53 8.00E−03 1.53 5.00E−03 1.75 0.45 0.515 1.53 8.00E−03 1.53 5.00E−03 1.75 0.45 0.550 1.53 8.00E−03 1.53 6.00E−03 1.75 0.44 0.633 1.53 8.00E−03 1.53 6.00E−03 1.75 0.43 0.694 1.53 8.00E−03 1.53 7.00E−03 1.75 0.43 0.860 1.52 8.00E−03 1.52 1.20E−02 1.75 0.43 1.536 1.4 8.00E−03 1.51 2.30E−02 1.77 0.46 2.250 1.22 9.00E−03 1.42 1.00E−02 1.81 0.50 3.750 1.27 1.10E−02 1.452 4.00E−03 1.90 0.57 R.sub.mean (μm) 0.50 0.005 0.0118 R.sub.std (σ) 2.99 2.99 2.00

[0030] (b) Obtaining spectral AOD derivatives s2: Based on the spectral AODs of the various aerosol types, NDAIs are used for calculation to derive particle-size distributions and complex refractive indices from first- and second-order derivatives of the spectral AODs of the various aerosol types for examination and to construct intrinsic parameters of the various aerosol types for calculating main components of aerosols and mixing ratios thereof.

[0031] According to traditional formula, a first-order derivative of spectral AOD of gap between λ.sub.1 and λ.sub.2 is figured out as shown in Eq.(1), which reflects the particle-size distribution as covering the influence of AOD yet unable to single out particle size information. For removing the influence of AOD, the present invention improves the first-order derivative, as shown in Eq.(2), which is defined as a normalized aerosol index. With the building of the normalized aerosol index, the affect of AOD on the particle-size distribution is greatly reduced, where the particle-size distributions of the DS, BB (black carbon), and AP (sulfate and nitrate) are clearly distinguished.

[00001]$\begin{array}{cc}\frac{\partial \tau }{\partial \lambda }\approx \nabla {\tau }_{\left({\lambda }_1,{\lambda }_2\right)}=\frac{{\tau }_{{\lambda }_1}-{\tau }_{{\lambda }_2}}{\mathrm{\Delta \lambda }}={\tau }_{{\lambda }_2}\times \left(1-A^{\alpha }\right)\times B& \left(1\right)\end{array}$$\begin{array}{cc}{\mathrm{NDAI}}_{\left({\lambda }_1,{\lambda }_2\right)}\equiv \nabla {\tau }_{\left({\lambda }_1,{\lambda }_2\right)}/{\tau }_{{\lambda }_{\mathrm{ref}}},& \left(2\right)\end{array}$

[0032] where Δλ=λ.sub.2−λ.sub.1, A=λ.sub.2/λ.sub.1 and B=1/(λ.sub.2−λ.sub.1) are constants of specific bands; λ is a wavelength (μm); α is an Ångstrom exponent (AE, related to particle-size distribution); ∇τ.sub.(λ.sub.1.sub.,λ.sub.2.sub.) is a first-order derivative of spectral AOD of gap between λ.sub.1 and λ.sub.2; τ.sub.λ.sub.ref is a normalization reference of various AOD size; and NDAI.sub.(λ.sub.1.sub.,λ.sub.2.sub.) is a spectral derivative of gap between λ.sub.1 and λ.sub.2 as being normalized with τ.sub.λ.sub.ref.

[0033] The second-order derivative of AOD spectrum (as shown in Eq.(3)) is related to the imaginary number part of the refractive index. After being normalized (as shown in Eq.4)), features of and differences between the various aerosol types on scattering and absorption are described to distinguish and identify the various aerosol types.

[00002]$\begin{array}{cc}\frac{{\partial }^2\tau }{\partial {\lambda }^2}\cong {\nabla }^2{\tau }_{\left({\lambda }_1,{\lambda }_2,{\lambda }_3\right)}=\frac{\left({\tau }_{{\lambda }_1}-2{\tau }_{{\lambda }_2}+{\tau }_{{\lambda }_3}\right)}{\left({\lambda }_1-{\lambda }_2\right)\left({\lambda }_2-{\lambda }_3\right)}& \left(3\right)\end{array}$$\begin{array}{cc}\frac{{\partial }^2\tau }{\partial {\lambda }^2}/{\tau }_{{\lambda }_{\mathrm{ref}}}\cong {\nabla }^2{\tau }_{\left({\lambda }_1,{\lambda }_2,{\lambda }_3\right)}/{\tau }_{{\lambda }_{\mathrm{ref}}}=\frac{\left({\tau }_{{\lambda }_1}-2{\tau }_{{\lambda }_2}+{\tau }_{{\lambda }_3}\right)/{\tau }_{{\lambda }_{\mathrm{ref}}}}{\left({\lambda }_1-{\lambda }_2\right)\left({\lambda }_2-{\lambda }_3\right)},& \left(4\right)\end{array}$

[0034] where τ.sub.λ.sub.ref is selected from τ.sub.λ.sub.1, τ.sub.λ.sub.2, and τ.sub.λ.sub.3 to bind a dynamic range to the spectral derivative of a various aerosol type.

[0035] For distinguishing AODs in a mixed aerosol of two main components comprising A-type component and B-type component, the change of AOD depends on the AOD fraction (fAOD) for each type. As shown in Eq.(5),

Δτ.sub.(λ.sub.1.sub.,λ.sub.2.sub.).sup.ABmixed=f.sub.AOD.sup.AΔτ.sub.(λ.sub.1.sub.,λ.sub.2.sub.).sup.A+f.sub.AOD.sup.BΔτ.sub.(λ.sub.1.sub.,λ.sub.2.sub.).sup.B  (5),

[0036] where f.sub.AOD.sup.A and f.sub.AOD.sup.B are the fAOD(NDAI) in the spectrum (λ.sub.1,λ.sub.2) of the mixed aerosol comprising the A-type component and B-type component; and f.sub.AOD.sup.A+f.sub.AOD.sup.B=1. Based on Eq.(2), Eq.(5) further derives Eq.(6) based on the normalized aerosol index.

NDAI.sub.(λ.sub.1.sub.,λ.sub.2.sub.).sup.ABmixed=∇τ.sub.(λ.sub.1.sub.,λ.sub.2.sub.).sup.ABmixed/τ.sub.ref=f.sub.AOD.sup.ANDAI.sub.(λ.sub.1.sub.,λ.sub.2.sub.).sup.A+f.sub.AOD.sup.BNDAI.sub.(λ.sub.1.sub.,λ.sub.2.sub.).sup.B  (6)

[0037] Eq.(6) is the theoretical basis for calculating the fraction ratios of the main components in the mixed aerosol based on the normalized aerosol index. With the coordination of the optical intrinsic parameters of the various aerosol types built with the normalized aerosol indices, specific ratios of the various aerosol types are obtained as shown in Eq.(7).

[00003]$\begin{array}{cc}{f}_{\mathrm{AOD}}^A=\frac{{\mathrm{NDAI}}_{\left({\lambda }_1,{\lambda }_2\right)}^{\mathrm{mean}-A}-{\mathrm{NDAI}}_{\left({\lambda }_1,{\lambda }_2\right)}^{\mathrm{ABmixed}}}{{\mathrm{NDAI}}_{\left({\lambda }_1,{\lambda }_2\right)}^{\mathrm{mean}-A}-{\mathrm{NDAI}}_{\left({\lambda }_1,{\lambda }_2\right)}^{\mathrm{mean}-B}};f_{\mathrm{AOD}}^B=1-f_{\mathrm{AOD}}^A,& \left(7\right)\end{array}$

[0038] where NDAI.sub.(λ.sub.1.sub.,λ.sub.2.sub.).sup.mean-A and NDAI.sub.(λ.sub.1.sub.,λ.sub.2.sub.).sup.mean-B are intrinsic parameters of A-type aerosol and B-type aerosol; and NDAI.sub.(λ.sub.1.sub.,λ.sub.2.sub.).sup.ABmixed is an intrinsic parameter of A/B mixed-type aerosol. Thus, a novel method of spectral AOD derivatives for partitioning type and loading is obtained.

[0039] The following states-of-use are only examples to understand the details and contents of the present invention, but not to limit the scope of patent of the present invention.

[0040] For actual measurement, the main observation data are the spectral AOD data obtained through long-term observation of the Aerosol Robotic Network (AERONET) observation stations distributed globally, which comprises main source areas of DS, BB (black carbon) and AP (sulfate and nitrate). As shown in Table 2, the control data set and verification data set obtained from the AERONET are used to identify aerosol type.

TABLE-US-00002 TABLE 2 Control data set (used to construct NDAI) AP DS BB August- April-May March-May September Beijing Chiang Mai Beijing (39 N, 116 E) (18 N, 98 E)  (39 N, 116 E) 2001-2012 2006-2012 2001-2012 Dalanzadgad Mukdahan Hong Kong (43 N, 104 E) (16 N, 104 E) (22 N, 114 E) 1997-2012 2003-2010 2005-2012 Solar Village Pimai Taihu (24 N, 46 E)  (15 N, 102 E) (31 N, 120 E) 1998-2012 2003-2008 2005-2012 Tamanrasset Taipei (22 N, 5 E)   (25 N, 121 E) 2006-2012 2002-2012 Validation data set (used to evaluate NDAI) AP DS BB August- April-May March-May September (2014-2016) (2014-2016) (2014-2016) Beijing Chiang Mai Beijing (39 N, 116 E) (18 N, 98 E)  (39 N, 116 E) La Laguna Doi Ang Khang Durban UKZN (28 N, 16 W)  (19 N, 99 E)  (30 S, 31 E)  XuZhou Luang Namtha Hong Kong (34 N, 117 E) (20 N, 101 E) (22 N, 114 E) Zinder Airport Maeson La Laguna (14 N, 9 E)   (19 N, 99 E)  (28 N, 16 W)  Mongu Inn Mongu Inn (15 S, 23 E)  (15 S, 23 E)  NhaTrang Taihu (12 N, 109 E) (31 N, 120E) Omkoi Taipei (17 N, 98 E)  (25 N, 121 E) Silpakorn Univ XuZhou (13 N, 100 E) (34 N, 117 E) Ubon Ratchathani (15 N, 104 E) Vientiane (17 N, 102 E)

[Experiment Result and Analysis]

Theoretical Spectral AOD Derivatives

[0041] The spectral distributions of different AODs at specific wavelengths (0.44 μm, 0.47 μm, 0.55 μm, 0.66 μm, 0.675 μm, 0.87 μm, and 1.02 μm) are simulated based on the 6S experimental data set using various aerosols (i.e. Table 1); and a Bezier curve method is used in FIG. 2 to connect discrete points. The circular data points in the figure are the spectrum distributions of DS aerosol, which shows a flat trend with a small increase in wavelength; and the square data points and triangular data points are the spectral distributions of AP and BB aerosols, respectively, which tends to decrease continuously as similar to traditional research results. As AE indices show, the spectral gradient is mainly related to particle-size information, which shows that the radius of DS aerosol is much larger than the radius of AP and BB aerosols. Although similar particle sizes may increase the difficulty of distinguishing AP from BB, their spectral gradients are still slightly different. Besides, it is worth noting that, even for the same type, the spectral gradient of AOD may also vary with the optical depth (AOD.sub.(0.55 μm)=0.4, 0.8, 1.2, 1.6 and 2.0).

[0042] As shown in FIG. 2, it may be difficult to distinguish AP and BB aerosols in the zero-order spectrums by simply simulating the detailed changes in the spectral AODs. However, spectral derivative can promote the identification of subtle changes in AODs caused by different scattering and absorption. Hence, spectral derivative is related to particle-size distribution and complex refractive index for enhancing aerosol intrinsic feature. Based on the data shown in FIG. 2, FIG. 3A describes first-order derivatives of AODs and ∇τ.sub.(λ.sub.1.sub.,λ.sub.2.sub.), which are the spectrum pairs of 0.44-0.55 μm, 0.55-0.675 μm, 0.675-0.87 μm, and 0.87-1.02 μm from DS, AP, and BB aerosols, respectively. Different curves with the same shape of data points represent AOD.sub.(0.55 μm) values of 0.4, 0.8, 1.2, 1.6, and 2.0, respectively. The values of DS aerosol become almost flat along wavelengths (i.e., they all tend to be zero). Obviously, except DS (due to the flat distribution of spectral AODs), the difference of ∇τ.sub.(λ.sub.1.sub.,λ.sub.2.sub.) between AP and BB aerosols becomes more obvious in the shorter wavelength spectrums following AOD changes and deviates. For the second-order derivatives of AOD before normalization, each of the three groups of DS, AP, and BB aerosols has three continuous spectral AODs (0.44-0.55-0.675 μm, 0.55-0.675-0.87 μm, and 0.675-0.87-1.02 μm) as shown in FIG. 3C. Although the difference between AP and BB aerosols can be further magnified in terms of ∇.sup.2τ, the value of the second-order derivative still depends on AOD size, which is similar to ∇τ.sub.(λ.sub.1.sub.,λ.sub.2.sub.) as shown in FIG. 3A. On acquiring intrinsic features, both the first- and second-order derivatives are important for normalization as shown in FIG. 3B and FIG. 3D. After normalization using AOD.sub.(0.44 μm), each type of curve under different AOD begins to merge into its own intrinsic spectrum. The optical intrinsic parameters of DS, AP, and BB aerosols effectively eliminate AOD effect.

[0043] According to the above simulation results, the first- and second-order derivatives of the unnormalized AODs are still affected by AOD size. But, as shown in FIG. 3A and FIG. 3C, a closely-overlapped line is obtained after normalization; and the intrinsic parameters of the particle-size distributions and scattering/absorption of various aerosols are clearly shown in FIG. 3B and FIG. 3D.

[0044] When different types of aerosols are mixed, the optical features are usually diverse. Thus, the first- and second-order derivatives are used to discuss the dynamic range caused by the mixing effect of DS, AP, and BB aerosols. As shown in FIG. 4A to FIG. 4F., the pre- and post-normalized first- and second-order derivatives (FIG. 4A and FIG. 4C vs. FIG. 4B and FIG. 4D) are compared for AOD (τ.sub.0.44 μm) based on the spectral AODs at 0.44 μm, 0.675 μm, and 0.87 μm. In the figures, the triangle-symbolized BB, the square-symbolized AP, and the circle-symbolized DS aerosols have their area positions proved by the first- and second-order derivatives before and after normalization (FIG. 4E and FIG. 4F). Therein, the dotted lines between the symbols represent the dynamic ranges of mixed DS-AP, DS-BB and AP-BB aerosols.

[0045] As shown in the results, FIG. 4B, FIG. 4D, and FIG. 4F show that the normalized first- and second-order derivatives do not change as following the changes in AODs, which highlights the importance of the normalization in acquiring intrinsic parameters of particle-size distribution and scattering/absorption. The measurement of spectral AOD derivatives is shown exactly as the first derivatives of the spectral AODs (∇.sub.(0.44,0.675)), where both 6S-model simulation (shown in diagram (a) of FIG. 5) and the AERONET measurement (shown in diagram (b) of FIG. 5) reveal the changes of particle size (α) and AOD (τ). Furthermore, the slope of each aerosol type in the theoretical simulation is more consistent with the in-situ measurement, which strongly supports the NDAI method proposed in the present invention.

[0046] Based on the data set used in diagram (b) of FIG. 5, FIG. 6 processes the filtration of AOD.sub.(0.44 μm)>0.8 under Single Scattering Albedo (SSA) and normalizes the result through AOD.sub.(0.87 μm). Besides, the mapping (dashed line) for normal distribution and the average value (black dot) are also indicated.

[0047] In the above results shown in the figures, the result of the first-order derivatives (particle-size distribution) of the ground observation data (AERONET) before and after normalization (FIG. 5 v.s. FIG. 6) strongly supports the result of the theoretical simulation (shown in FIG. 4A-FIG. 4F), where various aerosol types can be clearly distinguished (shown in FIG. 6).

[0048] FIG. 7 shows that the second-order derivatives (normalized ∇2τ) and first-order derivatives (normalized VT) of two-component aerosol mixtures (DS-BB, DS-AP, and AP-BB, denoted with “6S”) are obtained to process simulation through the mixed volume/density weights for the 6S model in 0.01 steps from 0.00 to 1.00. The ground observation data (AERONET) and the average value and standard deviation of the spectral derivatives of DS, AP, BB, DS-AP, DS-BB, and AP-BB aerosols are configured together.

[0049] FIG. 7 shows the comparison between the first- and second-order derivatives of the data of the ground observation (AERONET) and the theoretical simulation, where a good consistency is shown in the figure and the normalized first- and second-order derivatives of the spectral AODs constructed by the present invention have considerable feasibility and application value in identifying and quantifying aerosol types.

[0050] Regarding practical applications, the present invention often applies to a variety of mixed aerosols, where the component proportions of three global representative aerosols are constructed through theory, comprising DS, BB (black carbon), and AP (sulfate and nitrate), for practical observation applications. As with the result shown in FIG. 8A and FIG. 8B, the first-order derivative (VT) is the optical features of particle-size distribution at 0.44 μm and 0.675 μm calculated for the mixed weights of volume and density in 0.01 steps from 0.00 to 1.00. In contrast, the second-order derivatives (∇.sup.2τ) are the calculated optical features related to refractive indices at 0.44 μm, 0.675 μm, and 0.87 μm. It means that, with the mixed volume/density weights of the three-component mixtures (i.e. mixture of DS, AP, and, BB aerosols) in the 6S model, the normalized second-order derivatives (∇.sup.2τ, related to refractive index at 0.44 μm, 0.675 μm, and 0.87 μm) of AOD.sub.(0.47 μm) is used to process simulation for the first-order derivatives (VT, related to particle-size distribution at 0.44 μm and 0.675 μm) in 0.01 steps from 0.00 to 1.00. In this way, with the calculated first-order (X-axis) and second-order (Y-axis) derivatives of the spectral AODs, corresponding aerosol types and their proportions can be found according to the database of FIG. 8A and FIG. 8.B.

[0051] It is still a challenge to quantify the compositions of aerosols (atmospheric particulate matter) with the data obtained from satellite telemetry or ground observation. Based on multi-spectral AODs, particle-size distributions and refractive indices are derived by normalizing first- and second-order derivatives for processing quantitative calibration of main components. At first, according to the optical feature parameters of various aerosol types (DS, BB, and AP), a radiation theory (6S model) is applied to simulate the multi-spectral optical depth for each density, including those of mixed types. The intrinsic parameters of the aerosol types are figured out with the normalized derivative aerosol index (NDAI) constructed according to the present invention. The apparent differences between the features of aerosols are used to figure out the main components of any specific aerosol and its mixing ratio. A simulation result of the NDAIs of the various aerosol types derived through applying the theory proposed in the present invention is in good agreement with the ground observation data of AERONET. It shows that the NDAI constructed according to the present invention is quite practicable in the quantitative calibration of the main components of atmospheric aerosols.

[0052] Hence, the main contributions of the present invention are as follows:

[0053] 1. First- and second-order derivatives obtained through multi-spectral AOD normalization is applied for identifying and quantitatively distinguishing aerosol types.

[0054] 2. The potential of satellite applications is obtained for providing global or regional distributions of aerosol types.

[0055] 3. Information of the temporal and spatial distribution of SSA having very scarce global observation data can be provided.

[0056] To sum up, the present invention is a method of spectral AOD derivatives for partitioning type and loading, where NDAI is used to integrate data of theoretical simulation and actual observation for examining various aerosol types with the relationships of particle-size distributions and complex refractive indices together with first- and second-order derivatives of spectral AODs and constructing optical intrinsic parameters of DS, BB, and AP; and, thus, each single type of aerosol is identified and main components of each mixed aerosol are quantitatively distinguished.

[0057] The preferred embodiment herein disclosed is not intended to unnecessarily limit the scope of the invention. Therefore, simple modifications or variations belonging to the equivalent of the scope of the claims and the instructions disclosed herein for a patent are all within the scope of the present invention.