Traveling waves in the martian atmosphere from MGS TES Nadir data

We have characterized the annual behavior of martian atmospheric traveling waves in the MGS TES data set from the first two martian years of mapping. There is a high degree of repeatability between the two years. They are dominated by strong low zonal wavenumber waves with high amplitudes near the polar jets, strongest in late northern fall and early northern winter. The m=1 waves have amplitudes up to about 20 K, are vertically extended, and occasionally extend even into the tropics. Periods for m=1 range from 2.5 to 30 sols. Much weaker waves were identified in the south, with amplitudes less than about 3.5 K. Traveling waves with m=2 and m=3 are also seen, but their amplitudes are typically limited to less than 4 K, and are generally more confined near the surface. In the north, they are more evident in fall and spring rather than winter solstice, which is clearly dominated by m=1 waves. Some evidence of storm tracks has been identified in the data, with accentuated weather-related temperature perturbations near longitudes 200° to 320° E for both the southern and northern hemispheres near latitude ±65° at the surface. Some evidence was also found for a sharpening of longitudinal gradients into what may be frontal systems. EP flux divergences show the waves extracting energy from the zonal mean winds. When the m=1 waves were strongest, decelerations of the zonal jet of order 30 m/(s sol) were measured. Above 1 scale height, the waves extract energy from the jet predominately through barotropic processes, but their character is overall mixed barotropic/baroclinic. Inertial instabilities may exist at altitude on the equatorward flanks of the polar jets, and marginal stability extends through to the tropics. This may explain the coordination of the tropical behavior of the waves with that centered along the polar jet, consistent with the ideas expressed in Wilson et al. (2002, Geophys. Res. Lett. 29, #1684) and similar to those in Barnes et al. (1993, J. Geophys. Res. 98, 3125-3148). Throughout the year, there exist large regions with the meridional gradient of PV less than zero, but they are strongest near winter solstice. Poleward of the winter jet, the regions of instability reach the surface, equatorward they do not. These regions, satisfying a necessary criterion for instability, likely explain the genesis of the waves, and perhaps also their bimodal character between surface (faster waves) and altitude (slow m=1 waves).     

Forced waves in the martian atmosphere from MGS TES nadir data

We have analyzed the temperature retrievals from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) nadir spectra to yield latitude-height resolved maps of various atmospheric forced wave modes as a function of season for a full Mars year. Among the isolated wave modes is the zonal mean, time mean temperature, which we used to derive zonal mean zonal winds and stationary wave quasi-geostrophic indices of refraction, diagnostic of their propagation. The diurnal Kelvin wave was isolated in the data, with results roughly consistent with models (Wilson and Hamilton, 1996, J. Atmos. Sci. 33, 1290-1326). The s = 1 and s = 2 stationary waves were found to have significant amplitude in ducts extending up the winter polar jets, while the s = 3 stationary wave was found to be confined to near the surface. The s = 1 stationary wave was found to have little phase tilt with height during northern winter, but significant westward phase tilt with height in the southern winter. This indicates that the wave carries heat poleward, slightly more than that found in Barnes et al. (1996; J. Geophys. Res. 101, 12,753-12,776). The s = 1 stationary wave is likely the dominant mechanism for eddy meridional heat transport for the southern winter. We noted that the phase of the s = 2 stationary wave is nearly constant with time, but that the s = 1 stationary wave moved 90° of longitude from fall to winter and back in spring in the North. While interannual variability is not yet addressed, overall, these results provide the first comprehensive benchmark for forced waves in Mars’s atmosphere against which future atmospheric models of Mars can be compared.     

Observing the martian surface albedo pattern: Comparing the AEOS and TES data sets

