Search

Global dust storms have major implications for the past and present climate, geologic history, habitability, and future exploration of Mars. Yet their mysterious origins mean we remain unable to realistically simulate or predict them. We identify four key Knowledge Gaps and make four Recommendations to make progress in the next decade.

A Martian semiannual oscillation (SAO), similar to that in the Earth's tropical stratosphere, is evident in the Mars Analysis Correction Data Assimilation reanalysis dataset (MACDA) version 1.0, not only in the tropics, but also extending to higher latitudes. Unlike on Earth, the Martian SAO is found not always to reverse its zonal wind direction, but only manifests itself as a deceleration of the dominant wind at certain pressure levels and latitudes. Singular System Analysis (SSA) is further applied on the zonal-mean zonal wind in different latitude bands to reveal the characteristics of SAO phenomena at different latitudes. The second pair of principal components (PCs) is usually dominated by a SAO signal, though the SAO signal can be strong enough to manifest itself also in the first pair of PCs. An analysis of terms in the Transformed Eulerian Mean equation (TEM) is applied in the tropics to further elucidate the forcing processes driving the tendency of the zonal-mean zonal wind. The zonal-mean meridional advection is found to correlate strongly with the observed oscillations of zonal-mean zonal wind, and supplies the majority of the westward (retrograde) forcing in the SAO cycle. The forcing due to various non-zonal waves supplies forcing to the zonal-mean zonal wind that is nearly the opposite of the forcing due to meridional advection above ∼3 Pa altitude, but it also partly supports the SAO between 40 Pa and 3 Pa. Some distinctive features occurring during the period of the Mars year (MY) 25 global-scale dust storm (GDS) are also notable in our diagnostic results with substantially stronger values of eastward and westward momentum in the second half of MY 25 and stronger forcing due to vertical advection, transient waves and thermal tides.

We have reconstructed the climatology of airborne dust from Martian years (MY) 24 to 31 using multiple datasets of retrieved or estimated column optical depth. The datasets are based on observations of the Martian atmosphere from March 1999 to July 2013 by different orbiting instruments: the Thermal Emission Spectrometer (TES) on board Mars Global Surveyor, the Thermal Emission Imaging System (THEMIS) on board Mars Odyssey, and the Mars Climate Sounder (MCS) on board Mars Reconnaissance Orbiter (MRO). The procedure we have adopted consists in gridding the available retrievals of column dust optical depth (CDOD) from TES and THEMIS nadir observations, as well as the estimates of this quantity from MCS limb observations. Our gridding method calculates weighted averages on a regular but likely incomplete spatial grid, using an iterative procedure with weights in space, time, and retrieval uncertainty. The derived product consists of daily synoptic gridded maps of CDOD at a resolution of 6 degree longitude x 3 degree latitude for MY 24-26, and 6 degree longitude x 5 degree latitude for MY 27-31. We have statistically analyzed the gridded maps to present an overview of the dust climatology on Mars over eight years, specifically in relation to its intraseasonal and interannual variability. Finally, we have produced complete daily maps of CDOD by spatially interpolating the available incomplete gridded maps using a kriging method. These complete maps are used as dust scenarios in the Mars Climate Database (MCD) version 5, and should be useful for many other applications. The maps for the eight available Martian years are publicly available and distributed with open access, under Creative Commons Attribution-ShareAlike 3.0 Unported License. The current version and future updates can be downloaded from the MCD website at the Laboratoire de Meteorologie Dynamique: http://wwwmars. lmd.jussieu.fr/mars/dust_climatology

