NASA Contractor Report NASAICR-2008-2 14755 

Weather Research and Forecasting Model Wind 

Sensitivity Study at Edwards Air Force Base, CA 

... ->-.. -... --· .. ---- Leet1 R:·~Wats6fi-~-~'--"~· ',' 

William H. Bauman III 

Applied Meteorology Unit 

Kennedy Space Center, Florida 

December 2008 

https://ntrs.nasa.gov/search.jsp?R=20100036761 2018-05-17T12:02:24+00:00Z
NASA STI Program ... in Profile 

Since its founding, NASA has been dedicated 

to the advancement of aeronautics and space 

science. The NASA scientific and technical 

information (STI) program plays a key part in 

helping NASA maintain this important role. 

The NASA STI program operates under the 

auspices of the Agency Chief Information 

Officer. It collects, organizes, provides for 

archiving, and disseminates NASA's STI. The 

NASA STI program provides access to the NASA 

Aeronautics and Space Database and its public 

• CONFERENCE PUBLICATION. Collected 

papers from scientific and technical 

conferences, symposia, seminars, or other 

meetings sponsored or co-sponsored 

• 

by NASA. 

SPECIAL PUBLICATION. Scientific, 

technical, or historical information from 

NASA programs, projects, and missions, 

often concerned with subjects having 

substantial public interest. 

interface, the NASA Technical Report Server, • TECHNICAL TRANSLATION. Englishthus providing one of the largest collections of language translations of foreign scientific 

··aeronautIcaranaspacesCienceSTrli1·the·woila:~ ·········~~·anatechnlcarn1ateriarpertlnenTto·······-·~~~ ........ . Results are published in both non-NASA channels NASA's mission. 

and by NASA in the NASA STI Report Series, 

which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 

completed research or a major significant 

phase of research that present the results of 

NASA Programs and include extensive data 

or theoretical analysis. Includes compilations 

of significant scientific and technical data and 

information deemed to be of continuing 

reference value. NASA counterpart of peerreviewed formal professional papers but has 

less stringent limitations on manuscript length 

and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 

and technical findings that are preliminary or 

of specialized interest, e.g., quick release 

reports, working papers, and bibliographies 

that contain minimal annotation. Does not 

contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 

technical findings by NASA-sponsored 

contractors and grantees. 

Specialized services also include creating 

custom thesauri, building customized databases, 

and organizing and publishing research results. 

For more information about the NASA STI 

program, see the following: 

• Access the NASA STI program home page 

at~~,~~c~~~~~ 

• E-mail your question via the Internet to 

• Fax your question to the NASA STI Help 

Desk at (301) 621-0134 

• Phone the NASA STI Help Desk at 

(301) 621-0390 

• Write to: 

NASA STI Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076-1320 
NASA Contractor Report NASAICR-2008-214755 

Weather Research and Forecasting Model Wind 

Sensitivity Study at Edwards Air Force Base, CA 

_- _______ --_-,"~ ____ --_______ ~ ____ l..eela.R. __ WatsOl1»_~~" ___ - _-_-_-_~~~_~, ___ >_ __~-._~~_~_.~. ___ ~_~_ 

William H. Bauman III 

Applied Meteorology Unit 

Kennedy Space Center, Florida 

December 2008 
Acknowledgements 

The authors thank Mr. Brian Hoeth of the Spaceflight Meteorology Group for providing candidate wind 

cycling case days and wind tower data. 

Available from: 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076-1320 

(301) 621-0390 

This report is also available in electronic form at 
Executive Summary 

NASA prefers to land the space shuttle at Kennedy Space Center (KSC). When weather conditions violate 

Flight Rules at KSC, NASA will usually divert the shuttle landing to Edwards Air Force Base (EAFB) in Southern 

California. When this happens, forecasting surface winds at EAFB is a challenge for the Spaceflight Meteorology 

Group (SMG) forecasters due to the complex terrain that surrounds EAFB. One particular phenomenon identified by 

SMG that makes it difficult to forecast the EAFB surface winds is called "wind cycling". This occurs when wind 

speeds and directions oscillate among towers near the EAFB runway leading to a challenging deorbit bum forecast 

for shuttle landings. Accurately forecasting the winds with at least 2 hours lead time along the runway is crucial in 

supporting a landing decision. The large-scale numerical weather prediction models cannot properly resolve the 

wind field due to their coarse horizontal resolutions, so a properly tuned high-resolution mesoscale model is needed. 

The Weather Research and Forecasting (WRF) model meets this requirement. 

The Applied Meteorology Unit (AMU) was tasked to determine the skill of different WRF model 

configurations in forecasting wind cycling cases at EAFB. The AMU compared model output of wind speed and 

direction from two WRF model dynamical cores, two "hot-start" initialization options of the WRF model and three 

physics options. The various combinations of these options resulted in a total of 42 model runs for seven wind 

cycling case days and 14 model runs for two null case days. The initialization options are three-dimensional weather 

.. (lTl(lIysi~ .. §y~te!fl,; .. th.aLLntegr(l!(!!l1l1l!iple.Ill(!te~roI2gic<lI .. d.at<l ... §01l~(;.es.iI21Q.()Tl~~<lTl.aJy~i.~2y,e.L.!h~L.tI.§xr"s.<!2I!1Jl. !L9J ..... 

interest. These analysis systems allow mesoscale models to benefit from the addition of high-resolution data sources 

in their initial conditions. 

The goal of this work was to assess the different model options and to determine which configuration best 

predicted surface wind speed and direction at EAFB. To do so, the AMU compared the WRF model performance 

using both hot start initializations with both dynamical cores and compared model performance while varying the 

physics options. The model ou!put was compared to the wind tower data from the seven wind-cycling and two null 

case days. 

Based on both subjective and objective analyses comparing the model forecast data to observations, the AMU 

found the Advanced Research WRF core and Mellor-Yamada-Janjic planetary boundary layer scheme performed 

better than the other model configurations. The model did not produce a better wind forecast on null days vs. wind 

cycling days. However, the model was able to differentiate between wind cycling days and null cases which would 

provide added value to the SMG forecaster's shuttle deorbit bum landing forecast. 

4 
Table of Contents 

Executive Summary ...................................................................................................................................................... 4 

List of Figures ............................................................................................................................................................... 6 

1. Introduction ............................................................................................................................................................ 8 

1.1 Background Information .............................................................................................................................. 8 

1.2 Report Fonnat and Outline ........................................................................................................................ 10 

2. Data and Methodology ......................................................................................................................................... 11 

2.1 Wind Cycling Candidate Days .................................................................................................................. 11 

2.2 Model Core and Initialization Options ...................................................................................................... 13 

2.3 Data Ingest ................................................................................................................................................. 13 

2.4 Model Configuration ................................................................................................................................. 14 

3. Subjective Analysis .............................................................................................................................................. 16 

" ""3:+",, .................... ;;::":;;":::".";:".';;":;:7::"::":";"::;~:;"::;~.;"::;"."::;"::;"."::";"::.".::-."+6"" " 

3.2 Model Indications of Wind Cycling Events .............................................................................................. 18 

3.3 Subjective Analysis Summary ................................................................................................................... 20 

4. Objective Analysis ............................................................................................................................................... 21 

4.1 Forecast vs. Observed Wind Speed ........................................................................................................... 22 

4.2 Mean Error for Wind Speed and Direction ................................................................................................ 23 

4.3 Pearson Correlation Coefficient for Wind Speed and Direction ................................................................ 25 

4.4 Objective Analysis Summary .................................................................................................................... 27 

5. Conclusions .......................................................................................................................................................... 28 

5.1 Results ....................................................................................................................................................... 28 

5.2 Recommendations ..................................................................................................................................... 29 

References ................................................................................................................................................................... 30 

List of Acronyms ......................................................................................................................................................... 32 

5 
List of Figures 

Figure I. Location of EAFB in relation to surrounding mountainous terrain. Background image from MSN 

Live Search Maps .................................................................................................................................... S 

Figure 2. Wind tower locations on EAFB. The towers along the concrete and lakebed runways are indicated 

by arrows. Background image from MSN Live Search Maps ................................................................. 9 

Figure 3. Visible satellite image showing several mountain wave clouds on 30 Jan 200S near EAFB. The 

yellow circle indicates the location ofEAFB ........................................................................................... 9 

Figure 4. Wind direction (degrees) and wind speed (kts) for Tower 44 (top) and Tower 234 (bottom) on 14 

February 200S from 0000 to 12000 UTC. The green line is the wind direction, the red line is the 

peak wind and the blue line denotes the steady-state wind. Wind cycling occurs in the time period 

between the vertical dashed blue lines ................................................................................................... 12 

Figure 5. Locations of WSR-SSD sites (yellow) used in creating the high-resolution LAPS and ADAS 

analyses. The black circle shows the location ofEAFB ........................................................................ 14 

Figure 6. Maps of EAFB showing the locations of the three concrete runway towers (red dots) and the 

rng<!t:Jgri<!.P9int~.(Qli!c;Is .. dg1s1foLJheMYJ,t11.cl.yQ!1~()i.I~Bk .. ssh()lJ1~~.(ll) .. \llJ.c!QESr]~L.s.<.:b()me 

(b). The grid point closest to each corresponding tower and used in the analysis is shown within 

each green ellipse ................................................................................................................................... 16 

Figure 7. Graphs of a wind cycling case from Tower 224 on 7 June 200S. In each graph the observed 

steady-state wind speed is shown by the dark red line and the forecast steady-state wind speed is 

shown by the orange and blue lines representing various model configurations as identified by the 

graph legends. The light blue shaded box indicates the duration of the wind cycling event. ................ 17 

Figure S. Graphs of a null case from Tower 224 on 9 June 2007. In each graph the observed steady-state 

wind speed is shown by the dark red line and the forecast steady-state wind speed is shown by the 

orange and blue lines representing various model configurations as identified by the graph 

legends ................................................................................................................................................... IS 

Figure 9. Plots of ADAS ARW MYJ Eta model forecast steady-state winds for 7 June 200S in the EAFB 

area for four model output times of 0300 UTC (a), 0400 UTC (b), 0500 UTC (c) and 0600 UTC 

(d). The location of EAFB is indicated by the white circle. The approximate location of the 

concrete runway is shown by the red line. The blue line indicates the leading edge of westerly 

winds ...................................................................................................................................................... 19 

Figure 10. Plots of ADAS ARW MYJ Eta model forecast steady-state winds for 9 June 2007 in the EAFB 

area for four model output times of 0700 UTC (a), OSOO UTC (b), 0900 UTC (c) and 1000 UTC 

(d). The location of EAFB is indicated by the white circle. The approximate location of the 

concrete runway is shown by the red line. The blue line indicates the leading edge of westerly 

winds moving across the Tehachapi Mountains .................................................................................... 20 

Figure 11. Chart of the forecast mean wind speed from the 6 WRF model configurations vs. the observed 

mean wind speed from the EDW towers for the 12-hour forecast period for a) 22 Dec 2006, b) 4/5 

Mar 200S, c) 30 Ju1200S, and d) 9 Jun 2007 ......................................................................................... 23 

Figure 12. Chart of the mean error for wind speed from the 6 WRF model configurations for the 12-hour 

forecast period for a) 22 Dec 2006, b) 4/5 Mar 200S, c) 30 Ju1200S, and d) 9 Jun 2007 ...................... 24 

