## Info

Figure 9. SCAT analysis of ocean surface winds (m/s) at 40 m over the southern Asia—Pacific ocean for 0000 UTC May 2, 1994.

winds retrieved from the passive microwave radiometer of SSM/I on board the DMSP satellites, and from active microwave scat-terometer measurements on board SEASAT, NSCAT, ERS-1/2 and QuikSCAT satellites were discussed.

Of special import are the summaries of pre-implementation test results on the use of these various satellite ocean surface wind data at NCEP global data assimilation systems, based on which these surface wind data have been operationally implemented at various timelines during the last two decades. It is fair to state that from gross statistics based on many cases of forecasts, the impact of the satellite ocean surface winds on NWP short range forecasts is mostly positive, small but significant because the data are of a single level nature and only available at the ocean surface. Large impacts of the satellite ocean surface winds may only occur over some special synoptic situations and over the areas of severe weather storms when conventional data are lacking. This is well demonstrated in a special case investigation over the southwestern Pacific Ocean, where ERS-1 scatterometer winds are found to be most significant in identifying the storm center intensity and circulation.

### Acknowledgements

The author would like express his sincere appreciation to Dr Stephen Lord of NCEP and Dr A. Hollingsworth of ECMWF for permission to use their centers' forecast performance statistics on anomaly correlations as shown in Figs. 1 and 2 of this article. He is also grateful to the two anonymous reviewers for their constructive comments and suggestions, which improved the final version of the manuscript.

[Received 15 January 2006; Revised 3 May 2007; Accepted 2 June 2007.]

### References

Atlas, R., R. Hoffman, S. Leider, J. Sienkiewicz, T.-W. Yu, S. Bloom, E. Brin, J. Ardizzone, J. Terry, D. Bungato, and J. Jusem, 2001: The effects of marine winds from scatterometer data on weather analysis and forecasting. Bull. Amer. Meteorol. Soc., 82, 1965-1990.

Atlas, R., S. Bloom, J. Ardizzone, E. Brin, J. Terry, and T.-W. Yu, 2001: The impact of Quikscat on weather analysis and forecasting, Preprint paper presented at 18th Conference on Weather Analysis and Forecasting and 14th Conference on Numerical Weather Prediction, 30 July-2 August 2001, Fort Lauderdale, Florida, USA, pp. 172175.

Atlas, R., S. Bloom, R. Hoffman, E. Brin, J. Ardizzone, J. Terry, D. Bungato, and J. Jusem, 1999: Geophysical validation of NSCAT winds using atmospheric data and analysis. J. Geophys. Res., 104 (C5), 11405-11424.

Andrew, P. L., and R. S. Bell, 1998: Optimizing the United Kingdom Meteorological Office data assimilation for ERS-1 scatterometer results. Quart. J. Roy. Meteorol. Soc., 124, 1831-1860.

Baker, W., R. Atlas, E. Kalnay, M. Halem, P. Woicesyhn, S. Peteherych, and D. Edelmann,

1984: Large-scale analysis and forecast experiments wind data from Seasat-A scatterometer. J. Geophys. Res., 89 (D3), 4927-4936.

Courtier, P., E. Anerson, W. Heckley, J. Paileux, D. Vasiljevic, M. Hamrud, A. Holligsworth, F. Rabier, and M. Fisher, 1998: The ECMWF implementation of three-dimensional variational assimilation (3D-Var). I: Formulation. Quart. J. Roy. Meteor. Soc., 124, 1783-1807.

Cressman, G. P., 1959: An operational objective analysis scheme. Mon. Wea. Rev., 87, 367-374.

Derber, J., D. Parrish, and S. Lord, 1991: The new global operational analysis system at the National Meteorological Center. Wea. Forecast., 6, 538-547.

Duffy, D., and R. Atlas, 1986: The impact of Seasat-A scatterometer winds on the numerical prediction of the Queen Elizabeth II storm. J. Geophys. Res., 91 (C2), 2241-2248.

