# Presentation: 087725

## NCEP Regional ReanalysisFedor Mesinger1, Geoff DiMego2, Eugenia Kalnay3, Perry Shafran4, Wesley Ebisuzaki5, Dusan Jovic4, Jack Woollen4, Kenneth Mitchell2, Eric Rogers2, Michael Ek1, Yun Fan6, Robert Grumbine2, Wayne Higgins5, Hong Li3, Ying Lin2,  Geoff Manikin2, David Parrish2, and Wei Shi6 1NCEP/EMC and UCAR, 2NCEP/EMC, 3Univ. of MD, 4NCEP/EMC and SAIC/GSO, 5NCEP/CPC, 6NCEP/CPC and RSISNARR Users Workshop, San Diego, 11 January 2005

## Goal/ Motivation

- Create long-term set of consistent climate data on a regional scale for the North American domain
- Use of the GR to drive the Regional Reanalysis (RR) system
- Superior to NCEP/NCAR Global Reanalysis (GR) due to:
- Use of a regional model (the Eta Model);
- Advances in modeling and data assimilation since 1995,  especially:
- Precipitation assimilation
- Direct assimilation of radiances
- Land-surface model updates

## E.g.:

_**ETA / NOAH LAND-SURFACE MODEL UPGRADES**_**: ****24 July 01**

**- assimilation of hourly precipitation**
**  **-- hourly 4-km radar/gage analysis (Stage IV)
**- cold season processes**(*Koren et al 1999*)  -- patchy snow cover  -- frozen soil (new state variable)  -- snow density (new state variable)  **- bare soil evaporation refinements**
-   -- parameterize upper sfc crust cap on evap
**- soil heat flux**  -- new soil thermal conductivity
-     (*Peters-Lidard et al 1998*)  -- under snowpack (*Lunardini, 1981*)  -- vegetation reduction of thermal cond.    (*Peters-Lidard et al 1997*)
**- surface characterization****  **-- maximum snow albedo database
-     (*Robinson & Kukla 1985*)  -- dynamic thermal roughness length
-     refinements
**- vegetation**
-   -- deeper rooting depth in forests
-   -- canopy resistance refinements

- NOAH LSM tested in various land-model
- intercomparison projects, e.g., GSWP, PILPS
-  2a, 2c, 2d, 2e, Rhone, and (near-future) DMIP.

## NARR Support

**NOAA Office of Global Programs (OGP)**
- GAPP: GEWEX America Prediction Project
**NCEP (EMC and CPC)**
- Computational Resources / IT Staff
- Staff scientific and technical expertise
**Science Steering Committee**
- J. Roads, Chair
**NCDC, UCAR**
- public access to NARR database

## System Design

- Fully cycled 3-hr EDAS (Eta Data Assimilation System)
- Lateral boundary conditions supplied by Global Reanalysis 2
- Resolution: 32-km, 45 layers
- RR time period: 1979-2003 (continued later in near-real time, as in CDAS)
- •	Free forecasts out to 72 hr every 2.5 days, using GR2 forecast boundary conditions (“reforecasts”)

##       180 km                                 32 km

-       180 km                                 32 km

- (GR equivalent)

## The Domain:

## Data Used in Global Reanalysis and Regional Reanalysis

| Dataset | Details | Source |
| --- | --- | --- |
| Radiosondes | Temperature, winds, moisture | NCEP/NCAR Global Reanalysis (GR) |
| Dropsondes | Same as above | GR |
| Pibals | Wind | GR |
| Aircraft | Temp. and wind | GR |
| Surface | Pressure | GR |
| Geostationary satellite | Cloud drift winds | GR |

## Data Added or Improved Upon for Regional Reanalysis

| Dataset | Details | Source |
| --- | --- | --- |
| Precipitation | CONUS (with PRISM), Mexico, Canada, CMAP over oceans (<42.5°N) | NCEP/CPC,Canada, Mexico |
| TOVS-1B radiances | Temperature, precipitable water over ocean | NESDIS |
| NCEP Surface | Wind, moisture | GR |
| TDL Surface | Pressure, wind, moisture | NCAR |
| COADS (ships/buoys) | Pressure, wind, moisture | NCEP/EMC |
| Air Force Snow | Snow depth | COLA and NCEP/EMC |
| SST | 1-degree Reynolds, with Great Lakes SSTs | NCEP/EMC, GLERL |
| Sea and lake ice | Contains data on Canadian lakes, Great Lakes | NCEP/EMC, GLERL,
 Ice Services Canada |
| Tropical cyclones | Locations used for blocking of CMAP Precipitation | Lawrence Livermore National Laboratory |

## Climatologies

| Dataset | Details | Source |
| --- | --- | --- |
| Green vegetation fraction | Used to initialize the vegetation, daily | GR |
| Baseline snow-free albedo | Used to initialize the albedo, seasonal | GR |

## Fixed Fields

