Neural fields in the geosciences, 2025–2026
Report 18 · Compression, reconstruction and physics constraints across ocean, atmosphere, land and ice
Abstract¶
Neural fields now clearly win for gridless reconstruction from sparse, irregular observations (altimeter tracks, stations, CryoSat-2 point clouds, seismic and gravity inversion), and as differentiable state parameterizations for variational assimilation, where the network’s spectral bias stands in for a background covariance. However, for the archival compression of reanalysis they still lose to tuned error-bounded codecs, and a 2026 ERA5 head-to-head found that SZ3.1 leads at every error bound. Atmospheric-composition fields (methane, NO₂, XCO₂) have essentially no neural-field work. What GeoStack should do about it is on the neural-fields design page.
Key findings¶
Spherical handling has converged on a few recipes. Spherical-harmonic embeddings into a SIREN Rußwurm et al., 2024 are the baseline. However, work in 2025 showed that they fail on small or fragmented regions such as islands and coasts, and it proposed spherical wavelet encodings instead Cai & Balestriero, 2025. HEALPix hybrid grids and learnable spherical Voronoi partitions are the alternatives.
Compression claims need careful reading. Most INR compression papers compare against JPEG2000 or against vanilla INRs, and not against SZ3/ZFP at matched error bounds. When the comparison is done properly, SZ3 is hard to beat on smooth gridded fields.
Gridless reconstruction is the killer application. Neural fields consume point clouds without regridding, and they give analytic derivatives for geostrophy, elevation-change rates and PDE residuals.
Spectral bias is being reframed as a feature. In 4D-Var and in geophysical inversion, the MLP’s preference for low frequencies acts as an implicit regularizer, which replaces the hand-built background covariances or smoothness penalties.
Amortization across samples is thin in practice. Functa, CORAL and DINo-style modulations are mature on ML benchmarks, whereas the real geoscience deployments are few (PINROD for the ocean, and the multi-altitude climate INRs).
JAX code exists but is scattered. There is BayesNF (Google), DIFFICE-jax (ice shelves), PINNICLE through DeepXDE’s JAX backend, and DeepMind’s functa. However, no sphere-aware, geoscience-specific neural-field library exists in JAX.
Compression and storage reduction¶
Neural fields do not dominate archival compression in 2026. Their advantages are random-access continuous queries, analytic gradients, context conditioning and compression as an embedding.
| Work | Group, venue | Field represents | Architecture and encoding | Key result | Code |
|---|---|---|---|---|---|
| Compressing weather and climate data into neural networks (foundation) Huang & Hoefler, 2023 | Huang & Hoefler (ETH), ICLR 2023 | (lat, lon, level, time) → ERA5 variable | Coordinate MLP, Fourier features, quantized weights | 300× to 3,000×; beats SZ3 on weighted RMSE/MAE; 790× WeatherBench loader with < 2% RMSE rise | PyTorch |
| HiHa Xu & others, 2025 | Xu et al. (Jilin), Computers & Geosciences 2025 | Spherical coords + level → atmospheric variables | 3-D FFT band split; multi-scale INR for low, octree-local SIRENs for mid, sparse storage for high frequencies | Targets 204 variables on 37 levels; explicit fix for spectral bias | not confirmed |
| CA-INR Mostajeran et al., 2025 | Mostajeran et al. (Utah), Scientific Reports 2025 | (space, time, context) → E3SM surface temperature | MLP with topography and climatology as context inputs | Ratio 140 on benchmarks; 500 years of CESM at 200; context adds up to 11.46 dB PSNR | not confirmed |
| Multi-altitude climate INRs Abdul Qayyum et al., 2024 | Abdul Qayyum et al. (BNL / TAMU), MLSP 2024; NeurIPS 2025 CCAI | Low-res latent + coordinate → multi-altitude wind | Encoders + EDSR features + LIIF coordinate decoder | Compression, arbitrary-scale reconstruction and modality fusion in one model | not confirmed |
| Hybrid Neural Representations for Spherical Data Hybrid Neural Representations for Spherical Data, 2024 | ICML 2024, arXiv 2402.05965 | Sphere → ERA5 / CMB | HEALPix grid features + MLP decoder | Strong sphere-native baseline | — |
| Climplicit Dollinger et al., 2025 | Dollinger et al. (EcoVision UZH), arXiv 2504.05089 | (lon, lat, month) → CHELSA climatology embedding | Residual SIREN blocks | 16 MB model replaces 99.97% of CHELSA; ×3,500 less disk | GitHub |
| ImplicitTerrainV2 ImplicitTerrainV2, 2026 | arXiv 2026 | (x, y) → elevation | Wavelet-guided adaptive SIREN + quantization + entropy coding | Beats ZSTD, ZFP, LERC at 0.5–2 bpp; SZ3 still wins (80.45 dB at 1.96 bpp vs 65.97 dB at 1.23 bpp) | — |
| Can DNNs improve compression of very large scientific data? Can Deep Neural Networks Improve Compression of Very Large Scientific Data?, 2026 | arXiv 2026 | ERA5 via learned predictors | CRA5, GraphCast, Aurora as predictors vs SZ3.1, same backend | No learned predictor improves dataset-level compression; up to 9.6× per-variable gains | — |
For context outside INRs, DeepMind’s Neural Compression of Atmospheric States Mirowski & others, 2024 uses HEALPix with hyperprior and vector-quantized autoencoders to reach more than 1,000× at about one second per global state, which is the bar that INRs must beat. The 2025 review Lossy Neural Compression for Geospatial Analytics Gomes & others, 2025 calls the compression efforts for Earth-system models limited, and it says that the shortcomings of INRs must be fixed first.
So one should benchmark against SZ3/ZFP at a matched maximum error, and not just on PSNR.
Interpolation, gap-filling and continuous reconstruction¶
This is where neural fields are strongest, as point clouds and swaths go in and a continuous, differentiable field comes out.
Ocean¶
Neural fields for fast and scalable interpolation of geophysical ocean variables: Johnson et al. (MEOM/IGE), NeurIPS ML4PS 2022 (foundation) Johnson & others, 2022. It maps (lon, lat, t) to SSH from along-track altimetry with SIREN and Fourier-feature MLPs, and it matches optimal interpolation while scaling better. OceanBench SSH Johnson et al., 2023 is the benchmark harness.
Continuous field reconstruction from sparse observations with INRs: Luo et al., ICLR 2024 Luo & others, 2024. They use separable spatial and temporal basis functions, and they beat prior INRs on climate output and on 0.01° satellite SST.
SWOT noise. Zhang & Callies (JGR Oceans 2025) find that the KaRIn signal below 30 to 50 km correlates with the significant wave height, so aliased waves are the likely source. So SSH neural fields must now model anisotropic, wave-correlated noise, and not just gaps.
Atmosphere and near-surface weather¶
SpLIIF (del Valle Inclan Redondo et al., Recursive; arXiv 2510.20228) del Valle Inclan Redondo & others, 2025: it maps sparse JMA stations and topography to continuous sub-km temperature and wind, and it is up to 50% better than CorrDiff for temperature and 10 to 20% better for wind.
WindINR WindINR: Latent-State Implicit Neural Representation for Fast Local Wind Query and Correction in Complex Terrain, 2026: a latent-state INR for local wind queries and for observation-driven correction in complex terrain.
BayesNF (Saad et al., Google; Nature Communications 2024) Saad et al., 2024: it maps (space, time) to a full predictive distribution through hierarchical priors over a neural field with MAP/VI ensembles, and it beats sparse variational GPs, gradient boosting and GLMMs on tens to hundreds of thousands of points. It runs in JAX, on GPU/TPU.
Land and solid Earth¶
ImplicitTerrainV2 (above): DEMs as terrain with analytic gradients.
IntraSeismic and Bayesian seismic inversion with INRs (Romero, Heidrich, Ravasi, KAUST; JGR-MLC 2024, GJI 2025) Romero et al., 2025: they map (x, y, z) to elastic properties with a multigrid hash encoding and a tiny MLP inside a seismic forward operator, and they use variational inference for the uncertainty.
Implicit gravity inversion (Computers & Geosciences 2025), Xu & Heagy on why a neural-field parameterization helps DC-resistivity and seismic inversion, and a 1-D magnetotelluric INR (2026).
Ecological field reconstruction (2026): sinusoidal and Fourier INRs keep the fine structure, whereas ReLU models over-smooth.
For soil moisture and vegetation, we found no neural-field paper from 2025 or 2026.
Ice and cryosphere¶
INR for ice-sheet surface elevation (Naylor et al., ESA Φ-lab / DTU Space; SSRN 2025) Naylor et al., 2025: it maps (x, y, t) to the elevation of Petermann Glacier from CryoSat-2 points which cover 5% of the domain, and it is 83.66% better than its baseline over more than 20,000 km² and 12 years, although it is weak at the low-density rim. The code is not in JAX, and the paper is not yet peer reviewed.
Sea-ice thickness and Greenland bed mapping (DeepTopoNet, physics-guided residuals) still use gridded CNNs, and we found no INR paper on sea ice or firn.
Location encoders¶
Spherical harmonics + SIREN (Rußwurm et al., ICLR 2024) Rußwurm et al., 2024: an analytic SH basis up to degree L = 100, with no polar artefacts.
No Location Left Behind (Cai & Balestriero, ICLR 2025) Cai & Balestriero, 2025: the FAIR-Earth benchmark shows that SH encoders underperform on islands and coasts, and the authors propose spherical wavelets.
In 2026 came a spherical Voronoi encoder, SLED distilled encodings, and an Earth Embeddings survey which covers SatCLIP, RANGE and the hash-table LIANet.
Physics-constrained neural fields¶
The live frontier is in three places: 1) neural-field 4D-Var, 2) PINN inversion for the ice and the ocean, and 3) latent dynamics over INR codes.
| Work | Venue | Neural-field role | Physics | Result | Code |
|---|---|---|---|---|---|
| Neural field reparameterization for 4DVAR | Jaemin Oh, arXiv 2025 (v2 2026) Oh, 2025 | Full state as one field | 4D-Var cost + physics losses, parallel-in-time | Better initial conditions than classical 4D-Var on 2-D Kolmogorov and 3-D Taylor–Green; spectral bias replaces an explicit background covariance ; separable PINNs give large speedups; no ground-truth training data | — |
| StrAss-PINN: deep learning in the abyss | arXiv 2025 StrAss-PINN: Deep learning in the abyss, 2025 | SIREN-style field per ocean layer | 3-layer quasi-geostrophic | Deep flow from SWOT-like SSH, ARGO-like interior, sparse deep-ARGO | — |
| PINROD | Zhao, Luo et al. (BNL/LANL/ANL), ICLR 2025 CCAI Zhao et al., 2025 | DINo-style INR decoder + latent codes | Parameterized latent neural ODE | Resolution-independent baroclinic ocean surrogate | — |
| Deep learning the flow law of Antarctic ice shelves | Wang et al., Science 2025 Wang et al., 2025 | MLPs for velocity, thickness, viscosity | Shallow-shelf approximation | Grain-size-sensitive rheology in compression zones; anisotropic viscosity at suture zones | DIFFICE-jax (JAX) |
| PINNICLE v1.0 | Cheng, Krishna, Morlighem, GMD 2025 Cheng et al., 2025 | (x, y) → ice fields | SSA, MOLHO, mass conservation; forward and inverse | General ice-sheet PINN library | GitHub |
| Subglacial topography from two conservation laws | Krishna, Cheng, Morlighem, The Cryosphere 2026 Krishna et al., 2026 | PINN for thickness + basal shear | Mass conservation + SSA | Single-law PINNs insufficient in sparse regions; finds bed features absent from BedMachine | PINNICLE |
| PINN surface currents from drifters | Frontiers in Marine Science 2025 | (x, y, t) → velocity | Flow equations as soft constraints | Physics term essential for Eulerian fields from Lagrangian drifters | — |
The latent-dynamics foundations are DINo Yin & others, 2023, which pairs an INR decoder with a latent ODE, CORAL Serrano & others, 2023, which uses shift-modulated SIRENs with meta-learned codes for operator learning on irregular grids, and AROMA Serrano & others, 2024, which uses local neural fields with transformer latent dynamics and beats both on harder systems.
Cross-cutting trends and open problems¶
| Problem | Where the field is | Open gap |
|---|---|---|
| Sphere and geodesic coordinates | SH(λ, φ) → SIREN is standard; spherical wavelets, HEALPix hybrid grids, spherical Voronoi emerging | SH cost grows as L² and underserves small landmasses; no native multiresolution hash grid on the sphere |
| Multiscale and spectral bias | Explicit band split (HiHa), wavelet-guided local nets (ImplicitTerrainV2), spectral bias as regularizer (neural 4D-Var, inversion) | Fidelity tasks fight spectral bias; ill-posed inverse problems exploit it |
| Scalability | Separable nets (SPINN), hash grids (IntraSeismic), octree-local models (HiHa), minibatch VI (BayesNF) | Global, multi-decade, multi-variable fields are mostly out of reach for one INR |
| Uncertainty | BayesNF (priors + ensembles), Bayesian IntraSeismic (VI) | GP links (NTK, OI equivalence) understood but rarely used in deployed INRs |
| Generalization across samples | Functa, COIN++, CORAL modulations; applied: PINROD, multi-altitude LIIF, Climplicit | Amortized encoders trained on real observation data |
| Foundation models | Aurora and GraphCast decode to fixed grids; MLP heads on frozen Aurora latents are patch decoders | Coordinate-queried continuous decoders (LIIF/INR heads) on foundation latents |
| Atmospheric composition | Methane, NO₂, XCO₂ gap-filling uses CNN, LSTM, U-Net, trees, diffusion | No neural-field work found; averaging-kernel-aware observation operators untouched |
Software landscape¶
| Library | Framework | What it gives |
|---|---|---|
| google/bayesnf | JAX | Bayesian spatiotemporal neural fields, MAP/VI ensembles |
| DIFFICE-jax | JAX | Ice-shelf PINNs (SSA), viscosity inversion |
| PINNICLE | DeepXDE (TF / PyTorch / JAX) | Ice-sheet forward and inverse PINNs |
| DeepMind functa | JAX | Meta-learned latent modulations |
| MarcCoru | PyTorch | Spherical harmonics + SIREN location encoders |
| TorchSpatial | PyTorch | Location-encoder zoo |
| Climplicit | PyTorch | Pretrained climate implicit embeddings |
| CORAL, AROMA | PyTorch | INR operator learning and latent dynamics |
| IceSheetReconstruction | PyTorch + Nextflow | CryoSat-2 elevation INR pipeline |
The benchmarks are OceanBench (SSH), FAIR-Earth (stratified fairness of Earth INRs), and SDRBench-style comparisons against SZ3, ZFP and MGARD.
Caveats¶
Several key items are preprints or workshop papers (Naylor et al., PINROD, Oh 2025, ImplicitTerrainV2, the 2026 ERA5 comparison), so one should check their numbers against the final versions.
Compression ratios are not comparable across papers, because the datasets, the error metrics (PSNR vs maximum error) and whether the weights are quantized and entropy-coded all differ.
Relative gains (83.66%, “up to 50%”) depend on baselines which are not always strong.
The absence findings (no INR work on trace gases, soil moisture, sea ice or firn) come from a limited search.
Works whose title or authors we could not confirm from the sources opened for this review (Zhang & Callies 2025 on SWOT noise; Xu & Heagy 2025; the 2026 magnetotelluric and ecological INR studies; the 2025 drifter PINN; IntraSeismic, Implicit Gravity Inversion, DeepTopoNet) are linked inline rather than entered in
references.bib.
- Rußwurm, M., Klemmer, K., Rolf, E., Zbinden, R., & Tuia, D. (2024). Geographic Location Encoding with Spherical Harmonics and Sinusoidal Representation Networks. International Conference on Learning Representations (ICLR). https://github.com/MarcCoru/locationencoder
- Cai, & Balestriero, R. (2025). No Location Left Behind. International Conference on Learning Representations (ICLR). https://mlanthology.org/iclr/2025/cai2025iclr-location/
- Huang, L., & Hoefler, T. (2023). Compressing Multidimensional Weather and Climate Data into Neural Networks. International Conference on Learning Representations (ICLR). https://arxiv.org/abs/2210.12538
- Xu, & others. (2025). HiHa: Hierarchical Harmonic Decomposition to Implicit Neural Compression for Atmospheric Data. Computers & Geosciences. https://arxiv.org/abs/2411.06155
- Mostajeran, Pawar, Villarreal, & Faroughi. (2025). Context-aware implicit neural representations to compress Earth system model data. Scientific Reports, 15, 25932. https://www.nature.com/articles/s41598-025-11092-w
- Abdul Qayyum, Luo, Urban, Qian, & Yoon. (2024). Implicit Neural Representations for Simultaneous Reduction and Continuous Reconstruction of Multi-Altitude Climate Data. IEEE International Workshop on Machine Learning for Signal Processing (MLSP). https://arxiv.org/abs/2409.17367
- Hybrid Neural Representations for Spherical Data. (2024). https://arxiv.org/abs/2402.05965
- Dollinger, Robert, Plekhanova, Drees, & Wegner. (2025). Climplicit: Climatic Implicit Embeddings. https://arxiv.org/abs/2504.05089
- ImplicitTerrainV2. (2026). https://arxiv.org/html/2605.22556
- Can Deep Neural Networks Improve Compression of Very Large Scientific Data? (2026). https://arxiv.org/html/2606.14353
- Mirowski, P., & others. (2024). Neural Compression of Atmospheric States. https://arxiv.org/abs/2407.11666
- Gomes, & others. (2025). Lossy Neural Compression for Geospatial Analytics: A Review. https://arxiv.org/abs/2503.01505
- Johnson, J. E., & others. (2022). Neural Fields for Fast and Scalable Interpolation of Geophysical Ocean Variables. NeurIPS Workshop on Machine Learning and the Physical Sciences. https://arxiv.org/abs/2211.10444
- Johnson, J. E., Febvre, Q., Gorbunova, A., Metref, S., Ballarotta, M., Le Sommer, J., & Fablet, R. (2023). OceanBench: The Sea Surface Height Edition. Advances in Neural Information Processing Systems 36, Datasets and Benchmarks Track. https://proceedings.nips.cc/paper_files/paper/2023/file/f6ccbf94fa57c2ae372ece91b537574d-Paper-Datasets_and_Benchmarks.pdf
- Luo, & others. (2024). Continuous Field Reconstruction from Sparse Observations with Implicit Neural Networks. International Conference on Learning Representations (ICLR). https://arxiv.org/abs/2401.11611