The Modern Near-Surface Martian Climate: A Review of In-situ Meteorological Data from Viking to Curiosity

We analyze the complete set of in-situ meteorological data obtained from the Viking landers in the 1970s to today’s Curiosity rover to review our understanding of the modern near-surface climate of Mars, with focus on the dust, CO2 and H2O cycles and their impact on the radiative and thermodynamic conditions near the surface. In particular, we provide values of the highest confidence possible for atmospheric opacity, atmospheric pressure, near-surface air temperature, ground temperature, near-surface wind speed and direction, and near-surface air relative humidity and water vapor content. Then, we study the diurnal, seasonal and interannual variability of these quantities over a span of more than twenty Martian years. Finally, we propose measurements to improve our understanding of the Martian dust and H2O cycles, and discuss the potential for liquid water formation under Mars’ present day conditions and its implications for future Mars missions. Understanding the modern Martian climate is important to determine if Mars could have the conditions to support life and to prepare for future human exploration.

Map with the large asymmetry in topography and the landing sites of spacecraft that landed successfully on Mars and infrared radiation, regulating atmospheric and surface heating rates and thus influencing the Martian climate (Gierasch and Goody 1972;Kahn et al. 1992;Forget et al. 1999;Basu et al. 2004;Madeleine et al. 2011). The Martian regolith forms a loosely packed, porous medium, which is the source of the ubiquitous airborne dust and allows the exchange of volatiles with the atmosphere. The large spatial and temporal variability in atmospheric dust content creates complex feedbacks connecting dust lifting and transport with the evolving atmospheric circulations (Newman et al. 2002a(Newman et al. , 2002bBasu et al. 2004;Kahre et al. 2006;Madeleine et al. 2011). As an example of such feedback, Martian global dust storms occur every few years, enshrouding the entire planet with dust in just a few weeks. On a global scale, the most vigorous dust activity and the resulting largest dust opacity values occur during southern spring and summer (Martin 1986;Smith 2004Smith , 2008Montabone et al. 2015), when regional dust storms persist throughout the southern hemisphere and the northern baroclinic zone. Smaller regional dust events also persist throughout the northern spring and summer, but they are less efficient at increasing the overall atmospheric dust opacity than the southern spring and summer events (Cantor et al. 2001). The southern hemisphere surface is on average around 5 km higher than the northern hemisphere surface (e.g.,   (Fig. 1), and this asymmetry in topography alone would produce a more intense Hadley circulation in southern spring and summer (Joshi et al. 1995;Richardson and Wilson 2002). In combination with the asymmetry in solar insolation, which also favors this time of year, the result is a far stronger overturning circulation and stronger surface wind stresses (and hence dust lifting from the surface) during the so-called 'dust storm season,' which runs from L s ∼ 180-330°(see also Sect. 3.1).
Mars also has a strong seasonal H 2 O cycle by which water is exchanged between the surface and the atmosphere through various processes, and then redistributed in the atmosphere by atmospheric circulations (Jakosky and Farmer 1982;Jakosky 1985;Jakosky and Haberle 1992;Houben et al. 1997;Richardson and Wilson 2002;Mischna et al. 2003). Seasonal water reservoirs include water ice on or within the seasonal or residual polar caps, surface or subsurface ice in high-latitude regions, physically adsorbed and chemically-bound water within the near-surface regolith, and possibly surface and subsurface liquid water in the form of liquid saline solutions (brine) (Jakosky 1985). Satellite and telescopic observations of variations in the integrated H 2 O column abundance (Conrath et al. 1973;Jakosky and Farmer 1982;Sprague et al. 1996;Smith 2004;Fedorova et al. 2006;Fouchet et al. 2007;Melchiorri et al. 2006), in combination with numerical modeling, have been used to determine the relative importance of the various processes controlling the seasonal H 2 O cycle (Flasar and Goody 1976;Jakosky and Farmer 1982;Haberle and Jakosky 1990;Richardson and Wilson 2002;Böttger et al. 2005). The seasonal H 2 O cycle is important for the Martian climate because it leads to the formation and affects the stability of water ice deposits on the surface and of water ice clouds in the atmosphere (Jakosky and Haberle 1992;Richardson and Wilson 2002;Smith et al. 1997;Haberle et al. 1999;Montmessin et al. 2004;Madeleine et al. 2012).
The goal of this article is to analyze the complete set of in-situ meteorological measurements from the 1970s Viking landers to the current Curiosity rover in order to review our understanding of the modern climate of Mars, with focus on the dust, CO 2 and H 2 O cycles and their impact on the near-surface radiative and thermodynamic conditions. In particular, we examine ground-based measurements made by the Viking Lander 1 (VL1), Viking Lander 2 (VL2), Mars Pathfinder (MPF) lander and Sojourner rover, Mars Exploration Rover A (MER-A, Spirit) and B (MER-B, Opportunity), Phoenix lander (PHX) and Mars Science Laboratory (MSL, Curiosity) rover ( Fig. 1). We provide values of atmospheric opacity, atmospheric pressure, near-surface air temperature, ground temperature, near-surface wind speed and direction, and near-surface air relative humidity and water vapor content with the highest confidence possible and study the diurnal, seasonal and interannual variability of these quantities at the various landing sites over a span of more than twenty Martian years. Although satellite observations and modeling have contributed significantly to our understanding of the Martian climate, in this article we focus only on ground-based measurements.
Section 2 describes the characteristics and performance of the meteorological instruments onboard each spacecraft that landed successfully on Mars. Section 3 shows the results of the data analysis, including diurnal, seasonal and interannual variations of atmospheric opacity (Sect. 3.1), atmospheric pressure (Sect. 3.2), near-surface air temperature (Sect. 3.3), ground temperature (Sect. 3.4), near-surface wind speed and direction (Sect. 3.5), and nearsurface relative humidity and water vapor content (Sect. 3.6). Section 4 describes the meteorological stations being developed for future missions to Mars (Sect. 4.1), proposes further measurements to improve our understanding of the Martian dust and H 2 O cycles (Sect. 4.2), and discuss the potential for liquid water formation under Mars' present day conditions and its implications for future Mars missions (Sect. 4.3). Finally, Sect. 5 summarizes the main results.

Meteorological Stations on Mars
Over the last 40 years, in-situ measurements made by meteorological instruments onboard various landers and rovers have been used to study the climate of Mars. Here, we review these instruments, including their goals, locations, sensors, performances, operations and main findings. Table 1 summarizes these instruments.

Viking Lander Missions
After several failures, the first successful meteorological observations from the Martian surface were made in 1976 by the two Viking landers (Chamberlain et al. 1976;Hess et al. 1977;Tillman et al. 1994). The main goals of these meteorological observations were to study The Viking Meteorology Instrument System (VMIS) onboard VL1 and VL2 consisted of sensors to measure atmospheric temperature, wind velocity and direction at 1.61 m above the ground, and pressure at 0.22 m above ground. The Viking landers measured temperature and wind in modules of 9 min at 0.25 or 0.125 Hz sampling rate, 20 min at 0.5 Hz sampling rate, or 39 min at 0.25 Hz sampling rate, spaced 1.5 h apart throughout the sol. In addition, operational modes consisting of extended periods of sampling at 1 Hz were used. Initially, four pressure measurements at 1 Hz were made at the start of every module, but this was later reduced to one sample every 17 min (Hess et al. 1977). The data acquisition frequency was reduced towards the end of the landers' lives. Furthermore, dust opacities at the two landing sites were determined from Sun images obtained by the cameras at a wavelength of 670 nm over a span of more than one Martian year (Colburn et al. 1989). These dust opacity measurements cover 328 out of the first 920 sols of VL1 and 250 out of 872 sols of VL2. The range, accuracy and resolution of the various VMIS sensors are summarized in Table 1.
Problems with VMIS onboard VL2 started on sol 25 when its temperature sensor had its first anomaly. From this sol onward, a VL2 engineering sensor was used for making temperature measurements, with a worst-case reduced accuracy of ±4°C (Hess et al. 1977;Sutton et al. 1978). Problems with VMIS onboard VL1 started on sol 45 when the heater of the wind direction sensor failed. This could be partially corrected until sol 376, when one of the two hot film wind sensors failed (Murphy et al. 1990).

Mars Pathfinder Mission
After a hiatus in in-situ surface measurements on Mars throughout most of the 1980s and the first half of the 1990s, the MPF lander re-started in-situ meteorological measurements in July of 1997 (Seiff et al. 1997;). As part of NASA's "faster, better, cheaper" strategy, MPF was a low-cost mission intended to operate for 30 sols. The goal of the MPF meteorological measurements was to study the variability of the Martian atmosphere at Ares Valley during the midsummer season on diurnal and seasonal time scales and compare the data with that from the VL1 obtained at the same time of the year. MPF landed at 19.3°N and operated for 84 sols.
The Atmosphere Structure Instrument/Meteorology package (ASI/MET) experiment had similar meteorological sensors to those of the Viking landers: a pressure sensor in the lander body at 0.1 m height above the surface, three atmospheric temperature sensors at 0.25, 0.5 and 1 m above the lander deck, which is at 0.27 m above the ground, a wind sensor at 1.1 m above the lander deck and three wind socks placed between the temperature sensors (see Table 1 for details on the instruments). ASI/MET measurements were made at a rate of 51 sessions per day, each 3 min long at a sampling rate of 0.25 Hz. Additionally, MPF conducted various 15 min and 1 h sessions at 1 Hz, and one complete 24 h session at 0.25 Hz. Dust opacities were obtained from direct images of the Sun acquired by the Imager for Mars Pathfinder (IMP). Smith and Lemmon (1999) provided several atmospheric opacity values per sol at four wavelengths (450, 670, 883 and 989 nm) on 70 out of the 83 sols of the mission.
The ASI/MET hot wire wind sensor experienced problems because it relied on the difference in response during adjacent periods with and without heating of the sensor. Unfortunately, the imposed temperature difference was comparable to (only a few times larger than) the natural temperature fluctuations of the daytime near surface air. In addition, nighttime temperatures were lower than those for which the sensor had been calibrated. Thus, the MPF wind data is not available in the NASA PDS. However, wind direction could be retrieved with an uncertainty of ∼15° ). In addition, eight sequences of measurements were made using the windsocks and several of these measurements were used to deduce the wind profile (Sullivan et al. 2000). The MPF pressure sensor also had calibration problems (Haberle et al. 1999). During pre-flight testing, the pressure sensor was tested at temperatures that were significantly lower than its design limits (220 K), and much lower than it experienced within the spacecraft on the Martian surface (265-300 K). It appears that this extreme stressing changed the sensor offset and caused an increase in the variation of the offset with temperature. It was not possible to verify the sensor gain in flight, but laboratory tests with the spare flight sensor suggested that the sensor stressing would not produce significant changes in gain.

Mars Exploration Rover missions
Similar to the Viking mission, the MER mission delivered two systems to the surface of Mars in January of 2004. The meteorology science goal of the MER mission was to provide high vertical resolution temperature profiles in the boundary layer, to measure dust and ice aerosol optical depth, and to determine aerosol properties from imaging (Squyres et al. 2003). MER-A (Spirit) landed at 14.6°S and operated for 2208 sols. MER-B (Opportunity) landed at 1.9°S and has been operating for more than 4638 sols as on Feb. 9, 2017.
The MER rovers do not have a dedicated set of environmental instruments (i.e., a meteorology station) like those of the previous missions. Instead, thermal infrared spectra of the Martian atmosphere taken by the Miniature Thermal Emission Spectrometer (Mini-TES) were used to determine surface temperature, near-surface air temperature at about 1.1 m, atmospheric temperatures at altitudes ranging from about 30 to 2000 m, and column-integrated values of dust optical depth and water vapor abundance . Mini-TES infrared spectra in 167 channels from 340 to 1997 cm −1 with a spectral resolution of 10 cm −1 were acquired at rates of 0.5 Hz. The Mini-TES spatial resolution is 20 mrad. Mini-TES typically conducted between one and four atmospheric observation sequences per sol. Nighttime observations were rare , with the vast majority of the measurements being made between 09:00 and 18:00 Local Mean Solar Time (LMST). Observations with Mini-TES pointing upward were made either at three fixed elevations scan in sets of 100 spectra each, or as a "stare" at a single elevation angle (typically 30°) in sets of 100-1275 spectra. Observations with Mini-TES pointing downward in sets of 10 spectra were used to measure the surface temperature and the near-surface air temperature. Even though Opportunity is still active, Mini-TES stopped operating when a globalscale dust storm obscured its mirrors in Mars Year (MY) 28. Ground and near-surface air temperature retrievals were made for the MER-A first 1124 sols and the MER-B first 1110 sols.
Atmospheric opacities at the MER landing sites are also obtained from direct images of the Sun at the wavelengths of 440 and 880 nm with the Pancam cameras . MER measurements have produced the most complete atmospheric opacity dataset available from the surface of Mars, with the MER-A dataset spanning more than 2200 sols, and the MER-B dataset so far spanning more than 4600 sols as of February 2017. The datasets presented here include results of measurements during 2030 sols for MER-A and 4010 sols for MER-B. Here, we use atmospheric opacity values obtained by Pancam instead of from Mini-TES measurements because the Pancam instrument provides a more complete coverage and are available in the NASA PDS.

Phoenix Mission
The Phoenix mission landed at 68.2°N in May of 2008 to study the history of water in the  Martian Arctic and to search for environments suitable for microbial life. It operated for 151  sols, exceeding the mission design life of 90 sols. PHX carried a meteorological station (MET) containing sensors to measure air temperature at heights of 1.25, 1.5 and 2 m above the surface, sensors to measure air pressure, and a "telltale" to measure wind at 2 m above the surface . The characteristics of the MET sensors are listed in Table 1. The MET instrument measured temperature and pressure at a rate of 0.5 Hz throughout the duration of the mission. The data could be processed on board to produce mean, standard deviation, minimum and maximum values within a 512 s window to reduce communication requirements. As part of the MET instrument, Phoenix also included a light detection and ranging (LIDAR) system to probe the vertical structure of the atmospheric boundary layer using pulsed laser light with wavelength of 532 nm backscattered from dust and ice particles from ∼200 m to ∼20 km above the ground (Whiteway et al. 2008). The optical extinction coefficient derived from the LIDAR signal can then be related to the ice-water content (IWC) profile to resolve the internal structure of water-ice clouds (Whiteway et al. 2009;Daerden et al. 2010;Dickinson et al. 2010;Moores et al. 2011). In addition, Phoenix's Thermal and Electrical Conductivity Probe (TECP) was used to measure subsurface and near-surface temperature, conductivity, permittivity and humidity (Hecht et al. 2008;Zent et al. 2009). The TECP consisted of a single electronics box, fitted with four needles, which could be inserted into the Martian soil, and was mounted on the Phoenix Robotic Arm. Atmospheric opacity measurements were obtained from direct images of the Sun by the Phoenix's Surface Stereo Imager (SSI) at the wavelengths of 451, 671, 887 and 991 nm. Opacities are available for 117 out of the 151 sols of the mission .
Wind data could only be obtained when the SSI was used to observe the telltale, a simple Kapton cylinder hanging from a Kevlar string in front of a mirror . A source of uncertainty in MET temperature measurements is lander-generated heat; this uncertainty is expected to be most severe on the lowest thermocouple for wind from 340°to 90°and from 240°through 310°, the directions of the solar panel. Another source of uncertainty in the temperature measurements is direct sunlight on the thermocouples, which can cause errors of up to 1 K .
Only the raw output of the TECP Relative Humidity (RH) sensor (unprocessed RH values) is available at the NASA PDS because of uncertainties in pre-flight instrument calibration. Values of temperature and RH covered in the pre-flight calibration (Zent et al. 2009) overlap only partially with the environmental conditions found at the Phoenix landing site ). The calibration function was revised in 2016 to correct for inaccuracies at the lowest temperatures observed at the Phoenix landing site (Zent et al. 2016), but the new processed RH values are not available in the PDS yet.

