Synergetic retrieval of terrestrial AOD from MODIS images of twin satellites Terra and Aqua

Aerosol optical depth (AOD) is one of the most important indicators of atmospheric pollution. It can be retrieved from satellite imagery using several established methods, such as the dark dense vegetation method and the deep blue algorithm. All of these methods require estimation of surface reflectance prior to retrieval, and are applicable to a certain pre-designated type of surface cover. Such limitations can be overcome by using a synergetic method of retrieval proposed in this study. This innovative method is based on the fact that the ratio K of surface reflectance at different angles/geometries is independent of wavelength as reported by Flowerdew and Haigh (1995). An atmospheric radiative transfer model was then established and resolved with the assistance of the ratio K obtained from two Moderate Resolution Imaging Spectroradiometer (MODIS) spectral bands acquired from the twin satellites of Terra and Aqua whose overpass is separated by three hours. This synergetic method of retrieval was tested with 20 pairs of MODIS images. The retrieved AOD was validated against the ground observed AOD at the Taihu station of the AErosol RObotic NETwork (AERONET). It is found that they are correlated with the observations at a coefficient of 0.828 at 0.47 μm and 0.921 at 0.66 μm wavelengths. The retrieved AOD has a mean relative error of 25.47% at 0.47 μm and 24.3% at 0.66 μm. Of the 20 samples, 15 and 17 fall within two standard error of the line based observed AOD data on the ground at the 0.47 μm and 0.66 μm, respectively. These results indicate that this synergetic method can be used to reliably retrieve AOD from the twin satellites MODIS images, namely Terra and Aqua. It is not necessary to determine surface reflectance first.     

An all-sky survey at 230 GHz by MLS on Aura

The Microwave Limb Sounder (MLS) instrument is a small satellite-borne radio telescope. Its purpose is to make limb-scanning measurements of atmospheric composition. One of the gases to which it is sensitive is carbon monoxide (CO), detected via the J = 2 → 1 rotational transition at 230 GHz. CO is present in molecular gas clouds in the Milky Way. Although it was not designed for the purpose, MLS can detect emissions from galactic CO, allowing a map of the 230 GHz radio sky to be constructed. We report the MLS measurements of galactic radio emission and discuss their effect on the atmospheric mission of MLS. The region of the Milky Way with emissions strong enough to significantly affect MLS observations of atmospheric CO is identified. Ground-based radio astronomers have been mapping the sky using CO emission for many years. However, the MLS data are the first such survey to be carried out from space. The MLS survey covers a larger area of the sky than any other 230 GHz survey, but no previously unknown gas clouds are observed.     

Evaluation of MODIS GPP over a complex ecosystem in East Asia: A case study at Gwangneung flux tower in Korea

Moderate Resolution Imaging Radiometer (MODIS) gross primary productivity (GPP) has been used widely to study the global carbon cycle associated with terrestrial ecosystems. The retrieval of the current MODIS productivity with a 1 × 1 km² resolution has limitations when presenting subgrid scale processes in terrestrial ecosystems, specifically when forests are located in mountainous areas, and shows heterogeneity in vegetation type due to intensive land use. Here, we evaluate MODIS GPP (MOD17) at Gwangneung deciduous forest KoFlux tower (deciduous forest; GDK) for 2006–2010 in Korea, where the forests comprise heterogeneous vegetation cover over complex terrain. The monthly MODIS GPP data overestimated the GDK measurements in a range of +15% to +34% and was fairly well correlated (R = 0.88) with the monthly variability at GDK during the growing season. In addition, the MODIS data partly represented the sharp GPP reduction during the Asian summer monsoon (June–September) when intensive precipitation considerably reduces solar radiation and disturbs the forest ecosystem. To examine the influence of subgrid scale heterogeneity on GPP estimates over the MODIS scale, the individual vegetation type and its area within a corresponding MODIS pixel were identified using a national forest type map (∼71-m spatial resolution), and the annual GPP in the same area as the MODIS pixel was estimated. This resulted in a slight reduction in the positive MODIS bias by ∼10%, with a high degree of uncertainty in the estimation. The MODIS discrepancy for GDK suggests further investigation is necessary to determine the MODIS errors associated with the site-specific aerodynamic and hydrological characteristics that are closely related to the mountainous topography. The accuracy of meteorological variables and the impact of the very cloudy conditions in East Asia also need to be assessed.     