Gandin, L. S., 1965: Objective Analysis of Meteorological Fields, Israel Program for Scientific Translation, Jerusalem, 242 pp.

Gemmill, W. H., T.-W. Yu, V. Krasnopolsky, C. Peters, and P. Woiceshyn, 1999: NCEP experience with "Real-time" ocean surface wind retrievals from satellites. NCEP Technical Note, OMB Contribution No. 164, 26 pp.

Goodberlet, M. A., C. T. Swift, and J. Wikerson, 1989: Remote sensing of ocean surface winds with the Special Sensor Microwave/Imager. J. Geophys. Res., 94, 547- 555.

Hoffman, R. N., 1993: A preliminary study of the impact of the ERS1 C-band scatterometer wind data on the ECMWF global data assimilation system. J. Geophys. Res., 98 (C6), 10233-10244.

Hoffman, R. N., and C. Grassotti, 1996: A technique for assimilating SSM/I observations of marine atmospheric storms: tests with ECMWF analysis. J. Appl. Meteorol., 35, 11771187.

Hollingsworth, A. (2004): Performance of ECMWF global forecasts in terms of 500 mb geopotential height anomaly correlation statistics. Personal communication.

Ingleby, N. B., and R. A. Blomley, 1991: A diagnostic study of the impact of SEASAT scatte-rometer winds on numerical weather prediction. Mon. Wea. Rev., 119, 84-103.

Krasnopolsky, V., L. Breaker, and W. Gemmill, 1995a: A neural network as a nonlinear transfer function model for retrieving surface wind speeds from the Special Microwave Imager. J. Geophys. Res., 100, 11033-11045.

Le DiMet, F., and O. Talagrand, 1986: Variational algorithms for analysis and assimilation of meteorological observations: Theoretical aspects. Tellus, 38A, 97-110.

Lord, S. (2004): Performance of NCEP global forecasts in terms of 500 mb geopotential height anomaly correlation statistics. Personal communication.

Parrish, D., and J. C. Derber, 1992: The National Meteorological Center's Spectral Statistical Interpolation analysis scheme. Mon. Wea. Rev., 120, 1747-1763.

Rabier, F., A. McNally, E. Anderson, P. Courtier, P. Unden, J. Eyer, A. Hollingsworth, and F. Bouttier, 1998: The ECMWF implementation of three-dimensional varaitional assimilation (3D-Var). II: Structure function. Quart J. Roy. Meteorol. Soc., 124, 1809-1829.

Rabier, F., J.-F. Mahfouf, M. Fisher, H. Haseler, A. Holligsworth, L. Isaksen, E. Klinker, S. Saarinen, C. Temperton, J.-N. Thepaut, P. Unden, and D. Vasiljevic, 1997: The ECMWF operational implementation of four-dimensional variational assimilation. ECMWF Research Development Tech. Memo. No. 240.

Sasaki, Y., 1970: Some basic formalisms in numerical variational analysis. Mon. Wea. Rev., 98, 875-883.

Sela, J. G., 1980: Spectral modeling at the National Meteorological Center. Mon. Wea. Rev., 108, 1279-1292.

Stoffelen, A., and G. J. Cats, 1991: The impact of Seasat-A scatterometer data on high-resoltuion analysis and forecasts: the development of the QE II storm. Mon. Wea. Rev., 119, 2794-2802.

Stoffelen, A., and D. Anderson, 1997: Ambiguity removal and assimilation of scatterometer data. Quart. J. Roy. Meteor. Soc., 123, 491-518.

Thepaut, J.-N., R. N. Hoffman, and P. Courtier, 1993: Interactions of dynamics and observations in a four-dimensional variational assimilation. Mon. Wea. Rev., 3393-3414.

Yu, T.-W., and R. D. McPherson (1984): Global data assimilation experiments using SESAT-A scatterometer wind data. Mon. Wea. Rev.

Yu, T.-W., 1987: A technique for deducing wind direction from satellite microwave measurements of wind speed. Mon. Wea. Rev., 115, 1929-1939.

Yu, T.-W., and D. Deaven, 1991: Use of SSM/I wind speed data at NMC's GDAS. Preprints,