High spatial resolution images of Mars were acquired with the Advanced Electro-Optical System (AEOS) 3.63-meter telescope at the Maui Space Surveillance System (MSSS) during both the 2001 and 2003 Mars apparitions. Comparisons are made of the surface albedo patterns obtained from these AEOS images to the surface albedo maps constructed from the Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) data taken during the same time periods. These comparisons demonstrate that the images provide albedo information in a limited area surrounding the sub-Earth point that is consistent with the TES-derived albedo field. Additionally, it is shown that by employing adaptive optics (AO), the typical ground-based observing season of Mars can be extended. This is the only known published AO data set of Mars with temporal coverage over an entire apparition. Changes in the surface albedo affect the local ground temperature, which impacts the depth of the planetary boundary layer (PBL) above the surface. Since it is the state of the PBL that controls surface/atmospheric interaction, albedo variations have the power to alter the amount of dust that is lifted. A one-dimensional radiative/convective version of the NASA Ames Mars General Circulation Model is used to demonstrate that the measured albedo variations can alter the daytime ground temperatures by as much as 5 K, which in turn alters the structure of the planetary boundary layer (PBL). Therefore, albedo changes are thermodynamically important, and the ability to characterize them, should orbital observations become unavailable, is a valuable capability.     

Paleo-pole positions from martian magnetic anomaly data

Magnetic component anomaly maps were made from five mapping cycles of the Mars Global Surveyor's magnetometer data. Our goal was to find and isolate positive and negative anomaly pairs which would indicate magnetization of a single source body. From these anomalies we could compute the direction of the magnetizing vector and subsequently the location of the magnetic pole existing at the time of magnetization. We found nine suitable anomaly pairs and from these we computed paleo-poles that were nearly equally divided between north, south and mid-latitudes. These results suggest that during the existence of the martian main magnetic field it experienced several reversals and excursions.     

Water vapor variability in the north polar region of Mars from Viking MAWD and MGS TES datasets

Atmospheric water vapor abundances in Mars’ north polar region (NPR, from 60° to 90°N) are mapped as function of latitude and longitude for spring and summer seasons, and their spatial, seasonal, and interannual variability is discussed. Water vapor data are from Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) and the Viking Orbiter (VO) Mars Atmospheric Water Detector (MAWD). The data cover three complete northern spring-summer seasons in 1977-1978, 2000-2001 and 2002-2003, and shorter periods of spring-summer seasons during 1975, 1999 and 2004. Long term interannual variability in the averaged NPR abundances may exist, with Viking MAWD observations showing twice as much water vapor during summer as the MGS TES observations more than 10 martian years (MY) later. While the averaged abundances are very similar in TES observations for the same season in different years, the spatial distributions in the early summer season do vary significantly year over year. Spatial and temporal variabilities increase between L_s ∼ 80-140°, which may be related to vapor sublimation from the North Polar Residual Cap (NPRC), or to changes in circulation. Spatial variability is observed on scales of ∼100 km and temporal variability is observed on scales of <10 sols during summer. During late spring the TES water vapor spatial distribution is seen to correlate with the low topography/low albedo region of northern Acidalia Planitia (270-360°E), and with the dust spatial distribution across the NPR during late spring-early summer. Non-uniform vertical distribution of water vapor, a regolith source or atmospheric circulation ‘pooling’ of water vapor from the NPRC into the topographic depression may be behind the correlation with low topography/low albedo. Sublimation winds carrying water vapor off the NPRC and lifting surface dust in the areas surrounding the NPRC may explain the correlation between the water vapor and dust spatial distributions. Correlation between water vapor and dust in MAWD data are only observed over low topography/low albedo area. Maximum water vapor abundances are observed at L_s = 105-115° and outside of the NPRC at 75-80°N; the TES data, however, do not extend over the NPRC and thus, this conclusion may be biased. Some water vapor appears to be released in plumes or ‘outbursts’ in the MAWD and TES datasets during late spring and early summer. We propose that the sublimation rate of ice varies across the NPRC with varying surface winds, giving rise to the observed ‘outbursts’ at some seasons.     

Water-ice clouds and dust in the north polar region of Mars using MGS TES data

