Layer Precipitable Water PGE07 NWCSAF    Table of contents 1.- Goal of LPW product  2.- LPW algorithm summary description  3.- List of inputs for LPW  4.- Coverage and resolution 5.- Description of LPW outputs 6.- Example of LPW visualisation Access to "User Manual for the PGE07 of the SAFNWC/MSG Scientific Part" for a more detailed description.

1.- Goal of LPW product

The LPW, developed within the SAF NWC context, aims to support nowcasting applications. The LPW allows measures of the precipitable water in different layers in cloud free areas to be obtained.

The central aim of the LPW is therefore to provide information on the water vapor contained in a vertical column of unit cross-section area in three layers in the troposphere:

• Boundary Layer (BL): 1013hPa-840 hPa
• Middle Layer (ML): 840hPa-437hPa
• High Layer (HL):  <437hPa

Besides, the product provides the total precipitable water parameter for validation purposes (TPW control). A conversion from mm to grey levels or digital count is necessary for writing BL, ML, HL and TPW control output in  8 bits for each one.

2.- LPW algorithm summary description

LPW's main algorithms have been devised as one neural network for each precipitable water layer.  It is possible that during future developments, other topology with different number of neurones in the hidden layer could present better performances. Therefore, the software has been designed to change the topology and weights only changing the name of the new neural networks in the configuration file. When the PGE07 is started, all the names of the files with the neural networks are read from the configuration file and the topology and weights are allocated and loaded.

The implementation of LPW algorithms are based mainly on simulated radiances, due to real satellite MSG data availability did not present enough time to use them during the first tuning and checking. 

The detailed description of these developments has been shown in different scientific documents during the DOP: Scientific Report from MTR on PGE07 Layer Precipitable Water, the scientific part of the Software User Manual, the EUM/STG/SWG/12/02/DOC/08 and two verbal presentations in the EUMETSAT Meteorological Satellite Data User's Conferences (1999&2001). And in two scientific documents during the IOP:  "Improvements in SEVIRI and MODIS LPW&SAI Neural Networks using 60L-SD and RTTOV-7. Applications to MODIS images" and "Scientific Report on First Checking and Tuning for MSG PGE07".

The CMa output is used to identify clear air pixels. Only if the pixel is labelled as clear air and the satellite zenith angle of this pixel lies within the configurable threshold the pre-processing and processing is calculated in this pixel.

The whole process includes:
 

• Pre-processing:
  • Smoothing of SEVIRI radiances (optional) performing a median o mean/average (configurable in the model configuration file) over a window centred at the processing pixel. Boundaries of the region are not considered.
  • Normalisation of SEVIRI radiances.
  • Correction of SEVIRI radiances to 45º MSG zenith angle. This step will use a dedicated neural network, whose inputs are the SEVIRI radiances and the satellite zenith angle.
• Processing:

  The calculus of each Layer precipitable water parameter will use a dedicated Neural Networks (NN). The NNs will be fed with corrected normalized IR SEVIRI radiance values (WV6.2, WV6.3, IR8.7, IR9.7, IR10.8, IR12.0 & IR13.4m m), and each one will provide as output a precipitable water parameter.
  For controlling purposes Total Precipitable Water (BL+ML+HL) is calculated.
  For land pixels the normalised altitude is also needed as input for the neural network and the neural networks weights that correspond to IR8.7 radiances are fixed to zero.

• Quality control:

  To check the product quality, different flags have been added to the output products.

  • Out of radiance flag: It indicates how many bands are out of range for product confidence purposes.
  • Number of cloudy pixels in the largest smoothing processing window: indicating the confidence of the smooth pre-processing step.
  • The quality flag: it is an internal control flag, This flag checks the difference between the TPW control and the sum of the three layers in each pixel.
  • The consistence flag: it shows the difference between TPWcontrol and PGE06(TPW) in sea pixels.
• Post-processing:

  The post-processing step will include:

  • Denormalisation of the precipitable water to mm, using the thresholds defined in the configuration file.
  • Smoothing of each layer when it is required in the model configuration file, given the proper parameters such as method to use (mean/average or median), length of the processing box and minimum number of cloud free pixels in the window to perform the smoothing. The process consists of making a mean/average or a median of the parameter over a window centred at the pixel.
  • Conversion of precipitable water for each layer from mm to counts.
  • Over cloudy pixel, a configurable IR channel is shown.

3.- List of inputs for LPW

• MSG radiance values: The following SEVIRI radiance values are needed at full IR spatial resolution:

R13.4 m	R12.0 m	R10.8m	R9.7m	R8.7m	R7.3m	R6.2m
Mandatory	Mandatory	Mandatory	Mandatory	Mandatory	Mandatory	Mandatory

The LPW system processing checks the availability of SEVIRI channels for each pixel; no results are produced for pixels where even only one channel is missing.
The SEVIRI channels are input by the user in HRIT format, and extracted on the processed region by SAFNWC software package.
 

• NWCSAF/MSG Products used as inputs:

• PGE01(Cma): The CMa is mandatory. After the CMa process is executed, the LPW system processing is computed.
• PGE06(TPW): The TPW is not mandatory. After the TPW process is executed, the LPW system processing is computed.

• Satellite zenith angles associated to selected region: This information is mandatory. It is computed by the LPW software itself, using a function available in the NWCLIB.
• Ancillary data sets: The following ancillary data, remapped onto satellite images, are mandatory:
  • Land/sea mask
  • Elevation mask

These ancillary data are available in the SAFNWC software package on MSG full disk in the default satellite projection at full IR resolution. They are extracted on the processed region by the LPW software itself.
Input Files: Region Configuration Files, Neuronal Network Files, PGE07 Model Configuration Files, are available in the SAFNWC software package, and is needed by the LPW software.

