Satellite derived aerosol optical depth climatology over Bangalore,India

Climatological aerosol optical depths (AOD) over Bangalore, India have been examined to bring out the temporal heterogeneity in columnar aerosol characteristics. AOD values at 550 nm derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor onboard NASA’s Terra and Aqua satellites, for the period of 2002–2011 have been analyzed (independently) for the purpose. Frequency distributions of the AOD values are examined to infer the monthly mean values. Monthly and seasonal variations of AOD are investigated in the light of regional synoptic meteorology. Climatological monthly and seasonal mean Terra and Aqua AOD values exhibited similar temporal variation patterns. Monthly mean AOD values increased from January, peaks during May and thereafter (except for a secondary peak during July) fall off to reach a minimum during December. Monsoon season recorded the highest climatological seasonal mean AOD, while winter season recorded the lowest. AOD values show an overall increasing trend on a yearly basis, which was found mainly due to sustained increase in the seasonal averaged AOD during summer. The results obtained in the present study are compared with that of the earlier studies over the same location and also with AOD over various other Indian locations. Finally, the radiative and climatic impacts are discussed.     

Seasonal variability of atmospheric aerosol over the North Indian region during 2005–2009

The Indo-Gangetic basin (IGB) extends 2000 km in length along NW–SE and has 400 km width, in the north the basin is bounded by towering Himalaya. High aerosol optical depth (AOD) is observed over the IGB throughout the year. The Himalaya restricts the transport of aerosols across Tibet and China. We have used ground based Kanpur and Gandhi College Aerosol Robotic Network (AERONET) stations and Multiangle Imaging SpectroRadiometer (MISR) and Moderate Resolution Imaging Spectroradiometer (MODIS) Terra level-3 AOD products for the years 2005–2009 to study the variability of aerosol over the Indo-Gangetic (IG) plains. An increase in both satellite-derived as well as ground observed aerosol loading during 2005–2009 has been found over major cities located in the IG plains. The correlation coefficients between AERONET and MISR data are found to be 0.70, 0.36 0.82, in contrast the correlation coefficients between AERONET and MODIS 0.49, 0.68, and 0.43, respectively during summer, winter and monsoon seasons. The AOD estimation using MISR is found to be close to AERONET data during summer and monsoon seasons, in contrast MODIS estimation is better during winter season.     

Monitoring spatio-temporal variations in aerosols and aerosol–cloud interactions over Pakistan using MODIS data

Clouds are important elements in climatic processes and interactions between aerosols and clouds are therefore a hot topic for scientific research. Aerosols show both spatial and temporal variations, which can lead to variations in the microphysics of clouds. In this research, we have examined the spatial and temporal variations in aerosol particles over Pakistan and the impact of these variations on various optical properties of clouds, using Moderate Resolution Imaging Spectroradiometer (MODIS) data from the Terra satellite. We used the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model for trajectory analysis to reveal the origins of air masses, with the aim of understanding these spatial and temporal variabilities in aerosol concentrations. We also documented seasonal variations in patterns of aerosol optical depth (AOD) over Pakistan, for which the highest values occur during the monsoon season (June–August). We then analyzed the relationships between AOD and four other cloud parameters, namely water vapour (WV), cloud fraction (CF), cloud top temperature (CTT) and cloud top pressure (CTP). Regional correlation maps and time series plots for aerosol (AOD) and cloud parameters were produced to provide a better understanding of aerosol–cloud interaction. The analyses showed strong positive correlations between AOD and WV for all of the eight cities investigated. The correlation between AOD and CF was positive for those cities where the air masses were predominantly humid, but negative for those cities where the air masses were relatively dry and carried a low aerosol abundance. These correlations were clearly dependent on the meteorological conditions for all of the eight cities investigated. Because of the observed AOD–CF relationship, the co-variation of AOD with CTP and CTT may be attributable to large-scale meteorological variations: AOD showed a positive correlation with CTP and CTT in northern regions of Pakistan and a negative correlation in southern regions.     

Validation of space-born Moderate Resolution Imaging Spectroradiometer remote sensors aerosol products using application of ground-based Multi-wavelength Radiometer

