Surface release of radiatively active particles, with high infrared- (IR-)to-visible extinction ratios, has been proposed as a method of warming Mars. However, to warm Mars using aerosols, particles released locally must disperse globally. Here we provide an initial reference study in a plume tracking, dry Martian atmospheric model to address this question. The winds that transport aerosols respond to the aerosol's IR forcing, implying strong radiative-dynamical feedbacks (RDF). We investigate RDF from surface release of two particle compositions: carbon (graphene) and metal (Al). Self-lofting helps particles rise and spread locally and regionally, and the Hadley cell strengthens under warming, aiding latitudinal mixing. Within our model, Mars RDF enable engineered-aerosol warming. Warming is slightly greater for three-dimensional vs. 1D-models and also depends on spectral resolution of radiative transfer. We assess implications for Mars warming. Many open atmospheric science questions remain, including the role of agglomeration, dry-deposition rate uncertainty, and modeling water cycle feedbacks.

One-third of Mars’ surface has shallow-buried H 2 O, but it is currently too cold for use by life. Proposals to warm Mars using greenhouse gases require a large mass of ingredients that are rare on Mars’ surface. However, we show here that artificial aerosols made from materials that are readily available at Mars—for example, conductive nanorods that are ~9 micrometers long—could warm Mars >5 × 10 3 time smore effectively than the best gases. Such nanoparticles forward-scatter sunlight and efficiently block upwelling thermal infrared. Like the natural dust of Mars, they are swept high into Mars’ atmosphere, allowing delivery from the near-surface. For a 10-year particle lifetime, two climate models indicate that sustained release at 30 liters per second would globally warm Mars by ≳30 kelvin and start to melt the ice. Therefore, if nanoparticles can be made at scale on (or delivered to) Mars, then the barrier to warming of Mars appears to be less high than previously thought.

Latitude and topography affect the amount of shadow cast on a landscape, which in turn can influence where water, ice or snow are stable. For martian climate models, whose output sometimes disagrees with observational data, we show that adding the ability to represent shadows generated by topography and understanding how they interact with relevant modelled variables can be important. We included a shadows scheme in a Mars mesoscale climate model and report on the impact that shadows had on surface temperature, and the resultant impact on one variable – surface water ice, and its accumulation or ablation in the study area of Lyot crater. Incorporating a scheme simulating shadows cast by topographical features had a measurable effect, increasing by up to 83% the amount of surface ice that the model predicted would be found in shadowed areas. The largest differences in ice content between simulations with and without shadows were found in those areas with the greatest amount of shadow.

We report the ubiquitous occurrence of nighttime temperature inversions in the tropical martian atmosphere during the dusty season, as observed by the Mars Climate Sounder. The inversions are linked to the occurrence of large-scale regional dust storms, with their strengths largely correlated to the strengths of the dust storms. Inversions strengthen between 2 am and 4 am, with the bases of the inversions getting cooler, and the tops of the inversions getting warmer. The inversions are strongest around Tharsis and Terra Sabaea, which are higher-elevation regions, suggesting they are forming due to a combination of topographically-excited tides and cloud radiative cooling. However, inversions are also observed over the flat plains, and are likely associated with stronger tides resulting from the increased dust abundance. These results highlight an important interplay between the dust distribution, water ice clouds and thermal tides.

SignificanceMars is cold today but once had lakes. We evaluated the water ice cloud greenhouse hypothesis for warming early Mars. Our results reconcile previous discrepant results by showing that the cloud greenhouse provides strong warming if the surface has patchy surface H2O but not if there is very extensive surface H2O. In our model, arid, warm, stable climates emerge with surface H2O (and low clouds) only at locations much colder than average surface temperature. At locations horizontally distant from the surface cold traps, clouds are found only at high altitudes, which maximizes cloud warming. As this scenario is consistent with geologic data that suggest a warm, arid early Mars climate, our results support the cloud greenhouse hypothesis for warming early Mars.

We analyze Mars Climate Sounder temperature and aerosol data in the tropics to study atmospheric tides and their relation to the dust and water ice distributions. Our results from data covering Mars years (MY) 29-35 reveal that MY29 has large amplitude non-migrating thermal tides during southern hemisphere spring equinox that are not observed at the same local time in any other year. It is the nighttime temperatures that are most perturbed compared to other years, with strong temperature inversions at 35-55 km altitude. Analysis of data at different local times reveals that the temperatures and water ice clouds at 03:45 am in MY29 more closely resemble those at 05:00 am in other years, suggesting there was a shift in the phase of the diurnal tide to an earlier local time. This phase shift, and the large amplitude non-migrating thermal tides, appear to be related to early dust activity. Two early dust storms occurred in MY29 around the time there was upwelling over the tropics, associated with the Hadley circulation, enabling the dust to be transported to higher altitudes where it has a larger radiative influence. As well as dust, water ice clouds are also found to influence the tidal structure. Due to the interaction of non-migrating tides, water ice clouds occur in two discrete longitudinal regions at night. The increased radiative cooling results in increased downwelling above the clouds, leading to increased adiabatic warming and a strengthening the temperature inversions.