Mars Science Laboratory Mission
The MSL Curiosity rover landed in Gale crater (4.6°S) in August of 2012 and as of February 2017, it has been operating for almost 1600 sols. The main meteorological goals of the MSL mission are to study the climate, water and CO 2 cycles, and atmospheric processes ranging from local to synoptic scales.
Curiosity's Rover Environmental Monitoring Station (REMS) consists of wind velocity, air temperature, ground temperature, and air humidity sensors at 1.6 m height on the rover mast and pressure and UV sensors at 1 m on the rover deck (Gómez-Elvira et al. 2012, 2014. The REMS sampling strategy consists of 5 min-long hourly samples at 1 Hz throughout the mission, with interspersed full hour sample periods at 1 Hz to cover every time of the sol over a period of a few sols. The cadence of full hour measurements is restricted by memory size and downlink data times and volume. The range, accuracy and resolution of the REMS sensors are listed in Table 1. Atmospheric dust opacities are obtained from direct images of the Sun at the wavelengths of 440 and 880 nm by the Mastcam instrument . In this manuscript, we analyze REMS and Mastcam measurements corresponding to the first 1526 sols of the MSL mission. Because the REMS temperature sensor is located only about 0.6 m above the rover deck, the heating of the rover by solar radiation and by the Radioisotope Thermoelectric Generator (RTG) may affect the air temperature measurements. In addition, the rover affects ground temperature measurements because the field of view of the Ground Temperature Sensor (GTS) is within the area of the ground heated by thermal radiation from the RTG. Indeed, the RTG can contribute up to 4 K to the measured ground temperature . The RH sensor is affected by heating produced by the sensor itself, which increases its temperature by up to ∼1.5 K minutes after it is turned on. Thus, the local RH measured by the sensor is closest to the actual RH of the atmosphere only during the first few seconds of operation, and values of the water vapor volume mixing ratio (VMR) derived from these measurements represent the VMR of the atmosphere. In long measurement sessions, the temperature of the RH sensor stabilizes after ∼30 min and then starts to follow changes in the atmospheric temperature, making derived VMR values again representative of atmospheric values. The data obtained during the transient heating period (∼30 min) may be retrieved by compensating for the heating effect; this requires additional analysis, which is planned for the future. Furthermore, the RH sensor response time is a logarithmic function of temperature, ranging from 0.1 s at 293 K to 700 s at 203 K, but the sensor adjusts to a changing RH while turned off and therefore provides correct values immediately upon being turned on (Gómez-Elvira et al. 2012). Recent laboratory measurements with a spare flight unit of the RH sensor in a low pressure CO 2 environment indicate that low raw values measured on Mars are actually the result of CO 2 effects on the Humicaps, supporting the use of the compensation method described by . According to the preliminary results of the laboratory measurements, the dynamic range of the sensor does not change in a low pressure CO 2 atmosphere .
MSL has two wind sensor booms mounted on the remote sensing mast (RSM), but after landing the side-rear-facing sensor was damaged and only two sensor boards on the front-facing boom have been working properly. Since the rover body (primarily the RSM) strongly disturbs wind reaching the remaining boom from the rear of the rover, the wind sensor can only be used to retrieve wind speeds and accurate wind directions when wind comes from the hemisphere in front of the rover. 'Rear' wind directions are retrieved with a much lower accuracy and wind speeds measurements are not reliable. Wind measurements are not available during the coldest periods of each sol (lasting between four and ten hours in local summer and winter, respectively) because of high electronic noise due to the performance of the Application-Specific Integrated Circuit (ASIC)-based Sensor Front-End, which transmits the wind data to the rover (Newman et al. 2017). However, 'wind characterization' campaigns-in which the rover is pointed in multiple directions over several sols-have been performed at various times and locations to provide the most complete wind information possible (Newman et al. 2017). Here we show results from one campaign only, although the full REMS wind speed and direction dataset is available in the PDS.

Results
In this section we analyze meteorological measurements from the Viking landers to the Curiosity rover to provide values of atmospheric opacity (Sect. 3.1), atmospheric pressure (Sect. 3.2), near-surface air temperature (Sect. 3.3), ground temperature (Sect. 3.4), nearsurface wind speed and direction (Sect. 3.5), and near-surface relative humidity and water content (Sect. 3.6) over a span of more than twenty Martian years.