Water-ice and dust optical depths in Mars’ north polar region are mapped as function of season, latitude and longitude, and their characteristics and variability on a geographic, seasonal, and interannual basis are discussed. We use water-ice and dust optical depth data provided by the Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES), covering nearly three northern spring and summer periods. We find that interannual variability exists in both the water ice and dust behavior, although there are trends that repeat year to year as well. The optical thickness of the north polar hood (NPH) exhibits interannually varying longitudinal structure, both during springtime recession and late-summer onset. We define the characteristics associated with the transition to and from the NPH and find that the disappearance occurs near L_s=75° and the reappearance near L_s=160-165°. We find that the late spring to early summer time frame is characterized by very low water-ice optical depths and enhanced dust activity, with a preference for lower water-ice and higher dust optical depths in the 0-90°W quadrant. We see possible evidence for stationary wavenumber 2 systems in a few of the maps examined.     

Effect of atmospheric dust loading on martian albedo measurements

This work is devoted to the analysis of the variation of albedo measured by orbiting instruments with atmospheric opacity on Mars. The study has been conduced by analysing Mars Global Surveyor Thermal Emission Spectrometer (MGS-TES) data from martian regions with different surface albedo.In support of these data, synthetic spectra with different surface albedo and atmospheric opacities have been computed, so that a comparison has been performed. The synthetic spectra have been retrieved by using two different grain sizes for suspended dust (0.5 and 1.2 μm), allowing a comparison between the two models and the observations.Using the DCI, a parameter describing the quantity of dust deposited on the surface, the effectiveness of the single scattering approximation has been tested for low atmospheric opacity by analysing the quality of the linear fit up to different atmospheric opacity.For more opaque conditions two kinds of fits have been applied to the data, linear and second-order degree polynomial. In this case, we found that the polynomial fit better describes the observations.The analysis of these data made it possible to notice a peculiar trend, already reported by Christensen (1988), of the albedo over Syrtis Major after the occurrence of dust storms, but, differently from that work, now the study of DCI together with atmospheric opacity and albedo allowed us to robustly confirm the hypothesis made by Christensen.Finally, the comparison between observations and synthetic spectra computed with models with different particles grain sizes indicates that dust particles of 0.5 μm diameter are the most effective to change the aerosol atmospheric opacity on Mars.     

Diffraction theory modeling of near-forward radio wave scattering from particle clusters

We investigate the use of diffraction theory as an alternative to solving the full electromagnetic problem of near-forward scattering by a collection of dielectric spheres of arbitrary spacing and orientation. We define an amplitude screen, which is constructed by projecting the shadows of a cluster of spheres onto a plane perpendicular to the wave vector of the incident radio wave. The far-field (Fraunhofer) diffraction pattern of the amplitude screen is then computed and compared with the Mie-theoretic result. We show that for EM scattering from a cluster of electrically large spheres—both singly-sized and belonging to a size distribution—there is excellent agreement between the exact Mie solution and its diffraction theory approximation when near-forward scattering is the angular range of interest. This excellent agreement holds over a broad range of particle separation and orientation configurations relative to the incidence direction. It is also achieved at a much reduced computational cost compared with an exact solution of the electromagnetic interaction problem. Fortified by these results, the authors have applied diffraction theory to the analysis of Cassini radio occultation data, thereby detecting fine-scale structure in Saturn’s rings [Thomson, F.S., Marouf, E.A., Tyler, G.L., French, R.G., Rappoport, N.J., 2007. Periodic microstructure in Saturn’s rings A and B. Geophys. Res. Lett. 34, L24203].     

The 3–5 MHz global reflectivity map of Mars by MARSIS/Mars Express: Implications for the current inventory of subsurface H2O