The measurements of aerosol optical properties were carried out during April 2006 to March 2011 over Mohal (31.9°N, 77.12°E) in the northwestern Indian Himalaya, using the application of ground-based Multi-wavelength Radiometer (MWR) and space-born Moderate Resolution Imaging Spectroradiometer (MODIS) remote sensors. The average (±standard deviation) values of aerosol optical depth (AOD) at 500 nm, Ångström exponent and turbidity coefficient during the entire measurement period were 0.25 ± 0.09, 1.15 ± 0.42 and 0.12 ± 0.06 respectively. About 86% AOD values retrieved from MODIS remote sensor were found within an uncertainty limit (Δτ_pλ = ±0.05 ± 0.15τ_pλ). In general, the MWR derived AOD values were higher than that of MODIS retrieval with absolute difference ∼0.02. During the entire period of measurement space-born MODIS remote sensor and ground-based MWR observation showed good correspondence with significant correlation coefficient ∼0.78 and root mean square difference ∼0.06. For daily observations the relative difference between these two estimates stood less than 9%. However, satellite-based and ground-based observation showed good correspondence, but further efforts still needed to eliminate systematic errors in the existing MODIS algorithm.     

Evaluation and utilization of CloudSat and CALIPSO data to analyze the impact of dust aerosol on the microphysical properties of cirrus over the Tibetan Plateau

The present study elucidates on the evaluation of two versions (V3 and V4.10) of vertical feature mask (VFM) and aerosol sub-types data derived from the Cloud-Aerosol LiDAR and Infrared Pathfinder Satellite Observations (CALIPSO), and its utilization to analyze the impact of dust aerosol on the microphysical properties of cirrus over the Tibetan Plateau (TP). In conjunction to the CALIPSO, we have also used the CloudSat data to study the same during the summer season for the years 2007–2010 over the study area 25–40°N and 75–100°E. Compared to V3 of CALIPSO, V4.10 was found to have undergone substantial changes in the code, algorithm, and data products. Intercomparison of both versions of data products in the selected grid between 30–31°N and 83–84°E within the study area during 2007–2017 revealed that the VFM and aerosol sub-types are in good agreement of ～95.27% and ～82.80%, respectively. Dusty cirrus is defined as the clouds mixed with dust aerosols or existing in dust aerosol conditions, while the pure cirrus is that in a dust-free environment. The obtained results illustrated that the various microphysical properties of cirrus, namely ice water content (IWC), ice water path (IWP), ice distribution width (IDW), ice effective radius (IER), and ice number concentration (INC) noticed a decrease of 17%, 18%, 4%, 19%, and 10%, respectively due to the existence of dust aerosol, consistent with the classical “Twomey effect” for liquid clouds. Moreover, the aerosol optical depth (AOD) showed moderate negative correlations between ?0.4 and ?0.6 with the microphysical characteristics of cirrus. As our future studies, in addition to the present work undertaken, we planned to gain knowledge and interested to explore the impact of a variety of aerosols apart from the dust aerosol on the microphysical properties of cirrus in different regions of China.     

Retrieval of aerosol size distribution from extinction measurements and verification with observations at a tropical station

Aerosol size distributions were retrieved by computing aerosol extinction parameters using extensive measurements of direct solar radiation made in the 0.4 and 0.6 μm wavelengths at Pune with Volz type sunphotometer during winter (November–February), pre-monsoon (March–May), monsoon (June–August) and post-monsoon (September–October) of 1980–1981. The computer aerosol size distributions are compared with the direct measurements made using Anderson eight-stage cascade impactor. There is agreement between the retrieved and measured size distributions. The retrieval method is simple and useful for intensive aerosol measurement programmes.     

Impact of emissions from anthropogenic sources on satellite-derived reflectance

The present study deals with the impact of extensive anthropogenic activities associated with festivities and agricultural crop residues burning in the Indo-Gangetic Plains (IGP) on satellite-derived reflectance during November 2007. Intense smoke plumes were observed in the IRS-P4 OCM satellite data over IGP associated with agricultural crop residue burning during the study period. Terra-MODIS AOD and CALIPSO LIDAR backscatter datasets were analysed over the region to understand the spatial and temporal variation of the aerosol properties. Ground-based measurements on aerosol optical properties and black carbon (BC) mass concentration were carried out during 7–14 November 2007 over urban region of Hyderabad, India. Top of the Atmosphere (TOA) reflectance estimated from IRS-P4 Ocean Color Monitor (OCM) data showed large variations due to anthropogenic activities associated with crop residue burning and fireworks. Atmospheric corrections to OCM satellite data using 6S radiative transfer code with inputs from ground and satellite measurements could account for the variations due to differential aerosol loading. Results are discussed in this paper.     