Figure 13. Chart of the mean error for wind direction from the 6 WRF model configurations for the 12-hour 

forecast period for a) 22 Dec 2006, b) 4/5 Mar 200S, c) 30 Ju1200S, and d) 9 Jun 2007 ...................... 25 

Figure 14. Chart showing the Pearson Correlation Coefficient for wind speed for the 12-hour forecasts for all 

wind cycling and null cases for the six model configurations ................................................................ 26 

Figure 15. Chart showing the Pearson Correlation Coefficient for wind direction for the 12-hour forecasts for 

all wind cycling and null cases for the six model configurations ........................................................... 27 

6 
List of Tables 

Table 1. List of the all wind cycling candidate days, null cases, start of the wind cycling event, and the end of 

the event. ................................................................................................................................................... I 1 

Table 2. List of WRF model run times for all wind cycling and null cases ............................................................. 15 

Table 3. List of the physics options used for each model run for both the ARW and NMM cores ......................... 15 

Table 4. Average mean error and bias (parenthesis) for the 12-hour forecast for each model configuration for 

the wind cycling and null cases. The smallest average mean error/bias is highlighted in red for each 

case day ..................................................................................................................................................... 24 

7 
1. Introduction 

Occasionally, the space shuttle must land at Edwards Air Force Base (EAFB) in Southern California when 

weather conditions at Kennedy Space Center (KSC) violate Flight Rules (FR). However, the complex terrain in and 

around EAFB makes forecasting surface winds a challenge for the Spaceflight Meteorology Group (SMG). In 

particular, "wind cycling" cases, in which the wind speeds and directions oscillate among towers near the EAFB 

runway, present a challenging forecast problem for shuttle landings. An accurate depiction of the winds along the 

runway is crucial in making the landing decision. Global and national scale numerical weather prediction models 

cannot properly resolve the wind field due to their coarse horizontal resolutions, so a properly tuned high-resolution 

mesoscale model is needed. The Weather Research and Forecasting (WRF) model meets this requirement. 

The goal of this work was to assess different WRF model options and to determine how well the model could 

predict surface wind speed and direction at EAFB and if one model configuration perfonned better than others. The 

data presented in this report shows many interesting phenomena that were not investigated in the interest of focusing 

on the SMG goal of model capability and time constraints to complete the work. Given the time constraints and the 

observation that both wind speed and direction indicated cycling during an event, The Applied Meteorology Unit 

(AMU) chose to focus on assessing model performance forecasting wind speed. 

The AMU investigated the WRF model's two dynamical cores: I) the Advanced Research WRF (ARW) and 2) 

.thecNon-hydr.o.st.aticMes.oscale.M{)det(NMM~<There~ar~alsoctw{)Optionsfor"a:'hot-start;'.initializationoftheW~~ 

model: I) the Local Analysis and Prediction System (LAPS; McGinley 1995) and 2) the Advanced Regional 

Prediction System (ARPS) Data Analysis System (ADAS; Brewster 1996). Both LAPS and ADAS are threedimensional weather analysis systems that integrate multiple meteorological data sources into one analysis over the 

user's domain of interest. These analysis systems allow mesoscale models to benefit from the addition of highresolution data sources in their initial conditions. In addition to model core and initialization options, there are many 

different model parameterization options within each core. Having a series of initialization options and WRF cores, 

as well as many options within each core, provides SMG with considerable flexibility as well as decision-making 

challenges. The goal of this study is to assess the different model configurations and to determine which 

configuration will best predict surface wind speed and direction at EAFB. To accomplish this, the AMU compared 

the WRF model performance using ADAS versus LAPS for the ARW and NMM cores and compared model 

perfonnance while varying the physics options. 

1.1 Background Information 

The location of EAFB (Figure I) is on the western edge of the Mojave Desert. It is surrounded by mountains on 

three sides: the San Bernardino Mountains to the southwest, the San Gabriel Mountains to the south, and the 

Tehachapi Mountains and southern Sierra Nevada to the west and northwest. The base is located next to Rogers Dry 

Lake whose hard surface acts as an extension to the EAFB runways. There are twelve wind towers located in and 

around the EAFB runways. Figure 2 shows the locations of the towers along the EAFB runway complex. 

Figure I. Location of EAFB in relation to surrounding mountainous 

terrain. Background image from MSN Live Search Maps. 

8 
Figure 2. Wind tower locations on EAFB. The towers along the concrete and lakebed 

runways are indicated by arrows. Background image from MSN Live Search Maps. 

The complex terrain located in and around EAFB creates the erratic behavior of the wind speed and direction 

seen on the EAFB wind towers. Wind cycling events occur when there is an oscillation in wind direction and/or 

wind speed among the wind towers in the network near the EAFB runway complex. During these cycling events, the 

wind speed and direction reported from the towers near the concrete runway (Towers 44, 220, 224) are noticeably 

different than that reported from towers near the lakebed runway (Towers 154, 230, 234). These events usually last 

from 90 minutes up to 4 hours or longer and most often occur when the prevailing wind is from the northwest or 

west-northwest. 

The wind cycling events at EAFB are thought to be caused by mountain waves created as winds from the 

northwest flow over the Tehachapi Mountains. Figure 3 shows several mountain wave clouds that developed near 

EAFB on 30 January 2008. On this day, the wind along the concrete runway fluctuated from westerly and relatively 

lighter to northwesterly and stronger. The winds along the lake bed were northwesterly and strong during the wind 

cycling event. 

9 

Figure 3. Visible satellite image showing several mountain wave clouds on 30 January 2008 

near EAFB. The yellow circle indicates the location ofEAFB. 
1.2 Report Format and Outline 

This report presents the findings from a one-year study of model sensitivities for predicting wind cycling events 

at EAFB. This analysis examined different model cores, initializations, and physics options to detennine the impact 

on model skill. Two local data analysis systems that ingested satellite data, radar data, and surface observations 

across southern California were used to initialize a mesoscale model. Model skill was assessed using both a 

subjective and objective verification that compared forecast to observed wind speed and direction from the local 

EAFB wind towers. Section 2 describes the data and methodology including an example of a wind cycling case at 

EAFB, the different model cores, model initialization options, and the data ingested. The results of the subjective 

and objective analyses are presented in Sections 3 and Section 4, respectively. Section 5 summarizes the report and 

provides recommendations on the best mesoscale model configuration for operational use based on this study. 

10 
2. Data and Methodology 

The important aspects ofthis study were the choice of wind cycling case days, the model configuration, and the 

data used to initialize the models. The candidate wind cycling days chosen occurred between December 2006 and 

July 2008. Comparisons were made between the ARW and NMM core initialized with ADAS and LAPS, as well as 

with different physics parameterizations. Both subjective and objective analyses were used to verify the forecasts. 

2.1 Wind Cycling Candidate Days 

The period of record for choosing wind cycling candidate days was December 2006 through July 2008. SMG 

provided the AMU with the dates of seven wind cycling cases. Six of the seven were strong wind cycling cases and 

one, the 30/31 July 2008 case, was marginal in that the fluctuations in wind speed and direction were smaller than 

for all other wind cycling cases. Two null, or non-wind cycling, days were also randomly chosen from the available 

data. Table I lists the candidate days, the start of the wind cycling event, and the end of the event. 

Table 1. List of all wind cycling candidate days, null cases, start of the wind 

cycling event, and the end of the event. 

Candidate Day Start of Event : 

I End of Event 

22 Dec 2006 0900 UTC 1800 UTC 

30 Jan 2008 1030 UTC 1330 UTC 

14 Feb 2008 0300 UTC 0700 UTC 

5 Mar 2008 0000 UTC 1000 UTC 

4 Jun 2008 1100 UTC 2000 UTC 

7 Jun 2008 0200 UTC 1300 UTC 

30/31 Jul 2008 1900 UTC 0100 UTC (31 Jul) 

9 Jun 2007 NULL CASE 

1 Nov 2007 NULL CASE 

