{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"Collapsed": "false"
},
"source": [
"# CAMS global atmospheric composition forecasts"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"```{hint} \n",
"Execute the notebook on the training platform >>\n",
"```"
]
},
{
"cell_type": "markdown",
"metadata": {
"Collapsed": "false"
},
"source": [
"This notebooks provides an introduction to the CAMS global atmospheric composition forecasts and shows you how the variable `Dust Aerosol Optical Depth` can be used to monitor and predict dust events.\n",
"\n",
"CAMS produces global forecasts for atmospheric composition for the next five days twice a day. The forecasts consist of more than 50 chemical species (e.g. ozone, nitrogen dioxide, carbon dioxide) and seven different types of aerosol (desert dust, sea salt, organic matter, black carbon, sulphate, nitrate and ammonium aerosol).\n",
"\n",
"The initial conditions of each forecast are obtained by combining a previous forecast with current satellite observations through a process called data assimilation. This best estimate of the state of the atmosphere at the initial forecast time step, called the analysis, provides a globally complete and consistent dataset allowing for estimates at locations where observation data coverage is low or for atmospheric pollutants for which no direct observations are available.\n",
"\n",
"The notebook examines the **Saharan Dust event** which occured over Europe in April 2024.\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"```{admonition} Basic facts\n",
"**Spatial resolution**: `0.4° x 0.4°`
\n",
"**Spatial coverage**: `Global`
\n",
"**Temporal resolution**: `1-hourly up to leadtime hour 120`
\n",
"**Temporal coverage**: `since 2015 to present`
\n",
"**Data format**: `GRIB` or `zipped NetCDF`\n",
"```"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"```{admonition} How to access the data\n",
"CAMS global atmospheric composition forecasts are available for download via the Copernicus Atmosphere Data Store (ADS). You will need to create an ADS account here.\n",
"\n",
"Data from the ADS can be downloaded in two ways:\n",
"* `manually` via the ADS web interface\n",
"* `programmatically` with a Python package called cdsapi (more information)\n",
"\n",
"```"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"
<xarray.Dataset> Size: 101MB\n", "Dimensions: (longitude: 900, latitude: 451, time: 31)\n", "Coordinates:\n", " * longitude (longitude) float32 4kB 0.0 0.4 0.8 1.2 ... 358.8 359.2 359.6\n", " * latitude (latitude) float32 2kB 90.0 89.6 89.2 88.8 ... -89.2 -89.6 -90.0\n", " * time (time) datetime64[ns] 248B 2024-04-07T12:00:00 ... 2024-04-11T...\n", "Data variables:\n", " duaod550 (time, latitude, longitude) float64 101MB ...\n", "Attributes:\n", " Conventions: CF-1.6\n", " history: 2024-05-30 10:00:01 GMT by grib_to_netcdf-2.25.1: /opt/ecmw...
<xarray.DataArray 'time' (time: 31)> Size: 248B\n", "array(['2024-04-07T12:00:00.000000000', '2024-04-07T15:00:00.000000000',\n", " '2024-04-07T18:00:00.000000000', '2024-04-07T21:00:00.000000000',\n", " '2024-04-08T00:00:00.000000000', '2024-04-08T03:00:00.000000000',\n", " '2024-04-08T06:00:00.000000000', '2024-04-08T09:00:00.000000000',\n", " '2024-04-08T12:00:00.000000000', '2024-04-08T15:00:00.000000000',\n", " '2024-04-08T18:00:00.000000000', '2024-04-08T21:00:00.000000000',\n", " '2024-04-09T00:00:00.000000000', '2024-04-09T03:00:00.000000000',\n", " '2024-04-09T06:00:00.000000000', '2024-04-09T09:00:00.000000000',\n", " '2024-04-09T12:00:00.000000000', '2024-04-09T15:00:00.000000000',\n", " '2024-04-09T18:00:00.000000000', '2024-04-09T21:00:00.000000000',\n", " '2024-04-10T00:00:00.000000000', '2024-04-10T03:00:00.000000000',\n", " '2024-04-10T06:00:00.000000000', '2024-04-10T09:00:00.000000000',\n", " '2024-04-10T12:00:00.000000000', '2024-04-10T15:00:00.000000000',\n", " '2024-04-10T18:00:00.000000000', '2024-04-10T21:00:00.000000000',\n", " '2024-04-11T00:00:00.000000000', '2024-04-11T03:00:00.000000000',\n", " '2024-04-11T06:00:00.000000000'], dtype='datetime64[ns]')\n", "Coordinates:\n", " * time (time) datetime64[ns] 248B 2024-04-07T12:00:00 ... 2024-04-11T06...\n", "Attributes:\n", " long_name: time
<xarray.DataArray 'latitude' (latitude: 451)> Size: 2kB\n", "array([ 90. , 89.6, 89.2, ..., -89.2, -89.6, -90. ], dtype=float32)\n", "Coordinates:\n", " * latitude (latitude) float32 2kB 90.0 89.6 89.2 88.8 ... -89.2 -89.6 -90.0\n", "Attributes:\n", " units: degrees_north\n", " long_name: latitude
<xarray.DataArray 'longitude' (longitude: 900)> Size: 4kB\n", "array([ 0. , 0.4, 0.8, ..., 358.8, 359.2, 359.6], dtype=float32)\n", "Coordinates:\n", " * longitude (longitude) float32 4kB 0.0 0.4 0.8 1.2 ... 358.8 359.2 359.6\n", "Attributes:\n", " units: degrees_east\n", " long_name: longitude
<xarray.Dataset> Size: 101MB\n", "Dimensions: (latitude: 451, time: 31, longitude: 900)\n", "Coordinates:\n", " * latitude (latitude) float32 2kB 90.0 89.6 89.2 88.8 ... -89.2 -89.6 -90.0\n", " * time (time) datetime64[ns] 248B 2024-04-07T12:00:00 ... 2024-04-11T...\n", " * longitude (longitude) float32 4kB -180.0 -179.6 -179.2 ... 179.2 179.6\n", "Data variables:\n", " duaod550 (time, latitude, longitude) float64 101MB ...\n", "Attributes:\n", " Conventions: CF-1.6\n", " history: 2024-05-30 10:00:01 GMT by grib_to_netcdf-2.25.1: /opt/ecmw...
<xarray.DataArray 'duaod550' (time: 31, latitude: 451, longitude: 900)> Size: 101MB\n", "[12582900 values with dtype=float64]\n", "Coordinates:\n", " * latitude (latitude) float32 2kB 90.0 89.6 89.2 88.8 ... -89.2 -89.6 -90.0\n", " * time (time) datetime64[ns] 248B 2024-04-07T12:00:00 ... 2024-04-11T...\n", " * longitude (longitude) float32 4kB -180.0 -179.6 -179.2 ... 179.2 179.6\n", "Attributes:\n", " units: ~\n", " long_name: Dust Aerosol Optical Depth at 550nm
<xarray.DataArray 'duaod550' (time: 31, latitude: 107, longitude: 162)> Size: 4MB\n", "array([[[5.16533464e-03, 5.26098837e-03, 5.30881523e-03, ...,\n", " 1.57831991e-03, 1.53049305e-03, 1.43483932e-03],\n", " [5.45229582e-03, 5.40446896e-03, 5.30881523e-03, ...,\n", " 1.53049305e-03, 1.48266619e-03, 1.43483932e-03],\n", " [5.11750778e-03, 5.06968092e-03, 4.92620033e-03, ...,\n", " 1.48266619e-03, 1.43483932e-03, 1.38701246e-03],\n", " ...,\n", " [4.78302063e-04, 6.69609515e-04, 9.08743831e-04, ...,\n", " 2.67304371e-01, 3.71901720e-01, 4.89221015e-01],\n", " [1.00439756e-03, 1.24353187e-03, 1.38701246e-03, ...,\n", " 9.00101897e-02, 1.21049824e-01, 1.63280944e-01],\n", " [1.67397364e-03, 1.81745423e-03, 1.91310795e-03, ...,\n", " 9.02493240e-02, 9.81407564e-02, 1.46206754e-01]],\n", "\n", " [[4.83054660e-03, 4.92620033e-03, 4.92620033e-03, ...,\n", " 1.05222442e-03, 1.00439756e-03, 9.56570694e-04],\n", " [5.26098837e-03, 5.16533464e-03, 5.06968092e-03, ...,\n", " 1.43483932e-03, 1.38701246e-03, 1.33918560e-03],\n", " [5.26098837e-03, 5.11750778e-03, 4.97402719e-03, ...,\n", " 1.67397364e-03, 1.62614678e-03, 1.62614678e-03],\n", "...\n", " [1.16601926e-01, 1.83750841e-01, 2.56543327e-01, ...,\n", " 1.92933599e-01, 2.09625174e-01, 2.27416767e-01],\n", " [1.77820310e-01, 2.46930127e-01, 2.94804817e-01, ...,\n", " 2.11825210e-01, 2.22634081e-01, 2.32582068e-01],\n", " [2.30382033e-01, 2.82895928e-01, 3.25653144e-01, ...,\n", " 2.13020881e-01, 2.21534063e-01, 2.31721185e-01]],\n", "\n", " [[1.53049305e-03, 1.57831991e-03, 1.72180050e-03, ...,\n", " 2.96529894e-03, 2.77399149e-03, 2.53485717e-03],\n", " [1.57831991e-03, 1.57831991e-03, 1.62614678e-03, ...,\n", " 7.26971661e-03, 6.55231367e-03, 5.78708386e-03],\n", " [1.62614678e-03, 1.62614678e-03, 1.62614678e-03, ...,\n", " 1.10480388e-02, 1.16697880e-02, 1.08089045e-02],\n", " ...,\n", " [2.81221988e-01, 2.90356919e-01, 2.91743898e-01, ...,\n", " 1.58067816e-01, 1.68876687e-01, 1.80976883e-01],\n", " [2.85048137e-01, 2.82130699e-01, 3.00878829e-01, ...,\n", " 1.66533171e-01, 1.77246388e-01, 1.88724835e-01],\n", " [2.83135063e-01, 2.82369833e-01, 3.10826816e-01, ...,\n", " 1.77581176e-01, 1.89346584e-01, 1.99772840e-01]]])\n", "Coordinates:\n", " * latitude (latitude) float32 428B 70.8 70.4 70.0 69.6 ... 29.2 28.8 28.4\n", " * time (time) datetime64[ns] 248B 2024-04-07T12:00:00 ... 2024-04-11T...\n", " * longitude (longitude) float32 648B -21.6 -21.2 -20.8 ... 42.0 42.4 42.8\n", "Attributes:\n", " units: ~\n", " long_name: Dust Aerosol Optical Depth at 550nm