- Vegetation type [index, 1-13]
- Soil type [index, 1-9]
- Surface slope type [index]
- Snow-free albedo [%]
- Maximum snow albedo [%]
- Surface roughness [m]
- Soil column bottom temp. [K]
- Number of root zone soil layers [non-dim]

## It was a lot of work /                        There were problems to solve:

- Acquire, or create (improved reprocessing) a variety of datasets, e.g.,
- TOVS-1B data for 1979-2002 time period;
- Air Force snow depth dataset;
- High (RR) resolution sea and lake ice, and sea and lake surface temperatures;
- Precipitation data (Mexico, Canada, CMAP), processing (PRISM/”Mountain Mapper”), disaggregation to hourly;
- Improvements/ refinements of various components of EDAS, e.g.,
- Precipitation assimilation (additional regions, blending, ...)
- Upgrade 3D-Var to run using the satellite bias correction for numerous additional satellites;

## ... a lot of work (cont’d)

- •	Presented/ discussed RR-related work at numerous conferences
- AMS Annual Meeting in Orlando, FL, in January 2002
- GEWEX Mississippi River Hydrology & Climate Conference in New Orleans, LA, in May 2002
- AMS Annual Meeting in Long Beach, CA, in February 2003
- Climate Diagnostics Workshop in Reno, NV, in October 2003
- AMS Annual Meeting in Seattle, WA, in January 2004
- Numerous contacts with potential users;
- Completed 24 years of RR production in just over 3 months, running 4 streams, on all of the previously mainframe NCEP IBM ASP

## We’ve learned a few things along the way, e.g.,

- •	Assimilation of 2-m land surface station air temperatures, with the system we have:
- Harmful in the sense of making the first guess considerably worse, throughout the troposphere (!)

## “Full” Assimilation

- No 2-m T assimilation

- T

- T

- Wind

- Wind

- July 1998

## More things learned (disappointments ?):

- •  No benefit identified from direct assimilation of radiances 	(the Eta Model top, 25 mb, too low ?);
- •  While RR’s 3 vs GR’s 6 hr analysis frequency resulted in
- higher RR time resolution, it did not increase accuracy (3 hr first guess not closer to observations than the 6 hr first guess would have been)

## More things learned (disappointments ?):

- •  No benefit identified from direct assimilation of radiances 	(the Eta Model top, 25 mb, too low ?);
- •  While RR’s 3 vs GR’s 6 hr analysis frequency resulted in
- higher RR time resolution, it did not increase accuracy (3 hr first guess not closer to observations than the 6 hr first guess would have been)
- But what matters most is the result, thus

## Results: Upper-Air

- Compared both GR and RR against fits to raobs;
- Root-mean-square (RMS) analysis fits significantly better for temperatures and winds for RR;
- Wind improvement greatest in the upper troposphere (!), especially in winter;
- First guess (3-hr forecast, pre-3D-Var) temperatures not always as favorable for RR compared to GR, winds improved;
- Relative humidity improved for RR for both analysis and first guess (not shown)

## Analysis, RMS

- January

- July

- Wind

- Wind

- T

- T

- 25-yr averages

## Analysis, Bias

- January

- July

- Wind

- Wind

- T

- T

## Q, Analysis 

- January

- July

- RMS

- Bias

- RMS

- Bias

## First Guess, RMS

- January

- July

- Wind

- Wind

- T

- T

## First Guess, Bias

- January

- July

- Wind

- Wind

- T

- T

## Q, First Guess

- January

- July

- RMS

- Bias

- RMS

- Bias

## Results: Near Surface

- •	Analysis (not shown): Comparison against the GR2 not possible.  RMS and biases look reasonable;
- •	First guess, 1997:  Surface temperature RMS improved both in winter and in summer;
- •	1998: Surface temperatures RMS favorable for RR in both winter and summer.  RR biases closer to zero and little diurnal variation problem in summer;
- 10-m winds: RMS improved over GR greatly in winter, a little bit in summer;
- Slow wind biases improved in RR: just a little in summer, visibly in winter

## July 1997

- January 1997

- 2-m  Temperature, 1st Guess

- RMS

- Bias

- RMS

- Bias

## July 1998

- January 1998

- 2-m  Temperature, 1st Guess

- RMS

- Bias

- RMS

- Bias

## 10-m Wind, 1st Guess

- January 1997

- July 1997

- RMS

- Bias

- RMS

- Bias

## 10-m Wind, 1st Guess

- January 1998

- July 1998

- RMS

- Bias

- RMS

- Bias

## Results: Precipitation

- Several sources of precipitation
- ConUS data with PRISM (Mountain Mapper) to improve orographic effects
- Canada
- Mexico
- CMAP (combination of satellite and gauge data) over oceans; CMAP is blocked:
-  Near central areas of hurricanes (7.5 by 7.5 deg)
-  Observed precipitation > 100 mm/day
-  A 15-degree “blending belt” between 27.5 and 42.5 N, with no CMAP north of 42.5 N

