The glacial–interglacial paradox: the relation between the global means of sea level and sea surface temperature

Sea level variability during the Quarternary is simulated using a stochastic climate model, and a sensitivity relation for the change in net oceanic evaporation due to a change in sea surface temperature. In the application of this relation, it is assumed that the greater part of the change in net oceanic evaporation causes changes in the land ice storage, rather than being directly returned to the ocean by rivers. The analysis suggests that the observed sea level changes can be interpreted as due to the transfer of heat to the deep ocean from the surface mixed layer, arising from random radiation perturbations of the same variance as would give rise to the interannual variability of the global temperature series. The paradox is that glacial conditions (increase in ice storage) are favoured by positive (temperate) sea surface temperature anomalies, and interglacial conditions (decrease in ice storage) by negative (temperate) sea surface temperature anomalies. The evolution of both these regimes, which are inherently unstable, appears to be controlled by the deep water formation process, while albedo feedback is of minor importance. Fluvial feedback, (in which as the ice storage increases the fluvial inflow decreases), however, is found to be an important process, and a small sensitivity of river inflow to storage is consistent with forcing by random variability or by astronomical forcing. A simple analytical model incorporating the key processes of oceanic evaporation and fluvial feedback is presented. The analysis points to the importance of an accurate river model for climate system modelling.     

Thermal migration of water on the Galilean satellites

We have modeled the thermal migration of water on the Galilean satellites under the assumption of ballistic molecular trajectories. We find that water migrating owing to solar radiation on an ice-covered satellite will build up in the temperate latitudes, in general not reaching the poles. As much as 50 m of ice may have been lost by this process from the equatorial regions of Europa over the age of the solar system. The disappearance of patches of ice—for instance, the bright rays surrounding some impact craters—from the equatorial regions of Ganymede and Callisto may approach a value (the irreversible evaporation rate) three orders of magnitude larger than the net equatorial loss rate for ice-covered Europa. The presence of water ice pole caps on Ganymede extending to the latitudes at which thermal migration becomes important suggests that some process distributed an extensive, thin covering of water on the satellite, and that the equatorial regions were subsequently cleared by the thermal process.     

Generation of liquid water on Mars through the melting of a dusty snowpack

The possibility that snowmelt could have provided liquid water for valley network formation early in the history of Mars is investigated using an optical-thermal model developed for dusty snowpacks at temperate latitudes. The heating of the postulated snow is assumed to be driven primarily by the absorption of solar radiation during clear sky conditions. Radiative heating rates are predicted as a function of depth and shown to be sensitive to the dust concentration and the size of the ice grains while the thermal conductivity is controlled by temperature, atmospheric pressure, and bulk density. Rates of metamorphism indicate that fresh fine-grained snow on Mars would evolve into moderately coarse snow during a single summer season. Results from global climate models are used to constrain the mean-annual surface temperatures for snow and the atmospheric exchange terms in the surface energy balance. Mean-annual temperatures within Martian snowpacks fail to reach the melting point for all atmospheric pressures below 1000 mbar despite a predicted temperature enhancement beneath the surface of the snowpacks. When seasonal and diurnal variations in the incident solar flux are included in the model, melting occurs at midday during the summer for a wide range of snow types and atmospheric pressures if the dust levels in the snow exceed 100 ppmw (parts per million by weight). The optimum dust concentration appears to be about 1000 ppmw. With this dust load, melting can occur in the upper few centimeters of a dense coarse-grained snow at atmospheric pressures as low as 7 mbar. Snowpack thickness and the thermal conductivity of the underlying substrate determine whether the generated snow-melt can penetrate to the snowpack base, survive basal ice formation, and subsequently become available for runoff. Under favorable conditions, liquid water becomes available for runoff at atmospheric pressures as low as 30 to 100 mbar if the substrate is composed of regolith, as is expected in the ancient cratered terrain of Mars.     

MAIC-2, a latitudinal model for the Martian surface temperature, atmospheric water transport and surface glaciation