Abstract We reanalyse the global distribution of gullies in order to provide a set of observational constraints that models of gully formation must explain. We validate our results derived from the global data with four detailed case studies. We show that the availability of steep slopes is an essential factor to consider when assessing the spatial distribution and abundance of gullies. When the availability of steep slopes is taken into account, it reveals, with a few exceptions, that gullies are found almost uniformly across the whole 30°–90° latitude band. Our analysis also reveals that massive ice deposits are anti-correlated with gullies, and that the undulations in the equatorwards limits of the gully distribution could be explained by longitudinal variations in maximum surface temperatures (controlled by variations in surface properties, including thermal inertia and albedo). We find a sharp transition in both hemispheres between pole-facing gullies, which extend from 30° to 40°, to a more mixed, but dominantly equator-facing orientation of gullies polewards of 40°. We have no definitive explanation for this transition but, based on previous studies, we suggest it could be linked to the availability of near-surface ice deposits.

Detailed modelling of the Martian atmosphere is completed for every spacecraft designed to land on the planet’s surface. This provides the most complete picture of the environment that the descending module will be entering and travelling through, and facilitates planning of the Entry, Descent and Landing (EDL) phase of the mission. The selected resolution of an atmospheric model can impact the results of the experiments performed. The complexities of atmospheric modelling also require models of different scales to best represent the behaviour of different scale atmospheric phenomena. Comparisons between multiple model results and in situ data are crucial for improving future environmental predictions for missions landing on Mars. This work describes how changes in model scale and resolution (horizontal and vertical) can impact experimental results, using as a case study the selected landing site of the European Space Agency (ESA) Schiaparelli module. Schiaparelli was part of ESA’s ExoMars 2016 mission; the module descended through the Martian atmosphere on 19th October 2016. Experiments were completed that encompassed the period of Schiaparelli’s descent, using both a global-scale and a mesoscale model. The global model used in this work is the UK version of the LMD (Laboratoire de Meteorologie Dynamique) Mars Glob-al Circulation Model (“the MGCM”), a 3D multi-level spectral model of the Martian atmosphere up to an altitude of ~100 km [1]. The mesoscale model used in this work is the LMD Martian Mesoscale Model (MMM) [2]; in these experiments an altitude of ~50 km was modelled in the mesoscale. Multiple resolution experiments were completed using the MGCM; results range from a ‘low’ resolution ~5° latitude x ~5° longitude (a resolution typically used for Martian climate modelling) to a ‘high’ resolution ~1° lat x ~1° lon. The vertical dimension is modelled using a set number of vertical layers; in these experiments the number of vertical layers selected was between 23 and 100. Experiments were run for a simulated year, starting from initial conditions based upon prior atmospheric observations, thus providing an independent prediction of conditions through the period of this case study. The MMM experiments were com-pleted in a set of nested resolutions, ranging from the outer, lowest resolution results at 63 km x 63 km, to the inner, highest resolution results at 7 km x 7 km. MMM experiments were completed using 60 vertical layers. Previous comparisons of global-scale and meso-scale modelling have focused on areas containing small-scale topographical variation that is not present in the global scale models. This work considers the relatively flat topography of the Schiaparelli site – a location that is more representative of the majority of historical Martian landing sites than areas that contain severe, small-scale topographical variation. Initial analysis has focused on constructing vertical profiles from the model output at both experimental scales, following preliminary information on the descent trajectory of the Schiaparelli module. An example comparison of atmospheric profiles constructed from MGCM results at different model resolutions. The plot displays atmospheric temperature obtained from experiments completed at different vertical resolutions: 23 and 100 vertical levels. There is a good match between the results, with a root mean square deviation (RMSD) of 9.83 K be-tween the results for the full height of the profiles; the RMSD reduces to 2.05 K when considering only the lowest ~10 km of the profiles (approximately one scale height). A comparison of vertical atmospheric temperature profiles from MGCM and MMM results. While the trend in the results is similar, the results differ by ~10 K between the models through most of the profile, down to a height of ~3 km above the surface. Between 50 and 3 km above the surface the RMSD of the profiles is 9.79 K; below 3 km (down to the lowest MGCM model lay-er) the match is closer, with an RMSD of 2.59 K. Further comparisons have been completed between the MGCM and MMM results, such as wind speed and direction, including consideration of the wider topographical and atmospheric context of Schiaparelli’s landing site and EDL period. These results show that, for the region considered within this case study, changing the horizontal or verti-cal resolution used in MGCM experiments does not greatly impact the results obtained. Similarly, the MMM results do not vary more than ~4 K with chang-ing horizontal resolution. In both cases, lower resolu-tions results (which are quicker and less computationally expensive to complete) are a good approximation of higher resolution results. Additionally, the similarity of the trends seen in the results from the different scale models suggests that global-scale model results are a reasonable approximation for mesoscale model results, for a number of potential landing locations on Mars. The module successfully transmitted some data that was captured during its descent, primarily from engineering sensors; this data includes the module's trajectory and attitude during the mission’s EDL phase. The ExoMars AMELIA (Atmospheric Mars Entry and Landing Investigations and Analysis) team aim to use the data returned by Schiaparelli during descent, combined with dynamic modelling of the module's motion, to reconstruct atmospheric profiles of density, pressure, temperature and wind speed [3]. Upon the release of the Schiaparelli data, the results from both the MGCM and MMM experiments will be compared with the data, supporting the work of the AMELIA team. References: [1] Forget et al. (1999) JGR, 104, E10. [2] Spiga et al. (2009) JGR, 114, E2. [3] Ferri et al. (2012) 9th International Planetary Probe Workshop.

