Modeling System

Descriptions of the system hardware, and model and graphical software are provided below for our daily numerical model runs for the Los Angeles basin.

Model Configuration

The Weather Research and Forecasting (WRF) Model is a next-generation mesocale numerical weather prediction system designed to serve both operational forecasting and atmospheric research needs. It features multiple dynamical cores, a 3-dimensional variational (3DVAR) data assimilation system, and a software architecture allowing for computational parallelism and system extensibility. WRF is suitable for a broad spectrum of applications across scales ranging from meters to thousands of kilometers.

We use the Eulerian mass coordinate core which has been dubbed the Advanced Research WRF (ARW). We use the ARW with and without 3DVAR to produce two daily 36-hour forecasts over California and adjoining areas. The runs where we use observations and 3DVAR in combination with the ARW to produce a forecast are referred to as with data assimilation and those runs where we do not use 3DVAR are referred to as without data assimilation. Unless otherwise noted the following configuration information applies to both the with- and without- data-assimilation runs.

  • Configuration
    • Terrain derived from 30 sec (0.9 km) terrain database
    • Kain-Fitsch Cumulus Parameterization
    • Long and shortwave radiation scheme with cloud radiative cooling (Rapid Radiative Transfer Model & Dudhia schemes)
    • Mixed phase cloud micro-physics (WRF Single-Moment 5-class scheme)
    • Surface layer uses ETA similarity
    • Noah Land Surface Model
    • Mellor-Yamada-Janjic planetary boundary layer (PBL) scheme
    • 37 vertical (half sigma) levels with the top of the model at 100 mb
    • WRF V2.2


Model Dimensions
WRF Domain Graphic
Approx:
Latmin
Latmax
Lonmin
Lonmax
Domain 1 28.649 40.278 -126.087 -109.913
Domain 2 31.913 35.755 -121.364 -114.636
Model Pressure Layers
Model Domain
Model Vertical Profile



Data Assimilation Model - 3DVar

Theory. The 3DVar (3- Dimensional VARiational) data assimilation system is developed by NCAR. We selected it because of its ability to assimilate a wide variety of observations, especially those that are not direct measures of the model state variables (e.g., satellite data). 3DVar observation operators relate the values of the model state variables at the analysis time to observed quantities. Observations are categorized by type, each with its own error statistics. The goal is to minimize the difference between the analysis and observations and a prior estimate of the model state (background). The cost function is shown below. The analysis is a “weighted fit” of all sources of information and the “optimal analysis” is that which minimizes the cost-function:


x = a vector of the model variables at a given time
y = Hx where H is the “observation operator”
O = Observation (instrumental) error
F = Representivity (observation operator) error
B = Background (previous forecast from NAM or WRF) error



Initialization Data Sets

Lateral Boundary Conditions: North American Meso Model (NAM; 40 Km grid spacing)

Water Surface: Navy 1/4 degree Sea Surface Temperatures for with data assimilation runs; NAM Sea Surface Temperatures for without data assimilation runs

Land Surface: NAM (40 Km grid spacing)

Atmospheric Initial Conditions: NAM (40 Km grid spacing) for without data assimilation runs; 3DVAR output for with data assimilation runs.


Background/First Guess Fields (with data assmilation)

    NAM (40 Km grid spacing) 00Z cycle
    WRF forecast (5 Km and 15 Km grid spacing) 06Z and 12Z cycles

Observations (with data assimilation)

  • Surface Data (Click here to view plot)
    • Conventional NWS and DoD surface reports (including ships, buoys)
    • SCAQMD Meteorological observations
    • Bureau of Land Management Remote Automated Weather Obs (RAWS)
    • Buoys from the National Data Buoy Center (NDBC)
  • Profile Data (Click here to view plot)
    • Boundary Layer Profilers (Wind)
    • Aircraft Reports (Aireps)
    • Radiosondes
    • GPS-Met Total Precipitable Water (GPS-Met)
  • Satellite Data (Click here to view plot)
    • GOES Cloud Drift Winds
    • DMSP SSM/I winds and Total Precipitable Water
    • QuikScat Wind speed and direction


Concept of Operations with Data Assimilation






Concept of Operations without Data Assimilation





Postprocessing

The Aerospace real-time system uses the {RIP version4.0} (Read/Interpolate/Plot) visualization program with NCAR Graphics for the all post-processing products. The code is being continuously improved to meet the needs of Aerospace users. Currently the set of visual products includes plots of 3 hour precipitation accumulations, ground/sea-surface temperatures, 2 meter surface air temperatures, low-level relative humidity, winds and absolute vorticity and sounding plots along with other standard parameters used in weather forecasting and atmospheric assessment applications.

Computing Environment

Realtime Modeling Runs
Apple XServe Cluster Environment
27 Xserve G5 Dual 2.3 GHz compute nodes
Each node has 4GB SDRAM
Mac OS X Server 10.4.11

Post-processing
Single Apple XServe G5 Dual 2.3 GHz
4GB SDRAM
Mac OS X Server 10.4.11
and
Sun Fire V490
Solaris 10
8GB RAM

Storage
Overland Ultamus RAID 1200
9TB capacity