Direct radiative forcing from black carbon aerosols over urban environment

There is growing evidence that the earth’s climate is changing and will likely continue to change in the future. It is still debated whether these changes are due to natural variability of the climate system or a result of increases in the concentration of greenhouse gases in the atmosphere. Black carbon (BC) has become the subject of interest for a variety of reasons. BC aerosol may cause environmental as well as harmful health effects in densely inhabited regions. BC is a strong absorber of radiation in the visible and near-infrared part of the spectrum, where most of the solar energy is distributed. Black carbon is emitted into the atmosphere as a byproduct of all combustion processes, viz., vegetation burning, industrial effluents and motor vehicle exhausts, etc. In this paper, we present results from our measurements on black carbon aerosols, total aerosol mass concentration and aerosol optical depth over an urban environment namely Hyderabad during January to May, 2003. Diurnal variations of BC indicate high BC concentrations during 6:00–9:00 and 19:00–23:00 h. Weekday variations of BC concentrations increase gradually from Monday to Wednesday and gradually decrease from Thursday to Sunday. Analysis of traffic density along with meteorological parameters suggests that the primary determinant for BC concentration levels and patterns is traffic density. Seasonal variations of BC suggest that the BC concentrations are high during dry season compared to rainy season due to the scavenging by air. The fraction of BC to total mass concentration has been observed to be 7% during January to May. BC showed positive correlation with total mass concentration and aerosol optical depth at 500 nm. Radiative transfer calculations suggests that during January to May, diurnal averaged aerosol forcing at the surface is −33 Wm² and at the top of the atmosphere (TOA) above 100 km it is observed to be +9 Wm⁻². The results have been discussed in detail in the paper.     

A surface reflectance correction model to improve the retrieval of MISR aerosol optical depth supported by MODIS data

Combined use of different satellite sensors are known to improve retrievals of aerosol optical depth (AOD). In this study, we propose a new method for retrieving Multi-angle Imaging SpectroRadiometer (MISR) AOD data supported by Moderate Resolution Imaging Spectroradiometer (MODIS) data in Jiangsu Province, China, over the period of 2016–2018 using MODIS L1B, bidirectional reflectance distribution function (BRDF), MISR 1B2T, and ground-measured AOD data. This method is based on the surface reflectance determined by the MODIS V5.2 algorithm. Through the observation angle and spectral conversion between different sensors, the MISR AOD can be obtained. The correlation coefficient (R) and root-mean-square error (RMSE) between the retrieved MISR and ground-measured AOD data varied between different seasons. The accuracy of the MISR AOD retrieval was notably improved after correcting the MISR surface reflectance. Therefore, the method proposed in this study is feasible for the retrieval of MISR AOD supported by MODIS data, and will be applicable to atmospheric environmental monitoring over large areas in the future.     

Effects of CME and CIR induced geomagnetic storms on low-latitude ionization over Indian longitudes in terms of neutral dynamics

This paper presents the response of the ionosphere during the intense geomagnetic storms of October 12–20, 2016 and May 26–31, 2017 which occurred during the declining phase of the solar cycle 24. Total Electron Content (TEC) from GPS measured at Indore, Calcutta and Siliguri having geomagnetic dips varying from 32.23°N, 32°N and 39.49°N respectively and at the International GNSS Service (IGS) stations at Lucknow (beyond anomaly crest), Hyderabad (between geomagnetic equator and northern crest of EIA) and Bangalore (near magnetic equator) in the Indian longitude zone have been used for the storms. Prominent peaks in diurnal maximum in excess of 20–45 TECU over the quiet time values were observed during the October 2016 storm at Lucknow, Indore, Hyderabad, Bangalore and 10–20 TECU for the May 2017 storm at Siliguri, Indore, Calcutta and Hyderabad. The GUVI images onboard TIMED spacecraft that measures the thermospheric O/N₂ ratio, showed high values (O/N₂ ratio of about 0.7) on October 16 when positive storm effects were observed compared to the other days during the storm period. The observed features have been explained in terms of the O/N₂ ratio increase in the equatorial thermosphere, CIR-induced High Speed Solar Wind (HSSW) event for the October 2016 storm. The TEC enhancement has also been explained in terms of the Auroral Electrojet (AE), neutral wind values obtained from the Horizontal Wind Model (HWM14) and equatorial electrojet strength from magnetometer data for both October 2016 and May 2017 storms. These results are one of the first to be reported from the Indian longitude sector on influence of CME- and CIR-driven geomagnetic storms on TEC during the declining phase of solar cycle 24.     

