Project

General

Profile

Conservative remapping of emission fields

Added by Diego Jimenez 4 months ago

Hi CDOers

I have a database of methane emissions (mass per time per surface area) that is in a regular lat-lon grid. For my project, I would like to remap these emissions into a similarly-sized ICON grid. The net emission should be at least approximately conserved, so I use conservative remapping (remapcon or remapcon2). However, when I compute the net emissions, I concluded that the conservative remappings were not very conservative.

Since emissions are localised, the original fields should have very sharp gradients. I suspect that both conservative remappings are having a hard time with the strong gradients. Therefore, I have one question:

Is there any (soon-to-be) alternative to remap fields with sharp gradients in CDO?

Cheers,

Diego


Replies (4)

RE: Conservative remapping of emission fields - Added by Estanislao Gavilan 4 months ago

Hi Diego,

I am not sure if they are developing new remap tools. If you were aware of any good remap algorithm, you could post requesting the implementation in future versions of CDO.

Kind regards,

Estanislao

RE: Conservative remapping of emission fields - Added by Ralf Mueller 4 months ago

hi!

I believe that conservative remapping should be able to work with strong gradients. but there are corner cases, esp. when source and target cells are very close to each other. For further investigation the input and output data would be very helpful for debugging and verification.

thx in advance + merry Christmas
ralf

RE: Conservative remapping of emission fields - Added by Diego Jimenez 4 months ago

Hi Estanislao and Ralf,

thanks for your answers. I have right now no specific example to send but the main problem for anthropogenic emissions is that one usually has a discontinuous field that should be non-negative, but after the remapping one obtains new extrema outside of the original range (due to the Gibbs phenomenon), leading to negative emissions or increased emissions and, therefore, no mass flux conservation between the original and the remapped field. I have issued a ticket, giving some more explanation and two references in the literature, one showing the problem and evaluating a solution and another describing the algorithm for the solution. Here, I give the references:

Explanation of the problem and evaluation of a solution: Marsico and Ullrich (2022)

The algorithm of the solution: Bradley et al. (2019)

Cheers,

Diego

PS. This problem is quite important in my project: We would like to verify emission estimates by taking the existing emission datasets as priors in a simulation.

RE: Conservative remapping of emission fields - Added by Estanislao Gavilan 4 months ago

hi Diego,

I have experienced some of those problems while remapping river runoff. Are those negative points random pixels? If that is the case you could filter that data with python, and then linearly interpolate as I did.

Kind regards,

Estanislao

    (1-4/4)