We extracted the surface echo power from 2 years of MARSIS measurements. The retrieved values are calibrated to compensate for changes in the distance of the spacecraft to the surface and for the attenuation of the signal by the ionosphere. The results are used to build the first global map of surface echo power at 3–5 MHz. The surface echo power variations are primarily caused by kilometer-scale surface roughness. Then, we derive the values of dielectric constant of the shallow subsurface materials by normalizing the surface echo power map using a simulation of MARSIS signal from the MOLA topography. As a result, we obtain a map that characterizes the dielectric properties of the materials down to a few decameters below the surface. Dielectric properties vary with latitude, with high values in mid-latitudes belts (20–40°) and lower values at both equatorial and high latitudes. From the comparison of MARSIS reflectivity map to GRS observations, we conclude that the reflectivity decrease observed poleward of 50–60° corresponds to the onset of water-ice occurrence within the regolith. Assuming homogenous ground composition and texture at the scale of the MARSIS resolution cell, our inferred volume of ground water ice is of 10⁶ km³, equivalent to a polar cap. Low reflectivity areas are also observed in equatorial regions. From radar studies alone, equatorial low dielectric constant values could have different interpretations but the correlation with GRS hydrogen distribution rather points toward a water-related explanation.     

Broadband submillimeter measurements of the full Moon center brightness temperature and application to a lunar eclipse

We report on observations of the full Moon brightness temperature covering the frequency range of 300-950 GHz, and also on observations of the lunar eclipse of July 16, 2000, though only covering the frequency range of 165-365 GHz due to poor atmospheric transmission at higher frequencies. All observations were performed from the summit of Mauna Kea (HI) using a Fourier Transform Spectrometer mounted on the Caltech Submillimeter Observatory and supplemented by measurements of the atmospheric opacity using a 183 GHz Water Vapor Monitor. The telescope was pointed to the center of the lunar disk (with a footprint of ∼45-15 km on the Moon at 300 through 900 GHz). In order to obtain the correct values of the Moon brightness temperatures at all frequencies we carefully corrected for the atmospheric absorption, which varies across the submillimeter domain. This correction is fully described. The measured pre-eclipse brightness temperature is around 337 K in the 165-365 GHz range. This temperature slightly increases with frequency to reach ∼353 K at 950 GHz, according to previous broader band data. The magnitude of the temperature drop observed during the eclipse at 265 GHz (central frequency of the band covered) was about ∼70 K, in very good agreement with previous millimeter-wave measurements of other lunar eclipses. We detected, in addition, a clear frequency trend in the temperature drop that has been compared to a thermal and microwave emission model of the lunar regolith, with the result of a good match of the relative flux drop at different frequencies between model and measurements.     

Mineralogical constraints on the high-silica martian surface component observed by TES

The Thermal Emission Spectrometer (TES) has observed a high-silica material in the dark regions of Mars that is spectrally similar to obsidian glass and may have a volcanic origin. An alternate interpretation is that the spectrally amorphous material consists of clay minerals or some other secondary material, formed by chemical alteration of surface rocks. The regions where this material is observed (e.g., Acidalia Planitia) have relatively high spectral contrast, suggesting that the high-silica material exists as coarse particulates, indurated soils or cements, within rocks, or as indurated coatings on rock surfaces. The geologic interpretation of this spectral result has major implications for understanding magmatic evolution and weathering processes on Mars. One of the complications in interpreting spectral observations of glasses and clay minerals is that both are structurally and compositionally complex. In this study, we perform a detailed spectroscopic analysis of indurated smectite clay minerals and relate their thermal emission spectral features to structural and crystal chemical properties. We examine the spectral similarities and differences between smectite clay minerals and obsidian glass from a structural-chemical perspective, and make further mineralogical interpretations from previous TES results. The results suggest that neither smectite clays nor any clay mineral with similar structural and chemical properties can adequately explain TES observations of high-silica materials in some martian dark regions. If the spectrally amorphous materials observed by TES do represent an alteration product, then these materials are likely to be poorly crystalline aluminosilicates. While all clay minerals have Si/O ratios ?0.4, the position of the emissivity minimum at Mars suggests a Si/O ratio of 0.4-0.5. The spectral observation could be explained by the existence of a silica-rich alteration product, such as Al- or Fe-bearing opal, an intimate physical mixture of relatively pure silica and other aluminosilicates (such as clay minerals or clay precursors), or certain zeolites. The chemical alteration of basaltic rocks on Mars to phyllosilicate-poor, silica-rich alteration products provides a geologically reasonable and consistent explanation for the global TES surface mineralogical results.     