4.- Coverage and resolution

The software to extract Layer Precipitable Water (LPW) from MSG SEVIRI Imagery over MSG N (Europe, North Africa and adjacent seas) has been designed within the EUMETSAT SAF NWC.  The validity of this product for other areas will not be guaranteed. However, the processing of the  full MSG SEVIRI image will not be excluded  (in the case of a large area the process time increase and the 15 minutes may not be enough).  The selected sub area must be a rectangular area defined in satellite projection of numlin x numele SEVIRI IR pixels.

This product will be derived every 15 minutes on the SEVIRI IR pixel horizontal resolution, i.e. the current resolution will be given by the local pixel size, which depends on the latitude and longitude of the pixel. At the sub-satellite point, the resolution will be 3 km, whereas over Central  Europe the resolution will be approximately 5 km.
 

5.-Description of LPW outputs
 

• Description of the conversion from mm to brightness values

The units of precipitable water are mm. To allow for a resolution better than 1 mm in BL, ML, HL and TPW values it is necessary to use float values or to convert from mm to counts. The conversion to count is more efficient because the final output size is smaller.

Each LPW parameter output is written in 8 bits (0:255). The first 8 values are reserved values, the next 120 values are precipitable water information and the last 128 are the values of the configurable IR band degrade to 7 bits.

The maximum and minimum used in the LPW parameters are the following 

Parameter	Minimum (mm)	Maximum (mm)	Range(mm)
BL	0	35	35
ML	0	45	45
HL	0	8	8
TPW	0	70	70

There are 120 grey levels for each LPW parameter, therefore the scales and offsets for the parameters are the following:

Parameter	Scale	Offset
BL	35/119	-8*35/119
ML	45/119	-8*45/119
HL	8/119	-8*8/119
TPW	70/119	-8*70/119

 These parameters appear in the HDF files.

 

• Description of LPW output format

The format will be SAFNWC, the content of the LPW is HDF-5 format:

• 8 bits to show the BL:     BL clear air pixels.     Configurable IR channel cloudy pixels.

• 8 bits to show the ML:    ML clear air pixels.   Configurable IR channel cloudy pixels.

• 8 bits to show the HL:     HL clear air pixels.     Configurable IR channel cloudy pixels.

• 8 bits to show the TPWcontrol: TPWcontrol clear air pixels. Configurable IR channel cloudy pixels.

• 1 to describe clear sky flag:

   0     Air clear pixel. LPW is calculated
   1     Cloudy pixel. IR value is given

• 2 bits for Number of out of range radiances: Number of bands whose radiances are out of a range specified in the configuration file. The ranges of the radiances for the seven different channels are specified in the configuration file. Possible values are:

  0     All bands within the range
  1     One band out of range
  2     Between 2 and  4 bands out of range
  3     Between 5 and  7 bands out of range

• 2 bits to label number of cloudy pixels in the biggest window: This label presents information about the neighbours of the processed pixel depending on the cloudy pixels in the biggest window used to perform the spatial pre-processing of bands. The window width is a configurable parameter (normally: none, 3x3 or 5x5 pixels). Possible values are:

  0     [0, 1/8] cloudy pixels.   (If window 3x3: [0, 1] pixel     &     if window 5x5: [0, 3] pixels)
  1     (1/8, 1/4] cloudy pixels (If window 3x3:    [2] pixels     &     if window 5x5: [4, 6] pixels)
  2    (1/4, 1/2] cloudy pixels  (If window 3x3: [3, 4] pixels    &     if window 5x5: [7, 12] pixels)
  3    (1/2, 1] cloudy pixels     (If window 3x3: [5, 8] pixels    &     if window 5x5: [13, 24] pixels)

• 2 bits the quality flag: This flag checks the difference between the TPW control and the sum of the three layers in each pixel. Depend of the % of the differences.

  0     High quality [0,10]%
  1     Medium quality (10,30]%
  2     Low quality (30,50]%
  3     Bad quality (<50%) or cloudy pixels or MSG zenith angle > threshold.

• 2 bits the consistence sea flag This flag only is applied on sea pixels. It shows the difference between TPW control and PGE06(TPW). Depend of the % of the differences.

  0     High consistence (0,10]%
  1     Medium consistence (10,30]%
  2     Low consistence (30,50]%
  3     Bad quality (<50 %) or cloudy pixels or PGE06 not available or MSG zenith angle > threshold.

   

6.- Example of LPW visualisation

We present here as an example images corresponding to the same study case (13-September-2003,00Z) that we presented in the last EUMETSAT user conference (Martínez et al., 2003). The Layer Precipitable Water images have been derived with version 1.0. The images are in polar Stereographic projection with 13 km resolution, centred at (44ºN, 0ºW).
 

Click on thumbnail for full-sized these images
Figure 1: (Layer Precipitable Water output)(1) Boundary Layer (2) Middle Layer (3)  High Layer  (4) TPW control

• Description of the example

The products of LPW have been generated from SEVIRI using the statistical retrieval method coded in version 1.0 SAFNWC package. The products displayed are the four layers of atmospheric precipitable water vapour values: boundary Layer precipitable water vapor value (1), the middle layer precipitable water vapor value (2), the high layer precipitable water vapor value (3), and the total atmospheric precipitable water vapor value(4). (Precipitable water is the amount of liquid water, in millimeters, if all the atmospheric water vapor in the column is condensed). The value is color-coded with browns being the driest and reds being the moistest. The color-coded is similar to the CIMSS TPW color-coded. Clouds are represented as a grey color (see the color bar and the mm values at the bottom of the images).