Figure 4 shows a time series of wind speed and direction at Towers 44 and 234 from a wind cycling event that 

occurred on 14 February 2008. Inspection of Figure 4 reveals a wave-like behavior in the time series between 0300 

and 0700 UTC (denoted by the blue vertical lines in Figure 4) for wind speed and direction at Tower 44 (top), which 

is located near the concrete runway. During this time there was an oscillation between southwesterly and 

northwesterly wind, as well as a 5 to 15 kt change in wind speed, that occurred approximately every half hour. The 

wind at Tower 234, near the lakebed runway, did not exhibit this wave-like behavior. 

II 
-;;;- 350 <II 300 <II 

... 

tID 250 <II 

::!:. 200 c: 150 0 

:;:; 100 v 

<II 

... SO 

0 0 

Time (UTC) 

Steady Peak 

-;;;- 40 .... 

o 12 30 

-

"g. 20. <II 

~ 10 ·r------' 

o ~~~~~~r~'~rm_~"~~!~~m~IT~~~~~~~~m~rr~~~~~~~~~~ 

Time (UTC) 

-;;;- 350 <II 300 <II 

... 

tID 250 <II 

::!:. 200 c: 150 0 

:;:; 100 v 

<II 

... SO 

0 0 

Time (UTC) 

-;;;-40 

.... 

o 

c: 30 

=-

-g 20 +-~-. 

<II 

~ 10 ·t·······_·-

o ~~~m=~¥~~~=~~·m*~rn~~~~~mr~~rn~~~~=~=~m~~=mrnrn~=rnrn~ 

Time (UTC) 

Figure 4. Wind direction (degrees) and wind speed (kts) for Tower 44 (top) and Tower 234 (bottom) 

on 14 February 2008 from 0000 to 12000 UTe. The green line is the wind direction, the red line is the 

peak wind and the blue line denotes the steady-state wind. Wind cycling occurs in the time period 

between the vertical dashed blue lines. 

12 
2.2 Model Core and Initialization Options 

For this study, the AMU employed both the ARW and NMM cores of the WRF model. The ARW core was 

developed primarily at the National Center for Atmospheric Research (NCAR) and is a fully compressible, nonhydrostatic mesoscale model with a hydrostatic option. It consists of a mass-based hydrostatic pressure terrain 

following coordinate, Arakawa C-grid staggering for the horizontal grid, time-split integration using a third order 

Runge-Kutta scheme with a small step for acoustic and gravity wave modes, and up to sixth order advection options 

in the horizontal and vertical (Skamarock et al. 2005). There are also full physics options for microphysics, 

planetary boundary layer, cumulus parameterization, radiation, and land surface schemes (Skamarock et al. 2005). 

The NMM core is also a fully compressible, non-hydrostatic mesoscale model with a hydrostatic option (Janjic et al. 

200 I, Janjic 2003a,b). It consists of a hybrid sigma-pressure, terrain following vertical coordinate, Arakawa E-grid, 

a forward-backward time integration scheme, a second order advection option in the horizontal and vertical, and 

conservation of energy and enstrophy (Janjic 1984). Most physics packages available to the ARW have not been 

tested with the NMM and therefore, the physics options for the NMM are more limited than for the ARW. 

A "hot-start" initialization of the WRF model can be made using LAPS or ADAS. LAPS is a data assimilation 

tool that uses numerous meteorological observations, such as satellite data, radar data, and surface observations, to 

generate a three-dimensional representation of the atmospheric forcing fields, such as wind speed and direction, 

.. ~llrface .. temper(lture.al1dJ~r.~ssllre,.~~I.ati,:e.hlllT1i.ditx,.pr~cipitati()Il .. (ln~cl()lldc(),:~~.(~c<:Jil1lete!.<lI' .. L~~I; .. ~!~ers .. 

1995; Albers et al. 1996; Birkenheuer 1999; McGinley 1995). It includes a wind analysis and a three-dimensional 

cloud analysis, which are needed for the LAPS hot-start initialization. The LAPS cloud analysis is designed to create 

consistency with all data and the typical meteorology of clouds by combining data from infrared and visible satellite 

data, three-dimensional LAPS radar reflectivity derived from the full-volume radar data, and the LAPS threedimensional temperatures (Albers et al. 1996). Fields derived from the cloud analysis include cloud liquid water, 

cloud type, cloud droplet size, and icing severity (Albers et al. 1996). However, in the latest versions of LAPS, the 

option to initialize precipitation was turned off. Mr. Chris Anderson of the National Oceanic and Atmospheric 

Administration Earth System Research Laboratory Global Systems Division (GSD) indicated that at a l-km grid 

spacing, initializing WRF with hydrometeors could create a large outflow that would degrade the forecasts. 

Consequently, GSD zeroed out the precipitation particles in the initialization, but left the cloud particles. He 

recommended keeping the precipitation initialization turned off if it was not an essential component of the research. 

Therefore, the AMU decided not initialize to precipitation since this study was not focused on convective 

forecasting. 

The ADAS system, developed by the University of Oklahoma, has two main components. The first is a Bratseth 

objective analysis scheme that evaluates pressure, wind, potential temperature, and specific humidity. The second 

component is a three-dimensional cloud analysis scheme that is used for the hot-start initialization (Zhang et al. 

1998). The ADAS cloud analysis is based on the LAPS cloud analysis with some modifications (Case et al. 2002). It 

uses surface observations of cloud cover and height, satellite data, and radar data to determine the cloud cover, cloud 

liquid and ice water, cloud type, rain/snow/hail mixing ratios, icing severity, in-cloud vertical velocity, cloud base 

and top, and cloud ceiling (Case et al. 2002; Zhang et al. 1998; Brewster 2002). 

2.3 Data Ingest 

Data ingested by the model through either the LAPS or ADAS analysis packages included Level II Weather 

Surveillance Radar-1988 Doppler (WSR-88D) data from the Las Vegas, NV (KESX), Yuma, AZ (KYUX), Santa 

Ana Mountains, CA (KSOX), San Joaquin Valley, CA (KHNX), San Diego, CA (KNKX), and Los Angeles, CA 

(KVTX) radars (Figure 5), Geostationary Operational Environmental Satellites (GOES) visible and infrared satellite 

imagery, Meteorological Assimilation Data Ingest System (MADIS; data, and local EAFB 

wind tower data. The Level II WSR-88D data contained full volume scans of reflectivity at a resolution of 10 by 1 

km, radial velocity at 10 by 0.25 km, and spectrum width data at a 10 by 0.25 km (Fulton et al. 1998). These data 

were available every 4 to 6 minutes. The GOES-12 visible imagery was available at a 1 km horizontal resolution 

every IS minutes, and the infrared imagery was available at a 4 km horizontal resolution also every 15 minutes. 

Both visible and infrared imagery provided brightness temperatures to the analysis packages. Each initialization 

analysis used a different background model to create the first-guess background fields. The choice of background 

data for both LAPS and ADAS was based on the highest resolution model data currently supported in the latest 

versions of both analysis packages. The latest version of LAPS allows the use of the Rapid Update Cycle (RUC) 20-

km model as background data, while ADAS supports the North American Model (NAM) 12-km data. 

13 
Figure 5. Locations of WSR-88D sites (yellow) used in creating the high-resolution LAPS and ADAS analyses. 

The black circle shows the location of EAFB. 

A suite of scripts written by the AMU formed the core for initializing the WRF model with ADAS output. The 

analysis package ingested satellite data in Man computer Interactive Data Access System (McIDAS) Area format 15 

minutes prior to model initialization time; raw, full volume radar data within 10 minutes of the model initialization 

time; and surface data from the hour of and hour prior to model initialization. The largest task in configuring LAPS 

was working with the ingest code. The ingest code can only be used with raw data that have the same configuration 

and format as GSD's raw data, which is the Network Common Data Form (NetCDF). Therefore, all data files were 

converted to NetCDF format to be used within LAPS. Software to convert the WSR-88D data to NetCDF format 

was obtained from GSD. Scripts were written to convert raw wind tower data into a format ingestible by both LAPS 

and ADAS. Converting satellite data in MclDAS Area format to NetCDF required several steps. The data were first 

ported to the local Meteorological Interactive Data Display System (MIDDS) system where they were remapped to 

the Lambert Conformal projection using the IMGREMAP command. Next, the remapped data were run through a 

program that converted them to NetCDF format. This program is called AreaToNetCDF and is available from the 

Space Science and Engineering Center (SSEC) at the University of Wisconsin. All of the reformatted data files were 

then ingested into LAPS to create an initialization field for the model. 

2.4 Model Configuration 

Each model simulation was run at a l-km horizontal grid spacing over EAFB and adjacent mountainous areas 

with 51 irregularly spaced, vertical sigma levels. Each run was integrated 12 hours and began one to three hours 

before the onset of the observed wind cycling event. Table 2 lists the WRF model run times for each of the 

candidate days. Several factors were considered in detennining the grid spacing of the WRF runs: the average 

distance between the wind towers, the customer requirements, and the computational requirements to run a highresolution model. Grasso (2000) notes that the theoretical minimum resolvable wavelength by a model is twice the 

horizontal grid spacing, but Walters (2000) suggests that as many as 10 grid points may be required for a more 

reasonable representation of the true solution. The AMU detem1ined that a 2-km or less horizontal grid spacing 

would be sufficient to resolve the wind cycling events given the average distance of 8 km between the wind towers 

on the runway (44, 220, and 224) and the lakebed (154,230, and 234) and taking into account the spacing of all the 

wind towers. However, customer requirements indicated that a 1.3-km or less grid spacing was desired. To 

determine the highest resolution that could feasibly be run in real-time, three different runs using 1.33, 1, and 0.66-

km grid spacing were made to determine the optimal resolution for the WRF forecasts. The l-km run slightly 

14 
outperformed the l.3-km run with respect to wind speed and direction forecasts. The model was then run with a 

0.66-km grid spacing, but it took more than twice as long to run as compared to the I-km WRF forecast. Since this 

is not practical for real-time operations, the 0.66-km WRF run was not considered and the I-km grid spacing was 

used. 

Table 2. List ofWRF model run times for all wind cycling and null cases. 

Candidate Day Model Start Time Model End Time 

22 Dec 2006 0600 UTC 1800 UTC 

30 Jan 2008 0900 UTC 2100 UTC 

14 Feb 2008 0000 UTC i 

i 1200 UTC 

4/S Mar 2008 2100 UTC (4 Mar) 0900 UTC (S Mar) 

4 Jun 2008 0900 UTC 2100 UTC 

7 Jun 2008 0000 UTC 1200 UTC 

30/31 Jul 2008 1800 UTC (30 Jul) 0600 UTC (31 Jul) 

9 Jun 2007 0600 UTC 1800 UTC 

-------- -- -------lNO-v2007---- <. <'.»-'~ >----- 0900 I JTe' . ____ 2100 UTe' 

Boundary conditions for the WRF runs were obtained from the NAM model with a horizontal grid spacing of 

12-km. The NAM model domain is on a Lambert conformal projection with a horizontal domain size of 614 by 239 

grid points and 42 vertical levels. It produces an 84-hr forecast every six hours, or four times per day. This model 

was chosen for boundary conditions as it has the best resolution of all available datasets. 

Each candidate day consisted of six runs using different model configurations, for a total of 42 runs for the wind 

cycling case days and 14 for the null cases. The different configurations included: 

• LAPS-ARW with the Yonsei University planetary boundary layer (PBL) scheme and NCAR/Penn State 

Mesoscale Model Version S (MMS) similarity surface layer scheme (LAPSARW _ Yonsei), 

• LAPS-ARW with the Mellor-Yamada-Janjic PBL scheme and ETA similarity surface layer scheme 

(LAPSARW _MYJ), 

• LAPS-NMM with the National Center for Environmental Prediction (NCEP) Global Forecast Systems 

(GFS) PBL scheme and NCEP GFS surface layer scheme (LAPSNMM_GFS), 

• LAPS-NMM with the Mellor-Yamada-Janjic PBL scheme and ETA similarity surface layer scheme 

(LAPSNMM_MYJ), 

• ADAS-ARW with the Yonsei University PBL scheme and MMS similarity surface layer scheme 

(ADASARW _ Yonsei), and 

• ADAS-ARW with the Mellor-Yamada-Janjic PBL scheme and ETA similarity surface layer scheme 

(ADASARW _ MYJ). 

All other physics parameters were the same for each ARW model run. The NMM runs did not use the same physics 

as the AR W runs since the ultimate goa) of the task was to see if wind cycling could be simulated. Instead, the best 

and most widely tested choice of physics found in the literature for the NMM core were used. Table 3 lists the 

physics options used in both the ARW and NMM runs that were held constant. 

Table 3. List of the physics options used for each model run for both the ARW and NMM cores. 

ARW NMM 

Microphysical scheme Lin et al. (1983) Ferrier (1994) 

Cumulus Parameterization None None 

Land surface option Noah Land Surface Model (Chen Noah Land Surface Model (Chen and 

and Dudhia 2001) Dudhia 2001) 

Shortwave radiation scheme Goddard (Chou and Suarez 1994) GFDL (Lacis and Hansen 1974) 

Longwave radiation scheme RR TM (Mlawer et al. 1997) GFDL (Fels and Schwarzkopf 1975; 

Schwarzkopf and Fels 1985, 1991) 

IS 
3. Subjective Analysis 

The AMU completed a subjective analysis of the WRF forecasts for all candidate days. The goals of the 

subjective analysis were to determine 

• If the model was able to predict the timing and/or magnitude of the wind cycling events at the concrete 

runway towers by comparing the observed wind speed to the forecast wind speed and 

• If the model could provide the forecasters with an indication of whether or not a wind cycling event was 

likely to occur by assessing the model forecasts on wind cycling days and null case days. 

3.1 Model Performance at the Concrete Runway Towers 

Using the Grid Analysis and Display System software, the AMU extracted the model forecast winds for every 

model run at each grid point nearest to each of the three concrete runway towers (44, 224 and 220 as shown in 

Figure 2) for comparison. The three tower locations and nearest corresponding grid point used in the comparison are 

shown in Figure 6. The model grid was identical for the MYJ and Yonsei PBL schemes and is shown in Figure 6a, 

while the model grid was slightly different for the GFS PBL scheme and is shown in Figure 6b. 

Figure 6. Maps of EAFB showing the locations of the three concrete runway towers (red dots) and the 

model grid points (black dots) for the MY] and Yonsei PBL schemes (a) and GFS PBL scheme (b). The 

grid point closest to each corresponding tower used in the analysis is shown within each green ellipse. 

The AMU plotted graphs of the observed steady-state wind speed at all three towers and overlaid the forecast 

steady-state wind speed from each of the model configurations for 42 wind cycling and 14 null case model runs. 

There are three graphs per tower per day that display the output from the different model configurations. These 

graphs included the following configurations: 

• LAPS ARW 

o Yonsei PBL scheme and MM5 similarity surface layer scheme 

o MY] PBL scheme and Eta similarity surface layer scheme 