Low-order dynamical behavior in the martian atmosphere: Diagnosis of general circulation model results

The hypothesis of a low dimensional martian climate attractor is investigated by the application of the proper orthogonal decomposition (POD) to a simulation of martian atmospheric circulation using the UK Mars general circulation model (UK-MGCM). In this article we focus on a time series of the interval between autumn and winter in the northern hemisphere, when baroclinic activity is intense. The POD is a statistical technique that allows the attribution of total energy (TE) to particular structures embedded in the UK-MGCM time-evolving circulation. These structures are called empirical orthogonal functions (EOFs). Ordering the EOFs according to their associated energy content, we were able to determine the necessary number to account for a chosen amount of atmospheric TE. We show that for Mars a large fraction of TE is explained by just a few EOFs (with 90% TE in 23 EOFs), which apparently support the initial hypothesis. We also show that the resulting EOFs represent classical types of atmospheric motion, such as thermal tides and transient waves. Thus, POD is shown to be an efficient method for the identification of different classes of atmospheric modes. It also provides insight into the non-linear interaction of these modes.     

Spectral imaging of martian water ice clouds and their diurnal behavior during the 1999 aphelion season (Ls = 130°)

We report high-spectral-resolution (λ/δλ = 800-2300) near-infrared mapping observations of Mars at L_s = 130° (April 1999), which were obtained by drift-scanning the cryogenic long-slit spectrometer at the KPNO 2.2-m telescope across the disk. Data were reformatted into calibrated spectral image cubes (x,y,λ) spanning 2.19 to 4.12 μm, which distinguish atmospheric CO₂ features, solar lines, and surface and aerosol features. Maps of relative band depth between 3.0 and 3.5 μm trace water ice clouds and show the diurnal evolution of features in the persistent northern summer aphelion cloud belt, which was mapped contemporaneously but at fixed local time by the Mars Global Surveyor Thermal Emission Spectrometer (MGS/TES). Cloud optical depth, particle sizes, and ice aerosol content were estimated using a two-stream, single-layer scattering model, with Mie coefficients derived from recently published ice optical constants, followed by a linear spectral deconvolution process. A comparison of data and model spectra shows evaporating nighttime clouds in the morning followed by afternoon growth of a prominent orographic cloud feature on the west flank of Elysium Mons. Cloud optical depth at 3.2 μm evolved to 0.28 ± 0.13 and ice aerosol column abundance to 0.9 ± 0.3 pr μm in the afternoon. Column abundances as large as 0.17 pr μm were retrieved in nonorographic clouds within the aphelion cloud band around midday. These clouds exhibit a modest decline in optical depth during the afternoon. Results show that ice particle radii from <2 μm to >4 μm exist in both cloud types. However, large particles dominate the spectra, consistent with recent MGS/TES emission phase function measurements of aphelion cloud aerosol properties.     

Retrieval algorithm for atmospheric dust properties from Mars Global Surveyor Thermal Emission Spectrometer data during global dust storm 2001A