The Mars Atmosphere-Ice Coupler MAIC-2 is a simple, latitudinal model, which consists of a set of parameterisations for the surface temperature, the atmospheric water transport and the surface mass balance (condensation minus evaporation) of water ice. It is driven directly by the orbital parameters obliquity, eccentricity and solar longitude (L_s) of perihelion. Surface temperature is described by the Local Insolation Temperature (LIT) scheme, which uses a daily and latitude-dependent radiation balance. The evaporation rate of water is calculated by an expression for free convection, driven by density differences between water vapor and ambient air, the condensation rate follows from the assumption that any water vapour which exceeds the local saturation pressure condenses instantly, and atmospheric transport of water vapour is approximated by instantaneous mixing. Glacial flow of ice deposits is neglected. Simulations with constant orbital parameters show that low obliquities favour deposition of ice in high latitudes and vice versa. A transient scenario driven by a computed history of orbital parameters over the last 10 million years produces essentially monotonically growing polar ice deposits during the most recent 4 million years, and a very good agreement with the observed present-day polar layered deposits. The thick polar deposits sometimes continue in thin ice deposits which extend far into the mid latitudes, which confirms the idea of “ice ages” at high obliquity.     

The effects of reduced land fraction and solar forcing on the general circulation: results from the NCAR CCM

Land fraction and the solar energy at the top of the atmosphere (solar constant) may have been significantly lower early in Earth's history. It is likely that both of these factors played some important role in the climate of the early earth. The climate changes associated with a global ocean(i.e. no continents) and reduced solar constant are examined with a general circulation model and compared with the present-day climate simulation. The general circulation model used in the study is the NCAR CCM with a swamp ocean surface. First, all land points are removed in the model and then the solar constant is reduced by 10% for this global ocean case.Results indicate that a 4 K increase in air temperature occurs with global ocean simulation compared to the control. When solar constant is reduced by 10% under global ocean conditions a 23 K decrease in air temperature is noted. The global ocean warms much of the troposphere and stratosphere, while a reduction in the solar constant cools the troposphere and stratosphere. The largest cooling occurs near the surface with the lower solar constant.Global mean values of evaporation, water vapor amounts, absorbed solar radiation and the downward longwave radiation are increased under global ocean conditions, while all are reduced when the solar constant is lowered. The global ocean simulation produces sea ice only in the highest latitudes. A frozen planet does not occur when the solar constant is reduced—rather, the ice line settles near 30° of latitude. It is near this latitude that transient eddies transport large amounts of sensible heat across the ice line acting as a negative feedback under lower solar constant conditions keeping sea ice from migrating to even lower latitudes.Clouds, under lower solar forcing, also act as a negative feedback because they are reduced in higher latitudes with colder atmospheric temperatures allowing additional solar radiation to reach the surface. The overall effect of clouds in the global ocean is to act as a positive feedback because they are slightly reduced thereby allowing additional solar radiation to reach the surface and increase the warming caused by the removal of land. The relevance of the results to the “Faint-Young Sun Paradox” indicates that reduced land fraction and solar forcing affect dynamics, heat transport, and clouds. Therefore the associated feedbacks should be taken into account in order to understand their roles in resolving the “Faint-Young Sun Paradox”.     

Evaporation of ice in planetary atmospheres: Ice-covered rivers on mars

The evaporation rate of water ice on the surface of a planet with an atmosphere involves an equilibrium between solar heating and radiative and evaporative cooling of the ice layer. The thickness of the ice is governed principally by the solar flux which penetrates the ice layer and then is conducted back to the surface. These calculations differ from those of Lingenfelter et al. [(1968) Science161, 266–269] for putative lunar channels in including the effect of the atmosphere. Evaporation from the surface is governed by two physical phenomena: wind and free convection. In the former case, water vapor diffuses from the surface of the ice through a lamonar boundary layer and then is carried away by eddy diffusion above, provided by the wind. The latter case, in the absence of wind, is similar, except that the eddy diffusion is caused by the lower density of water vapor than the Martian atmosphere. For mean Martian insolations the evaporation rate above the ice is ～ 10^?8 g cm^?2 sec^?1. Thus, even under present Martian conditions a flowing channel of liquid water will be covered with ice which evaporates sufficiently slowly that the water below can flow for hundreds of kilometers even with quite modest discharges. Evaporation rates are calculated for a wide range of frictional velocities, atmospheric pressures, and insolations and it seems clear that at least some subset of observed Martian channels may have formed as ice-choked rivers. Typical equilibrium thicknesses of such ice covers are ～ 10 to 30 m; typical surface temperatures are 210 to 235°K. Ice-covered channels or lakes on Mars today may be of substantial biological interest. Ice is a sufficiently poor conductor of heat that sunlight which penetrates it can cause melting to a depth of several meters or more. Because the obliquity of Mars can vary up to some 35°, the increased polar heating at such times seems able to cause subsurface melting of the ice caps to a depth which corresponds to the observed lamina thickness and may be responsible for the morphology of these polar features.     