Detection of aerosol pollution sources during sandstorms in Northwestern China using remote sensed and model simulated data

The present paper has used a comprehensive approach to study atmosphere pollution sources including the study of vertical distribution characteristics, the epicenters of occurrence and transport of atmospheric aerosol in North-West China under intensive dust storm registered in all cities of the region in April 2014. To achieve this goal, the remote sensing data using Moderate Resolution Imaging Spectroradiometer satellite (MODIS) as well as model-simulated data, were used, which facilitate tracking the sources, routes, and spatial extent of dust storms. The results of the study have shown strong territory pollution with aerosol during sandstorm. According to ground-based air quality monitoring stations data, concentrations of PM₁₀ and PM_2.5 exceeded 400?μg/m³ and 150?μg/m³, respectively, the ratio PM_2.5/PM₁₀ being within the range of 0.123–0.661. According to MODIS/Terra Collection 6 Level-2 aerosol products data and the Deep Blue algorithm data, the aerosol optical depth (AOD) at 550?nm in the pollution epicenter was within 0.75–1. The vertical distribution of aerosols indicates that the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) 532?nm total attenuates backscatter coefficient ranges from 0.01 to 0.0001?km^?1?×?sr^?1 with the distribution of the main types of aerosols in the troposphere of the region within 0–12.5?km, where the most severe aerosol contamination is observed in the lower troposphere (at 3–6?km). According to satellite sounding and model-simulated data, the sources of pollution are the deserted regions of Northern and Northwestern China.     

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.     

Determination of aerosol optical depth from ground measurements

The methods used to determine the aerosol optical depth as a function of wavelength are briefly described and discussed. Some results from the operational network of the World Meteorological Organization and other, more research oriented studies, are reviewed and critically analysed to assess the reliability and accuracy of such determinations and their value as ground truth measurements for space applications.     

Temporal evolution of the EIA along 95°E as obtained from GNSS TEC measurements and SAMI3 model

The total electron content (TEC) derived from GNSS measurements at a trans-hemispheric meridional chain of ground stations around 95°E longitude are used to study the quiet time inter-hemispheric structure and dynamics of the equatorial ionization anomaly (EIA) during the period March 2015 to February 2016. The stations are Dibrugarh (27.5°N, 95°E, 43° dip), Kohima (25.6°N, 94.1°E, 39° dip), Aizawl (23.7°N, 92.8°E, 36° dip), Port Blair (11.63°N, 92.71°E, 9° dip) and Cocos Islands (12.2°S, 96.8°E, 43° dip). The observation shows that the northern crest of the EIA lies in the south of 23°N (Aizawl) in all seasons but recedes further south towards the equator during December solstice. The largest poleward expansion of the northern (southern) EIA is observed in the March equinox (December solstice). The equinoctial and hemispherical asymmetry of TEC is noted. The winter anomaly is observed in the northern hemisphere but not in the southern hemisphere. The highest midday TEC over any station is observed in the March equinox. The TEC in southern summer (December solstice) is significantly higher than that in the northern summer (June solstice). The observed northern EIA contracts equatorward in the postsunset period of solstice but the southern EIA persists late into the midnight in the December solstice. The asymmetry may be attributed to the different geographic location of the magnetically conjugate stations. The SAMI3 simulations broadly capture the EIA structure and the inter-hemispheric asymmetry during solstices. The difference between observations and the SAMI3 is higher in March equinox and December solstice. The higher E?×?B vertical drift in the 90–100°E sector and the large geographic-geomagnetic offset in observing stations may have contributed to the observed differences.     

