Introduction: Atmospheric dust is ubiquitous on Mars, and as a result of its absorption and scattering of radiation, is the key driver of the martian circulation. Accurately representing the complex spatial and temporal distribution of dust is therefore crucial for understanding Mars’ atmospheric dynamics. In particular, the vertical representation of the dust distribution in Mars’ atmosphere has been shown to have a significant effect on results from modelling and assimilation [1,2,3]. With the goal of more accurately representing this distribution, the assimilation of dust vertical information is a valuable technique which is being increasingly explored [4,5]. However, it brings with it its own challenges and methodological questions to be explored.
Model and assimilation details: We use the LMD-UK Mars Global Circulation Model (MGCM), which solves the meteorological primitive equations of fluid dynamics, radiative and other parameterised physics to calculate the state of the martian atmosphere [6,7]. The UK version of the MGCM possesses a spectral dynamical core and semi-Lagrangian advection scheme [8], and is a collaboration between the Laboratoire de Météorologie Dynamique, The Open University, the University of Oxford, and the Instituto de Astrofisica de Andalucia. The model was run using a range of spectral and vertical resolutions, the latter spaced logarithmically. The assimilation scheme used was a modified version of the Analysis Correction scheme developed at the Met Office [9], adapted for use on Mars [10]. This method has the advantage of being computationally in-expensive, and its use of repeated insertion, weighted over a time window of about six hours, helps counter the issue of relaxation of the atmospheric state – an especially significant problem given the low thermal inertia of Mars’ atmosphere.
Retrievals: The retrievals used in this study are from the Mars Climate Sounder (MCS) instrument aboard the Mars Reconnaissance Orbiter (MRO) [11], which now has amassed over five full martian years’ worth of data. For this study, the assimilated MCS variables were temperature and dust profiles. Temperature profiles extend from the surface to approximately 100 km, and dust profiles from as low as 10 km above the surface up to a maximum height of approximately 50 km. Retrieval of dust profiles allows MCS to observe the complex vertical dust structure in the atmosphere. The retrieval version used is 5.2, a re-processing using updated 2D geometry [12]. This results in improved retrievals, especially in the polar regions.
While not used in this study, the NOMAD instrument aboard ExoMars TGO will soon provide another high-volume source of dust profiles alongside MCS [13], and should return observations with an even higher vertical resolution.
Discussion: The assimilation of MCS dust profiles poses unique technical challenges, but presents the opportunity of representing Mars’ vertical dust distribution with unprecedented spatial and temporal accuracy within a GCM. Some outstanding questions for further experimentation and discussion include:
What are the optimal spatial and, in particular, vertical model resolutions for assimilation of this data?
Can dust profile assimilation aid in forecasting? Previous indirect assimilation of vertical dust via its MCS temperature signature has yielded a forecast time of 10 sols [5]; how dependent is this on the assimilation scheme and the choice of assimilating variables?
How should we approach the bimodal nature of MCS local times? Should we give higher weighting to nightside dust observations, which tend to have better vertical coverage due to reduced scattering? And how much can we validly infer from the high day-night variability seen in MCS dust profiles?
What are the best heuristics for filtering spurious opacities which could disrupt the assimilation, for example due to CO2 ice or surface reflectance [16])?
What are the optimal ways of dealing with spatial and temporal gaps in the dataset?
How can we best represent the dust distribution beyond the range of MCS, especially in the lowest 5-10 km of the atmosphere?
What are the advantages and disadvantages of directly assimilating the dust field vs indirectly up-dating the dust field via its temperature signature, as seen in Fig. 1?
Dust profile assimilation has been used to track individual dust storm events [4]; what can this tell us about storm formation and evolution, and can it be used for storm forecasting?
How can we best constrain and validate the column optical depths of MCS dust profiles?
Some ways forward regarding these questions will be explored, including comparative reanalyses and validation against different orbital datasets. Comparisons against MCS and other retrievals (such as NOMAD) should provide insight into the advantages of various in-model representations of features such as the dust distribution as well as the possible advantages or disadvantages of pruning the assimilated dataset. Meanwhile, alternate orbital or even ground-based sources of column opacity (such as Mars Express and MSL) could help better con-strain the distribution of dust not seen by MCS and offer clues how best to proceed in periods when MCS data is missing or limited. Some results of intercomparisons will be presented with the aim of fostering a more general discussion on MCS assimilation techniques.
References: [1] Lewis, S. R. et al., Icarus 192 (2), 327-347, 2007. [2] Rogberg, P. et al., QJRMS 136, 1614-1635, 2010. [3] Greybush, S. J. et al., JGR. 117, E11008, 2012. [4] Ruan, T., DPhil Thesis, 2015. [5] Navarro T. et al., Earth and Space Sci., 2017. [6] Forget, F. et al., JGR 104, 24155-24175, 1999. [7] Madeleine, J.-B. et al., JGR (Planets) 116, E11010, 2011. [8] Newman, C. E. et al., JGR 107, 5123, 2002. [9] Lorenc, A. C. et al., QJRMS 117, 59-89, 1991. [10] Lewis, S. R. et al., Icarus 192, 327-347, 2007. [11] McCleese, D. J. et al., J. Geophys. Res. 115, E12016, 2010. [12] Kleinböhl, A. et al., J. Quant. Spectrosc. Radiat. Transfer 187, 511-522, 2017. [13] Patel, M. R. et al., Appl. Opt. 56 (10), 2771-2782, 2017. [14] Navarro, T. et al., Geophys. Res. Lett. 41, 6620-6626, 2014. [15] Streeter, P. M. et al., 6th Intl. Workshop on the Mars Atmosphere, 2017. [16] Kleinböhl, A. et al., Icarus 261, 118-121, 2015