For many years, several teams around the world have developed GCMs (General Circulation Model or Global Climate Model) to simulate the environment on Mars. The GCM developed at the Laboratoire de Météorologie Dynamique in collaboration with several teams in Europe (LATMOS, France, University of Oxford, The Open University, the Instituto de Astrofisica de Andalucia), and with the support of ESA and CNES is currently used for many applications. Its outputs have also regularly been compiled to build a Mars Climate Database, a freely available tool useful for the scientific and engineering communities. The Mars Climate Database (MCD) has over the years been distributed to more than 150 teams around the world. Following the recent improvements in the GCM, a new series of reference simulations have been run and compiled into a new version (version5.1) of the Mars Climate Database, released in the first half of 2014. To summarize, MCD v5.1 provides: - Climatologies over a series of dust scenarios: standard year, cold (ie: low dust), warm (ie: dusty atmosphere) and dust storm, all topped by various cases of Extreme UV solar inputs (low, mean or maximum). These scenarios differ from those of previous versions of the MCD (version 4.x) as they have been derived from home-made, instrument-derived (TES, THEMIS, MCS, MERs), dust climatology of the last 8 Martian years. - Mean values and statistics of main meteorological variables (atmospheric temperature, density, pressure and winds), as well as surface pressure and temperature, CO2 ice cover, thermal and solar radiative fluxes, dust column opacity and mixing ratio, [H20] vapor and ice columns, concentrations of many species: [CO], [O2], [O], [N2], [H2], [O3], ... - A high resolution mode which combines high resolution (32 pixel/degree) MOLA topography records and Viking Lander 1 pressure records with raw lower resolution GCM results to yield, within the restriction of

Large-scale planetary waves are diagnosed from an analysis of profiles retrieved from the Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft during its scientific mapping phase. The analysis is conducted by assimilating thermal profiles and total dust opacity retrievals into a Mars global circulation model. Transient waves are largest throughout the northern hemisphere autumn, winter and spring period and almost absent during the summer. The southern hemisphere exhibits generally weaker transient wave behaviour. A striking feature of the low-altitude transient waves in the analysis is that they show a broad subsidiary minimum in amplitude centred on the winter solstice, a period when the thermal contrast between the summer hemisphere and the winter pole is strongest and baroclinic wave activity might be expected to be strong. This behaviour, here called the ‘solsticial pause,’ is present in every year of the analysis. This strong pause is under-represented in many independent model experiments, which tend to produce relatively uniform baroclinic wave activity throughout the winter. This paper documents and diagnoses the transient wave solsticial pause found in the analysis; a companion paper investigates the origin of the phenomenon in a series of model experiments.

A meteorological data assimilation system has been developed recently for analyzing measurements of temperature and dust opacity on Mars and has been successfully applied in several studies (e.g. Montabone et al. 2005, Lewis et al. 2007) to study various atmospheric phenomena. A more sophisticated data assimilation system, now with full dust transport incorporated, is becoming available to represent more accurately and realistically the physical transport of dust.

Ever-increasing numbers of atmospheric observations from orbiting spacecraft, and increasingly sophisticated numerical atmospheric models, have recently permitted data assimilation techniques to be applied to planets beyond the Earth. Mars is the first extra-terrestrial planet for which reanalyses of the atmospheric state are now available. The Thermal Emission Spectrometer (TES) on board NASA’s Mars Global Surveyor (MGS) has produced an extensive atmospheric data set during its scientific mapping phase between 1999 and 2004. Nadir thermal profiles for the atmosphere below about 40 km altitude, and total dust and water ice opacities, have been retrieved from TES spectra, covering almost three complete Martian seasonal cycles (each seasonal cycle on Mars corresponds to 668.6 mean solar days, and the Martian mean solar day is about 24 hours and 40 minutes). Note that dust on Mars plays a key role in the weather and climate, mainly through its strong absorption of short wave radiation with a short radiative relaxation timescale of 1-2 days. Assimilating dust opacities correctly is, therefore, particularly important for atmospheric data assimilation on the Red Planet. TES retrieved observations have been analysed by assimilation into a Mars general circulation model (MGCM), making use of a sequential procedure known as the Analysis Correction scheme, a form of successive corrections method which has proved simple and robust under Martian conditions, even during the less-than-ideal MGS aerobraking period. The MGCM used at the University of Oxford and at The Open University consists of a spectral dynamical solver and a tracer transport scheme developed in the UK. Its package of state-of-the-art physical parameterization routines is shared with the LMD-MarsGCM, developed by the Laboratoire de Météorologie Dynamique in Paris (France). One limitation of TES is that relatively few limb profiles are available, compared to nadir soundings. Our MGS/TES reanalysis, therefore, doe

Introduction: There is considerable evidence that shows that the Martian atmosphere behaves in a more regular fashion than its terrestrial counterpart [1, 2, 3, 4]. This evidence leads to the hypothesis of theMartian climate attractor being of a relatively low dimension, which, in turn, would imply the possibility of describing the state of the atmosphere by means of a relatively few degrees of freedom. We explore this hypothesis by assuming that the atmospheric total energy (TE), i.e. the sum of kinetic energy and total potential energy (gravitational potential energy plus internal energy), is confined in a few coherent structures which dynamically interact nonlinearly with each other.