Latitudinal distribution of temporary liquid cryobrines on Mars

Simulations of the surface temperature and atmospheric humidity with a modern Mars climate model (MCD) and with Phoenix data are used to study the conditions for a liquefaction of brines as a function of latitude and season. The results show that, in the presence of appropriate salts, liquid cryobrines can in course of the diurnal cycle temporarily evolve at high latitudes on Mars’ current climate. The conditions for the liquefaction of “Mars-relevant” cryobrines and time and duration of their stability during the diurnal cycle are calculated for northern spring and for the Phoenix landing site.     

Temporary liquid water in upper snow/ice sub-surfaces on Mars?

It is investigated whether conditions for melting can be temporarily created in the upper sub-surface parts of snow/ice-packs on Mars at subzero surface temperatures by means of the solid-state greenhouse effect, as occurs in snow- and ice-covered regions on Earth. The conditions for this possible temporary melting are quantitatively described for bolometric albedo values A = 0.8 and A = 0.2, and with model parameters typical for the thermo-physical conditions at snow/ice sites on the surface of present Mars. It is demonstrated by numerical modelling that there are several sets of parameters which will lead to development of layers of liquid water just below the top surface of snow- and ice-packs on Mars. This at least partial liquefaction occurs repetitively (e.g. diurnally, seasonally), and can in some cases lead to liquid water persisting through the night-time in the summer season. This liquid water can form in sufficient amounts to be relevant for macroscopic physical (rheology, erosion), for chemical, and eventually also for biological processes. The creation of temporary pockets of sub-surface water by this effect requires pre-existing snow or ice cover, and thus is more likely to take place at high latitudes, since the present deposits of snow/ice can mainly be found there. Possible rheologic and related erosion consequences of the appearance of liquid sub-surface water in martian snow/ice-packs are discussed in view of current observations of recent rheologic processes.     

The effects of the martian regolith on GCM water cycle simulations

This paper describes General Circulation Model (GCM) simulations of the martian water cycle focusing on the effects of an adsorbing regolith. We describe the 10-layer regolith model used in this study which has been adapted from the 1-D model developed by Zent, A.P., Haberle, R.M., Houben, H.C., Jakosky, B.M. [1993. A coupled subsurface-boundary layer model of water on Mars. J. Geophys. Res. 98 (E2), 3319-3337, February]. Even with a 30-min timestep and taking into account the effect of surface water ice, our fully implicit scheme compares well with the results obtained by Zent, A.P., Haberle, R.M., Houben, H.C., Jakosky, B.M. [1993. A coupled subsurface-boundary layer model of water on Mars. J. Geophys. Res. 98 (E2), 3319-3337, February]. This means, however, that the regolith is not able to reproduce the diurnal variations in column water vapour abundance of up to a factor of 2-3 as seen in some observations, with only about 10% of the atmospheric water vapour column exchanging with the subsurface on a daily basis. In 3-D simulations we find that the regolith adsorbs water preferentially in high latitudes. This is especially true in the northern hemisphere, where perennial subsurface water ice builds up poleward of 60° N at depths which are comparable to the Odyssey observations. Much less ice forms in the southern high latitudes, which suggests that the water ice currently present in the martian subsurface is not stable under present conditions and is slowly subliming and being deposited in the northern hemisphere. When initialising the model with an Odyssey-like subsurface water ice distribution the model is capable of forcing the simulated water cycle from an arbitrary state close to the Mars Global Surveyor Thermal Emission Spectrometer observations. Without the actions of the adsorbing regolith the equilibrated water cycle is found to be a factor of 2-4 too wet. The process by which this occurs is by adsorption of water during northern hemisphere summer in northern mid and high latitudes where it remains locked in until northern spring when the seasonal CO₂ ice cap retreats. At this time the water diffuses out of the regolith in response to increased temperature and is returned to the residual water ice cap by eddie transport.     

The initial responses of hot liquid water released under low atmospheric pressures: Experimental insights