Seasonal variations of GPS derived TEC at three different latitudes of the southern hemisphere during geomagnetic storms

Data from the archive of the International GNSS Services (IGS) were used to study the seasonal variations of Total Electron Content (TEC) over three stations located at different latitudes in the southern hemisphere during the geomagnetic storms of 11 January, 6 April, 8 June, and 13 October 2000, representing storms that occurred in summer, autumn equinox, winter and spring equinox, respectively. The percentage TEC deviation with respect to reference values differs substantially from season to season. A strong seasonal anomaly and clear equinoctial asymmetry in TEC response to the storms were observed. Weak and short-lived positive TEC deviations as well as strong and long-lasting negative trends were observed in summer storm during the main and recovery phases respectively over the high and low latitudes whereas in winter storm, the highest positive TEC deviations was recorded during the main phase over the entire latitudes. TEC enhancement dominated all the stations during the autumn (March) equinox storm while TEC depletion was majorly observed during the spring (September) equinox. All these variations find their explanations in the thermospheric composition change and circulation. Future work with direct or modeled measurement of atomic Oxygen to molecular Nitrogen ratio (O/N₂), large number of storms and other possible factors such as variations in storm’s intensity and local time dependence of the storm onset is expected to validate the observations in this study.     

The longitudinal and hemispherical variation of the IRI’s foF2 estimates at ±40° dip angle around 95°E and 130°E sector under fluctuating solar flux conditions

The deviation of the IRI estimates of the monthly mean foF2 in the low mid latitude of 95°E–130°E longitude sector is investigated using simultaneous ground measurements at four stations during 2010–2014. The stations form two conjugate pairs of the same geo-magnetic latitude at two fixed longitudes enabling direct longitudinal and hemispheric comparison. The temporal, spatial, seasonal and solar activity variations of the deviations are discussed with reference to the longitudinal density variation in the transition region between low and midlatitudes. Cases of underestimation/overestimation as well as good estimate are noted. Underestimation (overestimation) in the daytime and overestimation (underestimation) in the nighttime of 95°E (130°E) are common. The longitudinal difference in the measurements suggests negative (positive) foF2 gradient from west to east in daytime (nighttime). In contrast, the IRI predicts flatter or increasing longitudinal profiles from 95°E to 130°E. The local time and longitudinal variation of the IRI deviations can be attributed to the combined role of the longitudinal EIA structure as well as midlatitude zonal wind-magnetic declination effect. The station/season independent deviations relate the role of solar activity representation in the IRI. These deviations may be attributed to the weak IRI response to rapid solar flux fluctuations.     

Seasonal,annual, and interannual variability in MLT quasi-two-day waves over the low-latitude region Kolhapur (16.8°N; 74.2°E)

This study presents the quasi-two-day wave (Q2DW) characteristics of the mesosphere and lower thermosphere (MLT) region obtained by taking hourly mean values of horizontal wind velocities for 4? years (August 2013–July 2017) through continuous measurements using a medium-frequency (MF) radar (operating frequency – 1.98?MHz) located at the low-latitude Indian station Kolhapur (16.8°N; 74.2°E). The MF radar located at Kolhapur was upgraded in 2013, and these results of Q2DW have been reported for the first time after upgrading. The present study investigated variability in seasonal, annual, interannual, and solar indices of Q2DWs traveling in zonal (EW) and meridional (NS) components in the MLT region. The Q2DW activity is observed to be stronger during austral summer (January–February) (EW?=?～5?m/s and NS?=?～8–10?m/s) than during boreal summer (June–July) (EW = ～5 m/s and NS = ～6–8?m/s). The Q2DW amplitudes are larger in the meridional component than in the zonal one. A strong semiannual oscillation (SAO) has been observed in Q2DWs, with peak during January–February and June–July. In addition, small enhancement is seen in meridional Q2DW in October (～5–6?m/s). It is observed that the entire spectrum (40–60?h) measured between 86 and 94?km contributes to the SAO amplitudes during January–February and June–July, whereas the waves measured between 42?h and 52?h contribute to enhancement in October similar to that reported elsewhere. In general, the Q2DW amplitude shows large interannual variability. The easterlies developed in the global circulation model in Northern hemisphere during May intensify up to around summer solstice. Q2DW activity peaks during westerly shear zone and intensifies with time at a lower thermospheric altitude (above 90?km). Small positive correlations (r?=?0.2 for sunspot number and r?=?0.1 for 10.7?cm solar flux) have been observed between Q2DW amplitudes and solar activity.     