From the introduction: Given the quantity of data expected from current and forthcoming spacecraft missions to Mars, it is now possible to use data assimilation as a means of atmospheric analysis for the first time for a planet other than the Earth. Several groups have described plans to develop assimilation schemes for Mars [Banfield et al., 1995; Houben, 1999; Lewis and Read, 1995; Lewis et al., 1996, 1997; Zhang et al., 2001]. Data assimilation is a technique for the analysis of atmospheric observations which combines currently valid information with prior knowledge from previous observations and dynamical and physical constraints, via the use of a numerical model. Despite the number of new missions, observations of the atmosphere of Mars in the near future are still likely to be sparse when compared to those of the Earth, perhaps comprising one orbiter and a few surface stations at best at any one time. Data assimilation is useful as a means to extract the maximum information from such observations, both by a form of interpolation in space and time using model constraints and by the combination of information from different observations, e.g. temperature profiles and surface pressure measurements which may be irregularly distributed. The procedure can produce a dynamically consistent set of meteorological fields and can be used directly to test and to refine an atmospheric model against observations.

Our Global Climate Model (GCM) of the Martian atmosphere is the result of twenty years of ongoing collaboration between our teams and has matured to the point of enabling to study the main cycles (dust, CO2, water) of present-day and past Martian climates. At the 2014 scientific assembly, we will report on the latest developments and improvements of our GCM, and also present the latest version of the Mars Climate Database (version 5.1) that is derived from GCM outputs, along with comparisons with available measurements (from TES, MCS, Viking, Phoenix, Curiosity, etc.).

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.

We calculate a Lorenz energy budget for the Martian atmosphere from reanalysis derived from Mars Global Surveyor data for Mars years 24–27. We present global, annual mean energy and conversion rates per unit area and per unit mass and compare these to Earth data. The directions of the energy conversion terms for Mars are similar to Earth, with the exception of the barotropic conversion between zonal and eddy kinetic energy reservoirs. Further, seasonal and hemispheric decomposition reveals a strong conversion between zonal energy reservoirs over the year, but these balance each other out in global and annual mean. On separating the diurnal timescale, the contribution to the conversion terms and eddy kinetic energy for diurnal and shorter timescales in many cases (especially during planet-encircling dust storms) exceeds the contribution of longer timescales. This suggests that thermal tides have a significant effect on the generation of eddy kinetic energy.

A Martian semiannual oscillation (SAO), similar to that in the Earth's tropical stratosphere, is evident in the Mars Analysis Correction Data Assimilation reanalysis dataset (MACDA) version 1.0, not only in the tropics, but also extending to higher latitudes. Unlike on Earth, the Martian SAO is found not always to reverse its zonal wind direction, but only manifests itself as a deceleration of the dominant wind at certain pressure levels and latitudes. Singular System Analysis (SSA) is further applied on the zonal-mean zonal wind in different latitude bands to reveal the characteristics of SAO phenomena at different latitudes. The second pair of principal components (PCs) is usually dominated by a SAO signal, though the SAO signal can be strong enough to manifest itself also in the first pair of PCs. An analysis of terms in the Transformed Eulerian Mean equation (TEM) is applied in the tropics to further elucidate the forcing processes driving the tendency of the zonal-mean zonal wind. The zonal-mean meridional advection is found to correlate strongly with the observed oscillations of zonal-mean zonal wind, and supplies the majority of the westward (retrograde) forcing in the SAO cycle. The forcing due to various non-zonal waves supplies forcing to the zonal-mean zonal wind that is nearly the opposite of the forcing due to meridional advection above ∼3 Pa altitude, but it also partly supports the SAO between 40 Pa and 3 Pa. Some distinctive features occurring during the period of the Mars year (MY) 25 global-scale dust storm (GDS) are also notable in our diagnostic results with substantially stronger values of eastward and westward momentum in the second half of MY 25 and stronger forcing due to vertical advection, transient waves and thermal tides.

