SMOS BUFR to GRIB converter (SMOS_BUFR2GRIB)

Introduction

Fig1: SMOS Brightness Temperature (K) at 30 (+/-1) degrees of incidence angle and at XX polarisation, on 02 February 2010, plotted from the GRIB file obtained by convertion with the SMOS_BUFR2GRIB software.

For offline forward modelling and assimilation studies, as well as for bias correction developments, requesting the use of long time series, a light format and a noise filterd product are requested.
To this end, ECMWF developed a convertion software to get SMOS brightness temperature products in GRIB format on a grid compatible with that used for the land surface model and microwave emission model, at all resolutions. In addition the software includes angular binning and observation spatial average which allow to filter out the noise in the data.
The SMOS_BUFR2GRIB converter uses the SMOS BUFR file information content only. Since the BUFRfile contains a lot of information on the altitude, water fraction etc, there is no need to use auxilliary data for the conversion.

Description of the SMOS_BUFR2GRIB software

The SMOS_BUFR2GRIB software is suitable to be used on Linux or Unix platforms. It converts SMOS NRT BUFR files to GRIB files, either on a regular Latitue-longitude (LL) grid or a Gaussian Reduced (GR) grid. The software works in two steps:

• First a metview programme converts BUFR to GEOPOINT and then GEOPOINT to GRIB on a regular Latitude-Longitude Grid.

• Second, an interpolation software (scin) is used to interpolate the regular Latitude-Longitude (LL) GRIB file on a regular Gaussian grid if requested. Scin is available at ECMWF, so at the moment users outside ECMWF can not use it. The SMOS_BUFR2GRIB proposes the option to skip the step2 (conveting the BFR files to a regular LL grid only).

Fig2: SMOS_BUFR2GRIB Software

The control script is smos_bufr2grib.sh ; When editing the script, the user has to define the input file name and the requested incidence angle(s), angular binning angle, as well as the output grid type (regular Latitude-Longitude or Gaussian reduced) and the resolution. The user also has to choose if he would like to converted only the Alias Free Field of View data (AFOV), or the entire FOV. It is recommended to use the alias free field of view to ensure good data quality. The script launches the metview macro smos_bufr2gribLL.met

In the default smos_bufr2gribLL.met programme only land data are converted to GRIB. To convert the data over ocean only or on both ocean and land, change the the water fraction custom values, line 82 of smos_bufr2gribLL.met.

To run the converter place it in the directory where the smos file is (or give the full path/name of the input file) and type: sh smos_bufr2grib.sh

Requested tools:
- Linux or Unix computer
- Metview
- SCIN: ECMWF New interpolation software. Not yet available outside ECMWF, so skip this part if not at ECMWF by puting grid_resol="regular" in smos_bufr2grib.sh

Note on the SMOS NRT BUFR file content and the paramater identification numbers: radiometric accuracy: 012167, snapshot accuracy: 012166 (Corbella), polarization: 002099 (0h, 1V), snapshot overallquality: 033028 (must be =1 but is always=3,so far (Feb 2012)), SMOS info flag: 025174, Faraday angle: 025083, TB: 12080. For more information see the SMOS BUFR specification document.

Download the Software

Metview programme: smos_bufr2gribLL.met

Readme file: README

Examples

Example1:

The two brightness temperature output 0.25x0.25 GRIB files are each 0.208 Mb. The entire output (ten files) is 2Mb:
smosbufr_Faxx_0.25_angle30.grib: Faraday angle for XX polarised data
smosbufr_Fayy_0.25_angle30.grib: Faraday angle for YY polarised data
smosbufr_flagfov5_xx_0.25_angle30.grib: SMOS information flag, bit5 (AFOV), for XX polarised data
smosbufr_flagfov5_yy_0.25_angle30.grib: SMOS inforamtion flag, bit5 (AFOV), for YY polarised data
smosbufr_snpaxx_0.25_angle30.grib: Snapshot accuracy for XX polarised data
smosbufr_snpayy_0.25_angle30.grib: Snapshot accuracy for YY polarised data
smosbufr_snpqxx_0.25_angle30.grib: Snapshot quality for XX polarised data
smosbufr_snpqyy_0.25_angle30.grib: Snapshot quality for YY polarised data
smosbufr_tbxx_0.25_angle30.grib: Brightness Temperature for XX polarised data
smosbufr_tbyy_0.25_angle30.grib: Brightness Temperature for YY polarised data

The following figure illustrates the obtained GRIB 0.25x0.25 TB at XX polarisation:

Fig3: SMOS brightness temperature (K) at XX polatisation on 07 February 2010, on a regular 0.25 lat-lon grid, at global scale (top) and over Europe (bottom).

Example2:

Convert SMOS BUFR data for incidence angles of 30 (+/-1) degrees, on a Gaussian T511 grid. Do do so, select the following options in smos_bufr2grib.sh :
bin_angle="1"
reg_grid_resol="0.25"
grid="gaussian"
gr_resol="256" # grid N256 corresponds to T511 (see Gaussian grid )
inputfile="smosbufr"
afov="yes"
for AA in 30

The Gaussian Reduced format output is 0.076 Mb for the TB files, and it is 0.64 Mb for the 10 output files at 30 (+/-1) degrees incidence angle.
GR_N256_smosbufr_Faxx_0.25_angle30.grib: Faraday angle for XX polarised data
GR_N256_smosbufr_Fayy_0.25_angle30.grib: Faraday angle for YY polarised data
GR_N256_smosbufr_flagfov5_xx_0.25_angle30.grib: SMOS information flag, bit5 (AFOV), for XX polarised data
GR_N256_smosbufr_flagfov5_yy_0.25_angle30.grib: SMOS information flag, bit5 (AFOV), for YY polarised data
GR_N256_smosbufr_snpaxx_0.25_angle30.grib: Snapshot accuracy for XX polarised data
GR_N256_smosbufr_snpayy_0.25_angle30.grib: Snapshot accuracy for YY polarised data
GR_N256_smosbufr_snpqxx_0.25_angle30.grib: Snapshot quality for XX polarised data
GR_N256_smosbufr_snpqyy_0.25_angle30.grib: Snapshot quality for YY polarised data
GR_N256_smosbufr_tbxx_0.25_angle30.grib: Brightness Temperature for XX polarised data
GR_N256_smosbufr_tbyy_0.25_angle30.grib: Brightness Temperature for YY polarised data

Fig4: SMOS brightness temperature (K) at XX polatisation on 07 February 2010, on a gaussian Reduced T511 grid, over Europe.

Comments: