WaveDINO: Learning-Based Atmospheric Correction of Unwrapped InSAR Interferograms

1COMET, School of Computer Science, University of Bristol, UK,

2COMET, School of Earth Science, University of Bristol, UK

Abstract

Interferometric Synthetic Aperture Radar (InSAR) enables effective monitoring of volcanic deformation; however, the observed signals are often corrupted by atmospheric phase delays, seasonal surface changes, and decorrelation effects. Existing atmospheric correction methods, such as numerical weather model–based methods, can reduce these effects but do not consistently remove atmospheric artefacts and may introduce residual biases.

To address these limitations, we propose a novel learning-based method for denoising unwrapped InSAR interferograms, using a hybrid training strategy that combines physically motivated synthetic deformation with real atmospheric noise. Specifically, we introduce WaveDINO, a wavelet-based multi-scale denoising framework conditioned on frozen DINOv3 foundation-model features and terrain information. Training uses synthetic magma-source deformation superimposed on short-term interferograms to expose the network to realistic atmospheric statistics while retaining known ground truth.

Performance is evaluated on both controlled synthetic data and long-term real interferograms from Laguna del Maule (Chile) and Campi Flegrei (Italy), with independent GNSS measurements used for validation. WaveDINO consistently outperforms competing models, improving agreement with GNSS measurements and reducing mean GNSS misfit by approximately 3% and 19% at the two sites respectively, while surpassing weather-model-based corrections.

At a Glance

2
Volcanic test sites (Chile & Italy)
~38k
Synthetic deformation maps for training
~6.8k
Real long-term test interferograms
24
GNSS stations used for validation

Sentinel-1 IW interferograms are processed automatically by the LiCSAR pipeline and referenced via LICSBAS for the COMET Volcano Portal. For each volcano we analyse a 50×50 km region of interest at 50 m pixel spacing (500×500 pixels), spanning 622 epochs at Laguna del Maule (Oct 2014 – Aug 2024) and 1008 epochs at Campi Flegrei (Oct 2014 – Sep 2025). Independent ground truth is provided by 5 GNSS stations at Laguna del Maule and 19 GNSS stations at Campi Flegrei.

Method

WaveDINO is a wavelet-based multi-scale denoiser that explicitly targets the spatially coherent, non-stationary statistics of tropospheric phase delay. The input unwrapped phase is decomposed into multi-resolution sub-bands so that broad atmospheric components can be modelled primarily in the low-frequency pathway, while fine deformation-related structure is preserved through high-frequency skip connections. Encoder channels follow 32 → 64 → 128, with the decoder mirrored.

DINOv3 Cross-Attention (DCA)
Fuses frozen DINOv3 ConvNeXt-Base features with the model's low-frequency representation via cross-attention with a low-frequency attention bias, FiLM modulation and a learnable gated residual. Provides global semantic context that disambiguates large-scale atmosphere from real deformation.
Terrain Modulation Block (TMB)
DEM-based feature-wise modulation (FiLM) using elevation and its horizontal/vertical gradients, with a learnable gate that adapts terrain influence in coastal and low-relief regions where turbulence dominates.
Phase Ramp Head (PRH)
A lightweight head that predicts a bounded first-order planar trend (a·x + b·y + c) from globally pooled decoder features, capturing dominant scene-wide bias so the convolutional backbone can focus on residual structure.
Hybrid Training Strategy
Synthetic magma-source deformation (1–4 Mogi sources, log₁₀ΔV ∈ [5, 7.9], depth 1.5–8 km) is superimposed on short-baseline (≤18 day) real interferograms, exposing the network to realistic atmosphere while keeping a known target.

Study Areas

Laguna del Maule, Chile

Silicic volcanic system at 70.492°W, 36.058°S in the southern Andes. Rugged high-altitude topography with seasonal snow cover — conditions that promote strong stratified tropospheric delays and coherence loss. Persistent uplift exceeding several centimetres per year suggests significant shallow magma accumulation. Tracks 18A and 83D, 622 epochs (Oct 2014 – Aug 2024).

Campi Flegrei, Italy

Densely monitored caldera west of Naples (14.139°E, 40.827°N) in a low-relief, coastal and urbanised environment. Persistent hydrothermal and magmatic unrest (bradyseism) produces spatially heterogeneous atmospheric artefacts. Tracks 22D, 44A and 124D, 1008 epochs (Oct 2014 – Sep 2025).

(Left) Laguna del Maule, Chile. (Right) Campi Flegrei, Italy.

Results on Synthetic Data

We hold out synthetic deformation maps and superimpose them on real short-baseline interferograms to obtain a controlled test set with known ground truth. WaveDINO achieves the best perceptual (LPIPS), structural (SSIM) and pixel-wise (MSE, L3, ABS) scores at both volcanoes, indicating strong atmospheric attenuation without oversmoothing the deformation signal.