We calculate a Lorenz energy budget for the Martian atmosphere from reanalysis derived from Mars Global Surveyor data for Mars years 24–27. We present global, annual mean energy and conversion rates per unit area and per unit mass and compare these to Earth data. The directions of the energy conversion terms for Mars are similar to Earth, with the exception of the barotropic conversion between zonal and eddy kinetic energy reservoirs. Further, seasonal and hemispheric decomposition reveals a strong conversion between zonal energy reservoirs over the year, but these balance each other out in global and annual mean. On separating the diurnal timescale, the contribution to the conversion terms and eddy kinetic energy for diurnal and shorter timescales in many cases (especially during planet-encircling dust storms) exceeds the contribution of longer timescales. This suggests that thermal tides have a significant effect on the generation of eddy kinetic energy.

For many years, several teams around the world have developed GCMs (General Circulation Model or Global Climate Model) to simulate the environment on Mars. The GCM developed at the Laboratoire de Météorologie Dynamique in collaboration with several teams in Europe (LATMOS, France, University of Oxford, The Open University, the Instituto de Astrofisica de Andalucia), and with the support of ESA and CNES is currently used for many applications. Its outputs have also regularly been compiled to build a Mars Climate Database, a freely available tool useful for the scientific and engineering communities. The Mars Climate Database (MCD) has over the years been distributed to more than 150 teams around the world. Following the recent improvements in the GCM, a new series of reference simulations have been run and compiled into a new version (version5.1) of the Mars Climate Database, released in the first half of 2014. To summarize, MCD v5.1 provides: - Climatologies over a series of dust scenarios: standard year, cold (ie: low dust), warm (ie: dusty atmosphere) and dust storm, all topped by various cases of Extreme UV solar inputs (low, mean or maximum). These scenarios differ from those of previous versions of the MCD (version 4.x) as they have been derived from home-made, instrument-derived (TES, THEMIS, MCS, MERs), dust climatology of the last 8 Martian years. - Mean values and statistics of main meteorological variables (atmospheric temperature, density, pressure and winds), as well as surface pressure and temperature, CO2 ice cover, thermal and solar radiative fluxes, dust column opacity and mixing ratio, [H20] vapor and ice columns, concentrations of many species: [CO], [O2], [O], [N2], [H2], [O3], ... - A high resolution mode which combines high resolution (32 pixel/degree) MOLA topography records and Viking Lander 1 pressure records with raw lower resolution GCM results to yield, within the restriction of

Our Global Climate Model (GCM) of the Martian atmosphere is the result of twenty years of ongoing collaboration between our teams and has matured to the point of enabling to study the main cycles (dust, CO2, water) of present-day and past Martian climates. At the 2014 scientific assembly, we will report on the latest developments and improvements of our GCM, and also present the latest version of the Mars Climate Database (version 5.1) that is derived from GCM outputs, along with comparisons with available measurements (from TES, MCS, Viking, Phoenix, Curiosity, etc.).

We have reconstructed the climatology of airborne dust from Martian years (MY) 24 to 31 using multiple datasets of retrieved or estimated column optical depth. The datasets are based on observations of the Martian atmosphere from March 1999 to July 2013 by different orbiting instruments: the Thermal Emission Spectrometer (TES) on board Mars Global Surveyor, the Thermal Emission Imaging System (THEMIS) on board Mars Odyssey, and the Mars Climate Sounder (MCS) on board Mars Reconnaissance Orbiter (MRO). The procedure we have adopted consists in gridding the available retrievals of column dust optical depth (CDOD) from TES and THEMIS nadir observations, as well as the estimates of this quantity from MCS limb observations. Our gridding method calculates weighted averages on a regular but likely incomplete spatial grid, using an iterative procedure with weights in space, time, and retrieval uncertainty. The derived product consists of daily synoptic gridded maps of CDOD at a resolution of 6 degree longitude x 3 degree latitude for MY 24-26, and 6 degree longitude x 5 degree latitude for MY 27-31. We have statistically analyzed the gridded maps to present an overview of the dust climatology on Mars over eight years, specifically in relation to its intraseasonal and interannual variability. Finally, we have produced complete daily maps of CDOD by spatially interpolating the available incomplete gridded maps using a kriging method. These complete maps are used as dust scenarios in the Mars Climate Database (MCD) version 5, and should be useful for many other applications. The maps for the eight available Martian years are publicly available and distributed with open access, under Creative Commons Attribution-ShareAlike 3.0 Unported License. The current version and future updates can be downloaded from the MCD website at the Laboratoire de Meteorologie Dynamique: http://wwwmars. lmd.jussieu.fr/mars/dust_climatology

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.

