Effects of regolith properties on the Martian subsurface water distribution using a global climate model

In today’s extremely dry Mars, water vapor “adsorption” on regolith grains is thought to play crucial roles in subsurface water retention and water vapor exchange with the atmosphere (Fanale & Cannon, 1971; Zent et al., 1993, 1995, 2001; Böttger et al., 2005; Savijärvi et al., 2016, 2020). Global models that explicitly account for water diffusion in the shallow subsurface and calculate subsurface water distribution have assumed globally uniform regolith properties to simplify assumptions (Böttger et al., 2005; Schorghofer & Aharonson, 2005; Steele et al., 2017). However, Pommerol et al. (2009) examined the adsorption efficiency of six samples similar to the Martian regolith and found that the samples with smaller grain sizes store more adsorbed water due to their larger specific surface areas. Therefore, we have newly implemented a regolith scheme in a Mars Global Climate Model (MGCM), considering regolith properties like grain size, porosity, and the specific surface area. The grain size distribution was obtained from the empirical equation as a function of thermal conductivity (Presley & Christensen, 1997). The distributions of porosity and the specific surface area are also determined, referring to the laboratory experiments of Sizemore & Mellon (2008). Our results clarify that regolith grains with large specific surface areas in the northern low and mid-latitudes and the southern high latitudes, which have high adsorption coefficients, affect water storage. Subsurface water in the northern low and mid-latitudes exists up to 0.5–1wt% as adsorbed water. Regolith with high adsorption properties makes the depth of subsurface ice shallower in the southern high latitudes. Pore ice accumulates in regions poleward of 50°N and 50°S and the west of Elysium Mons and Olympus Mons, which is consistent with previous simulations. Also, with a homogeneous specific surface area, seasonal increases in pore ice were calculated at a depth of about 60 cm in mid-latitudes with low thermal inertia and high atmospheric water vapor content, but with the specific surface area map, the seasonal increases were not demonstrated. This study suggests that adsorption properties influence subsurface water dynamics, emphasizing the importance of considering inhomogeneous regolith properties in models of subsurface water distributions and the atmospheric water cycle including the regolith.