The Effect of Ground Ice Redistribution on the Martian Paleo-CO2 Cycle, David et al. 2024

Dataset to produce plots for JGR paper "The Effect of Ground Ice Redistribution on the Martian Paleo-CO2 Cycle"

1. Figure 2: "Ground ice distributions – MONS & modern MSIM" was produced using "icetable_gcm_GCM1182.mat" (ice table for Eq. case), "depth_static_H2OonMinDist_230723.mat" (ice table for MONS case), "InfH2OonCaps_grid_220623.mat" (perennial water ice cap distribution), "depth_data.mat" (depth vector) and "dimension_vars.mat" (different dimension vectors: Ls, Time, Longitude, Latitude, etc.).

2. Figure 3: "Pressure with MONS, modern MSIM ground ice and Viking" was produced using "gcm1116_xpsol125.mat" (pressure with MONS GI, default TI case), "gcm1115_xpsol125.mat" (pressure with MONS GI, high TI case) "gcm1184_xpsol140.mat" (pressure with Eq. GI, default TI case), "gcm1182_xpsol140.mat" (pressure with Eq. GI, high TI case), "VL1_lmd.txt" (VL1 pressure) & "dimension_vars.mat".

3. Figures 4-9: "Seasonal caps & GI areas at different obliquities", "Full grid maps of CO2 mass & time at obliquity 25", "Full grid maps of CO2 mass & time at obliquity 45", "CO2 ice and surface temperature at a low TI regions", " Surface pressure at different obliquities" and "Surface temperature at obliquities 25 and 45" were produced using "dimension_vars.mat", "co2_GCM1117-1121.mat" (CO2 from GCM simulations with Eq. GI, e = 0, obliquity = 25, 30, 35, 40, 45 respectively), "icetable_gcm_GCM1117-1121.mat" (Eq. GI ice table), "co2_GCM1050-1054.mat" (CO2 from GCM simulations with MONS GI, e = 0, obliquity = 25, 30, 35, 40, 45 respectively) and "depth_static_H2OonMinDist_230723.mat".

4. Figures 10-12: "Contours of northern seasonal cap as function of orbital parameters", "Contours of southern seasonal cap as function of orbital parameters" and "Interpolation of surface pressure for past 20 MA" were produced using "CO2_Orbit_e0-012_eq/st.mat" (CO2 mass and pressure for e = 0, 0.04, 0.08, 0.12; Eq./MONS GI cases, respectively).

5. Figure SM 1: "CO2 mass as a function of MSIM-GCM iteration" was produced using "run1050-1054.nc", "run1117-1021.nc", "run1185-1189.nc" (MONS/Eq. 1/Eq. 2 full GCM results for e = 0, obliquity = 25, 30, 35, 40, 45, respectively).