Skip to main content
Have a personal or library account? Click to login
Ensemble-based observation impact estimates using the NCEP GFS Cover

Ensemble-based observation impact estimates using the NCEP GFS

Open Access
|Dec 2013

Figures & Tables

Table 1

Observation types assimilated in the experiment

Type of data Description 800 hPa~400–800 hPa125–400 hPa ~125 hPaTotal (1000 s)Aircraftu, v, and T observations from the aircrafts37.253.981.70.0172.7RadiosondeRadiosonde observations (u, v, T, q and Ps)9.513.213.818.454.8Satellite_WindAtmospheric motion vectors (u and v) from geostationary satellites31.21.862.80.195.8GPSROGPS radio occultation1.49.417.867.095.5Land-surfacePs observations from land-surface stations53.41.00.00.054.4Marine-surfaceSurface u, v, T, q and Ps observations from the buoys and ships23.10.00.00.023.1MODIS_WindAtmospheric motion vectors (u and v) from MODIS0.826.19.20.036.1ASCAT_Windu and v observations from ASCAT scatterometer over ocean16.70.00.00.016.7PIBALu and v observations from pilot balloons0.80.60.40.21.9NEXRAD_Windu and v observations from the NEXRAD (radar)6.56.60.40.013.5Profiler_Windu and v observations from the wind profilers1.25.14.90.611.8DropsondeFlight-level reconnaissance and dropsonde (u, v, T and q)0.10.10.10.00.2WINDSAT_Windu and v observations from ASCAT scatterometer over ocean (super observation)0.80.00.00.00.8TCVitalPseudo surface pressure observations at tropical cyclone storm centre0.00060.00.00.00.0006OzoneOzone retrievals from satellite radiances0.91.81.89.914.4AMSUASatellite microwave sounder radiances (from five satellites)144.248.6123.8269.2585.8IASISatellite infrared hyperspectral sounder radiances140.5116.6258.5627.41143.0Aqua_AIRSSatellite infrared hyperspectral sounder radiances59.890.0195.1285.4630.2ATMSSatellite microwave sounder radiances (from Suomi-NPP)41.426.129.352.8149.6HIRSSatellite infrared radiances (from two satellites)13.018.040.847.1118.9MHSSatellite microwave sounder radiances (from three satellites)31.729.35.80.066.8GOESGOES infrared sounder radiances (GOES13 and 15)10.38.613.97.039.8SEVIRISEVIRI clear sky radiances0.22.73.20.06.1

[i] The third to sixth columns show the average number of observations assimilated on each vertical layer in one analysis (in thousands). Peak pressure of weighting function is used as vertical position of satellite radiance observation. u, v, T, q and Ps represents u and v wind components, temperature, specific humidity and surface pressure, respectively.

Fig. 1

Time series of the total forecast error reduction of each estimate (unit: J kg−1). Black, red and blue lines show the actual forecast error reduction verified against the own analysis, estimated error reduction from the EnKF-based method with fixed localisation (fixed) and with moving localisation (advected). Numbers on upper left corner show the correlation and RMSE of each estimate to the actual forecast error reduction.

Fig. 2

Estimated average 24-hour forecast error reduction contributed from each observation types (moist total energy, J kg−1). (a) represents the total error reduction and (b) represents error reduction per observation.

Fig. 3

Same as Fig. 2 but with the dry total energy norm (J kg−1)

Fig. 4

Estimated AIRS satellite radiance observation impacts classified by channel with the moist total energy (red for channels sensitive below 300 hPa and magenta above 300 hPa, J kg−1), and only the moist term of the total energy (blue for channels sensitive below 300 hPa and aqua above 300 hPa, J kg−1). Average estimated forecast error reduction from a single observation is shown. Vertical bars represent the 95% confidence interval of the average values.

Fig. 5

Estimated average observation impacts of a) radiosonde and b) aircraft classified by observed level (moist total energy, J kg−1). Average estimated forecast error reduction by a single observation is shown.

Fig. 6

