I am currently working on an alternative to the classical model-driven data assimilation. The proposed technique is completely data-driven and is based on a catalog of historical data. The basic idea is to combine analog forecasting and classical ensemble data assimilation. Presentation of the method is given in this chapter of the book Machine Learning and Data Mining Approaches to Climate Science. This book grew out of the workshop on Climate Informatics held in Boulder, Colorado. The Python code for the Analog Data Assimilation (AnDA) is available at: https://github.com/ptandeo/AnDA.
Uncertainty quantification in data assimilation:
The Ensemble Kalman Filter is a widely used data assimilation method. This algorithm is a stochastic approach based on Monte Carlo simulations and very useful to account for nonlinear dynamics and high dimensional problems. In this article of the Quarterly Journal of the Royal Meteorological Society, I proposed a statistical algorithm (Expectation-Maximization) to estimate the error covariance matrices in such nonlinear state-space models used in EnKF. I also presented an application for the estimation of physical parameters for subgrid-scale processes. The Python code for the Covariance Estimation in Data Assimilation (CEDA) is available at: https://github.com/ptandeo/CEDA.
Another part of my research is dedicated to the synergy in ocean remote sensing, studying the statistical relationships between SST and SSH mesoscale fields. In the Agulhas current, we identified 4 hidden relationships corresponding to different responses of geostrophic currents to surface temperature local variations. Results are published in the journal Transactions on Geoscience and Remote Sensing. Contact me (pierre.tandeo@imt-atlantique) if you want to get the Matlab version of the latent class regression.
Then, using the 4 dynamical modes identified by the latent class regression, I generated surface currents from single SST snapshots. Results are presented in this proceeding of the ICIP conference (10% best paper award).
I am also interested in the synergy between ocean colour, sea surface temperature and altimetry in the Malvinas current. See this presentation at the COSPAR meeting for more details.
Spatio-temporal variability and satellite interpolation:
During my Ph.D. at IFREMER, I was working on the interpolation of satellite data using statistical approaches. First, I learned the SST temporal variability using infrared METOP-AVHRR data. Results are published in the journal Stochastic Environmental Research and Risk Assessment.
Then using geostatistical tools on the same satellite data, I learned the SST spatial variability using anisotropic structures represented as ellipsoids to model the spatial correlation lengths and angles. Results are published in the journal Remote Sensing of Environment.
Finally, using the spatial and temporal correlation lengths, we apply a Kalman filter and smoother to interpolate raw SST satellite data (with missing values due to clouds). Results are published in the Copernicus Marine Environment Monitoring Service (at http://marine.copernicus.eu/, SST_NWS_SST_L4_REP_OBSERVATIONS_010_023 product).
Correction of atmospheric conditions:
Another part of my work is dedicated to the correction of SST data using drifting buoys. The idea is to fit a statistical model between delta SST (difference between satellite and in situ) and the atmospheric conditions. Results for the AATSR sensor onboard Envisat are published in the journal Geoscience and Remote Sensing Letters.