Reference : Correcting Biases in a lower resolution global circulation model with data assimilation
Scientific congresses and symposiums : Poster
Physical, chemical, mathematical & earth Sciences : Earth sciences & physical geography
http://hdl.handle.net/2268/204824
Correcting Biases in a lower resolution global circulation model with data assimilation
English
Canter, Martin mailto [Université de Liège > Département d'astrophys., géophysique et océanographie (AGO) > GeoHydrodynamics and Environment Research (GHER) >]
Barth, Alexander mailto [Université de Liège > Département d'astrophys., géophysique et océanographie (AGO) > GeoHydrodynamics and Environment Research (GHER) >]
2016
No
No
EGU
17–22 April 2016
[en] Data assimilation
[en] With this work, we aim at developping a new method of bias correction using data assimilation. This method is
based on the stochastic forcing of a model to correct bias.

First, through a preliminary run, we estimate the bias of the model and its possible sources. Then, we establish a forcing term which is directly added inside the model’s equations. We create an ensemble of runs and consider the forcing term as a control variable during the assimilation of observations. We then use this analysed forcing term to correct the bias of the model. Since the forcing is added inside the model, it acts as a source term, unlike external forcings such as wind.

This procedure has been developed and successfully tested with a twin experiment on a Lorenz 95 model. It is currently being applied and tested on the sea ice ocean NEMO LIM model, which is used in the PredAntar project. NEMO LIM is a global and low resolution (2 degrees) coupled model (hydrodynamic model and sea ice model) with long time steps allowing simulations over several decades. Due to its low resolution, the model is subject to bias in area where strong currents are present. We aim at correcting this bias by using perturbed current fields from higher resolution models and randomly generated perturbations. The random perturbations need to be constrained in order to respect the physical properties of the ocean, and not create unwanted phenomena. To construct those random perturbations, we first create a random field with the Diva tool (Data-Interpolating Variational Analysis). Using a cost function, this tool penalizes abrupt variations in the field, while using a custom correlation length. It also decouples disconnected areas based on topography. Then, we filter the field to smoothen it and remove small scale variations. We use this field as a random stream function, and take its derivatives to get zonal and meridional velocity fields. We also constrain the stream function along the coasts in order not to have currents perpendicular to the coast. The randomly generated stochastic forcing are then directly injected into the NEMO LIM model’s equations in order to force the model at each timestep, and not only during the assimilation step.

Results from a twin experiment will be presented. This method is being applied to a real case, with observations on the sea surface height available from the mean dynamic topography of CNES (Centre national d’études spatiales). The model, the bias correction, and more extensive forcings, in particular with a three dimensional structure and a time-varying component, will also be presented.
GHER/MARE/AGO
Fonds de la Recherche Scientifique (Communauté française de Belgique) - F.R.S.-FNRS ; CECI
http://hdl.handle.net/2268/204824
http://meetingorganizer.copernicus.org/EGU2016/EGU2016-14835.pdf

File(s) associated to this reference

Fulltext file(s):

FileCommentaryVersionSizeAccess
Open access
EGU2016-14835.pdfPublisher postprint33.56 kBView/Open

Bookmark and Share SFX Query

All documents in ORBi are protected by a user license.