1 / 94

An Overview of the NCEP Eta Model

An Overview of the NCEP Eta Model. COMET/UCAR SOO Symposium on NWP. 28 March 2000. Presented by Thomas Black. EMC/Mesoscale Modeling Branch. EDAS slides by Eric Rogers. Outline. Brief model description Eta Data Assimilation System Physics Examples of products/statistics Future. Domain.

mahlah
Download Presentation

An Overview of the NCEP Eta Model

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. An Overview of the NCEP Eta Model COMET/UCAR SOO Symposium on NWP 28 March 2000 Presented by Thomas Black EMC/Mesoscale Modeling Branch EDAS slides by Eric Rogers

  2. Outline • Brief model description • Eta Data Assimilation System • Physics • Examples of products/statistics • Future

  3. Domain • Semi-staggered Arakawa E grid • 32km horizontal resolution • 45 vertical eta layers • Silhouette step topography

  4. Eta Domains Past & Present

  5. 32km Domain Topography w/ Water Points

  6. 32km CONUS Topography w/ Water Points

  7. ground ground Sigma and Eta Coordinates P  P  MSL is small At point P: At point P: is small

  8. Eta Coordinate Reference heights and temperatures taken from the standard atmosphere P=  LM-3 PMSL LM-3 LM-2 P=  LM-2 PMSL Z= ZREF at P = LM-2  PSMSL P=  LM-1 PMSL LM-1 Z= ZREF at P = LM-1  PSMSL Z=0 LM=1 P=PMSL Mean Sea Level

  9. Eta Model 45-Layer Distribution 25 hPa 27 hPa 29 29 27 26 26 24 23 250 hPa 24 26 30 31 32 32 32 33 500 hPa 33 33 33 33 32 31 29 700 hPa 27 23 21 19 18 18 17 850 hPa 16 14 13 13 13 13 12 12 12 1000 hPa 12 11 8 6 5 2

  10. H H H V V d d H H V V V H H H V V H H V V V H H H V V The Semi-Staggered E Grid  H mass point  V velocity point constant transformed latitude constant transformed longitude

  11. Eta Data Assimilation System GOAL : Produce best possible initial conditions for the Eta Model forecast* • KEY COMPONENTS - State of the art analysis (variational) - Consistency between assimilating and forecast model (resolution, physics, dynamics) - Intelligent selection and use of observations * NOT necessarily the same as fitting all the observations exactly

  12. What is 3D-VAR? • An analysis technique that attempts to minimize analysis error • Takes background (forecast) and observation error into account • Variational method allows use of “non-traditional” data sources, such as GOES precipitable water

  13. - Stream function - Temperature - Potential function - Specific humidity - Surface pressure - Geopotential height ETA 3DVAR ANALYSIS (Parrish et al. 1996 NWP Preprint Volume) • Loosely patterned after NCEP global SSI analysis • Analysis variables: • More adaptable than OI for using new data types (e.g., NEXRAD radial velocities used in Eta-10 runs during 1996 Olympics)

  14. 3D-VAR vs. OI

  15. EDAS Original Configuration Eta-48 fcst 00Z/12Z Eta-29 fcst 03Z/15Z

  16. WHY DO CYCLING? • Initial conditions more consistent with forecast model • Less spinup of divergence, cloud, precipitation, and TKE • More accurate representation of soil moisture

  17. Observations Used By ETA 3DVAR • Upper air data - Rawinsonde height/temperature/wind/moisture - Dropwindsondes - Wind Profilers - NESDIS thickness retrievals from polar orbiting satellites (oceans only) - VAD winds from NEXRAD - Aircraft (conventional and ACARS) winds/temps - Satellite cloud drift winds - SSM/I and GOES precipitable water retrievals - Synthetic tropical cyclone data • Surface data - Surface land wind/temperature/moisture - Ships and buoys - SSM/I oceanic surface winds

  18. DATA QUALITY CONTROL • CQC: Complex QC of raob height/temps (baseline, hydrostatic, lapse rate, radiation correction, etc.) • ACQC : Quality control of conventional aircraft data (remove duplicates, track checks, create “superobs”) • SDMEDIT : NCEP Senior Duty Meteorologist can flag all or parts of suspect raobs • 3DVAR : Analysis performs gross check vs. first guess: - Temperature : +/- 15oC - Wind :+/- 25 ms-1 - RH : +/- 90% - Precipitable water : +/- 12 g/kg - Height : +/- 100 m

  19. Use of Surface Data: Eta OI vs. Eta 3DVAR Eta OI Analysis Eta 3DVAR Analysis

  20. New 3DVAR tested in July 1998 and showed improved fit to surface and raobs (especially moisture) • Re-tuned 3DVAR implemented on 3 November 1998 • We thought everything was OK….. BUT……..

  21. 24-H ACCUMULATED PRECIPITATION EQUITABLE THREAT SCORES: ALL FCSTS Solid = Eta Short Dash = NGM Long Dash = AVN/MRF 12/1/97 - 2/28/98 12/1/98 - 2/28/99 10-15% drop in Eta skill between 1997-98 and 1998-99

  22. Persistent synoptic error in Eta-32 during winter of 98-99: weaker and faster Eastern Pacific troughs/cyclones than observed Example : 48-h Forecasts valid 1200 UTC 17 March 1999

  23. PROBLEM 1: November 98 change degraded mass/wind balance in 3DVAR • If mass / wind balance well-behaved, positive height correction is coincident with center of anticyclonic wind correction 850 mb ANL-GUESS height/wind 80KM EDAS valid 00Z 3/15/99 • Note 10 degree longitude displacement between centers of wind and height correction • Problem is most severe in regions and at analysis times without widespread raob data but with large amounts of wind or mass only data (e.g., satellite winds)

  24. SOLUTION : Improve geostrophic coupling of mass/wind analysis corrections in 3DVAR (5/99) Original 3DVAR analysis Improved 3DVAR analysis Note: Improved height/wind coupling near Aleutians

  25. PROBLEM 2: Horizontal/vertical correlations too narrow : observation had VERY limited impact on analysis away from its level SOLUTION: Expand the influence of the observations (5/99) • One observation test : Insert one height observation 10 m greater than first guess at 200, 500, 900 mb and measure impact in horizontal/vertical

  26. Original 32-km 3DVAR Improved 32-km 3DVAR 200 mb 900 mb

  27. Performance of new 3DVAR : 3 December 1998 to 16 January 1999 test at 80 km resolution 24-h accumulated precipitation threat scores: All forecasts Dashed = Modified 3DVAR Solid = Operational 3DVAR Equitable Threat Score Threshold (in)

  28. The Forecast Model Split Explicit Integration: Dynamics • Fundamental prognostic variables • T, u, v, q, Psfc, TKE, cloud water/ice • Inertial gravity wave adjustment • forward-backward scheme (Dt=90s) • Vertical advection • Euler-backward scheme • centered in space • piecewise linear for q

  29. Split Explicit Integration: Dynamics • Horizontal advection • modified Euler-backward scheme • Janjic advection in space • conservative, (nearly) shape-preserving scheme for H20 • upstream advection near boundaries

  30. Split Explicit Integration: Physics • Betts-Miller-Janjic convection • Mellor-Yamada level 2.5 turbulent exchange • GFDL radiation • explicit cloud water/ice prediction • 4-layer NOAH land surface package • D2 horizontal diffusion

  31. One-wayBoundary Conditions • 3-hour tendencies • 6 hour old AVN forecast used

  32. vertical advection 1 5 Runstream Schematic of Eta Model Integration gridscale cloud gridscale precip convection turbulence horizontal advection inertial grv wave adjustment 1 2 3 4 5 6 7 8 9 1 0 1 1 1 2 1 3 1 4 1 6 timestep t = 90 s Radiative temperature tendency updates Shortwave: 40 timesteps (1 hour) Longwave: 80 timesteps (2 hours)

  33. The Betts-Miller-Janji Convection Scheme in the Eta Model References: Betts, 1986 (QJRMS) Betts and Miller, 1986 (QJRMS) Janji, 1994 (MWR) • Deep (precipitating) convection • Temperature reference profile • Moisture reference profile • Convective adjustment • Modification for “precipitation efficiency” • Shallow (non-precipitating) convection • Temperature reference profile • Moisture reference profile

  34. Find the Deep Convective Clouds 1. For all ‘parcels’ within 0.2xPsfc mb of the ground find Psat andES 2. At each point, select parcel with the maximum ES 3. Given Psat, choose cloud base as the model level just below it 4. Adjust cloud base if needed: (a) at least 25 mb above middle of lowest layer (b) at least one model layer above lowest layer 5. Compute Tmad above cloud base using ES and P in lookup tables 6. Set cloud top at highest level where Tmad<T-T (currently T=0) 7. Gather all clouds at least 0.2xPsfc mb deep

  35. Betts-Miller Reference Temperature Profile

  36. 1. Define ‘deficit from saturation pressure’ (DSP) for cloud bottom, freezing level, and cloud top (larger DSPdrier state) Construction of 1st Guess Humidity Reference Profile For Deep Convection 2. Linearly interpolate DSP’s for values between these 3 levels 3. Define the reference humidity profile as qsat in each layer

  37. The Enthalpy Correction Modify the profiles to ensure enthalpy in the cloud column is conserved during adjustment Initially: Corrections: Currently the above procedure is repeated two times

  38. 1. Relax model profiles of T and q toward the reference profiles where relaxation time  equals 2400 s Final Deep Convective Adjustment 2. Convective rainfall amount is 3. At any point, deep adjustment is ignored and a “swap” to shallow convection occurs if: (a) S < 0 (b) precipitation is negative

  39. Deep Convective Adjustment of Temperature Upward Transport of Heat cloud top REFERENCE TEMPERATURE AMBIENT TEMPERATURE cloud base

  40. Modification for ‘Precipitation Efficiency’  Examination of Eta integrations shows: As convective precipitation increases, entropy changes decrease. So defineprecipitation efficiencyas the ratio: Numerator: Q arising from entropy change Denominator: Q arising from precipitation (H=0)

  41. Thus, larger E less mature system smaller E more mature system USE E TO MODERATE HEAVY RAIN IN LONG-LIVED MATURE SYSTEMS (A) Modify the humidity reference profile (B) Modify the relaxation time 

  42. Humidity Reference Profile Limits q cloud top HUMIDITY REFERENCE PROFILE E=0.2 E=1 MOISTPROFILELIMIT p DRY PROFILE LIMIT cloud base DSP’s vary between DRY (fast) and MOIST (slow) limits

  43. Top -1875 Pa Freezing -5875 Pa Bottom -3875 Pa Values of DSP Limits Dry limits MOIST DSP limits equal 0.85 times DRY limits. Att=0, DSP’s are set to the DRY limits.

  44. Modification of Relaxation Time Multiply the standard change due to adjustment by some quantity F which is a function of the precipitation efficiency OR For simplicity assume F is linear. Then empirically: 0.7 <F< 1.0 for 0.2 <E< 1.0

  45. Find the Shallow Convective Clouds 1. Find the tops of the “swapped” clouds (a) set a preliminary top at pbot - 0.2xPsfc mb (pbot > 450 mb) (b) reset top to level where maximum (RH)/p occurs 2. Gather all clouds that are: (a) greater than 10 mb deep (b) less than 0.2xPsfc mb deep (c) at least two model layers deep

  46. (cpT p) = 0 Correct Tref assuming Construction of Temperature Reference Profile For Shallow Convection cloud top Mixing Line REFERENCE TEMPERATURE cloud base

More Related