Table I. Image-quality metrics on synthetic data computed over the full image.

Bold and underline denote the best and second-best results in each column.

Method Laguna del Maule Campi Flegrei
MSE ↓LPIPS ↓SSIM ↑L3 ↓ABS ↓ MSE ↓LPIPS ↓SSIM ↑L3 ↓ABS ↓
UNet54.890.250.840.01094.9935.300.290.730.00883.37
Diffusion (DDPM)93.640.470.710.01406.5866.520.510.620.01174.48
SwinIR69.110.340.790.01255.7440.530.350.720.00973.83
Restormer33.770.270.850.00903.9530.420.460.640.00813.34
WaveMamba44.130.160.880.00994.5546.290.250.760.00943.70
WaveDINO (Ours)30.120.120.910.00803.8126.440.090.920.00632.87

MSE in mm2; L3 and ABS in mm. SSIM is unitless.

Table II (top). MSE on synthetic data, binned by ground-truth deformation magnitude.

Lower is better. Bold = best, underline = second-best per column.

Method Laguna del Maule Campi Flegrei
<1010–2020–30>30Mean <1010–2020–30>30Mean
UNet27.7664.3986.76120.0274.7329.9258.5773.39112.7968.67
Diffusion (DDPM)39.9489.13138.63247.14128.7163.22103.25132.91204.29125.92
SwinIR28.3171.76108.32169.3894.4436.0261.2979.05105.5670.48
Restormer17.8135.6646.4866.8241.6922.9851.2973.0599.4161.68
WaveMamba22.2749.0364.4786.3955.5420.8841.5156.3079.9149.65
WaveDINO (Ours)19.1035.5541.8054.1237.6415.6226.7737.6155.2933.82

All values in mm. Bins refer to ground-truth deformation magnitude in mm.

Validation on Real Interferograms (GNSS)

For real-data evaluation, we project GNSS time series into the satellite line-of-sight per track and compute MSE against the corrected InSAR phase, binned by deformation magnitude. Up to 5,000 long-term interferograms per volcano are used, including intervals with weak deformation, to avoid subjective selection bias.

Learning-based denoisers substantially improve agreement with independent GNSS over both raw unwrapped and GACOS-corrected interferograms. WaveDINO yields the lowest mean error at both sites, with especially large gains at Campi Flegrei (mean MSE 173 mm vs 215 mm for the best baseline).

Table II (bottom). MSE on real interferograms vs GNSS, binned by deformation magnitude.

Lower is better. Bold = best, underline = second-best per column. Unwrapped Interf. = raw input; GACOS = weather-model correction.

Method Laguna del Maule Campi Flegrei
<1010–2020–30>30Mean <1010–2020–30>30Mean
Unwrapped Interf.2703844041484636167125143448221
GACOS39352383936671355191143165464241
UNet471883593188946136189212640294
Diffusion (DDPM)4418242331699551102063381495537
SwinIR3720541933601005113162158425215
Restormer38189397331298490147182437214
WaveMamba431933652812853104178175487236
WaveDINO (Ours)4819435527208294298157417173

All values in mm. Bins refer to GNSS-derived deformation magnitude in mm.

Ablation Study

Removing any single component degrades performance on Laguna del Maule synthetic data. DEM conditioning is the single most influential component (V2: MSE 37.57 vs 30.12 for the full model), followed by the phase-ramp head (V3), water-feature input (V1) and DINOv3 conditioning (V4) — indicating complementary benefits from terrain, water-aware inputs, ramp modelling and global semantic context.

Table III. Ablation study on Laguna del Maule synthetic data.

Variant Water DEM Ramp DINO MSE ↓
V133.51
V237.57
V334.08
V432.39
Full (Ours)30.12

Acknowledgements

This work is supported by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (MAST; Grant No. 101003173) and the NERC Centre for the Observation and Modelling of Earthquakes, Volcanoes and Tectonics (COMET). The authors thank L. Córdova and the Observatorio Volcanológico de los Andes del Sur (OVDAS), SERNAGEOMIN, for providing GNSS time series from Laguna del Maule. Campi Flegrei GNSS data are provided by Istituto Nazionale di Geofisica e Vulcanologia — Osservatorio Vesuviano (INGV-OV).

BibTeX

@article{popescu2026wavedino,
  title   = {WaveDINO: Learning-Based Atmospheric Correction of Unwrapped InSAR Interferograms Validated by GNSS: Results at Laguna del Maule and Campi Flegrei Volcanoes},
  author  = {Popescu, Robert and Biggs, Juliet and Zhu, Tianyuan and Anantrasirichai, Nantheera},
  journal = {arXiv preprint},
  year    = {2026}
}