Retrieving global sources of aerosols from MODIS observations by inverting GOCART model

International audienceKnowledge of the global distribution of tropospheric aerosols is important for studying the effects of aerosols on global climate. Chemical transport models rely on archived meteorological fields, accounting for aerosol sources, transport and removal processes can simulate the global distribution of atmospheric aerosols. However, the accuracy of global aerosol modeling is limited. Uncertainty in location and strength of aerosol emission sources is a major factor in limiting modeling accuracy. This paper describes an effort to develop an algorithm for retrieving global sources of aerosol from satellite observations by inverting the GOCART aerosol transport model. To optimize inversion algorithm performance, the inversion was formulated as a generalized multi-term least-squares-type fitting. This concept uses the principles of statistical optimization and unites diverse retrieval techniques into a single flexible inversion procedure. It is particularly useful for choosing and refining a priori constraints in the retrieval algorithm. For example, it is demonstrated that a priori limitations on the partial derivatives of retrieved characteristics, which are widely used in atmospheric remote sensing, can also be useful in inverse modeling for constraining time and space variability of the retrieved global aerosol emissions. The similarities and differences with the standard "Kalman filter" inverse modeling approach and the "Phillips-Tikhonov-Twomey" constrained inversion widely used in remote sensing are discussed. In order to retain the originally high space and time resolution of the global model in the inversion of a long record of observations, the algorithm was expressed using adjoint operators in a form convenient for practical development of the inversion from codes implementing forward model simulations. The inversion algorithm was implemented using the GOCART aerosol transport model. The numerical tests we conducted showed successful retrievals of global aerosol emissions with a 2°×2.5° resolution by inverting the GOCART output. For achieving satisfactory retrieval from satellite sensors such as MODIS, the emissions were assumed constant within the 24 h diurnal cycle and aerosol differences in chemical composition were neglected. Such additional assumptions were needed to constrain the inversion due to limitations of satellite temporal coverage and sensitivity to aerosol parameters. As a result, the algorithm was defined for the retrieval of emission sources of fine and coarse mode aerosols from the MODIS fine and coarse mode aerosol optical thickness data respectively. Numerical tests showed that such assumptions are justifiable, taking into account the accuracy of the model and observations and that it provides valuable retrievals of the location and the strength of the aerosol emissions. The algorithm was applied to MODIS observations during two weeks in August 2000. The global placement of fine mode aerosol sources retrieved from inverting MODIS observations was coherent with available independent knowledge. This was particularly encouraging since the inverse method did not use any a priori information about the sources and it was initialized under a "zero aerosol emission" assumption. The retrieval reproduced the instantaneous global MODIS observations with a standard deviation in fitting of aerosol optical thickness of ~0.04. The optical thickness during high aerosol loading events was reproduced with a standard deviation of ~48%. Applications of the algorithm for the retrieval of coarse mode aerosol emissions were less successful, mainly due to the currently existing lack of MODIS data over high reflectance desert dust sources. Possibilities for enhancing the global satellite data inversion by using diverse a priori constraints on the retrieval are demonstrated. The potential and limitations of applying our approach for the retrieval of global aerosol sources from aerosol remote sensing are discussed

Dust Optical Properties Over North Africa and Arabian Peninsula Derived from the AERONET Dataset

Dust optical properties over North Africa and the Arabian Peninsula are extracted from the quality assured multi-year datasets obtained at 14 sites of the Aerosol Robotic Network (AERONET). We select the data with (a) large aerosol optical depth (AOD >= 0.4 at 440 nm) and (b) small Angstrom exponent (A(sub ext)<= 0.2) for retaining high accuracy and reducing interference of non-dust aerosols. The result indicates that the major fraction of high aerosol optical depth days are dominated by dust over these sites even though it varies depending on location and time. We have found that the annual mean and standard deviation of single scattering albedo, asymmetry parameter, real refractive index, and imaginary refractive index for Saharan and Arabian desert dust is 0.944 +/- 0.005, 0.752 +/- 0.014, 1.498 +/- 0.032, and 0.0024 +/- 0.0034 at 550 nm wavelength, respectively. Dust aerosol selected by this method is less absorbing than the previously reported values over these sites. The weaker absorption of dust from this study is consistent with the studies using remote sensing techniques from satellite. These results can help to constrain uncertainties in estimating global dust shortwave radiative forcing

Retrieving global aerosol sources from satellites using inverse modeling