Experiments have been performed to simulate the shallow ascent and surface release of water and brines under low atmospheric pressure. Atmospheric pressure was treated as an independent variable and water temperature and vapor pressure were examined as a function of total pressure variation down to low pressures. The physical and thermal responses of water to reducing pressure were monitored with pressure transducers, temperature sensors and visible imaging. Data were obtained for pure water and for solutions with dissolved NaCl or CO₂. The experiments showed the pressure conditions under which the water remained liquid, underwent a rapid phase change to the gas state by boiling, and then solidified because of removal of latent heat. Liquid water is removed from phase equilibrium by decompression. Solid, liquid and gaseous water are present simultaneously, and not at the 611 Pa triple point, because dynamic interactions between the phases maintain unstable temperature gradients. After phase changes stop, the system reverts to equilibrium with its surroundings. Surface and shallow subsurface pressure conditions were simulated for Mars and the icy satellites of the outer Solar System. Freezing by evaporation in the absence of wind on Mars is shown to be unlikely for pure water at pressures greater than c. 670 Pa, and for saline solutions at pressures greater than c. 610 Pa. The physical nature of ice that forms depends on the salt content. Ice formed from saline water at pressures less than c. 610 Pa could be similar to terrestrial sea ice. Ice formed from pure water at pressures less than c. 100 Pa develops a low thermal conductivity and a ‘honeycomb’ structure created by sublimation. This ice could have a density as low as c. 450 kg m⁻³ and a thermal conductivity as low as 1.6 W m⁻¹ K⁻¹, and is highly reflective, more akin to snow than the clear ice from which it grew. The physical properties of ice formed from either pure or saline water at low pressures will act to reduce the surface temperature, and hence rate of sublimation, thereby prolonging the lifespan of any liquid water beneath.     

Effect of altered Arctic sea ice and Greenland ice sheet cover on the climate of the GENESIS general circulation model

Previous studies have examined the effect of reduced Arctic sea ice cover on the circulation of climate models. Generally, the response is restricted to high northern latitudes. Here we examine a variant on those simulations, specifying both reduced Arctic sea ice cover and no Greenland ice sheet. The GENESIS general circulation model is used in these experiments. As in earlier studies, we find the effect limited primarily to the high latitudes of the northern hemisphere, being greater in winter than in summer. New results reported herein involve: (1) in winter reduced Arctic ice cover has a significantly greater effect than reduced Greenland ice cover; (2) reduced ice cover had little effect on location of the winter freezing line over North America and Eurasia; (3) removal of ice caused a 30–50% increase in precipitation in high northern latitudes; however there were no significant effects elsewhere. This result does not support the hypothesis that past changes in Arctic ice cover were responsible for significant changes in area of tropical rainforests; (4) there is a peculiar surface pressure anomaly that extends into the high latitudes of the southern hemisphere. This anomaly may be a spurious artifact of the effect of the removed Greenland ice sheet on the spherical harmonic expansion terms in the model. These sensitivity experiments should serve as a useful frame of reference for future Pliocene simulations with a more complete set of altered boundary conditions.     

Metastability of Liquid Water on Mars

A simple model of local heat transport on Mars demonstrates that transient melting of ice may occur in depressions and gullies nearly anywhere on the planet where thin ice is illuminated by normal-incidence insolation. An experiment has been performed to confirm the model of evaporation rate at low pressure. Reduction of radiative cooling due to gully geometry is shown to be important. Since appropriate meteorological, topographic, and optical conditions may occur on slopes nearly anywhere on the planet, hydrological features such as gullies would likely form only where such ice might accumulate, notably in sheltered locations at high latitudes. It is suggested that cold-trapping of winter condensation could concentrate a sufficient amount of ice to allow seasonal melting in gullies.     

The dependence of surface tidal stress on the internal structure of Europa: The possibility of cracking of the icy shell

This study shows dependence of the surface diurnal tidal stress on the internal structure of Europa. Its purpose is to investigate possibility of cracking of the icy shell. The stress is evaluated under the plausible model of the internal structure constrained by the gravity field data. The possible effective stress at the sub-Jovian point decreases with thickening of the shell, while it does not depend on the core radius and the thickness of the H₂O layer. The range of this value is from 0.095 to 0.161 MPa, which does not exceed the tensile strength of ice. The stress required for the surface cracking would be mainly due to longer period deformations, especially non-synchronous rotation. And/or the actual strength of the ice at the surface would be smaller because of the preexisting cracks than that at a laboratory of the same temperature.     

The Mars water cycle