## Results: Precipitation (cont’d)

- Precipitation observations used to prescribe the latent heat profile in Eta model;
- Model uses given latent heat profile to simulate precipitation;
- Resulting precipitation pattern looks very much like the analyzed precipitation pattern in both summer and winter

## January 1997 Precipitation Results

## July 1997 Precipitation Results

## (A strong El Niño)

## July 1998 Precipitation Results

## Flood of 1993 vs Drought of 1988

- Shows impact of summer precipitation;
- Extreme years compared;
- Monthly mean June/July precip of 1993 minus monthly mean June/July precip of 1988;
- Shows the success of precipitation assimilation even for drought and flood years

## Flood of 1993 vs Drought of 1988

## R-CDAS

- Running in near real-time (similar to the current CDAS that runs with the Global Reanalysis system),
- Precipitation:
- CMORPH, a 1/2-hr high-resolution satellite-based dataset, replaces CMAP over oceans, used for 2003 and beyond;
- No Canadian precipitation; using EDAS precipitation instead
- Mexico and US (N-LDAS) precipitation in real-time;
- Data available after ~ 21 hr delay
- System, including data ingest script, created and tested
- Currently running ~ 4 days behind real time, responsibility of NCEP/CPC

## Our output files produced  (about 75 Tb total):

- Six sets of reanalysis files, at 3 hr intervals each:
- • “Restart files”: on “native” (model) grid, all variables needed to restart the model;
- •  On the “native” (E) grid, but interpolated to p levels:
- 		Analysis files,		First guess files;
- •  Interpolated in addition to Grid 221 (Lambert, “AWIPS”):
- 		Analysis files,		First guess files;
- •  “Merged” files (handout): Grid 221, a selection of (mostly) analysis variables, enhanced by, e.g., fluxes (“fcst”, from 1st guess), 3-h averages, and 3-hr accumulations (4 types of time validity)

## Output files produced, cont’d

- •  Input observations;
- •  Climatological and fixed fields used;
- •  Various plots;
- •  Special subsets:
- 	-  Land-surface subset - fields related to
- 		land-surface/hydrology modeling
- 	-  Climatological subset (CPC)
- •  “Free forecasts” (“reforecasts”): 3-day forecast 	every 2.5 days

## Output files produced, cont’d;                              more on the merged files (handout):

- •  Two main sets of the merged files produced; so-called merged-a and merged-b;
- •  Temporal frequency 3 hour (8 files per day, archived in a single tar file)
- •  Merged-a daily tar file is about 450 Mb, merged-b about 75 Mb;
- •  Horizontal grid standard NCEP grid 221 (Lambert conformal, 32 km).  Note that a sliver of 221 is sticking out of the RR domain/ variables undefined;
- •  In vertical: pressure levels, "hybrid levels” (lowest eta mid-layers), soil layers, boundary layer etc.;
- •  Fields in merged-a files mostly analysis; some 3-hr accumulations, 3-hr averages, and 3-hr forecasts;
- •  When the same filed is output as analysis and 3-hr forecast, 3-hr forecast is placed in merged-b file (grads problem);
- •  Some fields are interpolated to 221 using nearest-neighbor algorithm (labeled with asterisk in column three in the list)

## Archiving

- Two or more archiving centers/ two have some of the data:
- National Climatic Data Center (NCDC), merged files for all 25 years available;
- National Center for Atmospheric Research (NCAR), merged files being downloaded;
- San Diego Supercomputing Center (SDSC): contacts in progress;
- Perhaps also University of Maryland

## NARR on the Web

- _[http://wwwt.emc.ncep.noaa.gov/mmb/rreanl/index.html](http://wwwt.emc.ncep.noaa.gov/mmb/rreanl/index.html)_

## Summary (results)

- Long-term, consistent, high-resolution climate dataset for North America a significant improvement over earlier global reanalyses;
- Near-surface temperatures and winds closer to observations;
- Winds in the upper troposphere, in addition to land-surface hydrology over the ConUS, maybe the strongest areas of improvement;
- Improvements greater in winter than in summer;

## But also a few problems/ weaknesses,               that we mostly became aware of after the fact: 

- • Excessive strength of the summertime Gulf of California low-level jet,  :-(   NAME.  Not understood at this point;
- •  Precipitation over Canada: the number of gauge observations available may have been insufficient to do better than the model might have done;
- Precipitation over northern Atlantic not simulated well;
- Atlantic hurricanes apparently not an RR strong feature;
- •  . . . (TBD)

##   Even so, we are confident the NARR will for some years to come be the best/ most accurate North American weather and climate dataset, for numerous research and application purposes;

USER input most welcome !!  
Ready to help to the extent we are able to.

- USER input most welcome !!
- Ready to help to the extent we are able to.