• LAPS NMM 

o GFS PBL scheme and GFS similarity surface layer scheme 

o MY] PBL scheme and Eta similarity surface layer scheme 

• ADASARW 

o Yonsei PBL scheme and MM5 similarity surface layer scheme 

o MY] PBL scheme and Eta similarity surface layer scheme 

16 
An example of a wind cycling event at Tower 224 on 7 June 2008 is shown in Figure 7. The event began at 

0200 UTe and indicates a wave-like behavior in the time series that ended at 1300 UTe (denoted by the blue shaded 

box). During this time there were six oscillations in wind speed indicated by a change of 3 to 6 ms-' (6 to 12 kt) 

occurring approximately every 45 minutes. All of the model configurations correctly forecast the general trend of 

the wind speeds with the strongest occurring from 0000 to 0200 UTe then decreasing until 0800 UTe then 

increasing again until the end of the cycling event. Only the ADAS ARW configurations showed any indication of 

the wind speed cycling within the general trend. The ADAS ARW MY] Eta configuration indicated four to five 

oscillations whereas the ADAS ARW Yonsei MM5 configuration indicated three oscillations within the event 

timeline. The MYJ Eta configuration also had larger magnitudes in the wind speed oscillation compared to the 

Yonsei MM5. This result was consistent for the other two towers for this event. 

A review of graphs at Tower 224 from a null case on 9 June 2007 is shown in Figure 8. Wind cycling was not 

observed at any of the EAFB towers on this day. The model was run from 0600 to 1800 UTe in the six 

configurations identical to the cycling events. During this time, the observed winds indicated a general weakening 

trend that all model configurations forecast. The LAPS NMM configurations indicated a low bias while the LAPS 

ARW and ADAS ARW indicated no overall bias and produced a better forecast than the LAPS NMM. 

Tower 224 7 Jun 08 LAPS ARW a Tower 224 7 Jun 08 LAPS NMM b 

" 1-~, ---;:-;-- --~.-~--~<-~--.. ----------.---- - -- -- - - --------

o 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 IB 19 20 21 22 25 o 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 

TIme (UTe) Time (UTC) 

r···---------------------------------- ---- -----------::l 

I Tower 224 7 Jun 08 ADAS ARW C 

~ -Tower224 - MYJEta -Vonsei MMS I ,. - . , " . 12 ----r ~ 

: ~ 10 1 -

"0 

.. 8 .. 

Co 

'" .. 

r; 

.. . ~ 

o 1 :2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 

Time (UTC) 

Figure 7. Graphs of a wind cycling case from Tower 224 on 7 June 2008. In each graph the observed steady-state 

wind speed is shown by the dark red line and the forecast steady-state wind speed is shown by the orange and blue 

lines representing various model configurations as identified by the graph legends. The light blue shaded box 

indicates the duration of the wind cycling event. 

17 
Tower 224 9 Jun 07 LAPS ARW a -Tower224 ~MY Eta -YonseiMMS 

12 

! 8~·~~··~t··~·········································· .......................................................................... . 

".. 

"0 

" 

Tower 224 9 Jun 07 LAPS NMM b - Tower224 ""~MYJEta -GFSGFS 

·······················~··························· .......................................................................................... ~ ..................... . 

1O+.~ ......................................................................................................................................................................... . 

~ 6 ii' ... c·'H!J~:L··············· ........................................................................... i~··················· 

ill, E 4 -l·····················+·~1LC .. "'.Af~"" .. ~: .............................. ~ ................ ~; 

~ " 

o 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 

Time (UTe) 

Tower 224 9 Jun 07 ADAS ARW c ower224 ,=~~MY Eta -YonseiMM5 

12 

4 5 6 7 13 9 

Figure 8. Graphs of a null case from Tower 224 on 9 June 2007. In each graph the observed steady-state wind 

speed is shown by the dark red line and the forecast steady-state wind speed is shown by the orange and blue lines 

representing various model configurations as identified by the graph legends. 

3.2 Model Indications of Wind Cycling Events 

None of the model configurations were able to accurately forecast the timing and magnitude of the steady-state 

wind speed oscillations at model grid points nearest to the wind towers. However, the model can still provide 

valuable information to the forecasters to indicate whether a wind cycling event may occur. The AMU compared 