A model has been developed to test the hypothesis that the observed seasonal and latitudinal distribution of water on Mars is controlled by the sublimation and condensation of surface ice deposits in the arctic and antarctic, and the meridional transport of water vapor. Besides reproducing the observed water vapor distribution, the model correctly reproduces the presence of a large permanent ice cap in the arctic and not in the antarctic. No permanent ice reservoirs are predicted in the temperate or equatorial zones. Wintertime ice deposits in the arctic are shown to be the source of the large water vapor abundances observed in the arctic summertime, and the moderate water vapor abundances in the northern temperate region. Model calculations suggest that a year without dust storms results in very little change in the water vapor distribution. The current water distribution appears to be the equilibrium distribution for present atmospheric conditions.     

Diurnal behaviour of water on mars

We have developed a numerical model of the diurnal transport of water across the martian surface. The atmospheric boundary layer is modelled in terms of local radiative-convective processes. The radiative effects of ice fogs near the surface are included in the model. The diffusion of water in the ground is treated for the cases of adsorption and condensation.The model is applied to the diurnal variation of water vapour in the atmosphere as observed by Barker (1974a,b,; 1975). We can explain the morning rise in the amount of water vapour in terms of the evaporation of ground fogs. The evening decrease is compatible with our model if adsorption dominates in the soil. The average level of vapour concentration requires that the atmosphere above the boundary layer be relatively dry. The ground fogs persist until midmorning and should be observable. Some consequences of these conclusions are discussed.     

Thermal stability of water and hydroxyl on the surface of the Moon from temperature-programmed desorption measurements of lunar analog materials

The adsorption of molecular water onto lunar analog materials was investigated under ultra-high vacuum with the goal to better understand the thermal stability and evolution of water on the lunar surface. Temperature-programmed desorption (TPD) experiments show that lunar-analog basaltic-composition glass is hydrophobic, with water-water interactions dominating over surface chemisorption. This suggests that lunar agglutinates will tend not to adsorb water at temperatures above where water clusters and multilayer ice forms. The basalt JSC-1A lunar mare analog, which is a complex mixture of minerals and glass, adsorbs water above 180 K with an adsorption profile that extends to 400 K, showing evidence for a continuum of water adsorption sites. Bancroft albite adsorbs more water, more strongly, than JSC-1A, with a well-defined desorption peak near 225 K. This suggests that mineral surfaces will adsorb more water than mare or mature (glassy, agglutinate rich) surfaces and may explain the association of water with fresh feldspathic craters at high latitudes. The activation energies for the thermal desorption of water from these materials were determined, and along with values from the literature, used to model the grain-to-grain migration of water within the lunar regolith. These models suggest that a combination of recombinative desorption of hydroxyl along with molecular desorption of water and its subsequent migration within and out of the regolith may explain observed diurnal variations in the distribution of water and hydroxyl on the illuminated Moon.     

A model of temporal variations in ozone density in the Martian atmosphere

Temporal variations of the Martian ozone density profile at high latitudes have been calculated for an entire Martian year, taking into account the seasonal and diurnal variations in temperature, water vapor and solar radiation. A new technique facilitates the long-term model calculations, including diurnal variations. The result is in better agreement with MARINER 9 observations of the time and magnitude of the seasonal maximum than is the result of the previous seasonal model calculated for the diurnally averaged temperature, water vapor and solar radiation. The large scatter of the MARINER 9 data may be partly experimental, but the effect of surface condition, including the water vapor variability and the surface chemistry, may explain some of the dispersion of the observed data. The predicted diurnal variation is substantial except near solstices, and the nighttime total column density is generally larger than the daytime value. The magnitude of the day-and-night difference and the shape of the diurnal variation change markedly with season. The opposite temporal variation is predicted for ozone density between the upper and lower regions. The model predicts the production of a ozone layer at 35–50 km, which is consistent with observations at low latitudes by MARS-5. The observed ozone density may be explained, if the atmospheric temperature is as low as ～ 140 K or if the atmosphere is subsaturated. Effects of the simultaneous existence of an aerosol layer, also observed by MARS-5, are briefly discussed.     

On the Formation of Planetesimals: Radial Contraction of the Dust Layer Interacting with the Protoplanetary Disk Gas