9th Conference on Numerical Weather Prediction, Denver, CO, USA, American Meteorological Society, pp. 416-417.

Yu, T.-W., W. Gemmill, and J. Wollen, 1992: Use of SSM/I wind speed in the operational numerical weather prediction system at NMC. Abstracts, NMC/NESDIS/DOD Conference on DMSP Retrieval Products, Washington DC, USA.

Yu, T.-W., P. Woiceshyn, W. Gemmill, and C. Peters, 1994: Analysis and forecast experiments using ERS-1 scatterometer wind measurements. Preprint paper presented at 7th Conference on Satellite Meteorology and Oceanography of the American Meteorological Society, June 6-10, Monterey, California, USA, pp. 600601.

Yu, T.-W., and J. C. Derber, 1995: Assimilation experiments with ERS-1 winds: Part I — Use of backscatter measurements in the NCEP spectral statistical analysis system. NCEP/EMC/Ocean Modeling Branch Technical Note 116, 32 pp.

Yu, T.-W., 1995: Assimilation experiments with ERS-1 winds: Part II — Use of vector winds in NCEP spectral statistical analysis system. NCEP/EMC/Ocean Modeling Branch Technical Note 117, 29 pp.

Yu, T.-W., M. D. Iredell, and Y. Zhu, 1996: The impact of ERS-1 winds on NCEP numerical weather analysis and forecasts, Preprint paper presented at 11th Conference on Numerical Weather Prediction of the American Meteorological Society, August 19-23, 1996, Norfolk, Virginia, USA, pp. 276-277.

Yu, T.-W., M. Iredell, and D. Keyser, 1997: Global data assimilation and forecast experiments using SSM/I wind speed data derived from a neural network algorithm. Wea. Forecast., 12, 859-865.

Yu, T.-W., 2000: Assimilating satellite observations of ocean surface winds in the NCEP global data and regional atmospheric analysis systems. Proceedings of Third WMO Symposium on Assimilations of Observations in Meteorology and Oceanography, WMO/TD-No. 986, pp. 229-232.

Yu, T.-W., 2003: Operational use of Quikscat winds in NCEP GDAS. Preprint paper at 12th Conference on Satellite Meteorology of the American Meteorological Society, February 9-13, Long Beach, CA, USA.

### Purdue Atmospheric Models and Applications

Wen-Yih Sun*'**, Wu-Ron Hsu^**, Jiun-Dar Chern*, Shu-Hua Chen$'**, Ching-Chi Wut, Kate Jr-Shiuan Yang*, Kaosan Yeh*, Michael G. Bosilovich*, Patrick A. Haines^, Ki-Hong Min", Tae-Jin Oh*'**, B. T. MacCall*'**, Ahmet Yildirim*, Yi-Lin Chang*, Chiao-Zen Chang*

and Yi-Chiang Yut

* Purdue University, West Lafayette, Indiana, USA t National Taiwan University, Taipei, Taiwan *Goddard Space Fight Center/NASA, Greenbelt, Maryland, USA $ University of California, Davis, California, USA ^ Army Research Laboratory, White Sand, New Mexico, USA "Valparaiso University, Valparaiso, Indiana, USA **Taiwan Typhoon and Flood Research Institute, Hsin-chu, Taiwan [email protected]