Empirical model of cosmic ray spectrum in energy interval 1 MeV–100 GeV during 11-year solar cycle

The galactic cosmic rays (GCR) are the main ionization source at altitude of ∼3–35 km in the atmosphere. For high latitude anomalous cosmic ray (ACR) component has also a significant influence on the atmospheric ionization. We propose an empirical model for differential spectra D(E) of galactic and anomalous cosmic rays in energy interval 1 MeV–100 GeV during solar cycle. In the model data are used which cover three solar cycles: 20, 22 and 23. The LEAP87, IMAX92, CAPRICE94, AMS98 and BESS experimental spectra for protons and alpha particles are fitted to the proposed empirical model. The modulated GCR differential spectra are compared with force-field approximation to the one-dimensional transport equation and with solutions of two-dimensional cosmic ray transport equation. For experimental spectra, the calculation of the model parameters is performed by Levenberg–Marquardt algorithm, applied to the special case of least squares. Algorithm that combines the rapid local convergence of Newton–Raphson method with globally convergent method for non-linear systems of equations is applied for theoretically obtained differential spectra. The described programmes are realized in algorithmic language C++. The proposed model gives practical possibility for investigation of experimental data from measurements of galactic cosmic rays and their anomalous component.     

A simplified data assimilation method for reconstructing time-series MODIS NDVI data

The Normalized Difference Vegetation Index (NDVI) is an important vegetation index, widely applied in research on global environmental and climatic change. However, noise induced by cloud contamination and atmospheric variability impedes the analysis and application of NDVI data. In this work, a simplified data assimilation method is proposed to reconstruct high-quality time-series MODIS NDVI data. We extracted 16-Day L3 Global 1 km SIN Grid NDVI data sets for western China from MODIS vegetation index (VI) products (MOD13A2) for the period 2003–2006. NDVI data in the first three years (2003–2005) were used to generate the background field of NDVI based on a simple three-point smoothing technique, which captures annual features of vegetation change. NDVI data for 2006 were used to test our method. For every time step, the quality assurance (QA) flags of the MODIS VI products were adopted to empirically determine the weight between the background field and NDVI observations. Ultimately, more reliable NDVI data can be produced. The results indicate that the newly developed method is robust and effective in reconstructing high-quality MODIS NDVI time-series.     

Experiment for mapping land cover and it’s change in southeastern Sri Lanka utilizing 250 m resolution MODIS imageries

Utilizing freely available MODIS NDVI and Natural color imageries of 250 m spatial resolution produced by NASA, an experiment was made to map land-cover and its change with an emphasis on vegetation cover in southeastern Sri Lanka, which plays a vital role for control of green house gas. For the change detection purpose, 1987 land cover map made by present authors from Landsat MSS image and extensive ground truth survey data was used as the base map. The result of the experiment shows that MODIS data are useful to make a land cover map of 250 m spatial resolution for tropical areas with high cloud coverage like Sri Lanka. It was found that the forest cover decrease amounted as large as 21% in 19 years time span in southeastern Sri Lanka, the prominent forest region of the country. On the other hand homestead/vegetation and mixed vegetation/scrub dominant categories increased by 13.7% and 7.1%, respectively. These changes are considered due to a large clearance of forest areas for agriculture and building houses to accommodate increasing inhabitants.     

Ground-based millimeter-wave observations of water vapor emission (183 GHz) at Atacama,Chile

We report the first results of ground-based millimeter-wave measurements of 183 GHz atmospheric water vapor spectra from Atacama highland (4800 m alt.), Chile. The measurements were carried out in December 2005 by using a spectroscopic radiometer equipped with a superconductive heterodyne receiver. A conspicuous H₂O spectrum at 183 GHz was detected with an integration time of only 1.5 min, and this is the first high frequency-resolution H₂O spectrum at 183 GHz obtained in the southern subtropical region. The vertical profile of H₂O volume mixing ratio between 40 and 64 km were retrieved from the spectrum by using the modified optimal estimation method.     

