390 likes | 538 Views
Introduction to NWP. What is NWP?. Numerical Weather Prediction: Predicting the future state of the atmosphere by use of the mathematical equations that describe it. Equations of motion. 1. Momentum. Conservation Laws. Mass:. Energy:.
E N D
What is NWP? • Numerical Weather Prediction: Predicting the future state of the atmosphere by use of the mathematical equations that describe it.
Equations of motion 1. Momentum Conservation Laws Mass: Energy: Hydrostatic: eliminate acoustic, high-frequency buoyancy waves
“Physics” equations • model processes in atmosphere and at earth’s surface.
Equations • Equations can only be solved in modifed form, on super-computer. • First successful NWP was in 1950 with simplified (barotropic) model run on ENIAC.
Discretization: of fields • bring irregularly spaced data to regular grid
xi, at time t-∆t xi, at time t+∆t Discretization: of equations Simple example: 1D advection equation: approximated by differences in space and time:
GEM grids Global Regional LAM (east)
Discretization: Vertical Levels GEM: 80 levels LAM: 58 levels Terrain-following “eta” coordinate becomes level with increasing height
Parameterization e.g.: • convective cloud processes • formation of cloud droplets and precipitation • solar and infrared radiative transfer • surface-air exchanges of energy, momentum, and moisture
Use J(xi-1) to find xi Data assimilation concept Data Acquisition Error Statistics (Observation and Forecast) Data Quality Control Analysis (Spatial QC) Cost function J(x,..) First Guess xa MIN J(x) NWP model
Observations assimilated at CMC Type Variables Thinning radiosonde/dropsonde U, V, T, (T-Td), ps 28 levels Surface report T, (T-Td), ps, (U, V over water) 1 report / 6h Aircraft (BUFR, AIREP, AMDAR, ADS) U, V, T 1o x 1o x 50 hPa per time step Ocean Land AMSU-A 3-10 6-10 AMSU-B / MHS 2-5 3-4 ATOVS NOAA 15-16-17-18, AQUA 250 km x 250 km per time step Water vapor channel GOES 11-12 IM3 (6.7 m) 2o x 2o 3-hourly U,V (IR, WV, VI channels) 1.5o x 1.5o 11 layers, per time step AMV’s (METEOSAT 5-8, GOES 11-12, MTSAT-1R) ~180 km boxes 11 layers, per time step MODIS polar winds (Aqua, Terra) U,V Profiler (NOAA Network) U,V (750 m) Vertical hourly
e.g. satellite winds received at CMC(6 hr period, 1 single level wind vector per obs, ~380k obs, from GTS)
e.g. radiosonde data received at CMC(6 hr period, multi-level U,V,T,Es, ~600 stations, from GTS and MSC)
e.g. Buoy data received at CMC(6 hr period, multiple elements, ~12k reports, from GTS and MSC, ARGO also avail.)
Use J(xi-1) to find xi Data assimilation and forecast Data Acquisition Error Statistics (Observation and Forecast) Data Quality Control Analysis (Spatial QC) Cost function J(x,..) First Guess xa MIN J(x) NWP model Forecasts
CMC Assimilation & Forecast cycle (Only 12-hour Regional cycle shown)
Forecast Error • 10 day MLSP forecast is completely off: Approaching theoretical limit of 2 weeks for weather forecasts (Lorenz: “butterfly effect”)
Ensemble Prediction System Final states (forecasts) Initial states (analyses) True initial state True final state
Ensemble Prediction System observations random numbers perturbed observations trial field or data assimilation 6-h integration with model i perturbed trial field perturbed analysis random numbers and perturbed physics model physics medium-range forecast medium-range integration with model i
Summer SW Point forecasts Winter SW Aviation Air Quality Meteorological Fields Operational Charts EER Plain language Graphics Fcst Product Generation Systems NWP Fields Weather Elements Numerical Weather Prediction Models Verification seasons Long range fcst months Public fcsts EPS Error Feedback Statistics Weather Elements Contingency Tables Perfect Prog Direct Model Output MOS Analog Neural Net Fuzzy Logic Kalman Filters Diagnostic