From the introduction: Given the quantity of data expected from current and forthcoming spacecraft missions to Mars, it is now possible to use data assimilation as a means of atmospheric analysis for the first time for a planet other than the Earth. Several groups have described plans to develop assimilation schemes for Mars [Banfield et al., 1995; Houben, 1999; Lewis and Read, 1995; Lewis et al., 1996, 1997; Zhang et al., 2001]. Data assimilation is a technique for the analysis of atmospheric observations which combines currently valid information with prior knowledge from previous observations and dynamical and physical constraints, via the use of a numerical model. Despite the number of new missions, observations of the atmosphere of Mars in the near future are still likely to be sparse when compared to those of the Earth, perhaps comprising one orbiter and a few surface stations at best at any one time. Data assimilation is useful as a means to extract the maximum information from such observations, both by a form of interpolation in space and time using model constraints and by the combination of information from different observations, e.g. temperature profiles and surface pressure measurements which may be irregularly distributed. The procedure can produce a dynamically consistent set of meteorological fields and can be used directly to test and to refine an atmospheric model against observations.

A meteorological data assimilation system has been developed recently for analyzing measurements of temperature and dust opacity on Mars and has been successfully applied in several studies (e.g. Montabone et al. 2005, Lewis et al. 2007) to study various atmospheric phenomena. A more sophisticated data assimilation system, now with full dust transport incorporated, is becoming available to represent more accurately and realistically the physical transport of dust.

Ever-increasing numbers of atmospheric observations from orbiting spacecraft, and increasingly sophisticated numerical atmospheric models, have recently permitted data assimilation techniques to be applied to planets beyond the Earth. Mars is the first extra-terrestrial planet for which reanalyses of the atmospheric state are now available. The Thermal Emission Spectrometer (TES) on board NASA’s Mars Global Surveyor (MGS) has produced an extensive atmospheric data set during its scientific mapping phase between 1999 and 2004. Nadir thermal profiles for the atmosphere below about 40 km altitude, and total dust and water ice opacities, have been retrieved from TES spectra, covering almost three complete Martian seasonal cycles (each seasonal cycle on Mars corresponds to 668.6 mean solar days, and the Martian mean solar day is about 24 hours and 40 minutes). Note that dust on Mars plays a key role in the weather and climate, mainly through its strong absorption of short wave radiation with a short radiative relaxation timescale of 1-2 days. Assimilating dust opacities correctly is, therefore, particularly important for atmospheric data assimilation on the Red Planet. TES retrieved observations have been analysed by assimilation into a Mars general circulation model (MGCM), making use of a sequential procedure known as the Analysis Correction scheme, a form of successive corrections method which has proved simple and robust under Martian conditions, even during the less-than-ideal MGS aerobraking period. The MGCM used at the University of Oxford and at The Open University consists of a spectral dynamical solver and a tracer transport scheme developed in the UK. Its package of state-of-the-art physical parameterization routines is shared with the LMD-MarsGCM, developed by the Laboratoire de Météorologie Dynamique in Paris (France). One limitation of TES is that relatively few limb profiles are available, compared to nadir soundings. Our MGS/TES reanalysis, therefore, doe

Introduction: There is considerable evidence that shows that the Martian atmosphere behaves in a more regular fashion than its terrestrial counterpart [1, 2, 3, 4]. This evidence leads to the hypothesis of theMartian climate attractor being of a relatively low dimension, which, in turn, would imply the possibility of describing the state of the atmosphere by means of a relatively few degrees of freedom. We explore this hypothesis by assuming that the atmospheric total energy (TE), i.e. the sum of kinetic energy and total potential energy (gravitational potential energy plus internal energy), is confined in a few coherent structures which dynamically interact nonlinearly with each other.

Global dust storms have major implications for the past and present climate, geologic history, habitability, and future exploration of Mars. Yet their mysterious origins mean we remain unable to realistically simulate or predict them. We identify four key Knowledge Gaps and make four Recommendations to make progress in the next decade.