The thermospheric composition different responses to geomagnetic storm in the winter and summer hemisphere measured by “SZ” Atmospheric Composition Detectors

Three “SZ” Atmospheric Composition Detectors (ACDs) on board spacecraft “SZ-2”, “SZ-3” and “SZ-4” were launched on 10th January 2001, 26th March 2002 and 31st December 2002 separately. A large quantity of thermospheric composition data at the orbital altitude ranging from 330 to 362 km were collected from the in-situ measurement of ACDs. The spacecrafts’ lifetime was just in the second peak period of the 23rd solar cycle which includes two peaks and the solar activity value F_10.7 was from 89 to 228. During this period, several intense geomagnetic disturbances happened.     

An empirical model of electron density in low latitude at 600 km obtained by Hinotori satellite

An empirical model of electron density (Ne) was constructed by using the data obtained with an impedance probe on board Japanese Hinotori satellite. The satellite was in circular orbit of the height of 600 km with the inclination of 31 degrees from February 1981 to June 1982. The constructed model gives Ne at any local time with the time resolution of 90 min and between −25 and 25 degrees in magnetic latitude with its resolution of 5 degrees in the range of F_10.7 from 150 to 250 under the condition of K_p < 4. Spline interpolations are applied to the functions of day of year, geomagnetic latitude and solar local time, and linear interpolation is applied to the function of F_10.7. Longitude dependence of Ne is not taken into account. Our density model can reproduce solar local time variation of electron density at 600 km altitude better than current International Reference Ionosphere (IRI2001) model which overestimates Ne in night time and underestimates Ne in day time. Our density model together with electron temperature model which has been constructed before will enable more understanding of upper ionospheric phenomenon in the equatorial region.     

Low latitude nighttime ionospheric vertical E × B drifts at African region

The nighttime vertical E × B drifts velocities of the F2-region were inferred from the hourly hmF2 values obtained from ionosonde data over an African equatorial station, Ilorin (8.50^oN, 4.68^oE; dip lat. 2.95^o) during period of low solar activity. For each season, the plasma drift Vz is characterized by an evening upward enhancement, then by a downward reversal at 1900 LT till around 0000 LT, except for June solstice. This was explained using the Rayleigh–Taylor (R-T) instability mechanism. The occasional drift differences in Vz obtained by inferred and direct measurement over Ilorin and Jicamarca, respectively are reflective of the importance of chemistry and divergent transport system due to both the E region electric and magnetic fields instead of simple motions. The pre-reversal enhancement (PRE) magnitude is higher during the equinoctial months than the solsticial months over Jicamarca, highest during December solstice and the equinoctial months over Ilorin, suggesting the dominance of higher E × B fountain during equinoxes at both stations. The lowest PRE magnitude was in June solstice. The appearance of post-noon peak in NmF2 around 1700 LT is highest during the equinoctial months and lowest during the solsticial period. A general sharp drop in NmF2 around 1800 LT is distinct immediately after sunset, lowest during June solstice and highest in March equinox. Our result suggests that between 0930 and 2100 LT, the general theory that vertical drifts obtained by digisonde measurements only match the E × B drift if the F layer is higher than 300 km is reliable, but does not hold for the nighttime period of 2200–0600 LT under condition of solar minima. Hence, the condition may not be sufficient for the representation of vertical plasma drift at nighttime during solar minima. This assertion may still be tentative, as more equatorial stations needed to be studied for better confirmation.     

Tidal associated temperature disturbances observed at the middle atmosphere (30–65 km) by a Rayleigh lidar at 23°S

A lidar has been operated in São José dos Campos, Brazil (23.2°S, 45.8°W) since 1972, mainly dedicated to the study of mesospheric sodium at the 589 nm resonant line. The molecular Rayleigh scattering can also be used provided we limit the height to ∼75 km where the sodium scattering begins. Nevertheless, the weak signal obtained only permits the determination of density and temperature profiles by accumulating a large number of shots giving only nocturnal average profiles. Temporal variations in density and temperature on the scale of hours can however, be obtained by performing a superposed epoch analysis for a given time interval and covering a period of several days. In this way we obtained hourly mean profiles grouped by months, seasons and overall, with data acquired from 1993 to 2004. The difference between the hourly temperatures and the nocturnal means shows for some months, with enough data coverage, downward propagating structures that apparently have tidal origin. The seasonal averages show a recurrent feature with high temperatures before and low temperatures after midnight above 50 km. Some similarity is found with the GSWM model, but the observed temperature amplitudes are twice of that for the model.     

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.     