plots of forecast steady-state winds over a large part of the model domain for all of the wind cycling events and the 

null cases to see if the model could discern between the event days and null days. 

Figure 9 shows four time intervals of model forecast steady-state winds plotted in the EAFB area for a wind 

cycling event day, 7 June 2008, from the ADAS ARW MYJ Eta configuration. The figure shows a progression of 

changing wind speed and direction near the concrete runway beginning with westerly winds of 7 -10 kt at 0300 UTC 

(Figure 9a), then southwesterly at 10-15 kt at 0400 UTC (Figure 9b), changing to south-southeast at 10-12 kt at 0500 

UTC (Figure 9c) and finally southeast at 15-20 kt at 0600 UTC (Figure 9d). As the northwesterly winds move across 

the Tehachapi Mountains and southern extent of the Sierra Nevada Mountains, they also move through the 

Tehachapi Pass that lies between the northeast end of the Tehachapi Mountains and the southern end of the Sierra 

Nevada Mountains (see Figure 1). During wind cycling events, the winds moving through the Tehachapi Pass will 

extend all the way to EAFB causing the wind speed and direction changes as the winds oscillate through the pass. 

The blue line in Figure 9 shows the approximate leading edge of the northwesterly winds moving across the 

mountains and through the pass. The winds moving through the pass cause the line to bulge eastward towards EAFB 

over time and disrupt the existing wind flow in the region. It is interesting to note that the Tehachapi Pass Wind 

Farm, located just east and south of the Tehachapi Pass, is one of California's largest wind farms used for generating 

electricity (Edison International 2008). 

18 
Figure 9. Plots of ADAS ARW MYJ Eta model forecast steady-state winds for 7 June 2008 in the EAFB area for 

four model output times of 0300 UTe (a), 0400 UTe (b), 0500 UTe (c) and 0600 UTe (d). The location ofEAFB 

is indicated by the white circle. The approximate location of the concrete runway is shown by the red line. The 

blue line indicates the leading edge of westerly winds. 

Figure 10 shows four time intervals of model forecast steady-state winds plotted in the EAFB area for a null 

case, 9 June 2007 from the ADAS ARW MYJ Eta configuration. The figure shows a nearly constant southwest wind 

flow with speeds of 5-10 kt throughout the series of images from 0700 to 1000 UTe near the concrete runway. The 

approximate leading edge of the northwesterly winds moving across the mountains and through the pass (blue line) 

was semi-stationary throughout the entire 12-hr model run. The northwesterly wind flow over the mountains was 

weaker in this null case compared to the cycling case shown in Figure 9, and the model correctly forecast a null 

event at EAFB. 

19 
Figure 10. Plots of ADAS ARW MYJ Eta model forecast steady-state winds for 9 June 2007 in the EAFB area 

for four model output times of 0700 UTC (a), 0800 UTC (b), 0900 UTC (c) and 1000 UTC (d). The location of 

EAFB is indicated by the white circle. The approximate location of the concrete runway is shown by the red line. 

The blue line indicates the leading edge of westerly winds moving across the Tehachapi Mountains. 

3.3 Subjective Analysis Summary 

Based on the subjective analysis of all wind cycling and null case days, the ARW runs outperformed the NMM 

runs. The NMM core consistently under-predicted wind speeds by the largest margin and did not capture changes in 

wind direction as well as the AR W core. Changing the model core seemed to have the biggest impact on the 

forecasts, while changing the model physics seemed to have the least impact on the forecasts. Although it appears 

the MYJ Eta scheme may have captured the cycling a bit better than the other configurations within the ARW core. 

The model did not forecast wind speeds better for the null days vs. the wind cycling days, but it clearly 

differentiated between cycling days and null cases used in this task. This would help the forecaster discern when a 

wind cycling case would occur. 

20 
4. Objective Analysis 

For the objective analysis, observed wind speed was compared to forecast wind speed using the latest version of 

the Model Evaluation Tools (MET) software. This software was developed by the NCAR Developmental Testbed 

Center. It is a state-of-the-art suite of verification tools that uses output from the WRF model to compute standard 

verification scores comparing gridded model data to point or gridded observations. Each statistic computed for this 

task compared the gridded WRF model data, available every 15 minutes, to the observations from the 12 wind 

towers at EAFB. In addition, the observed wind direction from the 12 wind towers was compared to forecast wind 

direction from the WRF model data. However, the MET software does not currently support objective analysis of 

wind direction. Therefore, the forecast wind direction was manually puIIed from the model data and all statistics 

were computed using Microsoft Excel. 

Many output statistics are available within the MET software. This study looked at three of those statistics for 

wind speed: the forecast vs. the observed mean, the mean error, and the Pearson Correlation Coefficient (PCC); and 

two statistics for wind direction: the mean error and the PCe. The statistics compared all 12 towers combined to the 

corresponding locations in the model forecast output at a 15-minute interval. Towers that indicated wind cycling 

events were not separated from those that did not indicate wind cycling. Thus, the objective analysis is not a standalone analysis that indicates whether the wind cycling phenomena were captured, rather it only shows how well the 

.... 111 oci~ Ip~r[()rl11 e.d. ():v~rall. I~e .Sll b jective<ll1al)'~is il1 S ecti()n 3 .\v(l~used to ~eterJ11 il1e. if the .111()c1e,Ic.apture,c1the \Vin,c1 ...... . 

cycling events. However, comparison of the objective analyses from the null cases to the wind cycling cases may 

help to determine whether the model is adept at forecasting a mesoscale process, such as wind cycling. 

The mean forecast vs. observed wind speed comparison shows how w~1I the forecast wind speed corresponded 

to the observed wind speed. The mean forecast and observed wind speed, f and 0, respectively, were computed as 

follows: 

1 n 

0= I(o,) 

n ;=1 

where: 

n 12 for the number oftower locations, 

j, forecast wind speed at a tower location, and 

0, observed wind speed at a tower. 

The mean error, ME, is a measure of the overall bias ofthe wind speed or direction. A perfect forecast has ME 

o. It is defined as: 

1 n -

ME = ICt: -0,)= f-o 

n 1=1 

where n,j" and 0, are defined as above and: 

f average forecast wind speed/direction at 12 tower locations, and 

o average observed wind speed/direction at 12 towers. 

The PCC, r, measures the strength of the linear association between the forecast and observed wind speed or 

direction. It is defined as: 

21 
,=1 1=1 

where 

j; forecast wind speed/direction at all tower locations for one output period, 

0, = observed wind speed/direction at all tower locations for one output period, 

f = average forecast wind speed/direction for 12-hour forecast, and 

-

o = average observed wind speed/direction for 12-hour forecast. 

The pee can range between -1 and I; I indicates a perfect correlation, -I indicates a perfect negative correlation, 

and 0 indicates no correlation between the forecast and observations. Specifically, the pee for wind speed measures 

whether large values of forecast wind speed tend to be associated with large values of observed wind speed (positive 

""correlation);" whether "smalt values ufforecast" wlndspeed~tendto"mnrssuctatedwith·targe-vatuesofobserved·wind=="· 

speed or vice versa (negative correlation), or whether values of both variables are unrelated (correlation near 0). For 

wind direction, the pee measures the rotation of the wind with a perfect positive correlation indicating that the 

winds are shifting in the same direction with the same magnitude and a perfect negative correlation indicating that 

the winds are shifting in opposite directions with the same magnitude. 

4.1 Forecast vs. Observed Wind Speed 

All WRF model configurations under-predicted the wind speed throughout the forecasts for all wind cycling 

case days, except for the 30/31 July 2008 run, which was the marginal wind cycling case day. Figure II shows the 

forecast vs. observed mean wind speed for three of those wind cycling case days: 22 December 2006 (Figure II a), 

4/5 March 2008 (Figure lIb), and 30/31 July 2008 (Figure IIc), and one null case: 9 June 2007 (Figure lid). As 

can be seen from Figure II a and b, the only model configuration that did not under-predict the wind speed was the 

LAPSARW _MYJ on 4/5 March (Figure lIb). For 30 July (Figure II c), all model configurations over-predicted the 

wind speed, except for the LAPSARW_ Yonsei configuration. For the 9 June null case (Figure lId), the model 

slightly under-predicted the wind speed. 

22 
Forecast vs. Observed Mean a 22 Dec 2006 

14 

~12 "C 10 QI 

QI 

C. " LAPSARW_ Vonsei <1\ 

"C lAPSARW_MYJ C 

~ i ADASARW_Yonsei 

"C x ADASARW_MYJ QI 4 c: LAPSNMM_GFS QI 

.a '" LAPSNMM""MYJ 0 

0 4 10 12 14 

Forecast Wind Speed (m/s) 

Forecast vs. Observed Mean 

30/31 Jul 2008 

~ 14 

..§. 

"C 12 QI 

~10 <1\ 

"C 8 LAPSARW _MYJ C 

~ "C 

QI 

c: 4 ¥ LAPSNMM_GFS QI 

.a '" 2 LAPSNMM_MYJ 0 

0 4 10 12 14 

Forecast Wind Speed (m/s) 

14 

Vi 

112 

"C 

QI 10 QI 

C. 

<1\ B "C 

C 

~ "C 

QI 

c: 

QI 

II> 

.a 2 

0 

0 

~7 

..§. 

"C 6 QI 

QI 

C. 

<1\ 

"C 4 

~ "C 

QI 

c: 2 QI 

<II 

.a 1 0 

G 

Forecast vs. Observed Mean 

4/5 Mar 2008 

6 8 10 12 

Forecast Wind Speed (m/s) 

Forecast vs. Observed Mean 

9Jun 2007 

Forecast Wind Speed (m/s) 

b 

* lAPSARW ~ YOOSei 