Radial contraction of the dust layer in the midplane of a gas–dust protoplanetary disk that consists of large dust aggregates is modeled. Sizes of aggregates vary from centimeters to meters assuming the monodispersion of the layer. The highly nonlinear continuity equation for the solid phase of the dust layer is solved numerically. The purpose of the study is to identify the conditions under which the solid matter is accumulated in the layer, which contributes to the formation of planetesimals as a result of gravitational instability of the dust phase of the layer. We consider the collective interaction of the layer with the surrounding gas of the protoplanetary disk: shear stresses act on the gas in the dust layer that has a higher orbital velocity than the gas outside the layer, this leads to a loss of angular momentum and a radial drift of the layer. The stress magnitude is determined by the turbulent viscosity, which is represented as the sum of the α-viscosity associated with global turbulence in the disk and the viscosity associated with turbulence that is localized in a thin equatorial region comprising the dust layer and is caused by the Kelvin–Helmholtz instability. The evaporation of water ice and the continuity of the mass flux of the nonvolatile component on the ice line is also taken into account. It is shown that the accumulation of solid matter on either side of the ice line and in other regions of the disk is determined primarily by the ratio of the radii of dust aggregates on either side of the ice line. If after the ice evaporation the sizes (or density) of dust aggregates decrease by an order of magnitude or more, the density of the solid phase of the layer’s matter in the annular zone adjacent to the ice line from the inside increases sharply. If, however, the sizes of the aggregates on the inner side of the ice line are only a few times smaller than behind the ice line, then in the same zone there is a deficit of mass at the place of the modern asteroid belt. We have obtained constraints on the parameters at which the layer compaction is possible: the global turbulence viscosity parameter (α < 10^?5), the initial radial distribution of the surface density of the dust layer, and the distribution of the gas surface density in the disk. Restrictions on the surface density depend on the size of dust aggregates. It is shown that the timescale of radial contraction of a dust layer consisting of meter-sized bodies is two orders of magnitude and that of decimeter ones, an order of magnitude greater than the timescale of the radial drift of individual particles if there is no dust layer.     

Water in the upper martian surface at mid- and low-latitudes: presence, state, and consequences

Sublimation of water ice is more effective than evaporation of sorption water at the same temperature. Therefore, water in the form of ice must, over geologic time-scales, have left the upper martian surface (m-scale) at mid- and low-latitudes, leaving sorption water as a possible physical form of stable subsurface water. Adsorption water is “liquid-like” at these temperatures (in the sense of a 2D-liquid). This property is the reason for the specific importance of physisorbed water under martian conditions. It is shown that unfrozen adsorption water can cause numerous physical, chemical, and possibly also biological processes in the upper martian surface and may be responsible for a number of its properties.     

The influence of van der Waals forces on the state of water in the shallow subsurface of Mars

Microscopic liquid layers of water can evolve via adsorption on grain and mineral surfaces at and in the soil of the surface of Mars. The upper parts of these layers will start to freeze at temperatures clearly below the freezing point of bulk water (freezing point depression). A sandwich structure with layers of ice (top), liquid water (in between) and mineral surface (bottom) can evolve. The properties of the interfacial water (of adsorption water and premelted ice) on grain surfaces are described by a sandwich-model of a layer of liquid-like adsorption water between the adsorbing mineral surface layer and an upper ice layer. It is shown that the thickness or number of mono-layers of the interfacial water (of adsorption water and premelted ice) depends on temperature and atmospheric relative humidity. The derived equations for the sandwich model fit well to a known phenomenological relation between thickness of the liquid layer and relative humidity, and can be a tool to estimate or to determine for appropriate materials Hamaker's constant for van der Waals interactions on grains and in porous media. The curvature of grain surfaces is shown to have no remarkable effects for particles in the μm-range and larger. The application of these equations to thermo-physical conditions on Mars shows that the thickness of frost-layers, which can evolve over several hours on cooling surface parts of Mars, is typically of the order or a few tenths of one millimeter or less. This is in agreement with observations. Furthermore, an equation is derived, which relates the freezing point depression for van der Waals force governed interfacial water to the value of the Hamaker constant, to the latent heat of solidification, to the mass density of water ice, and to the thickness of the liquid-like layer. Again, this equation fits well to a known phenomenological relation between freezing point depression and thickness of the liquid-like layer. The derived equation shows that the lower limiting temperature of the liquid phase can reach about 180 K under martian conditions having an atmospheric water content of around 10 pr μm. An “Equilibrium Moisture Content” (EMC)/“Equilibrium Relative Humidity” (ERH) relation for the water content of martian soil has been derived, which relates, for equilibrium conditions, soil water content and atmospheric relative humidity. This relation indicates that the content of liquid interfacial water in the upper surface of Mars can reach up to 10% by weight and more in course of saturation during night hours, and it can be of about 2% by weight during the dry daytime hours.