We present a new parameter retrieval algorithm for Mars Global Surveyor Thermal Emission Spectrometer data. The algorithm uses Newtonian first-order sensitivity functions of the infrared spectrum in response to variations in physical parameters to fit a model spectrum to the data at 499, 1099, and 1301 cm⁻¹. The algorithm iteratively fits the model spectrum to data to simultaneously retrieve dust extinction optical depth, effective radius, and surface temperature. There are several sources of uncertainty in the results. The assumed dust vertical distribution can introduce errors in retrieved optical depth of a few tens of percent. The assumed dust optical constants can introduce errors in both optical depth and effective radius, although the systematic nature of these errors will not affect retrieval of trends in these parameters. The algorithm does not include the spectral signature of water ice, and hence data needs to be filtered against this parameter before the algorithm is applied. The algorithm also needs sufficient dust spectral signature, and hence surface-to-atmosphere temperature contrast, to successfully retrieve the parameters. After the application of data filters the algorithm is both relatively accurate and very fast, successfully retrieving parameters, as well as meaningful parameter variability and trends from tens of thousands of individual spectra on a global scale (Elteto, A., Toon, O.B. [2010]. Icarus, this issue). Our results for optical depth compare well with TES archive values when corrected by the single scattering albedo. Our results are on average 1–4 K higher in surface temperatures from the TES archive values, with greater differences at higher optical depths. Our retrieval of dust effective radii compare well with the retrievals of Wolff and Clancy (Wolff, M.J., Clancy, R.T. [2003]. J. Geophys. Res. 108 (E9), 5097) for the corresponding data selections from the same orbits.     

Radio occultation measurements and MGCM simulations of Kelvin waves on Mars

We have derived new results concerning thermal tides on Mars from a combination of radio occultation measurements and numerical simulations by a Mars General Circulation Model (MGCM). This investigation exploits a set of concurrent observations by Mars Express (MEX) and Mars Global Surveyor (MGS) in mid-2004, when the season on Mars was midspring in the northern hemisphere. The MEX occultations sampled the atmosphere near the evening terminator at latitudes ranging from 54° N to 15° S. The MGS occultations provided complementary coverage near the morning terminator at latitudes of 35° N and 71° S. The geopotential field derived from these measurements contains distinctive modulation caused by solar-asynchronous thermal tides. Through careful analysis of the combined observations, we characterized two prominent wave modes, obtaining direct solutions for some properties, such as the amplitude and phase, as well as constraints on others, such as the period, zonal wave number, and meridional structure. We supplemented these observations with MGCM simulations. After evaluating the performance of the MGCM against the measurements, we used the validated simulation to deduce the identity of the two tidal modes and to explore their behavior. One mode is a semidiurnal Kelvin wave with a zonal wave number of 2 (SK2), while the other is a diurnal Kelvin wave with a zonal wave number of 1 (DK1). Both modes are known to be close to resonance in the martian atmosphere. Our observations of the SK2 are more complete and less ambiguous than any previous measurement. The well-known DK1 is the dominant solar-asynchronous tide in the martian atmosphere, and our results confirm and extend previous observations by diverse instruments.     

The effects of topographically-controlled thermal tides in the martian upper atmosphere as seen by the MGS accelerometer

Mars Global Surveyor accelerometer observations of the martian upper atmosphere revealed large variations in density with longitude during northern hemisphere spring at altitudes of 130-160 km, all latitudes, and mid-afternoon local solar times (LSTs). This zonal structure is due to tides from the surface. The zonal structure is stable on timescales of weeks, decays with increasing altitude above 130 km, and is dominated by wave-3 (average amplitude 22% of mean density) and wave-2 (18%) harmonics. The phases of these harmonics are constant with both altitude and latitude, though their amplitudes change significantly with latitude. Near the South Pole, the phase of the wave-2 harmonic changes by 90° with a change of half a martian solar day while the wave-3 phase stays constant, suggesting diurnal and semidiurnal behaviour, respectively. We use a simple application of classical tidal theory to identify the dominant tidal modes and obtain results consistent with those of General Circulation Models. Our method is less rigorous, but simpler, than the General Circulation Models and hence complements them. Topography has a strong influence on the zonal structure.     

Continuous monitoring of nightside upper thermospheric mass densities in the martian southern hemisphere over 4 martian years using electron reflectometry