Average number of assimilated radiosonde observations (a) from 250 to 125 hPa, (b) from 800 to 600 hPa and aircraft observations (c) from 250 to 125 hPa and (d) from 800 to 600 hPa in each 5°×5° area for one analysis.

Fig. 7

Average impact (moist total energy, J kg−1) of a single radiosonde profile from the fixed land stations. Only the stations that have more than 20 profiles in the period are shown. Numbers 4220 (Egedesminde) and 4270 (Narssassuaq) indicate the location of the stations shown in Fig. 8.

Fig. 8

Comparison of the radiosonde observations from Narssassuaq (red line, 4270) and Egedesminde (blue line, 4220) showing (a) average impacts (J kg−1) of each observation element (solid: temperature, dashed: winds, dotted: humidity) on each pressure level by one profile and observation departure statistics (dashed: bias, solid: standard deviation) of (b) temperature (K) and (c) wind speed (m s−1).

Table 2

Width of longitude used for the local area on each latitude band

LatitudeWidth for longitude (°)60N–90N6050N–80N4040N–70N3030N–60N2020N–50N2010N–40N150N–30N1510S–20N15

[i] Only the values for Northern hemisphere are shown since it is symmetric in the Southern Hemisphere.

Table 3

List of local 24-hour forecast failure cases (initial time from 00 UTC, 8 January 2012, to 18 UTC 7 February 2012)

InitialAreaSizeRateNDenied observation (denied number/total number)Change (estimate)06 UTC JAN 1250N–80N 145E–175W1.991.365AMSUA ch4, 5, 6 (2735/125 063)−8.7% (−19.5%)00 UTC JAN 1630N–60N 20W–02.711.356GPSRO 600–950 hPa (50/4918)−0.6% (−4.3%)18 UTC JAN 2730S–0 105E–120E2.401.211AIRS (19 908/670 041)−0.2% (−6.0%)00 UTC JAN 3070S–40S 165E–165W2.001.256AMSUA ch1, 3, 4, 5, 15 (3822/164 934)−4.7% (−12.8%)06 UTC FEB 250N–80N 150W–110W3.011.225GPSRO 250–400 hPa, 600–850 hPa (407/13 092)−11.7% (−8.7%)06 UTC FEB 430N–60N 150W–130W1.811.263IASI (57 950/1 177 256), HIRS ch3, 4, 9, 11, 12, 14, 15 (785/73 419), Aircraft 950 hPa~, 125–600 hPa (5794/100 896)−25.5% (−81.6%)18 UTC FEB 660N–90N 40E–100E1.711.382MODIS_Wind (10 970/43 452)−28.4% (−47.7%)

[i] The third, fourth and fifth columns show the forecast error size normalized with the time-averaged error, the rate of the error size compared with the error of 30-hour forecast from the previous analysis, and the number of areas that meet the criteria. Sixth column shows the denied observation type based on the observation impact estimates, number of denied observations and total number of observations with this observation type assimilated in the original analysis. Finally, seventh column shows the change of the 24-hour local forecast error for observation denial experiment and corresponding estimated forecast error change in the parenthesis. The case shown in bold font is the example shown in Fig. 9.

Fig. 9

Twenty-four hour forecast error of 500 hPa geopotential height (unit: m, 18 UTC 6 February 2012 initial) from original analysis (left) and forecast change due to the removal of the MODIS polar wind observations in the data-denial experiment (middle: actual change and right: projection on the ensemble perturbations). Black contours show the analysis. Magenta cones show the target area of the observation impact estimate.

Language: English
Page range: 20038 - 20038
Submitted on: Nov 7, 2012
Accepted on: Aug 9, 2013
Published on: Dec 1, 2013
Published by: Stockholm University Press
In partnership with: Paradigm Publishing Services

© 2013 Yoichiro Ota, John C. Derber, Eugenia Kalnay, Takemasa Miyoshi, published by Stockholm University Press
This work is licensed under the Creative Commons Attribution 4.0 License.