Satellite evidence of harmful algal blooms and related oceanographic features in the Bohai Sea during autumn 1998

Harmful algal blooms (HABs) are truly global marine phenomena of increasing significance. Some HAB occurrences are different to observe because of their high spatial and temporal variability and their advection, once formed, by surface currents. A serious HAB occurred in the Bohai Sea during autumn 1998, causing the largest fisheries economic loss. The present study analyzes the formation, distribution, and advection of HAB using satellite SeaWiFS ocean color data and other oceanographic data. The results show that the bloom originated in the western coastal waters of the Bohai Sea in early September, and developed southeastward when sea surface temperature (SST) increased to 25–26 °C. The bloom with a high Chl-a concentration (6.5 mg m⁻³) in center portion covered an area of 60 × 65 km². At the end of September, the bloom decayed when SST decreased to 22–23 °C. The HAB may have been initiated by a combination of the river discharge nutrients in the west coastal waters and the increase of SST; afterwards it may have been transported eastward by the local circulation that was enhanced by northwesterly winds in late September and early October.     

Seasonal-latitudinal variations of the mesospheric ozone density

Quartz-UV occultation measurements by the satellite Interkosmos-16 have been used to calculate ozone densities at altitudes between 55 and 75 km for the period July 27 – October 28 of 1976. Although the densities agree quite well with the Krueger-Minzner-model below 65 km distinct seasonal-latitudinal variations have been found. During July and August latitudinal variations are more pronounced than in September and October with a slight maximum shifting from 5° S in July to 30 – 40° S in September. A comparison of different height levels shows a decreasing latitudinal variation for increasing altitude during July and August and rather modest variations for September and October.     

Global MGS TES data and Mars-GRAM validation

Mars Global Reference Atmospheric Model (Mars-GRAM 2001) is an engineering-level Mars atmosphere model widely used for many Mars mission applications. From 0 to 80 km, it is based on NASA Ames Mars General Circulation Model (MGCM), while above 80 km it is based on University of Michigan Mars Thermospheric General Circulation Model. Mars-GRAM 2001 and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA). Validation studies are described comparing Mars-GRAM with a global summary data set of Mars Global Surveyor Thermal Emission Spectrometer (TES) data. TES averages and standard deviations were assembled from binned TES data which covered surface to ∼40 km, over more than a full Mars year (February 1999–June 2001, just before start of a Mars global dust storm). TES data were binned in 10° × 10° latitude–longitude bins (36 longitude bins, centered at 5°–355°, by 18 latitude bins, centered at −85° to +85°), and 12 seasonal bins (based on 30° increments of L_s angle). Bin averages and standard deviations were assembled at 23 data levels (temperature at 21 pressure levels, plus surface temperature and surface pressure). Two time-of-day bins were used: local time near 2 or 14 h. Two dust optical depth bins were used: infrared optical depth, either less than or greater than 0.25 (which corresponds to visible optical depth less than or greater than about 0.5). For interests in aerocapture and precision entry and landing, comparisons focused on atmospheric density. TES densities versus height were computed from TES temperature versus pressure, using assumptions of perfect gas law and hydrostatics. Mars-GRAM validation studies used density ratio (TES/Mars-GRAM) evaluated at data bin center points in space and time. Observed average TES/Mars-GRAM density ratios were generally 1 ± 0.05, except at high altitudes (15–30 km, depending on season) and high latitudes (>45°N), or at most altitudes in the southern hemisphere at L_s ∼ 90° and 180°. Compared to TES averages for a given latitude and season, TES data had average density standard deviation about the mean of ∼2.5% for all data, or ∼1–4%, depending on time of day and dust optical depth. Average standard deviation of TES/Mars-GRAM density ratio was 8.9% for local time 2 h and 7.1% for local time 14 h. Thus standard deviation of observed TES/Mars-GRAM density ratio, evaluated at matching positions and times, is about three times the standard deviation of TES data about the TES mean value at a given position and season.