Time-Varying Causal Treatment for Quantifying the Causal Effect of Short-Term Variations on Arctic Sea Ice Dynamics
Akila Sampath, Vandana Janeja, Jianwu Wang
TL;DR
The paper addresses the challenge of estimating causal effects in complex Arctic ocean–ice systems, where short-term SSH variations influence sea ice thickness through physical processes and unobserved confounding. It introduces KGCM-VAE, a knowledge-guided variational autoencoder with a causal-adjacency decoder and MMD-based latent deconfounding to learn balanced, physically-consistent representations for time-varying treatments. The approach leverages hydrostatic and geostrophic physics to guide treatment generation and enforces a learned causal structure, yielding superior PEHE performance over state-of-the-art baselines and demonstrating robustness across treatment lags. The work advances data-driven causal understanding in polar climate dynamics, enabling more reliable counterfactual analysis and informing climate model refinement and prediction under real-world observational constraints.
Abstract
Quantifying the causal relationship between ice melt and freshwater distribution is critical, as these complex interactions manifest as regional fluctuations in sea surface height (SSH). Leveraging SSH as a proxy for sea ice dynamics enables improved understanding of the feedback mechanisms driving polar climate change and global sea-level rise. However, conventional deep learning models often struggle with reliable treatment effect estimation in spatiotemporal settings due to unobserved confounders and the absence of physical constraints. To address these challenges, we propose the Knowledge-Guided Causal Model Variational Autoencoder (KGCM-VAE) to quantify causal mechanisms between sea ice thickness and SSH. The proposed framework integrates a velocity modulation scheme in which smoothed velocity signals are dynamically amplified via a sigmoid function governed by SSH transitions to generate physically grounded causal treatments. In addition, the model incorporates Maximum Mean Discrepancy (MMD) to balance treated and control covariate distributions in the latent space, along with a causal adjacency-constrained decoder to ensure alignment with established physical structures. Experimental results on both synthetic and real-world Arctic datasets demonstrate that KGCM-VAE achieves superior PEHE compared to state-of-the-art benchmarks. Ablation studies further confirm the effectiveness of the approach, showing that the joint application of MMD and causal adjacency constraints yields a 1.88\% reduction in estimation error.
