CEDS Global Emissions

We have a data loader for CEDS (Community Emissions Data System) global gridded anthropogenic emissions, CEDS.

CEDS provides monthly global anthropogenic emissions at 0.5° × 0.5° resolution from 1750 to 2023 in units of kg m⁻² s⁻¹, with 8 anthropogenic sectors.

Reference: Hoesly, R. M., et al. (2018). Historical (1750–2014) anthropogenic emissions of reactive gases and aerosols from the Community Emissions Data System (CEDS). Geoscientific Model Development, 11, 369–408. https://doi.org/10.5194/gmd-11-369-2018

Equations

This is what its equation system looks like:

using EarthSciData, EarthSciMLBase
using ModelingToolkit, DynamicQuantities, DataFrames
using ModelingToolkit: t
using DynamicQuantities: dimension
using Dates

domain = DomainInfo(
    DateTime(2016, 5, 1), DateTime(2016, 5, 2);
    lonrange = deg2rad(-180):deg2rad(2.5):deg2rad(175),
    latrange = deg2rad(-85):deg2rad(2):deg2rad(85),
    levrange = 1:10,
    u_proto = zeros(Float32, 1, 1, 1, 1)
)

emis = CEDS(domain; species = ["SO2", "NOx", "CO"])

\[ \begin{align} \mathtt{SO2\_em\_anthro}\left( t \right) &= \mathtt{SO2\_em\_anthro\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right) \\ \mathtt{NOx\_em\_anthro}\left( t \right) &= \mathtt{NOx\_em\_anthro\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right) \\ \mathtt{CO\_em\_anthro}\left( t \right) &= \mathtt{CO\_em\_anthro\_itp}\left( t + \mathtt{t\_ref}, \mathtt{lon}, \mathtt{lat} \right) \end{align} \]

Variables

Here are the variables in tabular format:

function table(vars)
    DataFrame(
        :Name => [string(Symbolics.tosymbol(v, escape = false)) for v in vars],
        :Units => [dimension(ModelingToolkit.get_unit(v)) for v in vars],
        :Description => [ModelingToolkit.getdescription(v) for v in vars]
    )
end
table(unknowns(emis))
3×3 DataFrame
RowNameUnitsDescription
StringDimensio…String
1SO2_em_anthrom⁻² kg s⁻¹SO2 Anthropogenic Emissions
2NOx_em_anthrom⁻² kg s⁻¹NOx Anthropogenic Emissions
3CO_em_anthrom⁻² kg s⁻¹CO Anthropogenic Emissions

Parameters

Finally, here are the parameters in tabular format:

table(parameters(emis))
6×3 DataFrame
RowNameUnitsDescription
StringDimensio…String
1lonLongitude
2latLatitude
3t_refsReference time
4SO2_em_anthro_itpm⁻² kg s⁻¹Interpolated SO2_em_anthro
5NOx_em_anthro_itpm⁻² kg s⁻¹Interpolated NOx_em_anthro
6CO_em_anthro_itpm⁻² kg s⁻¹Interpolated CO_em_anthro