Modeling CO2 saturation distribution in eolian systems