International audienceUnderstanding aerosol effects on global climate requires knowing the global distribution of tropospheric aerosols. By accounting for aerosol sources, transports, and removal processes, chemical transport models simulate the global aerosol distribution using archived meteorological fields. We develop an algorithm for retrieving global aerosol sources from satellite observations of aerosol distribution by inverting the GOCART aerosol transport model. The inversion is based on a generalized, multi-term least-squares-type fitting, allowing flexible selection and refinement of a priori algorithm constraints. For example, limitations can be placed on retrieved quantity partial derivatives, to constrain global aerosol emission space and time variability in the results. Similarities and differences between commonly used inverse modeling and remote sensing techniques are analyzed. To retain the high space and time resolution of long-period, global observational records, the algorithm is expressed using adjoint operators. Successful global aerosol emission retrievals at 2°×2.5 resolution were obtained by inverting GOCART aerosol transport model output, assuming constant emissions over the diurnal cycle, and neglecting aerosol compositional differences. In addition, fine and coarse mode aerosol emission sources were inverted separately from MODIS fine and coarse mode aerosol optical thickness data, respectively. These assumptions are justified, based on observational coverage and accuracy limitations, producing valuable aerosol source locations and emission strengths. From two weeks of daily MODIS observations during August 2000, the global placement of fine mode aerosol sources agreed with available independent knowledge, even though the inverse method did not use any a priori information about aerosol sources, and was initialized with a "zero aerosol emission" assumption. Retrieving coarse mode aerosol emissions was less successful, mainly because MODIS aerosol data over highly reflecting desert dust sources is lacking. The broader implications of applying our approach are also discussed

Statistically Optimized Inversion Algorithm for Enhanced Retrieval of Aerosol Properties from Spectral Multi-Angle Polarimetric Satellite Observations

The proposed development is an attempt to enhance aerosol retrieval by emphasizing statistical optimization in inversion of advanced satellite observations. This optimization concept improves retrieval accuracy relying on the knowledge of measurement error distribution. Efficient application of such optimization requires pronounced data redundancy (excess of the measurements number over number of unknowns) that is not common in satellite observations. The POLDER imager on board the PARASOL microsatellite registers spectral polarimetric characteristics of the reflected atmospheric radiation at up to 16 viewing directions over each observed pixel. The completeness of such observations is notably higher than for most currently operating passive satellite aerosol sensors. This provides an opportunity for profound utilization of statistical optimization principles in satellite data inversion. The proposed retrieval scheme is designed as statistically optimized multi-variable fitting of all available angular observations obtained by the POLDER sensor in the window spectral channels where absorption by gas is minimal. The total number of such observations by PARASOL always exceeds a hundred over each pixel and the statistical optimization concept promises to be efficient even if the algorithm retrieves several tens of aerosol parameters. Based on this idea, the proposed algorithm uses a large number of unknowns and is aimed at retrieval of extended set of parameters affecting measured radiation

An Analysis of AERONET Aerosol Absorption Properties and Classifications Representative of Aerosol Source Regions

Partitioning of mineral dust, pollution, smoke, and mixtures using remote sensing techniques can help improve accuracy of satellite retrievals and assessments of the aerosol radiative impact on climate. Spectral aerosol optical depth (tau) and single scattering albedo (omega (sub 0) ) from Aerosol Robotic Network (AERONET) measurements are used to form absorption [i.e., omega (sub 0) and absorption Angstrom exponent (alpha(sub abs))] and size [i.e., extinction Angstrom exponent (alpha(sub ext)) and fine mode fraction of tau] relationships to infer dominant aerosol types. Using the long-term AERONET data set (1999-2010), 19 sites are grouped by aerosol type based on known source regions to: (1) determine the average omega (sub 0) and alpha(sub abs) at each site (expanding upon previous work); (2) perform a sensitivity study on alpha(sub abs) by varying the spectral omega (sub 0); and (3) test the ability of each absorption and size relationship to distinguish aerosol types. The spectral omega (sub 0) averages indicate slightly more aerosol absorption (i.e., a 0.0 < delta omega (sub 0) <= 0.02 decrease) than in previous work and optical mixtures of pollution and smoke with dust show stronger absorption than dust alone. Frequency distributions of alpha(sub abs) show significant overlap among aerosol type categories and at least 10% of the alpha(sub abs) retrievals in each category are below 1.0. Perturbing the spectral omega (sub 0) by +/- 0.03 induces significant alpha(sub abs) changes from the unperturbed value by at least approx. +/- 0.6 for Dust, approx. +/-0.2 for Mixed, and approx. +/-0.1 for Urban/Industrial and Biomass Burning. The omega (sub 0)440nm and alpha(sub ext) 440-870nm relationship shows the best separation among aerosol type clusters, providing a simple technique for determining aerosol type from surface- and future space-based instrumentation

Reduction of Aerosol Absorption in Beijing Since 2007 from MODIS and AERONET

