{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"Collapsed": "false"
},
"source": [
"# Copernicus Atmosphere Monitoring Service (CAMS) \n",
"# Global Fire Assimilation System (GFAS) data"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"```{hint} \n",
"Execute the notebook on the training platform >>\n",
"```"
]
},
{
"cell_type": "markdown",
"metadata": {
"Collapsed": "false",
"tags": []
},
"source": [
"The Copernicus Atmopshere Monitoring Service (CAMS) provides consistent and quality-controlled information related to air pollution and health and greenhouse gases. CAMS data consist of `global forecasts and analyses`, `global reanalyses`, `fire emissions` and `greenhouse gas flux inversions`. \n",
"\n",
"This notebook provides an introduction to the [CAMS Global Fire Assimilation System (GFAS) global fire emissions data](https://atmosphere.copernicus.eu/global-fire-emissions). The CAMS Global Fire Assimilation System (GFAS) assimilates fire radiative power (FRP) observations from satellite-based sensors to produce daily estimates of wildfire and biomass burning emissions. It also provides information about injection heights derived from fire observations and meteorological information from the operational weather forecasts of ECMWF.\n",
"\n",
"The GFAS data output includes spatially gridded `Fire Radiative Power (FRP)`, `dry matter burnt` and `biomass burning emissions` for a large set of chemical, greenhouse gas and aerosol species. Data are available globally on a regular latitude-longitude grid with horizontal resolution of 0.1 degrees from 2003 to present."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"```{admonition} Basic Facts\n",
"**Spatial resolution**: `~10km`
\n",
"**Spatial coverage**: `Global`
\n",
"**Time step**: `Daily data packaged in a monthly file`
\n",
"**Data availability**: `since 2003`\n",
"```"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"```{admonition} How to access the data\n",
"\n",
"CAMS GFAS data are available in either `GRIB` or `netCDF` format. They can be downloaded from the [Copernicus Atmosphere Data Store (ADS)](https://ads.atmosphere.copernicus.eu/cdsapp#!/dataset/cams-global-fire-emissions-gfas?tab=form). You will need to create an ADS account [here](https://ads.atmosphere.copernicus.eu/user/register). \n",
"\n",
"Get more information in the [CAMS GFAS data documentation](https://confluence.ecmwf.int/display/CKB/CAMS+global+biomass+burning+emissions+based+on+fire+radiative+power+%28GFAS%29%3A+data+documentation).\n",
"```"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"
<xarray.Dataset>\n", "Dimensions: (longitude: 3600, latitude: 1800, time: 31)\n", "Coordinates:\n", " * longitude (longitude) float32 0.05 0.15 0.25 0.35 ... 359.8 359.9 360.0\n", " * latitude (latitude) float32 89.95 89.85 89.75 ... -89.75 -89.85 -89.95\n", " * time (time) datetime64[ns] 2021-08-01 2021-08-02 ... 2021-08-31\n", "Data variables:\n", " frpfire (time, latitude, longitude) float32 ...\n", "Attributes:\n", " Conventions: CF-1.6\n", " history: 2022-02-15 11:56:28 GMT by grib_to_netcdf-2.24.0: grib_to_n...
array([5.0000e-02, 1.5000e-01, 2.5000e-01, ..., 3.5975e+02, 3.5985e+02,\n", " 3.5995e+02], dtype=float32)
array([ 89.95, 89.85, 89.75, ..., -89.75, -89.85, -89.95], dtype=float32)
array(['2021-08-01T00:00:00.000000000', '2021-08-02T00:00:00.000000000',\n", " '2021-08-03T00:00:00.000000000', '2021-08-04T00:00:00.000000000',\n", " '2021-08-05T00:00:00.000000000', '2021-08-06T00:00:00.000000000',\n", " '2021-08-07T00:00:00.000000000', '2021-08-08T00:00:00.000000000',\n", " '2021-08-09T00:00:00.000000000', '2021-08-10T00:00:00.000000000',\n", " '2021-08-11T00:00:00.000000000', '2021-08-12T00:00:00.000000000',\n", " '2021-08-13T00:00:00.000000000', '2021-08-14T00:00:00.000000000',\n", " '2021-08-15T00:00:00.000000000', '2021-08-16T00:00:00.000000000',\n", " '2021-08-17T00:00:00.000000000', '2021-08-18T00:00:00.000000000',\n", " '2021-08-19T00:00:00.000000000', '2021-08-20T00:00:00.000000000',\n", " '2021-08-21T00:00:00.000000000', '2021-08-22T00:00:00.000000000',\n", " '2021-08-23T00:00:00.000000000', '2021-08-24T00:00:00.000000000',\n", " '2021-08-25T00:00:00.000000000', '2021-08-26T00:00:00.000000000',\n", " '2021-08-27T00:00:00.000000000', '2021-08-28T00:00:00.000000000',\n", " '2021-08-29T00:00:00.000000000', '2021-08-30T00:00:00.000000000',\n", " '2021-08-31T00:00:00.000000000'], dtype='datetime64[ns]')
[200880000 values with dtype=float32]
<xarray.Dataset>\n", "Dimensions: (longitude: 3600, latitude: 1800, time: 31)\n", "Coordinates:\n", " * longitude (longitude) float32 -180.0 -179.9 -179.8 ... 179.8 179.9 180.0\n", " * latitude (latitude) float32 89.95 89.85 89.75 ... -89.75 -89.85 -89.95\n", " * time (time) datetime64[ns] 2021-08-01 2021-08-02 ... 2021-08-31\n", "Data variables:\n", " frpfire (time, latitude, longitude) float32 ...\n", "Attributes:\n", " Conventions: CF-1.6\n", " history: 2022-02-15 11:56:28 GMT by grib_to_netcdf-2.24.0: grib_to_n...
array([-179.95001, -179.85 , -179.75 , ..., 179.75 , 179.85 ,\n", " 179.95001], dtype=float32)
array([ 89.95, 89.85, 89.75, ..., -89.75, -89.85, -89.95], dtype=float32)
array(['2021-08-01T00:00:00.000000000', '2021-08-02T00:00:00.000000000',\n", " '2021-08-03T00:00:00.000000000', '2021-08-04T00:00:00.000000000',\n", " '2021-08-05T00:00:00.000000000', '2021-08-06T00:00:00.000000000',\n", " '2021-08-07T00:00:00.000000000', '2021-08-08T00:00:00.000000000',\n", " '2021-08-09T00:00:00.000000000', '2021-08-10T00:00:00.000000000',\n", " '2021-08-11T00:00:00.000000000', '2021-08-12T00:00:00.000000000',\n", " '2021-08-13T00:00:00.000000000', '2021-08-14T00:00:00.000000000',\n", " '2021-08-15T00:00:00.000000000', '2021-08-16T00:00:00.000000000',\n", " '2021-08-17T00:00:00.000000000', '2021-08-18T00:00:00.000000000',\n", " '2021-08-19T00:00:00.000000000', '2021-08-20T00:00:00.000000000',\n", " '2021-08-21T00:00:00.000000000', '2021-08-22T00:00:00.000000000',\n", " '2021-08-23T00:00:00.000000000', '2021-08-24T00:00:00.000000000',\n", " '2021-08-25T00:00:00.000000000', '2021-08-26T00:00:00.000000000',\n", " '2021-08-27T00:00:00.000000000', '2021-08-28T00:00:00.000000000',\n", " '2021-08-29T00:00:00.000000000', '2021-08-30T00:00:00.000000000',\n", " '2021-08-31T00:00:00.000000000'], dtype='datetime64[ns]')
[200880000 values with dtype=float32]
<xarray.DataArray 'longitude' (longitude: 3600)>\n", "array([-179.95001, -179.85 , -179.75 , ..., 179.75 , 179.85 ,\n", " 179.95001], dtype=float32)\n", "Coordinates:\n", " * longitude (longitude) float32 -180.0 -179.9 -179.8 ... 179.8 179.9 180.0
array([-179.95001, -179.85 , -179.75 , ..., 179.75 , 179.85 ,\n", " 179.95001], dtype=float32)
array([-179.95001, -179.85 , -179.75 , ..., 179.75 , 179.85 ,\n", " 179.95001], dtype=float32)
<xarray.DataArray 'frpfire' (time: 31, latitude: 1800, longitude: 3600)>\n", "[200880000 values with dtype=float32]\n", "Coordinates:\n", " * longitude (longitude) float32 -180.0 -179.9 -179.8 ... 179.8 179.9 180.0\n", " * latitude (latitude) float32 89.95 89.85 89.75 ... -89.75 -89.85 -89.95\n", " * time (time) datetime64[ns] 2021-08-01 2021-08-02 ... 2021-08-31\n", "Attributes:\n", " units: W m**-2\n", " long_name: Wildfire radiative power
[200880000 values with dtype=float32]
array([-179.95001, -179.85 , -179.75 , ..., 179.75 , 179.85 ,\n", " 179.95001], dtype=float32)
array([ 89.95, 89.85, 89.75, ..., -89.75, -89.85, -89.95], dtype=float32)
array(['2021-08-01T00:00:00.000000000', '2021-08-02T00:00:00.000000000',\n", " '2021-08-03T00:00:00.000000000', '2021-08-04T00:00:00.000000000',\n", " '2021-08-05T00:00:00.000000000', '2021-08-06T00:00:00.000000000',\n", " '2021-08-07T00:00:00.000000000', '2021-08-08T00:00:00.000000000',\n", " '2021-08-09T00:00:00.000000000', '2021-08-10T00:00:00.000000000',\n", " '2021-08-11T00:00:00.000000000', '2021-08-12T00:00:00.000000000',\n", " '2021-08-13T00:00:00.000000000', '2021-08-14T00:00:00.000000000',\n", " '2021-08-15T00:00:00.000000000', '2021-08-16T00:00:00.000000000',\n", " '2021-08-17T00:00:00.000000000', '2021-08-18T00:00:00.000000000',\n", " '2021-08-19T00:00:00.000000000', '2021-08-20T00:00:00.000000000',\n", " '2021-08-21T00:00:00.000000000', '2021-08-22T00:00:00.000000000',\n", " '2021-08-23T00:00:00.000000000', '2021-08-24T00:00:00.000000000',\n", " '2021-08-25T00:00:00.000000000', '2021-08-26T00:00:00.000000000',\n", " '2021-08-27T00:00:00.000000000', '2021-08-28T00:00:00.000000000',\n", " '2021-08-29T00:00:00.000000000', '2021-08-30T00:00:00.000000000',\n", " '2021-08-31T00:00:00.000000000'], dtype='datetime64[ns]')
<xarray.DataArray 'frpfire' (time: 31, latitude: 100, longitude: 200)>\n", "array([[[0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " ...,\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.]],\n", "\n", " [[0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " ...,\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.]],\n", "\n", " [[0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " ...,\n", "...\n", " ...,\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.]],\n", "\n", " [[0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " ...,\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.]],\n", "\n", " [[0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " ...,\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.]]], dtype=float32)\n", "Coordinates:\n", " * longitude (longitude) float32 10.05 10.15 10.25 10.35 ... 29.75 29.85 29.95\n", " * latitude (latitude) float32 44.95 44.85 44.75 44.65 ... 35.25 35.15 35.05\n", " * time (time) datetime64[ns] 2021-08-01 2021-08-02 ... 2021-08-31\n", "Attributes:\n", " units: W m**-2\n", " long_name: Wildfire radiative power
array([[[0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " ...,\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.]],\n", "\n", " [[0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " ...,\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.]],\n", "\n", " [[0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " ...,\n", "...\n", " ...,\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.]],\n", "\n", " [[0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " ...,\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.]],\n", "\n", " [[0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " ...,\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.],\n", " [0., 0., 0., ..., 0., 0., 0.]]], dtype=float32)
array([10.050003, 10.149994, 10.25 , 10.350006, 10.449997, 10.550003,\n", " 10.649994, 10.75 , 10.850006, 10.949997, 11.050003, 11.149994,\n", " 11.25 , 11.350006, 11.449997, 11.550003, 11.649994, 11.75 ,\n", " 11.850006, 11.949997, 12.050003, 12.149994, 12.25 , 12.350006,\n", " 12.449997, 12.550003, 12.649994, 12.75 , 12.850006, 12.949997,\n", " 13.050003, 13.149994, 13.25 , 13.350006, 13.449997, 13.550003,\n", " 13.649994, 13.75 , 13.850006, 13.949997, 14.050003, 14.149994,\n", " 14.25 , 14.350006, 14.449997, 14.550003, 14.649994, 14.75 ,\n", " 14.850006, 14.949997, 15.050003, 15.149994, 15.25 , 15.350006,\n", " 15.449997, 15.550003, 15.649994, 15.75 , 15.850006, 15.949997,\n", " 16.050003, 16.149994, 16.25 , 16.350006, 16.449997, 16.550003,\n", " 16.649994, 16.75 , 16.850006, 16.949997, 17.050003, 17.149994,\n", " 17.25 , 17.350006, 17.449997, 17.550003, 17.649994, 17.75 ,\n", " 17.850006, 17.949997, 18.050003, 18.149994, 18.25 , 18.350006,\n", " 18.449997, 18.550003, 18.649994, 18.75 , 18.850006, 18.949997,\n", " 19.050003, 19.149994, 19.25 , 19.350006, 19.449997, 19.550003,\n", " 19.649994, 19.75 , 19.850006, 19.949997, 20.050003, 20.149994,\n", " 20.25 , 20.350006, 20.449997, 20.550003, 20.649994, 20.75 ,\n", " 20.850006, 20.949997, 21.050003, 21.149994, 21.25 , 21.350006,\n", " 21.449997, 21.550003, 21.649994, 21.75 , 21.850006, 21.949997,\n", " 22.050003, 22.149994, 22.25 , 22.350006, 22.449997, 22.550003,\n", " 22.649994, 22.75 , 22.850006, 22.949997, 23.050003, 23.149994,\n", " 23.25 , 23.350006, 23.449997, 23.550003, 23.649994, 23.75 ,\n", " 23.850006, 23.949997, 24.050003, 24.149994, 24.25 , 24.350006,\n", " 24.449997, 24.550003, 24.649994, 24.75 , 24.850006, 24.949997,\n", " 25.050003, 25.149994, 25.25 , 25.350006, 25.449997, 25.550003,\n", " 25.649994, 25.75 , 25.850006, 25.949997, 26.050003, 26.149994,\n", " 26.25 , 26.350006, 26.449997, 26.550003, 26.649994, 26.75 ,\n", " 26.850006, 26.949997, 27.050003, 27.149994, 27.25 , 27.350006,\n", " 27.449997, 27.550003, 27.649994, 27.75 , 27.850006, 27.949997,\n", " 28.050003, 28.149994, 28.25 , 28.350006, 28.449997, 28.550003,\n", " 28.649994, 28.75 , 28.850006, 28.949997, 29.050003, 29.149994,\n", " 29.25 , 29.350006, 29.449997, 29.550003, 29.649994, 29.75 ,\n", " 29.850006, 29.949997], dtype=float32)
array([44.95, 44.85, 44.75, 44.65, 44.55, 44.45, 44.35, 44.25, 44.15, 44.05,\n", " 43.95, 43.85, 43.75, 43.65, 43.55, 43.45, 43.35, 43.25, 43.15, 43.05,\n", " 42.95, 42.85, 42.75, 42.65, 42.55, 42.45, 42.35, 42.25, 42.15, 42.05,\n", " 41.95, 41.85, 41.75, 41.65, 41.55, 41.45, 41.35, 41.25, 41.15, 41.05,\n", " 40.95, 40.85, 40.75, 40.65, 40.55, 40.45, 40.35, 40.25, 40.15, 40.05,\n", " 39.95, 39.85, 39.75, 39.65, 39.55, 39.45, 39.35, 39.25, 39.15, 39.05,\n", " 38.95, 38.85, 38.75, 38.65, 38.55, 38.45, 38.35, 38.25, 38.15, 38.05,\n", " 37.95, 37.85, 37.75, 37.65, 37.55, 37.45, 37.35, 37.25, 37.15, 37.05,\n", " 36.95, 36.85, 36.75, 36.65, 36.55, 36.45, 36.35, 36.25, 36.15, 36.05,\n", " 35.95, 35.85, 35.75, 35.65, 35.55, 35.45, 35.35, 35.25, 35.15, 35.05],\n", " dtype=float32)
array(['2021-08-01T00:00:00.000000000', '2021-08-02T00:00:00.000000000',\n", " '2021-08-03T00:00:00.000000000', '2021-08-04T00:00:00.000000000',\n", " '2021-08-05T00:00:00.000000000', '2021-08-06T00:00:00.000000000',\n", " '2021-08-07T00:00:00.000000000', '2021-08-08T00:00:00.000000000',\n", " '2021-08-09T00:00:00.000000000', '2021-08-10T00:00:00.000000000',\n", " '2021-08-11T00:00:00.000000000', '2021-08-12T00:00:00.000000000',\n", " '2021-08-13T00:00:00.000000000', '2021-08-14T00:00:00.000000000',\n", " '2021-08-15T00:00:00.000000000', '2021-08-16T00:00:00.000000000',\n", " '2021-08-17T00:00:00.000000000', '2021-08-18T00:00:00.000000000',\n", " '2021-08-19T00:00:00.000000000', '2021-08-20T00:00:00.000000000',\n", " '2021-08-21T00:00:00.000000000', '2021-08-22T00:00:00.000000000',\n", " '2021-08-23T00:00:00.000000000', '2021-08-24T00:00:00.000000000',\n", " '2021-08-25T00:00:00.000000000', '2021-08-26T00:00:00.000000000',\n", " '2021-08-27T00:00:00.000000000', '2021-08-28T00:00:00.000000000',\n", " '2021-08-29T00:00:00.000000000', '2021-08-30T00:00:00.000000000',\n", " '2021-08-31T00:00:00.000000000'], dtype='datetime64[ns]')