First all, my question was different from the question Converting grid resolution from degrees to kilometers.

In chemical transport modeling, the pollutant source inventory was an essential input which contain residential, industrial, agriculturial, vehicle emission in the area of interest.

The source inventory can be built independently or downloaded from the websites of several institutions which support some global/regional source information dataset.

For example:

  1. INTEX-B

    It contain different emission species like SO2, NOx, CO, NMVOC in the resolution of 0.5 degree x 0.5 degree

  2. Tracer-P

  3. REAS -> Regional Emission inventory in ASia
    Spatial resolution: 0.25 degree by 0.25 degree

  4. EDGAR The Emissions Database for Global Atmospheric Research
    spatial resolution: 0.1deg x0.1deg and 0.5deg x0.5deg

These source inventories are all gridded in geographic coordinate. But for regional/small scale research, the simulation grid was set in the unit of kilometers(For example, 36 km for regional simulation, 4 km for urban simulation).

Then, I want to re-arrange the source inventory into km x km cell. Due to different distance corresponding to the same degrees gap along the latitude, the degree-cell can't be transform into kilometer-cell directly.

If directly, the width of each cell

$$\Delta y=2\pi R_e\frac{\Delta lat}{360^\circ}$$

The length of each cell:

$$\Delta x=2\pi R_e\cos(lat)\frac{\Delta lon}{360^\circ}$$

So, every cell has its own size in kilometers. The uneven grid networks enhance the difficulty for simulation. I'm wondering to transform the inventory into uniform grid in kilometers(for example: 0.25 degree x 0.25 degree -> 36 km x 36 km)

How to achieve it? I'm familiar with Python language.
Any advice would be appreciate!

  • $\begingroup$ Which projection do you use for the 36 km x 36 km grid? $\endgroup$ Jul 5, 2016 at 11:00
  • 1
    $\begingroup$ Remark 1: The approach, which was described here, to approximate the grid cell size is based on the assumption of ideal spherical Earth. This is fine for a first order guess of the grid cell size but it should not be used when we want to know the 'real' size of the grid cell (because the Earth is not ideal sphere). $\endgroup$ Jul 6, 2016 at 14:18
  • 1
    $\begingroup$ Remark 2b: Preserving two metrics is generally problematic and not possible for the whole globe. However, a projection might approximately preserve these metrics in a predefined region. See for example the Lambert conic equal area projection. We, for example, use these type of projection for our chemistry transport model simulations. $\endgroup$ Jul 6, 2016 at 14:38
  • 1
    $\begingroup$ sounds reasonable :-) $\endgroup$ Jul 7, 2016 at 8:03
  • 1
    $\begingroup$ A colleague of mine does the emission preparation for our whole working group. He uses his own interpolation routines. I use cdo (cdo remap) when I do interpolations. Maybe we should switch to the chat, if you have more questions? $\endgroup$ Jul 7, 2016 at 9:25


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.