Behavior of the atomic oxygen 557.7 nm atmospheric emission in the solar cycle 23

We present the results of nightglow observation of the atomic oxygen 557.7 nm line emission in the solar cycle 23. We use the experimental data obtained at Geophysical observatory near Irkutsk (52°N, 103°E), Russia, for the 1997–2006 period. The 557.7 nm emission observations data are compared with atmospheric and solar parameters. We note a difference in correlation coefficients between the 557.7 nm emission intensity and the solar activity indices in different phases of the solar cycle. Airglow observation results are compared with the observational data obtained by other authors.     

Spiral eddies in the Baltic,Black and Caspian seas as seen by satellite radar data

Despite spiral eddies were first seen on the sea surface more than 40 years ago, there is still a lot of uncertainty concerning these eddies. The present paper is aimed to provide the comprehensive results on the occurrence and statistics of small-scale eddies in the three inner seas (the Baltic, Black and Caspian) using satellite synthetic aperture radar (SAR) images. The dataset used includes over 2000 medium resolution Envisat ASAR and ERS-2 SAR images obtained in 2009–2010 in the different parts of the seas mentioned. As a result of the analysis performed ∼14,000 vortical structures were detected. 71% of them were visualized due to surfactant films (“black” eddies), while 29% due to wave/current interactions (“white” eddies). Practically all the eddies detected were cyclonically rotating. Their diameter was within 1–20 km. Characteristic size of the “black” eddies in all the basins was discovered to be less than that of the “white” eddies. Characteristic eddy size for the Baltic, Black and Caspian seas proved to be strictly proportional to the values of the baroclinic Rossby radius of deformation typical for these basins. The “black” eddies did not demonstrate a significant connection with the basin- and meso-scale surface circulation of the seas. Most of the “white” eddies detected were attributed to the zones with the most intense drift currents, i.e. those along the western boundaries and (in the Baltic Sea only) in the elongated parts of the basin.     

Investigation of precursory phenomena in the ionosphere,atmosphere and groundwater before large earthquakes of M > 6.5

