Skip to content

This repo is work done for fun on trying to apply the parametric Kalman Filter (Pannekoucke et al. 2016) on the sphere for a transport equation.

Notifications You must be signed in to change notification settings

antoineperrot/spherical-advection

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

23 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

spherical-advection

This repo is work done for fun on trying to apply the parametric Kalman Filter (Pannekoucke et al. 2016) on the sphere for a transport equation.

Notebook 1: Semi-Lagrangian method for transport on $\mathbb{R}^2$

This notebook presents a numerical experiment where a semi-lagrangian method is employed to simulate the transport by the wind of an initial state.

observation: the wind field (black arrows) transports correctly the initial state, and forms some spirals around critical points.

Notebook 2: Semi-Lagrangian method for transport on $\mathbb{S}^2$

This time, the semi-lagrangian method is applied for a spherical domain. Unfortunately, it is a non-conservative method, meaning that additional treatment must be applied to conserve the total mass along one timestep.

observation: the initial state (a high concentration on the north pole) has been distorded by the flow.

Notebook 3: Generating isotropic gaussian random fields on the sphere using $\mathbf{B}^{1/2}$

In order to prepare future data assimilation experiments (with an EnKF), it is necessary to know how to generate isotropic gaussian random field (GRF). In this notebook, we do so by building the covariance matrix of an isotropic GRF (on a lat-lon grid), then computing its Cholesky decomposition to apply it to a centered normal vector. This method is numerically very expensive and therefore limited to very low-resolution grid.

Then, a sanity check indicates that the random fields generated are not perfectly isotropic (longer meridional length-scale)

Notebook 4: Generating isotropic GRF on the sphere with spherical harmonics

In this notebook we use a lemma that relies on spherical harmonics to generate isotropic GRF at high-resolution.

Then, a sanity chech confirms the isotropic character of the GRF.

Notebook 5: Plotting VLATcov modeled correlation function on the sphere

In this notebook we use the VLATcov covariance model to represent homogeneous correlation functions on the sphere.

Notebook 6: Enkf and PKF for Advection on the sphere

In this last notebook we compare the error statistics forecasted by an EnKF and the PKF for the transport on the sphere (to be continued...).

About

This repo is work done for fun on trying to apply the parametric Kalman Filter (Pannekoucke et al. 2016) on the sphere for a transport equation.

Topics

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published