LAPSARW_MYJ 

" ADASARW ~ Yonsel 

x ADASARW _MYJ 

x LAPSNMM_GFS 

LAPSNMM _MYJ 

d 

& LAPSARW~Yonsel 

LAPSARW_MYJ 

ADASARW ~ Yonsei 

ADASARW_MYJ 

x lAPSNMM_GFS 

lAPSNMM_MYJ 

Figure 11. Chart of the forecast mean wind speed from the six WRF model configurations vs. the observed mean 

wind speed from the EAFB towers for the 12-hour forecast period on a) 22 December 2006, b) 4/5 March 2008, c) 

30 July 2008, and d) 9 June 2007. 

4.2 Mean Error for Wind Speed and Direction 

As stated above, the mean error is a measure of the overall bias of the wind speed or direction. Figure 12 shows 

the mean error for wind speed for the same three wind cycling case days and one null case as in Figure 11. Again, it 

is apparent that all model configurations generally under-predicted the wind speed for both the 22 December (Figure 

12a) and 4/5 March (Figure 12b) cases. In both cases, the LAPSNMM_GFS was the worst performer, underpredicting the wind speed by 4 to 8 mls throughout the forecast. Figure 12a reveals that the only model 

configurations to over-predict wind speed during the forecast were both of the ADAS runs. By the second half of the 

forecast, all model configurations except the LAPSNMM_ GFS had nearly the same mean error values. Figure 12b 

shows that all of the ARW runs were clustered throughout the entire 12-hour forecast, indicating that they captured 

and missed all of the same fluctuations in the observed wind speed. The mean error for 30/31 July (Figure 12c) 

shows the same range in mean error values as for 22 December and 4/5 March, but reversed. That is, in general the 

model over-predicted the wind speeds throughout the forecast with the mean error falling in the positive range. 

Unlike the other two wind cycling cases, the NMM runs outperformed all ARW runs with the LAPSNMM_MYJ 

performing the best and the ARW MYJ runs performing the worst. For the null case (Figure 12d), the latter half of 

the forecast was more accurate than the beginning of the forecast. The mean error range for this case was smaller 

than for the wind cycling cases, but this was expected since the wind speed values were half of those for the wind 

cycling days. A comparison of the bias, given below, will give a more accurate comparison of the wind cycling days 

to the null cases. 

23 
10 

Mean Error - Wind Speed 

22 Dec 2006 

Forecast Hour 

-LAPSARW_Yonsei ~LAPSARW_MVJ =--ADASARW_Yonsei 

AOASARW.MYJ ·~···lAPSNMM.GFS -lAPSNMM.MYJ 

§ 

o 

Mean Error - Wind Speed 

30/31Jul 2008 

Forecast Hour 

-lAPSARW.Yonsei -LAPSARW.MYJ ~AOASARW.Yons.i 

ADASARW.MYJ --lAPSNMM.GFS -LAPSNMM.MYJ 

§ 

o 

Mean Error - Wind Speed 

4/5 Mar 2008 

o 

o 

C:\ 

Forecast Hour 

-lAPSARW.Yonsei ~LAPSARW.MYJ --ADASARW.Yonsei 

AOASARW.MYJ ~lAPSNMM.GFS -LAPSNMM.MYJ 

8 

8 § 

o 

Mean Error - Wind Speed 

9Jun 2007 

Forecast Hour 

-lAPSARW_'fonsei -lAPSARW_MYJ ~ADASARW_Yonsej 

ADASARW.MYJ --LAPSNMM.GFS -LAPSNMM_MYJ 

b 

d 

Figure 12. Chart of the mean error for wind speed from the six WRF model configurations for the 12-hour 

forecast period for a) 22 December 2006, b) 4/5 March 2008, c) 30 July 2008, and d) 9 June 2007. 

Table 4 lists the average mean error and bias for wind speed for the entire 12-hour forecast for each model 

configuration for each wind cycling case day and the two null cases. The smallest average mean erroribias is 

highlighted in red for each case day and the largest average mean erroribias is highlighted in blue. Overall, both the 

LAPSARW _MYJ and ADASARW MY] produced the smallest average mean errors and had the smallest bias of all 

configurations, followed by the ADASARW _ Yonsei configuration. Clearly the worst performer among all model 

configurations was the LAPSNMM_GFS which was the worst in seven out of the nine cases. 

Table 4. Average mean error and bias (parenthesis) for the I 2-hour forecast for each model configuration for the wind 

cycling and null cases. The smallest average mean error/bias is highlighted in red for each case day. 

Model Configurations 

Date LAPSARW LAPSARW ADASARW ADASARW LAPSNMM LAPSNMM 

Yonsei MYJ Yonsei MYJ GFS MYJ 

22 Dec 06 -1.96 (0.77) -2.28 (0.93) ·1.72 (1,03) -1.43 (1.09) 458 (053) -3.14 (0.75) 

30 Jan 08 -2.86 (0.79) ·1.81 (0.88) -3.43 (0.74) -2.25 (0.84) ·7,45 (0.39) -5.78 (0.54) 

14 Feb 08 -1.73 (0.86) -0.79 (0.97) -2.62 (0.71) -2.7 (0.68) -5.110,48) -3.25 (0.69) 

4/5 Mar 08 -1.23 (0.82) -0.6 (0.89) -1.71 (0.75) -1.19 (0.81) -4.96 (0.35) -3.77 (0.49) 

4 Jun 08 -3.52 (0.77) -2.57 (0.84) -2.5 (0.83) 1.I (0.95) ·7.44 (0,43) -5.69 (0.58) 

7 Jun 08 -I. I 7 (0.88) -0.75 (0.92) -0.94 (0.93) ·0.75 (0.94) -4.1 (0.54) . -2.33 (0.75) 

30/3 I Jul08 1.64 (1.23) 3.58 (1.5) 2.01 (1.28) 3.2 I (1.45) -1.68 (0.76) -0.2 (0.97) 

9 Jun 07 (null case) -0.58 (0.92) I -0.99 (0.78) 0.3 (0.99) -0.3 I (0.82) ·1.91 (0,47) -1.67 (0.53) 

I Nov 07 (null case) -0. I 8 (0.86) 0.170.12) -0.5 (066) 0.05 (1.05) -0.36 (0.77) -0.36 (0.77) 

24 
Figure 13 shows the mean error for wind direction for the same three wind cycling case days and one null case. 

In general, the model configurations that used the same cores but different physics give nearly the same mean error 

for all wind cycling and null cases. In fact, all model configurations seemed clustered throughout the forecast, 

indicating that they captured and missed all of the same shifts in wind direction. Overall, no one model 

configuration was the best or worst performer for a majority of the cases. The discrepancies between model 

configurations for mean error of wind direction were smaller than for the wind speed, indicating that the model may 

have been better at forecasting wind direction than wind speed. In the highlighted cases, the smallest mean error 

occurred on 30/31 July (Figure l3c), while the largest mean error occurred for the null case on 9 June (Figure 13d). 

50 .., 

<II 40 

~ 30 

<II 20 ::!:. 

c: 10 

0 0 .~ 

<II ·10 .. 

C ·20 

"0 

c: ·30 

§ ·40 

Mean Error - Wind Direction 

22 Dec 2006 

Forecast Hour 

-LAPSARW_Yonsei ~LAPSARW_MYJ -AOASARW_Yonsei 

ADASIIRW_MYJ ·""LAPSNMM_GfS -LAPSNMM_MYJ 

0 0 0 

0 s 0 0 

Mean Error - Wind Direction 

30/31Jul 2008 

0 

i2 

0 

Forecast Hour 

-LAPSARW_Yof1sei ~LAPSARW_MYJ ----AOASARW_Vonsei 

<5 

§ 

c 

.., 

<II 

i!! 

bD 

<II 

::!:. 

c: 

0 

.+:: 

... 

<II 

a 

"0 

c: 

§ ·zoo 

'" 

~ 

Mean Error - Wind Direction 

4/5 Mar 2008 

-LAPSARW_Yonsei -lAPSARW~MYJ --AOASARW_vonsei 

ADASARW_MYJ ~"'LAP$NMM_GfS -LAPSNMM_MYJ 

0 

0 

0 

0 

Mean Error - Wind Direction 

9Jun 2007 

0 

;;:: 

0 

Forecast Hour 

-lAPSARW_Yonsei --lAP$ARW _MYJ --==ADASARW_Yonsei 

ADASARW_MYJ ~LAPSNMM_GFS -LAPSNMM_MYJ 

d 

Figure 13. ehart of the mean error for wind direction from the 6 WRF model configurations for the 12-hour 

forecast period for a) 22 December 2006, b) 4/5 March 2008, c) 30 July 2008, and d) 9 June 2007. 

4.3 Pearson Correlation Coefficient for Wind Speed and Direction 

The pee can indicate whether the model configurations caught the overall trend of the wind speed. That is, it 

answers the question of whether the model winds fluctuated positively and negatively with the same magnitude as 

the observed winds. The closer the correlation coefficient is to I, the better the model was able to capture these 

trends. When comparing forecast vs. observed wind speed, only positive coefficients indicate any value in the model 

forecasts. 

Figure 14 shows the pee for wind speed from the 12-hour forecasts of all wind cycling and null cases 

simulated with the six model configurations. The most obvious feature is that on both 30 January 2008 and 4 June 

2008 each model configuration did a poor job at capturing the trends in the observed winds. The best model forecast 

days were 4/5 March 2008 and 30/31 July 2008. Both days had correlation coefficients of 0.6 or above for all model 

configurations. The best performer for 4/5 March was the LAPSNMM_MYJ and the LAPSNMM_GFS for 30/31 