An analysis of the time series of MODIS-based and AERONET aerosol records over Beijing reveals two distinct periods, before and after 2007. The MODIS data from both the Terra and Aqua satellites were processed with the new Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm. A comparison of MAIAC and AERONET AOT shows that whereas MAIAC consistently underestimated peak AOT values by 10-20% in the prior period, the bias mostly disappears after mid-2007. Independent analysis of the AERONET dataset reveals little or no change in the effective radii of the fine and coarse fractions and of the Angstrom exponent. At the same time, it shows an increasing trend in the single scattering albedo, by approx.0.02 in 9 years. As MAIAC was using the same aerosol model for the entire 2000-2010 period, the decrease in AOT bias after 2007 can be explained only by a corresponding decrease of aerosol absorption caused by a reduction in local black carbon emissions. The observed changes correlate in time with the Chinese government's broad measures to improve air quality in Beijing during preparations for the Summer Olympics of 2008

Column aerosol characterization in a semi‐arid region around Marrakech during the WATERMED 2003 campaign

In the framework of the WATERMED (WATer use efficiency in natural vegetation and agricultural areas by Remote sensing in the MEDiterranean basin) project, an experimental field campaign was carried out in a semi‐arid region near Marrakech, Morocco, during March 2003. This work focused on the columnar aerosol characterization from spectroradiometric and photometric measurements of direct solar irradiance and sky radiance at ground level. The results show a high dependence of the aerosol optical properties on the air masses present in the area. Two periods with different aerosol loads were observed in the campaign. In the first, Atlantic‐origin aerosols showed a marine particle accumulation mode combined with a dominant mode corresponding to coarse particles of local origin. In the second period the synoptic situation changed towards a Saharan+Atlantic origin, the Atlantic influence persisting at altitude while at lower levels the presence of air masses proceeding from the Saharan desert was observed. This dry air carried a large quantity of mineral particles that increased the area's turbidity, leading to a substantial increase in the volume of the whole distribution, a displacement of the accumulation mode radius towards higher values and an increasingly evident spread of the coarse particle mode, although the radius remained the same as that in the first period of the campaign

Climatological aspects of the optical properties of fine/coarse mode aerosol mixtures

Aerosol mixtures composed of coarse mode desert dust combined with fine mode combustion generated aerosols (from fossil fuel and biomass burning sources) were investigated at three locations that are in and/or downwind of major global aerosol emission source regions. Multiyear monitoring data at Aerosol Robotic Network sites in Beijing (central eastern China), Kanpur (Indo-Gangetic Plain, northern India), and Ilorin (Nigeria, Sudanian zone of West Africa) were utilized to study the climatological characteristics of aerosol optical properties. Multiyear climatological averages of spectral single scattering albedo (SSA) versus fine mode fraction (FMF) of aerosol optical depth at 675 nm at all three sites exhibited relatively linear trends up to similar to 50% FMF. This suggests the possibility that external linear mixing of both fine and coarse mode components (weighted by FMF) dominates the SSA variation, where the SSA of each component remains relatively constant for this range of FMF only. However, it is likely that a combination of other factors is also involved in determining the dynamics of SSA as a function of FMF, such as fine mode particles adhering to coarse mode dust. The spectral variation of the climatological averaged aerosol absorption optical depth (AAOD) was nearly linear in logarithmic coordinates over the wavelength range of 440-870 nm for both the Kanpur and Ilorin sites. However, at two sites in China (Beijing and Xianghe), a distinct nonlinearity in spectral AAOD in logarithmic space was observed, suggesting the possibility of anomalously strong absorption in coarse mode aerosols increasing the 870 nm AAOD

Optical Properties of Boreal Region Biomass Burning Aerosols in Central Alaska and Seasonal Variation of Aerosol Optical Depth at an Arctic Coastal Site

Long-term monitoring of aerosol optical properties at a boreal forest AERONET site in interior Alaska was performed from 1994 through 2008 (excluding winter). Large interannual variability was observed, with some years showing near background aerosol optical depth (AOD) levels (<0.1 at 500 nm) while 2004 and 2005 had August monthly means similar in magnitude to peak months at major tropical biomass burning regions. Single scattering albedo (omega (sub 0); 440 nm) at the boreal forest site ranged from approximately 0.91 to 0.99 with an average of approximately 0.96 for observations in 2004 and 2005. This suggests a significant amount of smoldering combustion of woody fuels and peat/soil layers that would result in relatively low black carbon mass fractions for smoke particles. The fine mode particle volume median radius during the heavy burning years was quite large, averaging approximately 0.17 micron at AOD(440 nm) = 0.1 and increasing to approximately 0.25 micron at AOD(440 nm) = 3.0. This large particle size for biomass burning aerosols results in a greater relative scattering component of extinction and, therefore, also contributes to higher omega (sub 0). Additionally, monitoring at an Arctic Ocean coastal site (Barrow, Alaska) suggested transport of smoke to the Arctic in summer resulting in individual events with much higher AOD than that occurring during typical spring Arctic haze. However, the springtime mean AOD(500 nm) is higher during late March through late May (approximately 0.150) than during summer months (approximately 0.085) at Barrow partly due to very few days with low background AOD levels in spring compared with many days with clean background conditions in summer