This article summarizes our research related to geofluid dynamics and numerical modeling. In order to have a better understanding of the motion in the atmosphere, we have been working on various forms of the Navier-Stokes equations, including the linearized and nonlinear systems as well as turbulence parametrization, cumulus parametrization, cloud physics, soil-snow parametrization, atmospheric chemistry, etc. We have also been working on numerical methods in order to solve the equations more accurately. The results show that many weather systems in the initial/growing stage can be qualitatively described by the linearized equations; on the other hand, many developed weather phenomena can be quantitatively reproduced by the nonlinear Purdue Regional Climate Model, when the observational data or reanalysis is used as the initial and lateral boundary conditions. The model can also reveal the detailed structure and physics involved, which sometimes can be misinterpreted by meteorologists according to the incomplete observations. However, it is also noted that systematic biases/errors can exist in the simulations and become difficult to correct. Those errors can be caused by the errors in the initial and boundary conditions, model physics and parametrizations, or inadequate equations or poor numerical methods. When the regional model is coupled with a GCM, it is required that both models should be accurate so as to produce meaningful results. In addition to the Purdue Regional Climate Model, we have presented the results obtained from the nonhydrostatic models, the one-dimensional cloud model, the turbulence-pollution model, the characteristic system of the shallow water equations, etc. Although the numerical model is the most important tool for studying weather and climate, more research should be done on data assimilation, the physics, the numerical method and the mathematic formulation in order to improve the accuracy of the models and have a better understanding of the weather and climate.

### 1. Introduction

The motion of the atmosphere and ocean can be represented by Navier-Stokes equations, which should be solved numerically. However, the evolution of a weather/climate system consists of motions on many different scales. At the beginning, some of the small disturbances may be described by the linearized equations. When disturbances grow, nonlinear equations become necessary. It is also noted that eigen values/vectors can be easily obtained and interpreted in a linearized system. Hence, for the past few decades, we have been working on linear instability and nonlinear numerical models to study meteorological phenomena, ranging from convection, turbulence, air pollution, cumulus clouds, cloud streets, symmetric instabilities, mountain waves, lee vortices, and land-sea breezes, to synoptic scale waves, barotropic instabilities, cyclones, fronts, and regional climate in East Asia and North America. In order to obtain accurate numerical results, we have also been developing aspects of the model such as the new diffusion equation; turbulence parametrization; snow-vegetation-soil and snow-sea ice packages; forward-backward, advection, and semi-Lagrangian schemes; the pressure gradient force; the multigrid method; the transport of dust and trace gases; and the interaction between aerosols and regional climate. The important components of our mesoscale model are shown in Fig. 1.

We developed the Purdue Regional Climate Model (PRCM). However, several mesoscale regional models exist, with varying approaches. They include the fifth-generation Penn State/ NCAR Mesoscale model (MM5; Grell et al, 1995) and the Weather Research and Forecast (WRF) model (Skamarock et al, 2005). Table 1 shows the important features that are used in the PRCM, MM5, and WRF.

A few topics will be briefly discussed here. Many of them require further study. The basic structure of the PRCM and applications will be discussed in Sec. 2, the nonhydrostatic models in Sec. 3, other topics in Sec. 4, and a summary is provided at the end.

2. Purdue Regional Climate Model

2.1. Basic equations

The PRCM is a hydrostatic primitive equation model that utilizes the terrain-following coordinate (ap) in the vertical direction. It has

Purdue Regional Climate Model

Snow-vegetation-soil Sea-ice-mixed layer ocean

( Sun and Wu 1992, Bosilovich and Sun 1995, Sun and Chern 1998,2005 )

Lin et al. (1983) Rutledge and Hobbs (1983) Chern (1994)

Kuo (1965, 1974) Anthes (1978) Molinari (1982) Sun and Haines (1996)

Purdue Regional Climate Model

Liou et al. (1988) Chou and Suarez (1994) Chouetal. (1001)

Chern (1994)

Cumulus Cloud

Dynamics

Hydrostatic Arakawa C Grid Forward-Backward scheme 4th-order advection Scheme Pressure gradient force (Sun and Hsu 1988, Sun 1993c, 1995)

Liou et al. (1988) Chou and Suarez (1994) Chouetal. (1001)

Chern (1994)

Cumulus Cloud

Dynamics

Hydrostatic Arakawa C Grid Forward-Backward scheme 4th-order advection Scheme Pressure gradient force (Sun and Hsu 1988, Sun 1993c, 1995)

Figure 1. The schematic diagram of Purdue Regional Climate Model.

Table 1. The characteristics of PRCM, WRF, and MM5. WRF has three different dynamical cores. The one that mentioned here is Advanced Research WRF (ARW).

## Post a comment