Most of Mars' ancient sedimentary rocks by volume are in wind‐eroded sedimentary mounds within impact craters and canyons, but the connections between mound form and wind erosion are unclear. We perform mesoscale simulations of different crater and mound morphologies to understand the formation of sedimentary mounds. As crater depth increases, slope winds produce increased erosion near the base of the crater wall, forming mounds. Peak erosion rates occur when the crater depth is ∼2 km. Mound evolution depends on the size of the host crater. In smaller craters mounds preferentially erode at the top, becoming more squat, while in larger craters mounds become steeper sided. This agrees with observations where smaller craters tend to have proportionally shorter mounds and larger craters have mounds encircled by moats. If a large‐scale sedimentary layer blankets a crater, then as the layer recedes across the crater it will erode more toward the edges of the crater, resulting in a crescent‐shaped moat. When a 160 km diameter mound‐hosting crater is subject to a prevailing wind, the surface wind stress is stronger on the leeward side than on the windward side. This results in the center of the mound appearing to “march upwind” over time and forming a “bat‐wing” shape, as is observed for Mount Sharp in Gale crater.

Most of Mars’ ancient sedimentary rocks by volume are in wind-eroded sedimentary mounds within impact craters and canyons, but the connections between mound form and wind erosion are unclear. We perform mesoscale simulations of different crater and mound morphologies to understand the formation of sedimentary mounds. As crater depth increases, slope winds produce increased erosion near the base of the crater wall, forming mounds. Peak erosion rates occur when the crater depth is ∼2 km. Mound evolution depends on the size of the host crater. In smaller craters mounds preferentially erode at the top, becoming more squat, while in larger craters mounds become steeper sided. This agrees with observations where smaller craters tend to have proportionally shorter mounds and larger craters have mounds encircled by moats. If a large-scale sedimentary layer blankets a crater, then as the layer recedes across the crater it will erode more toward the edges of the crater, resulting in a crescent-shaped moat. When a 160 km diameter mound-hosting crater is subject to a prevailing wind, the surface wind stress is stronger on the leeward side than on the windward side. This results in the center of the mound appearing to “march upwind” over time and forming a “bat-wing” shape, as is observed for Mount Sharp in Gale crater.

We investigate the exchange of water vapour between the regolith and atmosphere of Mars, and how it varies with different orbital parameters, atmospheric dust contents and surface water ice reservoirs. This is achieved through the coupling of a global circulation model (GCM) and a regolith diffusion model. GCM simulations are performed for hundreds of Mars years, with additional one-dimensional simulations performed for 50 kyr. At obliquities ε = 15◦ and 30◦, the thermal inertia and albedo of the regolith have more control on the subsurface water distribution than changes to the eccentricity or solar longitude of perihelion. At ε = 45◦, atmospheric water vapour abundances become much larger, allowing stable subsurface ice to form in the tropics and mid-latitudes. The circulation of the atmosphere is important in producing the subsurface water distribution, with increased water content in various locations due to vapour transport by topographically-steered flows and stationary waves. As these circulation patterns are due to topographic features, it is likely the same regions will also experience locally large amounts of subsurface water at different epochs. The dustiness of the atmosphere plays an important role in the distribution of subsurface water, with a dusty atmosphere resulting in a wetter water cycle and increased stability of subsurface ice deposits.

The formation of large dust storms on Mars is believed to be driven by dust lifting due to near-surface wind stress (NSWS). Accurately representing this dust lifting within Mars Global Circulation Models (MGCMs) is important in order to gain a full understanding of the Martian dust storm cycle and climatology. Model resolutions of ~ 5°x5° in latitude and longitude are often used in simulations of the Martian climate. This resolution represents large-scale weather patterns well, but small scale phenomena (e.g. near-surface winds) are not accurately depicted. Parameterisations of dust lifting by NSWS exist within several MGCMs, but few studies have explored in detail how the results of these parameterisations are affected by changing the resolution of the model. We use the LMD-UK MGCM to complete simulations across multiple model resolutions. Our experiments range from ‘low’ resolution ~5° lat x ~5° lon to ‘high’ resolution ~1° lat x ~1° lon. In experiments with fixed lifting parameters, we find that increasing the model’s horizontal resolution results in more dust being lifted; we also find that increasing the model’s number of vertical layers results in more dust being lifted. Geographical and temporal distributions of dust lifting are investigated, with the aim of developing a calibration scheme for this dust lifting across model resolutions. The scheme is verified through comparison with spacecraft observations of dust optical depths and dust storm locations.