Precursory phenomena in the ionosphere, atmosphere and groundwater before large earthquakes (M > 6.5) are extensively investigated toward the earthquake prediction. Upward tornado type seismic clouds occurred near the epicenter associated with strong LF-VLF radio noises from lightning discharges in the evening of January 9, 1995 [Yamada, T., Oike, K. On the increase of electromagnetic noises before and after the 1995 Hyogo-Ken Nanbu earthquake. In: Hayakawa M. (Ed.), Atmospheric and Ionospheric Electromagnetic Phenomena Associated with Earthquakes. TERRAPUB, Tokyo, pp. 417–427, 1999] and anomalous foEs increases up to 10 MHz were detected at Shigaraki, 90 km of the epicenter and at Kokubunji, 500 km east of the epicenter [Ondoh, T. Anomalous sporadic-E layers observed before M7.2 Hyogo-ken Nanbu earthquake; Terrestrial gas emanation model. Adv. Polar Upper Atmos. Res. 17, 96–108, 2003; Ondoh, T. Anomalous sporadic-E ionization before a great earthquake, Adv. Space Research 34, 1830–1835, 2004] associated with strong ELF noises from lightning discharges in the daytime on January 15, 1995 [Hata, M., Fujii, T., Takumi, I. EM precursor of large-scale earthquakes in Japan, in: Abstracts of International Workshop on Seismo Electromagnetics (IWSE 2005), Univ. Electro-Communications, Chofu, Tokyo, Japan, March 15–17, pp. 182–186, 2005] before the M7.2 Hyogoken–Nanbu earthquake of January 17, 1995. The anomalous foEs increases occurred at epicentral distances within 500 km that are the same as those of the terrestrial gas emanations along active faults before large earthquakes [King, C.-Y. Gas geochemistry applied to earthquake prediction: An overview. J. Geophys. Res. 91(B12), 12269–12281, 1986]. The anomalous foEs increases seem to be a seismic precursor because geomagnetic and solar conditions were very quiet all day on January 15,1995 and the normal foEs in Japanese winter is below 6 MHz. No significant pre-seismic geomagnetic field variation was detected at epicentral distance of 100 km before this earthquake [Ondoh, T., Hayakawa, M. Anomalous occurrence of sporadic-E layers before the Hyogoken–Nanbu earthquake, M7.2 of January 17, 1995. In: Hayakawa, M. (Ed.), Atmospheric and Ionospheric Electromagnetic Phenomena Associated with Earthquakes, TERRAPUB, Tokyo, pp. 629–639, 1999; Ondoh, T., Hayakawa, M. Seismo discharge model of anomalous sporadic E ionization before great earthquakes. In: Hayakawa, M., O.A. Molchanov, (Eds.), Seismo Electromagnetics: Lithosphere–Atmosphere–Ionosphere Couplings, TERRAPUB, Tokyo, pp. 385–390, 2002; Ondoh. T., Hayakawa, M. Synthetic study of precursory phenomena of the M7.2 Hyogo-ken Nanbu earthquake. Phys. Chem. Earth 31, 378–388, 2006]. The foF2 decrease and h’F increase occurred before the M7.8 Hokkaido Nansei-Oki earthquake of July 12,1993 in a geomagnetic quiet period [Ondoh, T. Ionospheric disturbances associated with great earthquake of Hokkaido southwest coast, Japan of July 12, 1993. Phys. Earth Planet. Interiors. 105, 261–269, 1998; Ondoh, T. Seismo ionospheric phenomena. Adv. Space Res. 26, 8, 1267–1272, 2000]. Characteristic phase changes at terminator times of Omega 10.2 kHz waves passing 70 km of the epicenter extended toward darker local times by 1 h for 3 days before this earthquake due to lowering of the wave reflection height or ion density increases in the D region [Hayakawa, M., Molchanov, O. A., Ondoh, T., Kawai, E. The precursory signature effect of the Kobe earthquake on VLF subionospheric signals. J. Commun. Res. La., 43, 00. 169–180, 1996]. The radon concentration in the atmosphere over Ashiya fault, Kobe [Yasuoka, Y., Shinogi, M. Anomaly in atmospheric radon concentration: a possible precursor of the 1995 Kobe, Japan, earthquake. Health Phys. 72(5), 759–761, 1997] and in the groundwater at 17 m well in Nishinomiya, Japan [Igarashi, G., Saeki, S., Takahata, N., Sumikawa, K., Tasaki, S., Sasaki, Y., Takahashi, M., Sano, Y. Ground-water radon anomaly before the Kobe earthquake in Japan. Science 269, 60–61, 1995] had gradually increased since 2 months before the M7.2 earthquake, increased suddenly in December 1994, and rapidly returned to the normal low level of October, 1994 [Yasuoka, Y., Shinogi, M. 1997. Anomaly in atmospheric radon concentration: a possible precursor of the 1995 Kobe. Japan, earthquake. Health Phys. 72(5), 759–761.]. Radon concentration changes in the groundwater before the M 7.0 Izu-Oshima-kinkai earthquake, Japan on January 14, 1978 [Wakita, H., Nakamura, Y., Notsu, K., Noguchi, M., Asada, T. 1980. Radon anomaly: a possible precursor of the 1978 Izu-Oshima-kinkai earthquake. Science 207, 882–883] and the M6.8 Chengkung earthquake, Taiwan on December 10, 2003 [Kuo, T., Fan, K., Chen, W., Kuochen, H., Han, Y., Wang, C., Chang, T., Lee, Y. Radon anomaly at the Antung Hot Spring before the Taiwan M6.8 Chengkung earthquake. Proceedings, Thirty-First Workshop on Geothermal Reservoir Engineering, Stanford University, Stanford, California, January 30–February 1, 2006, SGP-TR-179, 2006] are also investigated to find common features of the groundwater radon concentration changes before large earthquakes (M > 6.5) in comparison with those before the M7.2 Hyogoken–Nanbu earthquake. Groundwater radon concentrations before the 3 large earthquakes had shown common characteristic changes of gradually initial ones from the normal level since about 2 months before the earthquake onsets, rapid decreases down to the minimum, and quick increases up to the maximum at 7–20 days before the earthquake onsets, respectively. These are very useful characteristics of pre-seismic radon anomaly for the earthquake prediction or warning. Promising observations toward the earthquake prediction are also discussed.     