Atmospheric Opacity and Surface Radiative Environment
Here we analyze atmospheric opacity measurements taken at 670 nm by the VL cameras, 883 nm by the MPF/IMP, 880 nm by the MER/Pancam cameras, 887 nm by the PHX/SSI and 880 nm by the MSL/Mastcam. In combination with values of the atmospheric opacity, we provide estimated values of the surface solar insolation (defined as the solar irradiance integrated over one sol) at each landing site to show the impact of atmospheric opacity on the surface radiative environment. Figure 2 (left column) shows the atmospheric opacity values at the VL1, VL2, MPF, MER-A, MER-B, PHX and MSL landing sites as a function of solar longitude. Among these missions, only the MERs and MSL can be used to study the interannual variability of the atmospheric opacity (Figs. 2c, 2e and 2g), with VL1 and VL2 measurements barely overlapping in L s during MY 12 and MY 13 (Fig. 2a). During the aphelion season (L s = 0°-180°), the opacity values are low (generally below 1, which indicates that the direct radiation at the surface has been reduced to exp(−1/μ 0 ) times its value at the top of the atmosphere, where μ 0 is the cosine of the solar zenith angle assuming a plane-parallel atmosphere), and the interannual variability is small (mean standard deviation around 0.05 and 0.03 at MER-A and MER-B locations). On the contrary, the opacity values and their interannual variability are larger during the perihelion season (L s = 180°-360°) at each landing site (Fig. 2, left column). Differences in solar insolation between the perihelion and aphelion seasons produce on average higher temperatures during the perihelion season, enhancing the injection of dust into the atmosphere via wind stress and dust devils (Haberle et al. 1982;Newman et al. 2002aNewman et al. , 2002bKahre et al. 2006;Greeley et al. 2010). Satellite measurements of atmospheric opacity show similar aphelion-versus-perihelion behavior on a planet-wide scale (Smith 2004;Cantor et al. 2001;Montabone et al. 2015).
Global dust storms, which occur every few years, are of particular interest because of their large impact on the atmospheric circulation and surface radiation budget. So far, these global dust storms have been observed from the surface only by the VL1 and VL2 in MY 12, and by the MER-A and MER 2c and 2d). In both years, the peak atmospheric opacity was extraordinarily high (τ > 3) and occurred shortly after the perihelion, at L s ∼ 300°in MY 12 and at L s ∼ 280°in MY 28. Interestingly, in years when global dust storms did not occur, a local minimum in the atmospheric opacity occurred around this time (Figs. 2c, 2e and 2g). A global dust storm has not occurred since the Curiosity rover started operations at Gale crater (Fig. 2g). The MPF and PHX missions did not observe any major dust storms because they were short and ended prior to the beginning of the dusty season (L s = 188°and 148°, respectively).
Despite the large interannual variability observed at each landing site during the perihelion season, there are common aspects in the evolution of the atmospheric opacity during this season. At each landing site, two periods of enhanced dust activity were observed each year, the first period between L s = 210°and 240°, and the second between L s = 320°and 340° (Fig. 2, left column). A third period of enhanced dust activity is apparent at the MER-A landing site between L s = 150-170° (Fig. 2c). The values of atmospheric opacity at the MER-B and MSL sites during this period were relatively high but did not present clear peaks (Figs. 2e and 2g). Although dust opacity measurements were only taken during the first 920 sols of the VL1 mission (corresponding to MY 12 and 13), analysis of air temperature data suggests that major dust storms also occurred between 270°and 300°in MY 14, and at around 210°in MY 15 (Savijärvi 2014). Figure 2 (right column) shows the seasonal and interannual variability of the solar insolation at the TOA (solid lines) and at the surface (crosses) at each landing site. A comprehensive spectral radiative transfer model with the latest wavelength-dependent radiative properties of dust, water ice clouds and atmospheric gases is used to estimate the solar insolation at the surface, with values of atmospheric opacity measured at each landing site as inputs (Vicente-Retortillo et al. 2015). Although there are situations in which water ice clouds can significantly contribute to the atmospheric opacity, such as during the aphelion cloud belt or at the edges of the polar caps (Toigo and Richardson 2000;Smith 2004;Madeleine et al. 2012;, their effect on the solar insolation at the surface is generally minor compared to that of the dust. To quantify this effect, we have used the radiative transfer model described in Vicente-Retortillo et al. (2015) to simulate the daily solar insolation at L s = 90°(when opacity due to water ice clouds shows relatively high values), assuming that half of the opacity is caused by water ice clouds. The difference in the solar insolation at the various landing sites ranges from 3% to 5% with respect to values obtained when only dust is assumed to contribute to the atmospheric opacity.
Surface solar insolation departs from that at the TOA because the solar radiation is absorbed and scattered by atmospheric constituents such as dust particles. This is particularly noticeable during the dusty perihelion season, when abrupt changes in atmospheric dust content (Fig. 2, left column) cause strong intraseasonal and interannual variations in surface solar insolation (Fig. 2, right column). In contrast, during the aphelion season when dust activity declines, intraseasonal and interannual variations in solar insolation at each landing site are small, and the variation in surface solar insolation follows that at the TOA smoothly. However, water ice clouds can also partially contribute to variations in solar insolation during this season, which we are not considering here.
At the VL sites in the northern hemisphere, the decreasing trend of solar insolation at the TOA during the perihelion season (northern fall and winter; see solid lines in Fig. 2b) is further reduced at the surface by enhanced dust activity (Fig. 2a), while the increasing trend at the TOA during the aphelion season (northern spring and summer) coincides with the lowest annual values in atmospheric opacity, which further enhances the insolation at the surface. Consequently, the highest annual values of surface solar insolation at the VL1 and VL2 sites occur during the aphelion season at around L s ∼ 150°and 120°, respectively. This is around the time that the insolation at the TOA is the highest and atmospheric opacity values are the lowest. In contrast, at the MER and MSL sites in the southern hemisphere, the increasing trend in solar insolation at the TOA during the perihelion season (southern spring and summer; see solid lines in Figs. 2d, 2f and 2h) is counteracted at the surface by the increase in dust opacity (Figs. 2c, 2e and 2g). Therefore, the highest annual values of solar insolation at the MER and MSL sites occur each year at L s ∼ 200°and L s ∼ 300°( with the exception of MY 28 at both MER landing sites), coinciding with local minima in atmospheric opacity and relatively high values (but not annual maximum) of solar insolation at the TOA.

Atmospheric Pressure
Surface pressure variations depend on many processes. On seasonal time-scales, variations in the bulk atmospheric mass due to the condensation and sublimation of the sea- Interannual and seasonal evolution of the daily mean atmospheric pressure (top) and the diurnal amplitude (difference between the daily maximum and minimum atmospheric pressure; bottom) at the VL1 (green colors), VL2 (gray colors), MPF (black), PHX (purple) and MSL (brown-red-orange) landing sites sonal polar caps cause large pressure variations across the planet (Hess et al. 1980;Tillman et al. 1993). On daily to monthly time scales, various components of the global circulation such as baroclinic eddies (Barnes 1980(Barnes , 1981, tropical overturning circulations (Magalhães 1987;Haberle et al. 1993a;Forget et al. 1999;Read and Lewis 2004;Moudden and McConnell 2005), and stationary waves (Hollingsworth and Barnes 1996;Banfield et al. 2003;Hinson et al. 2003) are associated with changes in surface pressure. On diurnal time scales, the global scale thermal tides (Zurek 1976;Zurek and Leovy 1981;Wilson and Hamilton 1996;Guzewich et al. 2013) and localized slope flows (Blumsack et al. 1973;Tyler and Barnes 2013;Rafkin et al. 2016) are important. On time-scales ranging from a few seconds to minutes, convective dust devils also induce measurable changes in surface pressure Steakley and Murphy 2016). Next we show evidence of pressure variations caused by these systems. Figure 3 (top) shows the daily mean atmospheric surface pressure measured during the VL1, VL2, MPF, PHX, and MSL missions as a function of solar longitude. To calculate representative daily mean surface pressure values at each site, we only use measurements taken on sols that have full diurnal coverage, resulting in 2245 sols at the VL1 site, 1050 sols at the VL2 site, 15 sols at the MPF site, 143 sols at the PHX site and 1096 sols at the MSL site. Of these five data sets, only the two VLs and the MSL have sufficient coverage to allow assessment of the annual repeatability of the pressure cycle. From MY 12 to MY 15, the interannual variability at the VL sites during the aphelion season was generally small (of the order of a few Pa), while it increased significantly during the perihelion season (up to of tens of Pa). This is particularly noticeable at the VL2 site during the northern winter (L s ∼ 280°) in MY 12 (Fig. 3, top). The pressure increase was caused by a global dust storm (Fig. 2a), which intensified the tropical overturning circulation whose descending branch in the VL2 region caused such a pressure increase (Haberle et al. 1982;Wilson 1997). At the MSL landing site, a noticeable (∼15 Pa) decrease in the mean surface pressure is observed during the mission (Fig. 3, top). This pressure decrease is mostly accounted for by the increase in altitude as the Curiosity rover has been moving uphill, climbing from about −4520 m at the beginning of the operations to about −4336 m on sol 1526, corresponding to the last MSL sol analyzed in this article. The good agreement between the predicted hydrostatic pressure decrease and the measurements (de la  suggests that the drift in the pressure is largely explained by changes in the rover elevation. On seasonal time-scales, the variability of the daily mean surface pressure is large at each landing site, with two distinct minima and maxima per Martian year caused by the condensation and sublimation of the CO 2 polar caps (Fig. 3, top). During the spring season in the southern hemisphere, parts of the southern polar cap sublimate and the released CO 2 increases the atmospheric pressure with an annual maximum occurring at each landing site in late southern spring at L s ∼ 260°. Then, during the cold season in the southern hemisphere, the deposition of CO 2 into the southern polar cap causes an annual minimum pressure in late southern winter at L s ∼ 150°. The relative minimum in late southern summer (L s ∼ 345°) and relative maximum in late southern fall (L s ∼ 55°) at each landing site are caused by the deposition and sublimation of CO 2 into and from the northern polar cap.
The sol-to-sol variability at each landing site is relatively small during the aphelion season, but increases significantly during the perihelion season (northern autumn and winter) particularly at the VL2 site (Fig. 3, top). The larger sol-to-sol variability during the perihelion season is caused by baroclinic instability in autumn and winter, which generates perturbations in surface pressure that peak at the latitude of the polar jet stream (Barnes 1981). Since the VL2 site is closest to the high-latitude jet, the largest sol-to-sol variations occur at the VL2 site, followed by those at the VL1 (northern edge of the tropics) site, with the lowest variability in this sequence at the MSL site in the southern tropics. The six-toeight-sol oscillations observed by REMS are consistent with the periods observed at VL1 and VL2 (de la Torre-Juárez et al. 2014). Indeed, the consistent seasonal trend in sol-to-sol variability at VL2, VL1 and MSL is due to strong inter-hemispheric asymmetry in transient eddy activity forced by the inter-hemispheric asymmetry in topography (Barnes et al. 1993), with baroclinic activity strongly suppressed in southern high latitudes during southern autumn and winter (near aphelion). As the signature of northern storm activity traveling into the southern hemisphere has been seen in models and orbital imaging data (Wang et al. 2003), it is highly plausible that REMS is seeing their signature in Gale crater .
Differences in the absolute values of the surface pressure measured by the various missions cannot be fully accounted for by differences in the elevation of the landing sites (Table 1). For instance, the surface pressure at the VL2 site (−4.5 km) during late autumn and early northern winter is significantly higher than at the MSL site (−4.4 km) (Fig. 3, top). As explained by , this is due to the fact that temperatures are lower at the VL2 than at the MSL site for this season, and that elevation differences between the two hemispheres allow air to pool in the northern lowlands at the expense of air in the southern highlands. Thus, the pressure increases in the northern hemisphere and decreases in the southern hemisphere because of this orographic effect (see also Hourdin et al. 1993).
Figure 3 (bottom) shows the amplitude of the diurnal surface pressure variation as a function of solar longitude at each landing site. This amplitude (also called the diurnal atmospheric tide) is highly correlated to atmospheric dust opacity (Zurek and Leovy 1981;Bridger and Murphy 1998;Wilson and Hamilton 1996;, with an Fig. 4 Diurnal evolution of the atmospheric pressure at the VLs, MPF, PHX and MSL landing sites at L s ∼ 145°, when pressure measurements were taken by each spacecraft (see Fig. 3). The hours-long periods of constant pressure in the Viking data are a consequence of the low resolution of the digitization of the data (∼9 Pa) additional contribution from background seasonal cycles of the individual tide modes (e.g., the westward-propagating migrating diurnal tide and the eastward-propagating diurnal Kelvin wave). Therefore, the amplitude of the diurnal surface pressure variation is larger in the perihelion season (when dust activity is greater) than during the aphelion season (Fig. 3,  bottom). Indeed, the amplitude of the diurnal surface pressure variation is largest at the VL1 and VL2 landing sites in MY 12 around L s ∼ 280° (Fig. 3, bottom) when the atmospheric opacity is high because of a global dust storm (Fig. 2a).
The amplitude of the diurnal tide at the MSL site is larger than at the VL2 site throughout the year (Fig. 3, bottom), in spite of the similar elevation of the two sites. This occurs even during the aphelion season when dust opacity values are similar in both sites (Figs. 2a and 2g). This partly results from the fact that the thermal forcing of the diurnal tide is strongest at the equator, which also explains the small amplitude of the tides at the PHX site ). Its amplitude can be further enhanced by the interaction of eastward and westward propagating tidal modes (e.g., Wilson and Hamilton 1996). At the longitude of Gale crater this interaction is constructive and boosts the diurnal tidal amplitudes. However, the tides observed by Curiosity are much larger than those simulated by global models, even when accounting for these tidal interactions. As such, various researchers have suggested further enhancements at the Curiosity site by crater circulations Tyler and Barnes 2015), topographic variations in the zonal band (Rafkin et al. 2016), and a hydrostatic adjustment flow . Figure 4 shows the diurnal evolution of surface pressure at each landing site at L s = 145°, a seasonal date at which there is an overlap in data collected by all missions included in this study (Fig. 3, top). Except for the PHX mission, the daily maximum pressure value at each landing site occurs at around 08:00 LMST, with the daily minimum occurring between 16:00 and 18:00 LMST. The timing of these daily maxima and minima is strongly controlled by the superposition of the thermal tides. At the MSL site, the superposition of the diurnal and semidiurnal tidal components, whose phases are roughly consistent with those predicted by classical tidal theory (see , largely accounts for the timing of the observed maxima and minima. The fact that the other sites have the same timing suggests a similar superposition. However, interactions with the eastward propagating diurnal Kelvin mode, which is closest to resonance at this season (Zurek 1988), can complicate this simple picture.

Fig. 5
Interannual and seasonal evolution of the daily mean near-surface air temperature (top) and the diurnal amplitude (difference between the daily maximum and minimum near-surface air temperature; bottom) at the VL1 (green colors), VL2 (gray colors), MPF (black), PHX (purple) and MSL (brown-red-orange) landing sites. MPF and PHX values correspond to measurements taken at ∼1.1 m and ∼2 m above the surface (furthest distance from the body of the landers), while VL1, VL2 and MSL values correspond to measurements taken at ∼1.6 m

Near-Surface Air Temperature
In this section we analyze the seasonal and interannual variability of the near-surface air temperature at the VL, MPF, MER, PHX and MSL landing sites. Figure 5 (top) shows the daily mean air temperature at these sites as a function of solar longitude. Only measurements on sols with full diurnal coverage are used to calculate these daily mean air temperature values, resulting in 1201 sols at the VL1 site, 830 sols at the VL2 site, 15 sols at the MPF site, 143 sols at the PHX and 1323 sols at the MSL site. At the MPF and PHX sites, the daily mean air temperature is calculated at 1.27 and 2 m height above the surface, respectively, to minimize the effects of thermal contamination by the lander.
The interannual temperature variability at the VL1, VL2 and MSL sites is generally small during the aphelion season (of the order of a few Kelvin), while during the perihelion season it noticeably increases at the VLs sites (Fig. 5, top). This is particularly pronounced at the VL2 site between L s = 260°and 300°in MY 12, when daily mean air temperatures were consistently larger than in MY 13 coinciding with a global dust storm (Fig. 5 top, gray colors). Since much less change in daily mean temperature is seen at VL1 during the storm (see also Ryan and Henry 1979), the increase at VL2 may be the result of an increase the downward infrared fluxes at the VL2 site during the storm. Large (20-30 K) increases in atmospheric temperatures (Martin and Kieffer 1979) associated with strong adiabatic warming in the downwelling branch of the overturning circulation (e.g., Haberle et al. 1982;Wilson 1997) would increase infrared emission to the surface possibly strongly enough to compensate for any decrease in solar insolation (which is small at this latitude and season).
Consistent with the increase in the annual amplitude of surface solar insolation with latitude (Figs. 2b and 2h), the seasonal variability of the daily mean air temperature is larger at the VL2 site than at the VL1 and MSL sites (Fig. 5, top), with the annual amplitude decreasing from ∼60 K at the VL2 site to ∼30 K at the VL1 and ∼20 K at the near-equatorial MSL site. The minimum annual temperature at the VL sites occurs during northern winter, between L s = 270°and 300° (Fig. 5, top), when the surface solar insolation presents an annual minimum (Fig. 2b), while the maximum annual temperature occurs in northern summer between L s = 120°and 150°at the VL1 site, and between L s = 90°and 120°at the VL2 site, when the surface solar insolation presents an annual maximum at both landing sites (Fig. 2b). Similarly, the maximum annual temperature at the MSL site coincides with the annual maximum in surface solar insolation around L s ∼ 210°(southern spring), while the minimum annual value coincides with the annual minimum in surface solar insolation around L s ∼ 90° (Figs. 5, top and 2h).
The sol-to-sol variability at each landing site is relatively small during the aphelion season, increasing moderately during the perihelion season particularly at the VL2 site (Fig. 5,  top). As with the surface pressure, larger sol-to-sol variations during the perihelion season are caused by baroclinic eddies in the northern hemisphere during autumn and winter (Barnes 1980(Barnes , 1981. The daily mean air temperatures at the PHX site in Mars' Arctic were similar to those at the VL2 site between L s = 78°and 125° (Fig. 5, top), and larger than that at the VL1 site. The surface solar insolation at the PHX site is also similar to that at the VL2 and higher than at the VL1 site between L s = 78°and 125° (Figs. 2a and 2h). Similarly, daily mean air temperatures at the MPF site during most of the mission were slightly higher than those at the MSL site (Fig. 5, top), with the surface solar insolation at MPF site slightly higher than at the MSL site (Fig. 2h).
The diurnal amplitude of near-surface air temperature at the VL1, VL2, MPF, PHX and MSL sites is shown in Fig. 5 (bottom) as a function of solar longitude. The abrupt decrease in amplitude from ∼42 to ∼12 K at the VL1 site around L s ∼ 275°was caused by a dramatic increase in atmospheric opacity during the global dust storm occurred in MY 12 (Fig. 2a). At the VL2 site, the contemporaneous decrease in the diurnal amplitude of the air temperature was noticeable but much less pronounced than at the VL1 site because the surface solar insolation at the VL2 site is low at this time of the year (northern winter). Therefore, the increase in the atmospheric opacity does not decrease the diurnal amplitude of the air temperature significantly at this location. At the MSL site, the interannual and seasonal variability of the diurnal amplitude of near-surface air temperature is smaller than at the VL sites. This occurs because seasonal variability in solar insolation is small in the tropics and because there were no major dust storms during the observational period. The diurnal amplitude of the air temperature is similar at the MPF and PHX sites during the aphelion season when there is an overlap in date from the two missions (Fig. 5, bottom).
Observations with Mini-TES pointing downward were used to measure the ground and near-surface air temperature at ∼1.1 m for the MER-A first 1124 sols and the MER-B first 1110 sols. Since nighttime observations were rare at both landing sites, representative values of the daily mean air temperature at ∼1.1 m cannot be obtained. Instead, we show in Fig. 6 the diurnal evolution of ∼1.1 m air temperatures at the MER-A (top) and MER-B (bottom) sites as a function of the local mean solar time with variation in L s shown using color code for MY 27 (the only year with complete annual coverage). For comparison with southernhemispheric MSL air temperatures, we also show the diurnal evolution of the air temperature at the MSL site at L s = 210°(sol 103) and L s = 83°(sol 528) as solid lines, coinciding with the warmest and coldest periods at the MSL site (Fig. 5, top). At the MER-A site (14.6°S), the annual amplitude of surface solar insolation is larger than at the MSL (4.6°S) due to its location further from the equator, resulting in warmer daytime temperatures during the annual warmest season and colder daytime temperatures during the annual coldest season than at the MSL site (compare orange and green lines with solid circles of the same color in top panel of Fig. 6). It is interesting to note that temperatures in the range 200-220 K between 12:00-14:00 LMST (see isolated green dots corresponding to L s ∼ 171°in top Fig. 6), were retrieved during a regional dust storm that raised the atmospheric opacity values from 0.6 to 1.4 in the span of a few sols (see Fig. 2c, MY 27; L s ∼ 170°), causing these lower daytime temperatures. At the MER-B site (1.9°S), the values and annual range of the surface solar insolation are similar to those at the MSL site (Figs. 2f and 2h). Indeed, at the MER-B site the air temperature at ∼1.1 m stays approximately within the limits of MSL's warmest and coldest sols (Fig. 6, bottom).

Ground Temperature
The values and variability of ground temperature are mainly governed by the surface radiation budget, which depends on the time of year, location, atmospheric opacity, and the thermal and physical properties of the terrain. In this section, we analyze the diurnal, seasonal and interannual variability of ground temperature at the MER, PHX and MSL sites.
Among the MER, PHX and MSL missions, only MSL conducted full diurnal measurements of ground temperature. Figure 7 shows the daily mean, maximum, and minimum ground temperature during the first 1526 sols (MY 31-33) of the MSL mission as a function of solar longitude. To calculate these temperatures, we use only the REMS Ground Fig. 7 Interannual and seasonal evolution of the daily maximum, mean and minimum ground temperature at the MSL landing site. Only the MSL mission provides temporal coverage good enough to show interannual and seasonal evolution of T g Temperature Sensor (GTS) data with the highest confidence possible, that is, with the ASIC power supply in its operation range, the highest recalibration quality and with no shadows in the GTS field of view (see  for details on GTS uncertainties) from measurements on sols with full diurnal coverage (resulting in 413 sols).
The interannual variability of the daily mean ground temperature at the MSL site is small, of the order of a few Kelvin (Fig. 7). This occurs in spite of the measurements being acquired in different locations as the rover traversed Gale crater. Small variations could be caused by changes in the albedo of the surface along the traverse, which range from 0.1 on sandy terrains to 0.3 on mudstone . The seasonal variability of the daily mean ground temperature is also small at the MSL site, with an annual amplitude of about 20 K (Fig. 7). The minimum annual temperatures occur during late fall at L s ∼ 85° (Fig. 7), when the surface solar insolation is the smallest (Fig. 2h), whereas the maximum annual values occur during spring at L s ∼ 210°, when the surface solar insolation is the largest.
The diurnal amplitude of ground temperature at the MSL site is primarily governed by the type of terrain, and in particular by its thermal inertia (I ). Terrains with low thermal inertia experience higher daytime maximum temperatures and lower nighttime minimum temperatures than terrains with higher thermal inertia. As an example, the sudden decrease in the diurnal maximum ground temperature (and increase in the diurnal minimum temperature) which occurred at L s ∼ 220°in MY 31 (Fig. 7) coincided with Curiosity's traverse from a terrain with I ∼ 300 J m −2 K −1 s −1/2 to a terrain (Yellowknife Bay) with I ∼ 420 J m −2 K −1 s −1/2 . During this traverse (around sol 120), dust opacity values barely changed and thus neither did the surface solar insolation. Indeed,  used ground temperature measurements and modeling to calculate the thermal inertia at the MSL site, obtaining values ranging from about 170 to 550 J m −2 K −1 s −1/2 in the first 1250 sols of the mission. The differences in terrain traversed by the Curiosity likely account for the variability of the diurnal amplitude of ground temperature.
Changes in atmospheric opacity also regulate the amplitude of variations in ground temperature (e.g. Määttänen and Savijärvi 2004). For a given type of terrain, lower atmospheric opacity causes higher daytime maximum temperatures and lower nighttime minimum temperatures and vice versa (same as with thermal inertia). This occurs because aerosols like dust scatter and absorb solar radiation, reducing the insolation at the surface. At night, downward infrared emission from aerosols warms the ground. However, the impact of variations in atmospheric opacity on the amplitude of the ground temperature has been very small at Gale crater compared with that of changes in thermal inertia. For instance, the sudden increase in atmospheric optical depth in MY 31, from ∼0.6 at L s = 200°to ∼1.2 at L s ∼ 208° (  Fig. 2g) was not followed by a decrease in the daily maximum ground temperature (Fig. 7). Instead, the daily maximum ground temperature kept increasing until L s ∼ 220°, when the rover reached a terrain with much higher thermal inertia .
As for the near-surface air temperature discussed in Sect. 3.3, we show in Fig. 8 the diurnal evolution of ground temperature at the MER-A (top), MER-B (middle) and PHX (bottom) sites as a function of the local mean solar time with variation in L s during MY 27 shown as colors at the MER sites, and in MY 29 at the PHX site. For the sake of comparison, we also show in Fig. 8 the diurnal evolution of the ground temperature at the MSL site at L s = 210°(sol 103) and L s = 86°(sol 534) in solid lines, coinciding with the warmest and coldest periods at the MSL site (Fig. 7). MER-A's further distance from the equator and larger annual amplitude of surface solar insolation than those of the MSL site result in significantly lower daytime temperatures at L s ∼ 86°and slightly greater daytime temperatures at L s ∼ 210° (Fig. 8, top). At the MER-B site, where the surface solar insolation is comparable to that at the MSL, ground temperature values stay to a very good extent within the limits marked by MSL's warmest and coldest sols (Fig. 8, middle).
The PHX/TECP instrument was used to measure the soil temperature during 16 sols, with sporadic diurnal coverage on each of these sols. Since the coupling of the TECP needles with the regolith was insufficient (Zent et al. 2010), and the TECP needles are 15 mm in length, the measured temperature represents an integrated temperature over the top 15 mm of the subsurface rather than surface skin temperature of the regolith. Figure 8 (bottom) shows the complete set of TECP soil temperature measurements. The daytime maximum temperature of ∼259 K was achieved on sol 71, coinciding with the highest value of surface solar insolation during PHX operations at L s ∼ 107° (Fig. 2h). As the mission progressed and the surface solar insolation declined (Fig. 2h), temperatures quickly decreased, with daytime maximum values of only ∼230 K and nighttime minimum values as low as ∼180 K at L s ∼ 133°on sol 122 (the last sol with nighttime measurements). The local minimum between 13:00 and 14:00 LMST is caused by shadows of different parts of the PHX spacecraft on the surface.

Near-Surface Wind Speed and Direction
Of all quantities discussed here, wind speed and direction are most strongly influenced by the local and regional topography. In addition, they also depend on the surface radiation budget and large-scale variation in atmospheric opacity, via their impact on the global overturning circulation, atmospheric tides and baroclinic/barotropic waves, among other processes. The extent to which each dominates the cycles of wind depends on the instrument's locationand for topography, its precise location, since moving only a few meters can significantly impact the wind provided that the local topography changes dramatically (Newman et al. 2017). It is worth noting that the VLs, MPF and PHX all landed on smooth, gently sloping plains in the northern hemisphere, largely due to entry-descent-landing (EDL) safety constraints. In such regions, the influence of thermal tides and regional slopes would be expected to dominate the diurnal variation of wind, at least in periods without significant wave activity (such as local summer in the northern hemisphere). By contrast, the landing system used for MSL was more sophisticated and permitted landing inside a ∼155 km-diameter crater, hence the topography around MSL is more extreme than at any previous site and strongly dominates in most seasons. In this section we analyze the diurnal, seasonal, and interannual variability of wind speed and direction at the VL, PHX and MSL sites. Figure 9 shows the variation of the daily minimum and maximum wind speed as a function of season at the VL1 (up to sol 44, when the hot wire heater failed), VL2, and PHX sites. The daily mean wind speed is also shown for the VL sites. The VL2 dataset shows the most complete seasonal data set, and demonstrates the strong dependence of wind speed on global and regional rather than local forcing at this site. Wind speeds are lowest during northern hemisphere spring and summer (L s ∼ 10°to 190°), despite solar insolation at the VL2 site peaking during this period (Fig. 2b). The processes dominating wind speed variation at the relatively high northern latitude of VL2 (47.9°N) are twofold: (i) the global circulation and (ii) wave activity. The global circulation is much stronger in northern fall and winter due to the stronger heating around southern summer solstice (due to Mars's large eccentricity and the timing of perihelion), the topographic dichotomy (which tends to enhance a Hadley Minimum (dark blue), mean (green), and maximum (red) daily wind speed at the VL2 site at a function of L s in MY 12 and 13. Also shown is the minimum (purple), mean (cyan), and maximum (orange) daily wind speed at the VL1 site over the first 44 sols of its mission in MY 12 (see left hand side of the figure), and the minimum (light blue) and maximum (yellow) daily wind speed at the PHX site in MY 29. Note that the VL minimum and maximum values are calculated using the 1/25th of a sol binned dataset, thus do not correspond to the smallest and largest individual winds measurements. The PHX values shown do correspond to individual measurements; however, the range of local times covered varied over the period shown, thus these data contain inherent bias circulation with rising motion in the south), and the dust feedback (with the higher winds causing more dust lifting, which in turn causes more radiative heating, which in turn causes a stronger circulation). Baroclinic waves produced by the strong thermal contrasts around the northern seasonal polar cap's edge, which are strongest in northern fall through early spring, cause wind speeds to peak during northern fall and winter in both VL2 years. This also likely explains the gradual increase in wind speed throughout the even higher northern latitude PHX mission.
The second VL2 year presents a drop in wind speed in the middle of the fall/winter period, around perihelion (L s ∼ 251°) and solstice. This is known as the 'solsticial pause' (e.g., Lewis et al. 2016) and its cause is not well understood. By contrast, the same period in the first year of VL2 shows only a slight drop in wind speed before a rapid increase begins at about L s = 270°. In addition, the period from L s = 270°to L s = 320°during the first year of VL2 exhibits far less sol-to-sol variability than the rest of the fall/winter seasons. This is due to the presence of the large 1977b global dust storm over this exact period, as shown by the increase in VL1 and (to a lesser extent) VL2 opacities in Fig. 2a. While the impact of the earlier 1977a storm (at about L s ∼ 220°) is less clear, the dust-storm-induced enhancement and broadening of the global circulation during the 1977b storm clearly extends out to VL2 latitudes, producing a noticeable increase in daily mean, maximum and minimum wind speeds. Again, this demonstrates the global and regional rather than local scale response of winds to change in dust loading and hence radiative forcing. The local opacities at the VL2 site are slightly higher during the earlier storm, which would suggest it should have the greatest effect if local forcing dominates. Instead, it is the global scale changes in dust loading (and hence circulation) that dominate and produce a far greater impact due to the later storm, which had the higher dust opacities planet-wide (as shown by the peak in VL1 dust loading, for example). The reduction in sol-to-sol variability reflects the suppression of baroclinic wave activity at VL2 during the 1977b storm (Barnes 1980). Wind direction may vary through a full 360°over a sol, thus the mean value (or maximum and minimum) is usually difficult to interpret and potentially misleading. For this reason we only show wind direction as diurnal variation at selected times of year, rather than trying to show a mean value or range as a function of season. Figures 10 and 11 show respectively the wind speed and wind direction as a function of local time over several sols at one time of year at the VL1 site, at six times of year (including two examples of the same season in different years) at the VL2 site, at one time of year at the PHX site, and at the MSL site during the Bagnold Dunes wind characterization investigation. Figure 12 shows the same information, but now averaged over the sols shown in Figs. 10 and 11 for each site, and plotted in the form of hodographs. In general, wind speeds peak during the daytime at all landing sites and seasons displayed, likely because of mixing by convectively-forced turbulence. Exceptions are at the VL2 site at L s ∼ 195°(Figs. 10b and f), where peak wind speeds occur at ∼21:00 in some sols, and at L s ∼ 285°during the 1977b global storm (Fig. 10c), as discussed further below. The greatest VL2 wind speeds occur in the L s = 345°plot (Fig. 10d), late in the northern winter period discussed above, which contains some of the highest wind speeds measured on Mars. The lowest wind speeds occur in the VL2 L s = 105°plot (Fig. 10e), during northern summer.
The low interannual variability in both wind speed and direction during the equinoctial season (L s = 195°) at the VL2 site can be seen by comparing plots (b) and (f) in Figs. 10-12. Shortly after northern winter solstice at L s = 285°, however, there is significant year-to-year variability, as shown by comparing plots (c) and (g). The first year has stronger peak winds, a different time of peak wind occurrence (shortly after sunrise and sunset at ∼07:00 and ∼20:00, rather than at ∼15:00 as occurs in the second year) and also far more sol-to-sol repeatability. The wind direction is also more westerly than northerly in the first year than in the second. These major year-to-year difference are due to the presence of the global 1977b dust storm in the first year, which strengthened the large-scale circulation, altered the relative amplitude of different tidal modes, and reduced the contribution of high-latitude frontal winds to the pattern here. From sol to sol, the wind direction is most repeatable at L s = 105°at both VL sites ( Fig. 11a and 11e), but appears to be most variable from sol to sol at the PHX site at the same season (Fig. 11h), most likely due to the strong influence of baroclinic waves at this high latitude. Wind directions at the MSL site during the Bagnold Dunes wind characterization investigation are also rather repeatable from sol to sol, most likely due to the strong topographic control (described below). Figure 12 demonstrates that, for all landing sites and seasons other than MSL, the wind direction exhibits continuous turning over the majority of each sol. At the VL1 site there is predominantly anti-clockwise turning from 22:00 to 16:00, with southeasterlies turning to strong southwesterlies by late morning, then back to southeasterlies again. By contrast, there is predominantly clockwise turning at the VL2 site in all seasons shown. This is likely due to the different phasing of the tides relative to the background circulation and the stronger Coriolis effect at the higher latitude of VL2 (e.g. Haberle et al. 1993b). The wind turns clockwise through a full 360°in some seasons (L s ∼ 345°and 105°, Figs. 12d and e). MPF wind direction data are not shown here, but wind direction as a function of time of day for 30 sols following L s = 142°was presented in Fig. 6 of . The MPF winds also turned clockwise through a full 360°rather steadily and very repeatedly throughout these sols, which might suggest that they were largely controlled by tidally-induced flow of mass at this relatively flat, near-equatorial landing site. However, the phasing is inconsistent with this explanation, suggesting either that the tidal forcing in this region is more complex or that the flow is more affected by slope winds. Indeed, southerly or southeasterly winds appear to dominate overnight, which is consistent with downslope winds related to the regional topography. The situation at PHX is more complex but seems to show clockwise turning overall. By contrast, at MSL winds are dominated by topographic flows at the season shown (late southern fall), due to its location on the NW slope of Aeolis Mons, the large central peak inside Gale Crater. As demonstrated by Figs. 11 and 12i, there is a clear transition from upslope winds during the day (from the NW/N, ∼270 to 360°) to downslope winds during the night (from the SE, ∼135°), with only short periods of rapid turning in between (e.g. from ∼16:30-18:30).

Near-Surface Relative Humidity and Water Vapor Content
Of all spacecraft landed on Mars, only the PHX and MSL carried sensors to measure the near-surface relative humidity (Table 1) (Zent et al. 2009;Gómez-Elvira et al. 2012). For Among the full set of REMS RH measurements, only those taken during the first four seconds after the RHS has been turned on after at least ∼5 min of inactivity are considered here. After four seconds, the local RH values are affected by the heating effect of the sensor and the local RH measured by the sensor significantly departs from the actual RH of the atmosphere (Sect. 2.5). Reliable values include RH measurements taken during the nominal and the so-called high-resolution interval mode (HRIM), which consists of alternately switching the sensor on and off at periodic intervals to minimize heating and is only used on selected one to two hour observation blocks. As an example, Fig. 13 shows in gray the complete set of RH measurements made on sol 1137 (L s ∼ 56°), with RH values with the highest confidence highlighted in green (HRIM mode) and blue (nominal mode).
Using REMS RH measurements selected as described above, we show in the Fig. 14 (top) the interannual and seasonal variability of the daily maximum RH value during the first 1526 sols of the MSL mission. REMS-derived values of VMR, calculated as VMR = RH × e s (T )/P using the daily maximum RH value and simultaneous measurements of temperature and pressure, where e s is the saturation vapor pressure over ice, are shown in Fig. 14 (bottom). The daily maximum RH generally occurs between 04:00 and 06:00 LMST throughout the year, and thus VMR values shown in Fig. 14 (bottom) represent nighttime estimates. Daytime REMS-derived VMR values present large uncertainties because of the extremely low RH (with daytime values <5% between 10:00 and 18:000 LMST throughout The daily maximum RH is generally achieved between 04:00 and 06:00 LMST, with daytime values staying below 5% between 10:00 and 18:00 LMST throughout the mission. (Bottom) Interannual and seasonal evolution of the most accurate VMR value per sol at the MSL landing site. This value is obtained at the same time as the RH shown above the year) and thus VMR cannot be obtained reliably during the day .
The large seasonal variability of the daily maximum RH shown in Fig. 14 (top) is mainly controlled by near-surface temperatures and not by atmospheric water content. The highest annual values of about 70% occur in early winter (L s = 93°), coinciding with the minimum annual near-surface air and ground temperatures (Fig. 5, top and Fig. 7), while the lowest annual (daily maximum) RH values of ∼10% occur around late spring and early summer, coinciding with the highest annual near-surface air and ground temperatures. Interestingly, the seasonal evolution of the daily maximum RH is anti-correlated with that of the nearsurface nighttime VMR shown in Fig. 14 (bottom). The lowest annual values of VMR occur in late fall and winter, when the daily maximum RH is the highest, while the highest annual values of VMR occur in late winter, spring and early summer, when the daily maximum RH is the lowest.
The sol-to-sol variability in VMR seems to be larger in spring and summer than in late fall and winter (Fig. 14, bottom). This variability may be related to the uncertainty in VMR estimations rather than to actual variations in VMR . This occurs because the most accurate daily VMR values are obtained when the relative humidity is the highest. Therefore, VMR values obtained in spring and summer, when even the diurnal maximum RH is low, have greater uncertainties than VMR values obtained in late fall and winter, when the diurnal maximum RH is the highest.
The detection of surface frost or near-surface fog is important because it allows the independent estimation of the water vapor content if the temperature is known. Nighttime frost was detected at the PHX site by the Robotic Arm Camera from sol 80 (early summer, L s = 113°) until the end of the mission on sol 151 (late summer, L s = 149°) (Smith et al. 2009), coinciding with a gradual decrease in atmospheric pressure and ground and nearsurface air temperatures (Fig. 3 top, Fig. 5 top, Fig. 8 bottom). Also during the second half of the PHX mission, the LIDAR instrument observed two regular daily patterns of water-ice cloud formation: a surface-based cloud (ice fog) forming near midnight, and a 4-km-high cloud layer forming around 01:00 LMST (Whiteway et al. 2009;Dickinson et al. 2010;Moores et al. 2011). Using values of temperature and atmospheric pressure near the surface at midnight, and at 4 km height at 01:00 LMST, the evolution of the VMR at the PHX site can be inferred. Near the surface, we obtain an estimate for the near-surface VMR that decreases from ∼210 ppm on sol 80 to ∼15 ppm on sol 150. These values were obtained using MET midnight measurements of air temperature at 2 m above the surface and atmospheric pressure at the lander deck (200 K and 775 Pa on sol 80 and 185 K and 734 Pa on sol 150). At 4 km height, Dickinson et al. (2010) obtained an estimate for the VMR decreasing from ∼1300 ppm (6.8 mg/m 3 at −64.3°C and ∼500 Pa) on sol 99 to ∼575 ppm (2.9 mg/m 3 at −70.2°C and ∼470 Pa) on sol 150 ). The lower values and larger relative decrease in VMR near the surface highlight the key role played by the deposition and sublimation of H 2 O onto and from the ground in the seasonal H 2 O cycle at the PHX site. Additionally, the decrease in VMR estimated at 4 km height is consistent with that in the column-integrated water abundance at PHX site measured by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) instrument on the Mars Reconnaissance Orbiter (MRO) between sols 99 and 150 , indicating that the H 2 O seasonal cycle at the PHX site is mostly confined to the lowest 4 km of the atmosphere .
At the Viking Lander 2 site, a thin layer of frost was observed on the ground for about 250 sols, from late fall (L s = 230°) until early spring (L s = 16°) in the first year of the mission (Jones et al. 1979;Wall 1981). At this location, surface frost formed at night and persisted throughout the day. Since the VL missions did not carry a ground temperature sensor, an estimate for the water vapor VMR at the ground cannot be obtained. Furthermore, air temperatures measured at 1.6 m (Fig. 5, top) during this period, which ranged from about 160 to 190 K at night, cannot be used unambiguously to determine the near-surface water vapor VMR because near-surface fog was not detected. However, Ryan et al. (1982) estimated VMR values at both VL sites indirectly. By attributing an increase in the rate of nocturnal air temperature decrease to ice fog formation (thus representing the frost point temperature), they inferred values of the VMR between 01:00 and 03:00 LMST throughout the year at both sites. At the VL1, nocturnal VMR values ranged from ∼17 ppm in late autumn (L s ∼ 250°) to ∼200 ppm in mid summer (L s ∼ 130°), while at the VL2 VMR values ranged from <1 ppm in early winter (L s ∼ 280°) to ∼200 ppm in late spring (L s ∼ 60°).
Ground frost has not been observed in situ at lower latitudes such as those of the VL1, MER, MPF and MSL missions landing sites. However, possible visual evidence for the formation of frost has been found on the deck and calibration target of the MER-B on sol 257 (L s ∼ 105°) (Möhlmann 2008). Additionally, near-surface fog was imaged by the Mars Express High Resolution Stereo Camera in Valles Marineris (∼10°S) (Möhlmann et al. 2009). At the MSL site, indirect evidence for frost formation on a few sols during late fall and early winter in MY 32 has been provided from the analysis of relative humidity and ground temperature measurements , further supported by numerical modeling (Savijärvi et al. 2015;Steele et al. 2017).

Discussion
In this section we briefly describe the meteorological stations being developed for future in-situ measurements on Mars, propose measurements to improve our understanding of the Martian dust and H 2 O cycles, and discuss the potential for liquid water formation under Mars' present day conditions and its implications for future Mars missions.

Future Meteorological Stations
In the next few years, in-situ meteorological measurements at the surface of Mars will continue with the Temperature and Wind Sensors (TWINS) instrument onboard NASA's 2018 InSight mission (Velasco and Rodríguez-Manfredi 2015), the Meteorological (METEO) and the Habitability, Brine, Irradiation and Temperature (HABIT) instruments onboard the Ex-oMars 2020 surface platform (ESA 2016), and the Mars Environmental Dynamics Analyzer (MEDA) instrument onboard NASA's Mars 2020 mission ). In addition to extending the current record of meteorological measurements to other locations on Mars, HABIT will measure changes in the hydration state (including potential formation of liquid brine) of hygroscopic salts exposed to Martian air to shed light on the water cycle, while MEDA will measure the surface energy balance (upward and downward long-wave and short-wave radiation fluxes) and dust particle morphology to characterize the energy exchange with the atmosphere and the radiative properties of dust particles.

Improving Our Understanding of the Dust and H 2 O Cycles
The characterization of dust-lifting processes is important to improve our knowledge of the Martian dust cycle and the various feedback processes associated with it. The Dust Characterization, Risk Assessment and Environment Analyzer on the Martian Surface (DREAMS) instrument onboard the ExoMars 2016 surface lander (Bettanini et al. 2014) was designed to shed light on dust-lifting processes, but it was unfortunately lost during descent. Windblown sand, or 'saltation,' causes wind erosion and ejects dust aerosols into the Martian atmosphere through impact of saltating particles on the ground (Greeley and Iversen 1985;Almeida et al. 2008). Electrification of wind-blown sand, dust devils and dust storms could produce electric discharges (Eden and Vonnegut 1973) and reduce the wind stress required to lift particles from the surface of Mars (Kok and Renno 2006). In addition, dust lifting might have important implications for habitability because, along with the transport of minerals by the wind, aeolian processes provide a means for the exchange of nutrients and energy sources between the atmosphere and the shallow subsurface (Farrell et al. 2004). The growth of segregated ice and salt crystals in rock pores is a critical weathering process on Earth, and there is evidence that this process is currently active on Mars (Hallet 2006). To improve our understanding of dust-lifting processes, we propose measurements of the particle mass flux, the near surface electric field, and the depth of the saltation layer at heights ranging from 0.01 to 2 m above the ground Renno 2009a, 2009b;Esposito et al. 2016), in combination with measurements of the near surface wind speed, the wind shear, and the soil particle size distribution and particle cohesion (Kok and Renno 2006).
Determining the relative importance of the exchange of physically adsorbed water between the regolith and the atmosphere at different time scales is important to improve our understanding of the current H 2 O cycle. Although there is consensus that adsorbed water exchanges seasonally between the regolith and the atmosphere, the significance of such an exchange at diurnal time-scales remains unclear (Zent et al. 1993;Titov et al. 1995;. In particular, it is not known whether kinetics would allow for the exchange of adsorbed water on the hourly time scales necessary for tracking diurnal insolation (Titov et al. 1995;. Unfortunately, ground-based estimations of near-surface water VMR with full diurnal coverage are not available (Sect. 3.6). Processed PHX/TECP atmospheric RH values are not available in the NASA PDS (Sect. 2.4), while MSL/REMSderived daytime VMR values present large uncertainties and cannot be used to address the diurnal cycle (Savijärvi et al. 2015;. However, the strong correlation between REMS-derived nighttime VMR values and surface and near-surface temperature values at Gale crater can be used to shed some light on the exchange of water between the regolith and the atmosphere. The lowest annual values of VMR (Fig. 14, bottom) coincide with the lowest annual values of surface and near-surface air temperature (Fig. 5, top and Fig. 7), suggesting that the regolith plays a significant role in exchanging water vapor with the atmosphere at Gale crater. Two processes could account for such an exchange: adsorption/desorption of water vapor and deposition/sublimation of frost. Adsorption by the regolith is expected to be more intense at lower temperatures, favoring this mechanism as a driver for the exchange of water vapor between the regolith and the atmosphere at Gale. In addition,  show that diurnal cycles of RH measured at Gale can be more accurately simulated when a temperature dependent soil adsorption scheme is used. Potential nighttime frost formation in late fall and early winter Steele et al. 2017) might also explain the annual lowest values of near-surface nighttime VMR at this time of the year (Fig. 14,  bottom), although as of February 2017 frost has not been detected by the MSL Curiosity rover.
While the strong correlation between VMR and temperature is consistent with atmosphere-regolith interactions, large-scale seasonal circulation patterns have the potential to create seasonal differences in the near-surface VMR values. Rafkin et al. (2016) show that the air mass in Gale mixes well with external air during southern summer (L s ∼ 270°), coinciding with relatively high values of near-surface VMR (Fig. 14, bottom). On the contrary, the mixing between the air in and external to Gale is more limited in other seasons, and in particular during southern winter (L s ∼ 90°) when the near-surface VMR shows the lowest annual values. However, air entering Gale crater during southern summer (when the North Pole is a water trap) has a northerly component (Rafkin et al. 2016), and thus the air advected into Gale is not expected to be particularly moist. Therefore, soil-atmosphere interactions seem to play a prevalent role in exchanging water vapor at Gale.
To improve our understanding of the exchange of water between the regolith and the atmosphere, we propose measurements of near-surface water vapor content throughout the diurnal cycle. Such measurements could be made either by a relative humidity sensor with accuracy better than 1% during the daytime (so that daytime derived-VMR are reliable given typical daytime values of RH and air temperatures at Gale crater), or by a sensor capable of measuring VMR directly. Additionally, we propose simultaneous measurements of the vertical component of the wind speed at high frequency (>1 Hz) to calculate vertical turbulent transport of water vapor. Currently planned missions will not directly address dust lifting processes and the exchange of water vapor between the atmosphere and the regolith. Thus, we propose that missions past 2020 study these processes.

Liquid Water and the H 2 O Cycle
Liquid water at the surface and in the shallow subsurface of Mars is not required to account for current seasonal variations in the integrated H 2 O column abundance observed from orbit (Jakosky 1985;Richardson and Wilson 2002), which, in combination with modeling, have been used to determine the relative importance of the various processes controlling the seasonal H 2 O cycle. However, the potential presence of liquid water has important implications for the current habitability of Mars.
Mars has large reservoirs of water as ice and vapor, but pure liquid water is unlikely to be present, even temporarily, on the surface of Mars because it would quickly evaporate in the regions where it neither freezes nor boils (Ingersoll 1970;Haberle et al. 2001;Martínez and Renno 2013). However, the presence of water ice in the shallow subsurface midlatitude and polar regions Mitrofanov et al. 2002;Byrne et al. 2009), along with the detection of perchlorates in polar and equatorial soil Glavin et al. 2013;Ming et al. 2014) and of chloride-bearing deposits in the southern highlands at low and midlatitudes (Osterloo et al. 2008), is important because they can melt this ice at Mars' present-day environmental conditions and produce liquid saline water (brine) (Clark 1978;Brass 1980;Clark and Van Hart 1981;Haberle et al. 2001;Chevrier and Altheide 2008;Rennó et al. 2009;McEwen et al. 2011;Ojha et al. 2015).
Two mechanisms have been suggested for brine formation on Mars: its formation by salts absorbing atmospheric water vapor (deliquescence) when the relative humidity exceeds a threshold value known as the deliquescence relative humidity and the temperature is above the salts' eutectic value (Clark 1978;Rennó et al. 2009;Davila et al. 2010;Gough et al. 2011;Nuding et al. 2015;Martín-Torres et al. 2015;Nikolakakos and Whiteway 2015), and its formation by salts melting water ice when the temperature exceeds the eutectic value of salts in contact with water ice (Brass 1980;Clark and Van Hart 1981;Fairén et al. 2009;. Perchlorate salts, in particular Ca(ClO 4 ) 2 , Mg(ClO 4 ) 2 and Na(ClO 4 ) salts, are extremely relevant for brine formation via deliquescence and melting because they have low eutectic temperatures and were found in polar and equatorial regions Glavin et al. 2013), suggesting that they are distributed globally.
The necessary conditions for both mechanisms to occur are more likely to be met in the Martian polar regions than at low and equatorial latitudes. In the polar regions, brine formation by the melting of ice in contact with salts is plausible because frost or snow are seasonally deposited on saline soils (Whiteway et al. 2009;Martínez et al. 2012), water ice is seasonally present in the shallow subsurface (Paige 1992;Mellon and Jakosky 1993;Schorghofer and Aharonson 2005;Cull et al. 2010), and temperatures exceed the eutectic value of salts detected on Mars during a significant fraction of the sol (from minutes to hours depending on the location and time of the year) (Möhlmann 2011;Nuding et al. 2014;. Moreover, numerical modeling and laboratory experiments indicate that brine formation by deliquescence in the polar region is theoretically possible because simulated diurnal cycles of RH and temperature at the ground at the PHX site are compatible with deliquescence of Ca(ClO 4 ) 2 and Mg(ClO 4 ) 2 salts (Möhlmann 2011;Nuding et al. 2014;Nikolakakos and Whiteway 2015). However, other laboratory experiments indicate that deliquescence of perchlorates at the PHX site is not rapid enough to occur during the short periods of the day during which the environmental conditions are favorable for deliquescence because the kinetics of this reaction is too slow at Mars relevant values of temperature and water vapor pressure .
At equatorial and low latitudes, brine formation by the melting of ice is unlikely because water ice is not thermodynamically stable in the shallow subsurface at these latitudes (Mellon and Jakosky 1993; Mellon et al. 2004;Schorghofer and Aharonson 2005). Similarly, brine formation by deliquescence is less likely to occur in these regions than at higher latitudes. Analyses of simultaneous measurements of the highest confidence possible of RH and ground temperature at Gale crater (Fig. 15, gray points) indicate that the environmental Fig. 15 Stability diagram of Ca(ClO 4 ) 2 , Mg(ClO 4 ) 2 and NaClO 4 salts, with superimposed values of relative humidity and temperature at the ground colored in gray for the first 1526 sols of the MSL mission. The environmental conditions at the MSL landing site are unlikely to lead to brine formation of perchlorate salts via deliquescence because they are too dry and cold. For each salt, the colored thick-dashed line represents the deliquescence relative humidity at which the various salts form aqueous solutions. The horizontal blue line represents the eutectic temperature of Ca(ClO 4 ) 2 . The eutectic temperatures of Mg(ClO 4 ) 2 and NaClO 4 are also shown in orange and gray colors. For simplicity, transitions between different hydrations states are not shown. Note that the relative humidity at the ground has been calculated as [RH 1.6 m × e s (T 1.6 m )]/e s (T g ), where RH 1.6 m and T 1.6 m are the relative humidity and air temperature at 1.6 m measured by REMS and e s is the saturation vapor pressure over ice. Therefore, the water vapor pressure has been assumed constant in the first 1.6 m conditions are unlikely to lead to brine formation of perchlorates by deliquescence. Martín-Torres et al. (2015) showed that the relative humidity and ground temperature at Gale crater are compatible with deliquescence of Ca(ClO 4 ) 2 salt. However, RH values used in that study were higher than those obtained using the most accurate REMS RH sensor recalibration performed in August 2016. Also, the theoretical value for the eutectic temperature of Ca(ClO 4 ) 2 of 194.5 K assumed by Martín-Torres et al. (2015) is lower than ∼199 K, the value obtained experimentally in recent laboratory experiments (Marion et al. 2010;.
Testing hypotheses for brine formation on Mars is important to support future missions focusing on habitability and the search for life. Moreover, human exploration of Mars is expected to occur in the next few decades and in preparation for it, it is important to understand the resources available for sustaining humans on Mars. Therefore, a future mission to the polar regions could have a greater potential than a mission to equatorial and low latitudes to shed light on liquid water formation under present-day Martian environmental conditions.
A variety of instruments have been suggested to investigate brine formation on Mars. A Raman spectrometer can be used to clearly differentiate between ice and liquid brine ). One such instrument, the Raman Laser Spectrometer (RLS) will fly on the 2020 ExoMars Rover (Hutchinson et al. 2014) and another one in the same year in SHERLOC on Mars 2020 (Beegle et al. 2015). Although their main goal is to look for geological biosignatures, both have the capability to detect liquid brine. Brine on the Martian surface or in the shallow subsurface can also be detected by measuring soil conductivity (Heinz et al. 2016), microwave resonance frequencies (Zandonadi et al. 2015), or radar reflections. RIMFAX on the Mars 2020 rover (Hamran et al. 2015) and WISDOM on the ExoMars rover (Ciarletti et al. 2011) are future ground penetrating radars with the ability to investigate possible liquid water in the shallow Martian subsurface.

Summary and Conclusions
We have analyzed the complete set of in-situ meteorological data obtained from the Viking landers in the 1970s to today's Curiosity rover. We have processed these data to provide values of the highest confidence possible of atmospheric opacity, atmospheric pressure, nearsurface air temperature, ground temperature, near-surface wind speed and direction, and near-surface air relative humidity and water vapor content. Then we have studied the diurnal, seasonal and interannual variability of these quantities over a span of more than twenty Martian years to review our understanding of the modern near-surface climate of Mars.
The seasonal and latitudinal variation of solar insolation caused by the Martian obliquity and the large eccentricity of Mars' orbit drives the dust, CO 2 and H 2 O cycles, which, coupled to radiative and dynamic processes, regulate the modern climate of Mars. Dust is ubiquitous in the Martian atmosphere and interacts strongly with solar and thermal radiation, affecting weather and climate and playing a major role in regional and global circulations. In contrast to Earth, the spatial and temporal variability of the dust content of the Martian atmosphere is extremely large, and this effect is enhanced by the low mass of the Martian atmosphere. During the aphelion season, dust opacity values are low and the intraseasonal and interannual variability is small at each landing site (Fig. 2, left column). In contrast, during the perihelion season the opacity values and the intraseasonal and interannual variability are larger (Fig. 2, left column), regardless of the hemispheric location. Repeatable regional dust storms during the perihelion season between L s = 210°and 240°and between L s = 320°and 340°produce seasonal increases in the atmospheric opacity at each landing site every year, while sporadic global dust storms between L s = 280°and 300°produce larger increases in atmospheric opacity every few years (Fig. 2, left column). During global dust storms, dust opacity values are extraordinarily high and the surface environmental conditions are significantly altered, with dramatic increases in the diurnal amplitude of the surface pressure (Fig. 3, bottom, MY 12) and dramatic decreases in the diurnal amplitude of the near-surface air temperature (Fig. 5,bottom,MY 12). However, values of the daily mean near-surface air temperature are not significantly changed because lower daytime temperatures are counterbalanced by greater nighttime temperatures (Fig. 5,top,MY 12).
The seasonal deposition of about 30% of the mass of the atmosphere in the polar caps during the cold seasons and its sublimation during the warm seasons strongly affects surface pressures (Fig. 3, top). At each landing site, the surface pressure shows an absolute annual maximum at L s ∼ 260°(sublimation of the southern polar cap) and an absolute annual minimum at L s ∼ 150°(deposition onto the southern polar cap). Additionally, the surface pressure presents a relative maximum at L s ∼ 55°(sublimation of the northern polar cap) and a relative minimum at L s ∼ 345°(deposition onto the northern polar cap). The larger sol-to-sol variability during the perihelion season than during the aphelion season at each landing site is caused by baroclinic instability in northern autumn and winter, which generates perturbations in surface pressure that peak at the latitude of the polar jet stream. Therefore, the largest sol-to-sol variations occur at the VL2 site (closest mission to the polar jet), followed by those at the VL1 (northern edge of the tropics) site, with the lowest variability in this sequence at the MSL site in the southern tropics. The diurnal amplitude of the surface pressure, also called the diurnal atmospheric tide, increases with increasing atmospheric aerosol content and proximity to the equator. Therefore, it shows larger values in the perihelion season than during the aphelion season at each landing site, with the largest values observed at the MSL site due to its proximity to the equator (Fig. 3, bottom).
The near-surface thermal environment at each landing site is primarily governed by the solar insolation reaching the surface. At each landing site, the maximum and minimum annual near-surface air and ground temperature (Fig. 5, top and Fig. 7) coincide with the annual maximum and minimum in surface solar insolation (Fig. 2). Following stronger seasonal variations in surface solar insolation at higher latitudes (Fig. 2), the annual amplitude in near-surface air temperature at each site increases with increasing latitude (Fig. 5, top). As for the surface pressure, larger sol-to-sol variations in the near-surface air temperature during the perihelion season are caused by baroclinic eddies in the northern hemisphere during autumn and winter. The diurnal amplitude of ground temperature at the MSL site (the only mission with full diurnal coverage) is primarily governed by the type of terrain, and in particular by its thermal inertia. Changes in atmospheric opacity can also regulate the diurnal amplitude of ground temperature, but their impact has been smaller than that produced by changes in thermal inertia, at least for the moderately low opacity values measured at Gale as of February 2017.
Of all the variables examined here, wind speed and direction are most strongly influenced by local and regional topography. This is particularly true for the MSL rover, where daytime upslope/nighttime downslope winds dominate the results shown (Figs. 10-12i). They are typically least affected by local variations in atmospheric opacity, although regional and global scale changes may greatly affect winds via their impact on the global circulation, thermal tides, and wave activity. For example, the altered circulation 1977b global dust storm produced faster wind speeds (Fig. 9), less sol-to-sol variability, and a complete change in the diurnal variation of wind speed and direction measured at the VL2 landing site, compared to the same time the following year (cf. Figs. 10-12c and 12g). This was due to the much stronger global circulation, suppression of baroclinic wave activity, and altered amplitudes of tidal components.
Only the PHX and MSL carried sensors to measure the near-surface relative humidity. At the MSL site, the daily maximum RH occurs between 04:00 and 06:00 LMST (Fig. 13) throughout the year, coinciding with the coldest time period of the day. Daytime values are extremely low, staying below <5% between 10:00 and 18:00 LMST regardless of the season. The large seasonal variability of the daily maximum RH (Fig. 14, top) is controlled by near-surface temperatures. The highest annual values of about 70% occur in early winter (L s = 93°), coinciding with the minimum annual near-surface air and ground temperatures (Fig. 5, top and Fig. 7). Similarly, the lowest annual (daily maximum) RH values of ∼10% occur around late spring and early summer, coinciding with the highest annual near-surface air and ground temperatures. At the PHX site, processed values of RH are not available in the NASA PDS, but LIDAR and SSI observations show that the ground and near-surface air were saturated at nighttime throughout the second half of the mission.
The near-surface water content has not been directly measured by any rover or lander, but it can be derived indirectly either from simultaneous measurements of RH, temperature and pressure at the MSL site, or from the detection of surface frost or near-surface fog (provided that the temperature is known) at the PHX site. At the MSL site, daytime REMS-derived VMR values present large uncertainties due to extremely low RH, and only nighttime values can be derived reliably. The lowest annual values (<5 ppm) occur in late fall and winter (Fig. 14 bottom), when the daily maximum RH is the highest (Fig. 14, top) and surface and near-surface air temperatures are the coldest (Fig. 5, top and Fig. 7), while the highest annual values of VMR occur in late winter, spring and early summer (∼70 ppm), when the daily maximum RH is the lowest and surface and near-surface air temperatures are the warmest. The strong correlation between nighttime values of VMR and temperature suggests that the regolith plays a significant role in exchanging water vapor with the regolith, with adsorption/desorption as the most likely mechanism. At the PHX site, the LIDAR instrument observed two regular daily patterns of water-ice cloud formation from sol 80 (early summer, L s = 113°) until the end of the mission on sol 151 (late summer, L s = 149°): a surfacebased cloud forming near midnight, and a 4-km-high cloud forming at about 01:00 LMST. Assuming saturated conditions at both heights and times between sols 80 and 150, nighttime VMR values decreasing from ∼210 ppm to ∼15 ppm near the surface, and from ∼1300 ppm to ∼575 ppm at 4 km high are obtained. The lower values and stronger decrease in VMR near the surface highlight the key role played by the deposition and sublimation of H 2 O onto and from the ground in the seasonal H 2 O cycle at the PHX site.
The characterization of dust-lifting processes and the exchange of adsorbed water between the regolith and the atmosphere at different time scales are important to improve our knowledge of the Martian dust and H 2 O cycles, but currently planned missions will not directly address these processes. We have proposed in Sect. 5 key measurements to study these processes and shed light on both cycles. Furthermore, we have shown that the necessary conditions for brine formation to occur, either by deliquescence or melting, are more likely to be met in the Martian polar regions than at low and equatorial latitudes. Therefore, a future mission to the polar regions could have a greater potential than a mission to equatorial and low latitudes to shed light on liquid water formation under present-day Martian environmental conditions.