Details are presented of an improved technique to use atmospheric absorption of magnetically reflecting solar wind electrons to constrain neutral mass densities in the nightside martian upper thermosphere. The helical motion of electrons on converging magnetic field lines, through an extended neutral atmosphere, is modeled to enable prediction of loss cone pitch angle distributions measured by the Magnetometer/Electron Reflectometer (MAG/ER) experiment on Mars Global Surveyor at 400 km altitude. Over the small fraction of Mars' southern hemisphere (∼2.5%) where the permanent crustal magnetic fields are both open to the solar wind and sufficiently strong as to dominate the variable induced martian magnetotail field, spherical harmonic expansions of the crustal fields are used to prescribe the magnetic field along the electron's path, allowing least-squares fitting of measured loss cones, in order to solve for parameters describing the vertical neutral atmospheric mass density profile from 160 to 230 km. Results are presented of mass densities in the southern hemisphere at 2 a.m. LST at the mean altitude of greatest sensitivity, 180 km, continuously over four martian years. Seasonal variability in densities is largely explained by orbital and latitudinal changes in dayside insolation that impacts the nightside through the resulting thermospheric circulation. However, the physical processes behind repeatable rapid, late autumnal cooling at mid-latitudes and near-aphelion warming at equatorial latitudes is not fully clear. Southern winter polar warming is generally weak or nonexistent over several Mars years, in basic agreement with MGS and MRO accelerometer observations. The puzzling response of mid-latitude densities from 160° to 200° E to the 2001 global dust storm suggests unanticipated localized nightside upper thermospheric lateral and vertical circulation patterns may accompany such storms. The downturn of the 11-year cycle of solar EUV flux is likely responsible for lower aphelion densities in 2004 and 2006 (Mars years 27 and 28).     

A sensitive search for SO2 in the martian atmosphere: Implications for seepage and origin of methane

Mars was observed near the peak of the strongest SO₂ band at 1364-1373 cm⁻¹ with resolving power of 77,000 using the Texas Echelon Cross Echelle Spectrograph on the NASA Infrared Telescope Facility. The observation covered the Tharsis volcano region which may be preferable to search for SO₂. The spectrum shows absorption lines of three CO₂ isotopomers and three H₂O isotopomers. The water vapor abundance derived from the HDO lines assuming D/H = 5.5 times the terrestrial value is 12±1.0 pr. μm, in agreement with the simultaneous MGS/TES observations of 14 pr. μm at the latitudes (50° S to 10° N) of our observation. Summing of spectral intervals at the expected positions of sixteen SO₂ lines puts a 2σ upper limit on SO₂ of 1 ppb. SO₂ may be emitted into the martian atmosphere by seepage and is removed by three-body reactions with OH and O. The SO₂ lifetime, 2 years, is longer than the global mixing time 0.5 year, so SO₂ should be rather uniformly distributed across Mars. Seepage of SO₂ is less than 15,000 tons per year on Mars which is smaller than the volcanic production of SO₂ on the Earth by a factor of 700. Because CH₄/SO₂ is typically 10⁻⁴-10⁻³ in volcanic gases on the Earth, our results show seepage is unlikely to be the source of the recently discovered methane on Mars and therefore strengthen its biogenic origin.     

An ensemble Kalman filter data assimilation system for the martian atmosphere: Implementation and simulation experiments

The local ensemble transform Kalman filter (LETKF) is applied to the GFDL Mars general circulation model (MGCM) to demonstrate the potential benefit of an advanced data assimilation method. In perfect model (aka identical twin) experiments, simulated observations are used to assess the performance of the LETKF-MGCM system and to determine the dependence of the assimilation on observational data coverage. Temperature retrievals are simulated at locations that mirror the spatial distribution of the Thermal Emission Spectrometer (TES) retrievals from the Mars Global Surveyor (MGS). The LETKF converges quickly and substantially reduces the analysis and subsequent forecast errors in both temperature and velocity fields, even though only temperature observations are assimilated. The LETKF is also found to accurately estimate the magnitude of forecast uncertainties, notably those associated with the phase and amplitude of baroclinic waves along the boundary of the polar ice cap during Northern Hemisphere winter.