{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"tags": []
},
"source": [
"# Metop-B IASI - Total Column Carbon Monoxide - Level 2"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"```{hint} \n",
"Execute the notebook on the training platform >>\n",
"```"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The [Infrared Atmospheric Sounding Interferometer (IASI)](https://www.eumetsat.int/iasi) is an instrument onboard the Metop-B/C satellites. It was also previously aboard Metop-A. It provides information on the vertical structure of temperature and humidity as well as main atmospheric species.\n",
"\n",
"This notebook provides you an introduction to data from Metop-B IASI. This dataset is used as a proxy for data from [IASI-NG](https://www.eumetsat.int/eps-sg-iasi-ng), a passive infrared sounder which has the capability to measure the temperature and water vapour profiles of the Earth’s atmosphere. IASI Level 2 data can be downloaded from the [IASI portal](https://iasi.aeris-data.fr/).\n",
"\n",
"The event featured is the [August Complex fire](https://www.fire.ca.gov/incidents/2020/8/16/august-complex-includes-doe-fire/) in California, USA in 2020. This was the largest wildfire in CA history, spreading over 1,000,000 acres (over 4,000 sq km). The image shown in this notebook is taken from 11 September 2020."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"```{admonition} Basic Facts\n",
"**Spatial resolution**: `4 x 12-km coverage close to the centre of a 48 x 48 km2 cell (average sampling distance: 24 km)`
\n",
"**Spatial coverage**: `Near global`
\n",
"**Revisit time**: `less than one day`
\n",
"**Data availability**: `since 2007`\n",
"```"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"```{admonition} How to access the data\n",
"IASI Level 2 are disseminated in the `netCDF` format and can be downloaded via the [IASI portal](https://iasi.aeris-data.fr/).\n",
"```"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"
<xarray.Dataset>\n", "Dimensions: (time: 542455, nlayers: 19, npressures: 20)\n", "Coordinates:\n", " * time (time) float64 4.322e+08 ... 4.323e+08\n", "Dimensions without coordinates: nlayers, npressures\n", "Data variables: (12/21)\n", " time_string (time) |S16 ...\n", " time_in_day (time) float64 ...\n", " latitude (time) float32 ...\n", " longitude (time) float32 ...\n", " solar_zenith_angle (time) float32 ...\n", " satellite_zenith_angle (time) float32 ...\n", " ... ...\n", " CO_total_column (time) float32 ...\n", " CO_total_column_error (time) float32 ...\n", " CO_degrees_of_freedom (time) float32 ...\n", " air_partial_column_profile (time, nlayers) float32 ...\n", " atmosphere_pressure_grid (time, npressures) float32 ...\n", " averaging_kernel_matrix (time, nlayers, nlayers) float32 ...\n", "Attributes: (12/31)\n", " title: IASI/METOPB ULB-LATMOS carbon monoxide (CO) ...\n", " institution: ULB-LATMOS for algorithm development ; EUMET...\n", " product_version: 6.5.0\n", " history: 2020-09-12 06:32:03 (date of data extraction...\n", " summary: This dataset contains Level 2 carbon monoxid...\n", " source: EUMETSAT IASI Level 2 carbon monoxide (CO) d...\n", " ... ...\n", " sensor: IASI\n", " spatial_resolution: 12km at nadir\n", " creator_type: institution\n", " creator_name: ULB-LATMOS\n", " contact_email: contact form at http://iasi.aeris-data.fr/co...\n", " data_policy: see https://iasi.aeris-data.fr/data-use-policy/
array([4.321728e+08, 4.321728e+08, 4.321728e+08, ..., 4.322592e+08,\n", " 4.322592e+08, 4.322592e+08])
[542455 values with dtype=|S16]
[542455 values with dtype=float64]
[542455 values with dtype=float32]
[542455 values with dtype=float32]
[542455 values with dtype=float32]
[542455 values with dtype=float32]
[542455 values with dtype=int64]
[542455 values with dtype=int32]
[542455 values with dtype=int32]
[542455 values with dtype=int32]
[542455 values with dtype=int32]
[542455 values with dtype=float32]
[10306645 values with dtype=float32]
[10306645 values with dtype=float32]
[10306645 values with dtype=float32]
[542455 values with dtype=float32]
[542455 values with dtype=float32]
[542455 values with dtype=float32]
[10306645 values with dtype=float32]
[10849100 values with dtype=float32]
[195826255 values with dtype=float32]
<xarray.DataArray 'CO_total_column' (time: 542455)>\n", "[542455 values with dtype=float32]\n", "Coordinates:\n", " * time (time) float64 4.322e+08 4.322e+08 ... 4.323e+08 4.323e+08\n", "Attributes:\n", " units: mol m-2\n", " long_name: retrieved carbon ...\n", " standard_name: atmosphere_mole_c...\n", " ancillary_variables: CO_total_column_e...\n", " multiplication_factor_to_convert_to_molecules_per_cm2: 6.02214086e+19
[542455 values with dtype=float32]
array([4.321728e+08, 4.321728e+08, 4.321728e+08, ..., 4.322592e+08,\n", " 4.322592e+08, 4.322592e+08])
<xarray.DataArray 'CO_total_column' (ground_pixel: 542455)>\n", "array([0.03949106, 0.03213775, 0.04679925, ..., 0.03180185, 0.03522488,\n", " 0.03115139], dtype=float32)\n", "Coordinates:\n", " latitude (ground_pixel) float32 8.665 8.915 8.798 ... 84.91 84.73 85.22\n", " longitude (ground_pixel) float32 -47.95 -48.01 -48.49 ... -66.97 -54.55\n", "Dimensions without coordinates: ground_pixel\n", "Attributes:\n", " long_name: retrieved carbon monoxide total column in mole/m2\n", " units: mol m-2
array([0.03949106, 0.03213775, 0.04679925, ..., 0.03180185, 0.03522488,\n", " 0.03115139], dtype=float32)
array([ 8.664803, 8.915398, 8.798302, ..., 84.9109 , 84.7296 ,\n", " 85.2249 ], dtype=float32)
array([-47.949097, -48.0121 , -48.4937 , ..., -68.1499 , -66.9676 ,\n", " -54.545998], dtype=float32)
<xarray.DataArray 'retrieval_quality_flag' (time: 542455)>\n", "[542455 values with dtype=int32]\n", "Coordinates:\n", " * time (time) float64 4.322e+08 4.322e+08 ... 4.323e+08 4.323e+08\n", "Attributes:\n", " long_name: retrieval quality flag summarizing processing flags\n", " comment: = 2 for the most reliable pixels; = 1 for the valuable pixels...
[542455 values with dtype=int32]
array([4.321728e+08, 4.321728e+08, 4.321728e+08, ..., 4.322592e+08,\n", " 4.322592e+08, 4.322592e+08])
<xarray.DataArray 'retrieval_quality_flag' (ground_pixel: 542455)>\n", "array([1, 1, 1, ..., 2, 2, 2], dtype=int32)\n", "Coordinates:\n", " latitude (ground_pixel) float32 8.665 8.915 8.798 ... 84.91 84.73 85.22\n", " longitude (ground_pixel) float32 -47.95 -48.01 -48.49 ... -66.97 -54.55\n", "Dimensions without coordinates: ground_pixel\n", "Attributes:\n", " long_name: retrieval quality flag summarizing processing flags\n", " units: -
array([1, 1, 1, ..., 2, 2, 2], dtype=int32)
array([ 8.664803, 8.915398, 8.798302, ..., 84.9109 , 84.7296 ,\n", " 85.2249 ], dtype=float32)
array([-47.949097, -48.0121 , -48.4937 , ..., -68.1499 , -66.9676 ,\n", " -54.545998], dtype=float32)
<xarray.DataArray (ground_pixel: 392634)>\n", "array([0.02875579, 0.0296423 , 0.03025403, ..., 0.03180185, 0.03522488,\n", " 0.03115139], dtype=float32)\n", "Coordinates:\n", " latitude (ground_pixel) float32 9.989 10.37 10.54 ... 84.91 84.73 85.22\n", " longitude (ground_pixel) float32 -43.75 -41.56 -41.6 ... -66.97 -54.55\n", "Dimensions without coordinates: ground_pixel\n", "Attributes:\n", " long_name: retrieved carbon monoxide total column in mole/m2\n", " units: mol m-2
array([0.02875579, 0.0296423 , 0.03025403, ..., 0.03180185, 0.03522488,\n", " 0.03115139], dtype=float32)
array([ 9.989197, 10.371399, 10.535301, ..., 84.9109 , 84.7296 ,\n", " 85.2249 ], dtype=float32)
array([-43.745903, -41.5605 , -41.599 , ..., -68.1499 , -66.9676 ,\n", " -54.545998], dtype=float32)
<xarray.DataArray (ground_pixel: 392634)>\n", "array([1.7317139e+18, 1.7851012e+18, 1.8219404e+18, ..., 1.9151522e+18,\n", " 2.1212919e+18, 1.8759807e+18], dtype=float32)\n", "Coordinates:\n", " latitude (ground_pixel) float32 9.989 10.37 10.54 ... 84.91 84.73 85.22\n", " longitude (ground_pixel) float32 -43.75 -41.56 -41.6 ... -66.97 -54.55\n", "Dimensions without coordinates: ground_pixel
array([1.7317139e+18, 1.7851012e+18, 1.8219404e+18, ..., 1.9151522e+18,\n", " 2.1212919e+18, 1.8759807e+18], dtype=float32)
array([ 9.989197, 10.371399, 10.535301, ..., 84.9109 , 84.7296 ,\n", " 85.2249 ], dtype=float32)
array([-43.745903, -41.5605 , -41.599 , ..., -68.1499 , -66.9676 ,\n", " -54.545998], dtype=float32)