July; however, the bias for these two configurations was poor compared to the other configurations for the same day. 

This indicates that although the model forecast wind speeds were too low, the model configuration was able to 

25 
capture the fluctuations in wind speed maximums and minimums. Disregarding the PCCs for the two worst days (30 

January 2008 and 4 June 2008), the most consistently good performers were the LAPSAR W _ Yonsei and 

LAPSNMM _ MY J configurations. 

Pearson Correlation Coefficient - Wind Speed 

1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

-0.2 

-0.4 

22 Oec 06 30 Jan 08 14 Feb 08 Mar 08 4 Jun 08 7 lun 08 30/3lJul 9 Jun 07 1 Nov 07 

... Wind Cycling Cases 08 ~ '-Null Cases-+ 

• 

Figure 14. Chart showing the Pearson Correlation Coefficient for wind speed for the 12-hour forecasts for 

all wind cycling and null cases for the six model configurations. 

For wind direction, the PCC can indicate whether the model configurations caught the overall shifts in wind 

direction. That is, it answers the question of whether the model wind direction fluctuated positively and negatively 

with the same magnitude as the observed wind direction. The closer the correlation coefficient is to 1, the better able 

the model to capture the same magnitude and direction of the observed wind shift. A correlation coefficient closer to 

-1 indicates that the model was able to capture the magnitude of the wind shift, but in the opposite direction. 

Figure 15 shows the PCC of wind direction for the 12-hour forecasts for all wind cycling and null cases for the 

six model configurations. Unlike the PCC for the wind speed forecasts, each model configuration did a poor job at 

capturing any of the trends in the observed winds on 30/31 July 2008. In fact, the average PCC for this day for all 

model configurations combined is 0, or no correlation. The model configurations did well capturing the shifts in 

wind direction for the rest of the wind cycling case days. Overall, the model was better able to capture shifts in wind 

direction for the wind cycling days than for the null cases. The best model forecast day was 7 June 2008, which had 

a correlation coefficient of 0.6 or above for all model configurations. As was found with the mean error calculations, 

there is no indication of one model configuration as clearly the best or worst performer. 

26 
1.0 

0.8 

0.6 

0.4 

0.2 

0.0 

-0.4 

Pearson Correlation Coefficient - Wind Direction 

22 Dec 06 30 Jan 08 14 Feb 08 4/5 Mar 08 4 Jun 08 7 Jun 08 30/31 Jul 9 Jun 07 1 Nov 07 

08 

Wind Cycling Cases-------~~ .-Null Cases-. 

• LMr' :>""wv 

Figure 15. Chart showing the Pearson Correlation Coefficient for wind direction for the 12-hour forecasts 

for all wind cycling and null cases for the six model configurations. 

4.4 Objective Analysis Summary 

Based on the objective analysis of all wind cycling and null case days, it appears that overall the AR W runs 

outperformed the NMM runs. Although the NMM core was more successful in capturing the increases and decreases 

in wind speed, it did a poor job at capturing the magnitude ofthe wind and consistently under-predicted wind speeds 

by the largest margin. Changing the model core seemed to have the biggest impact on the forecasts, while changing 

the model physics seemed to have the least impact on the forecasts. It does not appear the model does any better 

forecasting wind speeds for the null days vs. the wind cycling days. 

27 
5. Conclusions 

Occasionally, the space shuttle must land at EAFB in Southern California when weather conditions at KSC 

violate Flight Rules. This can create a challenge for SMG as the complex terrain in and around EAFB makes 

forecasting surface winds difficult. In particular, "wind cycling" cases, in which the wind speeds and directions 

oscillate among towers near the EAFB runway, present a challenging forecast problem for shuttle landings. An 

accurate depiction of the winds along the runway is crucial in making the landing decision. Since global and national 

scale models cannot properly resolve the wind field due to their coarse horizontal resolutions a properly tuned highresolution mesoscale model, like WRF, is needed. 

The evaluation described in this report was designed to assess the skill of different WRF model configurations 

for forecasting wind cycling events at EAFB. This study was done in two steps: 1) Initializing both cores of the 

WRF model with either the LAPS or ADAS software, while also varying some of the physical parameters, and 2) 

Assessing the accuracy of model output by comparing surface wind forecasts to local EAFB tower observations. 

SMG identified seven wind cycling case days from December 2006 to July 2008. Six of the seven days were 

considered strong wind cycling cases, while one was a marginal case. The AMU also randomly chose two null 

cases. To compare the different model configurations, the WRF model was run for 12 hours at a I-km horizontal 

grid spacing over EAFB and the surrounding mountainous areas using data from the identified case days. Data 

lru~",",4e(lj,v Ute:~~mftil,~j.l1tTftHfl'htfle !la't:a~~i'tna+VSts vaekagesinelttdedLevel ~HWSR-8-8l}data,OOES visible ~ ~nd ~~ 

infrared satellite imagery, MADIS data, and local EAFB wind tower data. Each candidate day consisted of six runs 

using different model configurations, for a total of 42 runs for the wind cycling case days and 14 for the null cases. 

The different configurations included: 

• LAPS-ARW with the Yonsei University PBL scheme and MM5 similarity surface layer scheme, 

• LAPS-ARW with the Mellor-Yamada-Janjic PBL scheme and ETA similarity surface layer scheme, 

• LAPS-NMM with the NCEP GFS PBL scheme and NCEP GFS surface layer scheme, 

• LAPS-NMM with the Mellor-Yamada-Janjic PBL scheme and ETA similarity surface layer scheme, 

• ADAS-ARW with the Yonsei University PBL scheme and MM5 similarity surface layer scheme, and 

• ADAS-ARW with the Mellor-Yamada-Janjic PBL scheme and ETA similarity surface layer scheme. 

5.1 Results 

The major results from the evaluation show the following: 

• Subjective analysis, 

o Overall, the ARW runs outperformed the NMM runs. 

o Changing the model core seemed to have the biggest impact on the forecasts, while changing the 

model physics seemed to have the least impact on the forecasts. 

o Model configurations that used the MYJ PBL scheme seemed to slightly outperform other PBL 

schemes. 

o The model was able to differentiate between wind cycling days and null cases. 

• Objective analysis, 

o Overall the AR W runs outperformed the NMM runs. 

o The NMM core was more successful in matching the increasing and decreasing trends in wind speed 

with the observations, but did a poor job at capturing the magnitude of the wind and consistently 

under-predicted wind speeds by the largest margin. 

o As was found in the subjective analysis, changing the model core seemed to have the biggest impact 

on the forecasts, while changing the model physics seemed to have the least impact on the forecasts. 

o It does not appear that the model did any better forecasting wind speeds for the null days vs. the wind 

cycling days. 

o Model configurations that used the MYJ PBL scheme seemed to slightly outperform other PBL 

schemes. 

28 
In both the subjective and objective analyses, the AMU found the ARW core and MYJ PBL scheme performed 

better than the other model configurations. The model did not produce a better wind forecast on null days vs. wind 

cycling days. The model was able to differentiate between wind cycling days and null cases which would provide 

added value to the shuttle landing forecast. 

5.2 Recommendations 

Of the six model configurations tested, the AMU recommends either the ADAS ARW MYJ or LAPS ARW 

MY] configuration for operational use for predicting wind cycling events at EAFB for the following reasons: 

• Both configurations consistently had the lowest bias for wind speed compared to the other model 

configurations, 

• Both configurations best captured wind speed oscillations when compared to the observed wind speeds at 

the concrete runway towers, and 

• The NMM core produced wind speed forecasts well below those observed and did not capture the wind 

direction changes well on the cycling days. 

5.,1 . Euture.WodL. 

The goal of this work was to assess different WRF model options and to determine how well the model could 

predict surface wind speed and direction at EAFB and if one model configuration performed better than others. 

While that goal was met, the data indicated further investigations outside the scope of this task could include 

sophisticated signal analysis, rotating the model coordinate systems orthogonal and/or parallel to the prevailing wind 

flow and closely comparing the observed and forecast wind speed cycles to wind direction cycles which may lead to 

discoveries of why the winds cycle as they do and how that relates to model configuration and performance. 

29 
References 

Albers, S. c., 1995: The LAPS wind analysis. Wea. Forecasting, 10,342-352. 

Albers, S. C., 1. A. McGinley, D. L. Birkenheuer, and J. R. Smart, 1996: The Local Analysis and Prediction System 

(LAPS): Analyses of clouds, precipitation, and temperature. Wea. Forecasting, 11, 273-287. 

Birkenheuer, D., 1999: The effect of using digital satellite imagery in the LAPS moisture analysis. Wea. 

Forecasting, 14,782-788. 

Brewster, K., 1996: Implementation of a Bratseth analysis scheme including Doppler radar. Preprints, 15th Can! on 

Weather Analysis and Forecasting, Norfolk, VA, Amer. Meteor. Soc., 92-95. 

Brewster, K., 2002: Recent advances in the diabatic initialization of a non-hydrostatic numerical model. Preprints, 

