A new set of basis functions was constructed using the Rossby-Haurwitz waves, which are the eigenfunctions of nondivergent barotropic vorticity equations on the sphere. The basis functions were designed to be non-separable, that is, not factored into functions of either the longitude or the latitude. Due to this property, the nodal lines of the functions are aligned neither along with the meridian nor the parallel. The basis functions can be categorized into groups of which members have the same degree or the total wavenumber-like index on the sphere. The orthonormality of the basis functions were found to be close to the machine roundoffs, giving the error of O(10-15) or O(10-16) for double-precision computation (64 bit arithmetic). It was demonstrated through time-stepping procedure that the basis functions were also the eigenfunctions of the non-divergent barotropic vorticity equations. The projection of the basis functions was carried out onto the low-resolution geopotential field of Gaussian bell, and compared with the theory. The same projections were performed for the observed atmospheric-geopotential height field of 500 hPa surface to demonstrate decomposition into the fields that contain disturbance of certain range of horizontal scales. The usefulness of the new basis functions was thus addressed for application to the eigenmode analysis of the atmospheric motions on the global domain.

Spherical harmonics power spectrum of the geopotential field of Gaussian-bell type on the sphere was investigated using integral formula that is associated with Legendre polynomials. The geopotential field of Gaussian-bell type is defined as a function of sine of angular distance from the bell's center in order to guarantee the continuity on the global domain. Since the integral-formula associated with the Legendre polynomials was represented with infinite series of polynomial, an estimation method was developed to make the procedure computationally efficient while preserving the accuracy. The spherical harmonics power spectrum was shown to vary significantly depending on the scale parameter of the Gaussian bell. Due to the accurate procedure of the new method, the power (degree variance) spanning over orders that were far higher than machine roundoff was well explored. When the scale parameter (or width) of the Gaussian bell is large, the spectrum drops sharply with the total wavenumber. On the other hand, in case of small scale parameter the spectrum tends to be flat, showing very slow decaying with the total wavenumber. The accuracy of the new method was compared with theoretical values for various scale parameters. The new method was found advantageous over discrete numerical methods, such as Gaussian quadrature and Fourier method, in that it can produce the power spectrum with accuracy and computational efficiency for all range of total wavenumber. The results of present study help to determine the allowable maximum scale parameter of the geopotential field when a Gaussian-bell type is adopted as a localized function.

기울어진 자전축을 갖는 회전계에서, 일정한 각속도로 회전하는 동서풍이 있는 경우에 대해서 로스비하우어비츠 파동의 섹터모드(적도에 대한 반구 비대칭의 첫 번째 모드)와 균형을 이루는 지위고도장을 해석적으로 유도하였다. 균형장은 발산방정식으로부터 시간변화를 제거하고 라플라시안 연산자를 역산함으로써 구하였다. 역산은 비선형항의 계산과 포이슨 방정식의 해를 구하는 두 단계의 연산과정으로 이루어져 있다. 두 번째 단계에서, 구면조화함수로 표현되는 강제력의 항은 구면조화함수의 선형관계를 이용하였고, 그 이외의 항은 구면조화함수를 적분함으로써 구하였다. 균형장은 여섯 개의 동서파수 성분으로 표현됨이 드러났다. 본 연구에서 구한 균형장은 적도에 대하여 비대칭의 구조를 가지기 때문에, 대칭의 구조만을 가지는 것에 비하여 미분방정식의 수치해의 검종법으로서의 활용도가 높다. 일정한 각속도를 갖는 배경 동서풍이 지구의 자전각속도와 같거나 1/2에 해당하는 경우에는, 일부 동서파수 성분이 제거되는 것으로 나타났다. 이론적으로 구한 균형장은 정교한 수치모델을 통하여 구한 균형장과 거의 정확하게 같은 것으로 밝혀져, 이론적 해의 타당성이 입증되었다. 마지막으로, 로스비하우어비츠 파동의 섹터모드와 균형을 이루는 지위고도장의 안정성을 장기간시간적분을 통하여 살펴보았다.