Cosmic-ray electron spectrum above 100 GeV from PPB-BETS experiment in Antarctica

Cosmic-ray electrons have been observed in the energy region from 10 GeV to 1 TeV with the PPB-BETS by a long duration balloon flight using a Polar Patrol Balloon (PPB) in Antarctica. The observation was carried out for 13 days at an average altitude of 35 km in January 2004. The PPB-BETS detector is an imaging calorimeter composed of scintillating-fiber belts and plastic scintillators inserted between lead plates. In the study of cosmic-ray electrons, there have been some suggestions that high-energy electrons above 100 GeV are a powerful probe to identify nearby cosmic-ray sources and search for particle dark matter. In this paper, we present the energy spectrum of cosmic-ray electrons in the energy range from 100 GeV to 1 TeV at the top of atmosphere, and compare our spectrum with the results from other experiments.     

The OIV 1407.3 Å /1401.1 Å emission-line ratio in a plasma

Line ratio of O IV 1407.3 Å/1401.1 Å is calculated using mostly our own atomic and collisional data.     

Cancer risk above 1 Gy and the impact for space radiation protection

Analyses of the epidemiological data on the Japanese A-bomb survivors, who were exposed to γ-rays and neutrons, provide most current information on the dose–response of radiation-induced cancer. Since the dose span of main interest is usually between 0 and 1 Gy, for radiation protection purposes, the analysis of the A-bomb survivors is often focused on this range. However, estimates of cancer risk for doses larger than 1 Gy are becoming more important for long-term manned space missions. Therefore in this work, emphasis is placed on doses larger than 1 Gy with respect to radiation-induced solid cancer and leukemia mortality. The present analysis of the A-bomb survivors data was extended by including two extra high-dose categories and applying organ-averaged dose instead of the colon-weighted dose. In addition, since there are some recent indications for a high neutron dose contribution, the data were fitted separately for three different values for the relative biological effectiveness (RBE) of the neutrons (10, 35 and 100) and a variable RBE as a function of dose. The data were fitted using a linear and a linear-exponential dose–response relationship using a dose and dose-rate effectiveness factor (DDREF) of both one and two. The work presented here implies that the use of organ-averaged dose, a dose-dependent neutron RBE and the bending-over of the dose–response relationship for radiation-induced cancer could result in a reduction of radiation risk by around 50% above 1 Gy. This could impact radiation risk estimates for space crews on long-term mission above 500 days who might be exposed to doses above 1 Gy. The consequence of using a DDREF of one instead of two increases cancer risk by about 40% and would therefore balance the risk decrease described above.     

Is there an abnormal enhancement of atmospheric aerosol before the 2008 Wenchuan earthquake?

There are extensive reports of ionospheric disturbances before the great 2008 Wenchuan earthquake, which are possibly explained by seismogenic electric field hypotheses linked with the aerosols injected in atmosphere. This paper attempts to investigate the possible change of atmospheric aerosol optical depth (AOD) associated with this earthquake by using MODIS data from both Terra and Aqua satellites. The result shows a clear enhancement of AOD along the Longmenshan faults 7 days before the quake, which is 1 day and 4 days earlier than the reported negative and positive ionospheric disturbances, respectively, and is 1 day earlier than or quasi-synchronism with other reported atmospheric anomalies including air temperature, outgoing longwave radiation and relative humidity. Particularly, the spatial distribution of AOD enhancement is very local and it is correlated well with the active faults and surface ruptures. We suggest that this unique enhancement could be associated with the Lithosphere–Atmosphere–Ionosphere coupling process during the preparation of the Wenchuan earthquake.