15th Can[ on Numerical Weather Prediction and 21st Conf on Severe Local Storms, San Antonio, TX, Amer. 

Meteor. Soc., J6.3. 

Case, J. L., 1. Manobianco, T. D. Oram, T. Garner, P. F. Blottman, and S. M. Spratt, 2002: Local Data Integration 

........ Qyer.Easl·CeDtraIElQridallsinglhe.ARrS~DalaAnalysjsSYSlellLWgg .. EQ!:.e."Cl.s.til1g,Jl,J~2§.L 

Chen, F., and J. Dudhia, 2001: Coupling an advanced land-surface/ hydrology model with the Penn State/ NCAR 

MM5 modeling system. Part I: Model description and implementation. Man. Wea. Rev., 129, 569-585. 

Chou M.-D., and M. 1. Suarez, 1994: An efficient thermal infrared radiation parameterization for use in general 

circulation models. NASA Tech. Memo. 104606, 3, 85pp. 

Edison International, 2008: Southern California Edison Signs Largest Wind Energy Contract in U.S. Renewable 

Industry History. [Available online at =='-'.:..~"-,=c:.===~'J 

Fels, S. B., and M. D. Schwarzkopf, 1975: The simplified exchange approximation: A new method for radiative 

transfer calculations. J Atmos. Sci., 32, 1475-1488. 

Ferrier, B. S., 1994: A double-moment mUltiple-phase fourclass bulk ice scheme. Part I: Description. J Atmos Sc. i, 

51,249-280. 

Fulton, R. A., J. P. Breidenbach, D. Seo, D. A. Miller, T. O'Bannon, 1998: The WSR-88D rainfall algorithm. Wea. 

Forecasting, 13,377-395. 

Grasso, L. D., 2000: Comments on "The Differentiation between Grid Spacing and Resolution and Their 

Application to Numerical Modeling" Reply. BAMS, 81, 2479-2479. 

Janjic, Z. I., 1984: Non-linear advection schemes and energy cascade on semi-staggered grids. Man. Wea. Rev., 112, 

1234-1245. 

Janjic, Z. I., 2003a: A nonhydrostatic model based on a new approach. Meteor. Atmos. Phys., 82,271-285. 

Janjic, Z. I., 2003b: The NCEP WRF core and further development of its physical package. 5th Int. SRNWP 

Workshop on Non-Hydrostatic Modeling, Bad Orb, Germany, 27-29 October. 

Janjic, Z., I., J. P. Gerrity Jr., and S. Nickovic, 2001: An alternative approach to nonhydrostatic modeling. Man. 

Wea. Rev., 129, 1164-1178. 

Lacis, A. A., and J. E. Hansen, 1974: A parameterization for the absorption of solar radiation in the earth's 

atmosphere. J Atmos. Sci., 31,118-133. 

Lin, Y.-L., R. D. Farley, and H. D. Orville, 1983: Bulk parameterization of the snow field in a cloud model. J 

Climate Appl. Meteor., 22, 1065-1092. 

30 
McGinley, 1. A., 1995: Opportunities for high resolution data analysis, prediction, and product dissemination within 

the local weather office. 14th Con! on Weather Analysis and Forecasting, Dallas, TX, Amer. Meteor. Soc., 

478-485. 

McGinley, 1. A., S. C. Albers, and P. A. Stamus, 1991: Validation of a composite convective index as defined by a 

real-time local analysis system. Wea. Forecasting, 6,337-356. 

Mlawer, E. 1., S. J. Taubman, P. D. Brown, M. 1. Iacono, and S. A. Clough, 1997: Radiative transfer for 

inhomogeneous atmosphere: RRTM, a validated correlated-k model for the longwave. J Geophys. Res., 102 

(D14), 16663-16682. 

Schwarzkopf, M. D., and S. B. Fels, 1985: Improvements to the algorithm for computing C02 transmissivities and 

cooling rates. J Geophys. Res., 90, 541-550. 

Schwarzkopf, M. D., and S. B. Fels, 1991: The simplified exchange method revisited: An accurate, rapid method 

for computations of infrared cooling rates and fluxes. J Geophys. Res., 96, 9075-9096. 

Skamarock, W. c., 1. B. Klemp, 1. Dudhia, D. O. Gill, D. M. Barker, W. Wang, and J. G. Powers, 2005: A 

......... descriptioUQf theAdy.~ncedRe.se~J:!:;h . .wR.f Y e!,sjQI)) .. NC::;AR Te.ch.nicaINote.NC:AR[[N::1~§:-~STJ3.y .8.~pp ... 

Walters, M. K., 2000: Comments on "The Differentiation between Grid Spacing and Resolution and Their 

Application to Numerical Modeling" Reply. BAMS, 81, 2475-2477. 

Zhang, 1., F. H. Carr, and K. Brewster, 1998: ADAS cloud analysis. Preprints, 12th Conf on Numerical Weather 

Prediction, Phoenix, AZ, Amer. Meteor. Soc., 185-188. 

31 
List of Acronyms 

ADAS Advanced Regional Prediction System MM5 NCAR/Penn State Mesoscale Model 

(ARPS) Data Analysis System Version 5 

AMU Applied Meteorology Unit MSE Mean Square Error 

ARW Advanced Research WRF MYJ Mellor-Yamada-Janjic 

EAFB Edwards Air Force Base NAM North American Mesoscale model 

FR Flight Rules NCAR National Center for Atmospheric Research 

GFS Global Forecast Systems NCEP National Centers for Environmental 

GOES Geostationary Operational Environmental Prediction 

Satellites NetCDF Network Common Data Form 

GSD Global Systems Division NMM Non-hydrostatic Mesoscale Model 

KSC Kennedy Space Center PBL Planetary Boundary Layer 

kt Knots PCC Pearson Correlation Coefficient 

LAPS Local Analysis and Prediction System POR Period of Record 

MADIS Meteorological Assimilation Data Ingest RUC Rapid Update Cycle 

System SMG Spaceflight Meteorology Group 

McIDAS Man computer Interactive Data Access SSEC Space Science and Engineering Center System 

MET Model Evaluation Tools WRF Weather Research and Forecasting model 

MIDDS Meteorological Interactive Data Display WSR-88D Weather Surveillance Radar-l 988 Doppler 

System 

32 
NOTICE 

Mention of a copyrighted, trademarked or proprietary product, service, or document does not constitute endorsement 

thereof by the author, ENSCO Inc., the AMU, the National Aeronautics and Space Administration, or the United 

States Government. Any such mention is solely for the purpose of fully informing the reader of the resources used to 

conduct the work reported herein. 

33 
REPORT DOCUMENTATION PAGE Form Approved 

OMS No. 0704-0188 

The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing 

data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information, Send comments regarding this burden estimate or 

any other aspect of this collection of information, including suggestions for reducing this burden, to Department of Defense, Washington Headquarters Services, Directorate 

for Information Operations and Reports (0704-0188),1215 Jefferson Davis Highway, Suite 1204, Arlinwon, VA 22202-4302, Respondents should be aware that 

notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comp y with a collection of information if it does not display a currentiy 

~~id EA ~~'b~~~! T n~'1~r~1 RUN Y R (OUI Fe ORM T THE ABOVE ADDRESS. 

1. REPORT DATE (DD-MM-VYYY) 2. REPORT TYPE 3. DATES COVERED (From - To) 

30-12-2008 Final January 2008 - December 2008 

4. TITLE AND SUBTITLE 5a. CONTRACT NUMBER 

Weather Research and Forecasting Model Wind Sensitivity Study at NNK06MA70C 

Edwards Air Force Base, CA 5b. GRANT NUMBER 

5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 5d. PROJECT NUMBER 

Leela R. Watson 

William H. Bauman III 5e. TASK NUMBER 

5f. WORK UNIT NUMBER 

"',' r' " , 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION 

ENSCO, Inc. 

REPORT NUMBER 

1980 N. Atlantic Ave Suite 230 

Cocoa Beach, FL 32931 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING/MONITOR'S ACRONYM(S) 

NASA 

John F, Kennedy Space Center 

Code KT-C-H 11. SPONSORING/MONITORING 

Kennedy Space Center, FL 32899 REPORT NUMBER 

CR-2008-214755 

12, DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified, Unlimited 

13. SUPPLEMENTARY NOTES 

An electronic version can be found at http://science.ksc.nasa,gov/amu/final.html 

14. ABSTRACT 

NASA prefers to land the space shuttle at Kennedy Space Center (KSC). When weather conditions violate Flight Rules at KSC, NASA will 

usually divert the shuttle landing to Edwards Air Force Base (EAFB) in Southern California. But forecasting surface winds at EAFB is a 

challenge for the Spaceflight Meteorology Group (SMG) forecasters due to the complex terrain that surrounds EAFB, One particular phenomena 

identified by SMG is that makes it difficult to forecast the EAFB surface winds is called "wind cycling", This occurs when wind speeds and 

directions oscillate among towers near the EAFB runway leading to a challenging deorbit bum forecast for shuttle landings. The large-scale 

numerical weather prediction models cannot properly resolve the wind field due to their coarse horizontal resolutions, so a properly tuned 

high-resolution mesoscale model is needed. The Weather Research and Forecasting (WRF) model meets this requirement. The AMU assessed the 

different WRF model options to determine which configuration best predicted surface wind speed and direction at EAFB, To do so, the AMU 

compared the WRF model performance using two hot start initializations with the Advanced Research WRF and Non-hydrostatic Mesoscale 

Model dynamical cores and eompared model performance while varying the physics options, 

15. SUBJECT TERMS 

Surface Winds, Wind Speed, Wind Direction, Numerical Weather Prediction, NWP, Modeling, Weather Research and Forecasting 

model, WRF, Edwards AFB, Shuttle Landing, Deorbit Burn Forecast, Space Shuttle, Weather, Meteorology, wind cycling, Local 

Analysis and Prediction System (LAPS), Advanced Regional Prediction System (ADAS) 

16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF 18. NUMBER 19a. NAME OF RESPONSIBLE PERSON 

ABSTRACT OF 

a. REPORT b.ABSTRACT c. THIS PAGE PAGES Jennifer G. Wilson 

U UU UU UU 34 

19b. TELEPHONE NUMBER (Inc/ude area code) 

(321) 867-0824 

Standard Form 298 (Rev, 8-98) Prescnbed by ANSI Std, Z39-18