• Journal of the Korean earth science society
• /
• v.37 no.2
• /
• pp.107-116
• /
• 2016

Stability of the barotropic Rossby-Haurwitz wave is investigated using the numerical models on the global domain. The Rossby-Haurwitz wave under investigation is composed of the basic zonal flow of super-rotation and a finite amplitude spherical harmonic wave. The Rossby-Haurwitz wave is given as either steady or unsteady wave by adjusting the strength of the super-rotating zonal flow. Stability as well as the growth rate of the wave in the numerical simulation is determined by comparing the perturbation amplitude at two different time stages. Unstable modes of the Rossby-Haurwitz wave exhibited a horizontal structure composing of various zonal-wavenumber components. The vorticity perturbation for some modes showed a discontinuity around the area of weak flow, which was found robust regardless of the horizontal resolution of the model. Fourier finite element model was shown to generate the unstable mode in earlier stage of the time integration due to less accuracy compared to the spherical harmonic spectral model. Taking the overall accuracy of the models into consideration, the time by which the unstable mode begin to dominate over the spherical harmonic wave was estimated.

• Journal of the Korean earth science society
• /
• v.34 no.5
• /
• pp.393-401
• /
• 2013

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.

• Korean Journal of Remote Sensing
• /
• v.6 no.1
• /
• pp.25-37
• /
• 1990

Methods of geometric correction for the Advanced Very High Resolution Radiometer imagery of NOAA satellites were developed and applied to the software for image processing of meteorological satellite data. The software for finding the earth location of each scan position and the software for gridding on original imagery were dedigned. On the assumption of circular orbits and the spherical earth, the methods developed were sufficiently accurate in the purpose of most meteorological data analyses.

• Atmosphere
• /
• v.31 no.1
• /
• pp.113-141
• /
• 2021

The impacts of ice clouds on the energy budget of the Earth and their representation in climate models have been identified as important and unsolved problems. Ice clouds consist almost exclusively of non-spherical ice crystals with various shapes and sizes. To determine the influences of ice clouds on solar and infrared radiation as required for remote sensing retrievals and numerical models, knowledge of scattering and microphysical properties of ice crystals is required. A conventional method for representing the radiative properties of ice clouds in satellite retrieval algorithms and numerical models is to combine measured microphysical properties of ice crystals from field campaigns and pre-calculated single-scattering libraries of different shapes and sizes of ice crystals, which depend heavily on microphysical and scattering properties of ice crystals. However, large discrepancies between theoretical calculations and observations of the radiative properties of ice clouds have been reported. Electron microscopy images of ice crystals grown in laboratories and captured by balloons show varying degrees of complex morphologies in sub-micron (e.g., surface roughness) and super-micron (e.g., inhomogeneous internal and external structures) scales that may cause these discrepancies. In this study, the current idealized models representing morphologies of ice crystals and the corresponding numerical methods (e.g., geometric optics, discrete dipole approximation, T-matrix, etc.) to calculate the single-scattering properties of ice crystals are reviewed. Current problems and difficulties in the calculations of the single-scattering properties of atmospheric ice crystals are addressed in terms of cloud microphysics. Future directions to develop physically consistent ice-crystal models are also discussed.

• Proceedings of the KSEEG Conference
• /
• 2000.04b
• /
• pp.308-310
• /
• 2000

Sampling rates become inconsistent when spatial data in spherical coordinate are re-sampled with respect latitudinal or longitudinal degree for mathematical processes such as Fourier Transfrom, and this results in the distrtions of the processed data in the wavenmber domain. This distortions are more evident in the polar regions. An example is presented to show such distortions during the recovery process of free-air gravity anomalies from ERSI radar atimeter data in Russian Arctic Barents Sea, and a method is present to minimize the distortion using Lambect Conformal Conic map projection.

• Atmosphere
• /
• v.20 no.4
• /
• pp.483-494
• /
• 2010

The vertical profiles of the extinction coefficient, the backscatter coefficient, and the lidar ratio (i.e., extinction-to-backscattering ratio) for Asian dust and pollution aerosols are determined from Raman (inelastic) and elastic backscatter signals. The values of lidar ratios during two polluted days is found between 52 and 82 sr (July 22, 2009) and 40~60 sr (July 31, 2009) at 52 nm, with relatively low value of particle depolarization ratio (<5%) and high value of sun photometer-derived Angstrom exponent (> 1.2). However, lidar ratios between 25 and 40 sr are found during two Asian dust periods (October 20, 2009 and March 15, 2010), with 10~20% of particle depolarization ratio and the relatively low value of sun photometer-derived Angstrom exponent (< 0.39). The lidar ratio, particle depolarization ratio and color ratio are useful optical parameter to distinguish non-spherical coarse dust and spherical fine pollution aerosols. The comparison of aerosol extinction profiles determined from inelastic-backscatter signals by the Raman method and from elastic-backscatter signals by using the Fernald method with constant value of lidar ratio (50 sr) have shown that reliable aerosol extinction coefficients cannot be determined from elastic-backscatter signals alone, because the lidar ratio varies with aerosol types. A combined Raman and elastic backscatter lidar system can provide reliable information about the aerosol extinction profile and the aerosol lidar ratio.

• Journal of the Korean earth science society
• /
• v.32 no.6
• /
• pp.586-594
• /
• 2011

Gravity Recovery and Climate Experiment (GRACE), launched in April, 2002, makes it possible to monitor Earth's mass redistribution with its time-varying gravity observation. GRACE provides monthly gravity solutions as coefficients of spherical harmonics, and thus ones need to convert the gravity spectrum to gravity grids (or mass grids) via the spherical harmonics. GRACE gravity solutions, however, include spatial alias error as well as noise, which requires to suppress in order to enhance signal to noise ratio. In this study, we present the GRACE data processing procedures and introduce some applications of time-varying gravity, which are studies of terrestrial water storage changes, Antarctic and Greenland ice melting, and sea level rise. Satellite missions such as GRACE will continue up to early 2020, and they are expected to be an essential resource to understand the global climate changes.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2011.10a
• /
• pp.177-178
• /
• 2011

• Journal of the Korean earth science society
• /
• v.38 no.3
• /
• pp.173-181
• /
• 2017

The stability of the steady Rossby-Haurwitz wave (R-H wave) in the nondivergent barotropic model (NBM) on the sphere was investigated with the normal mode method. The linearized NBM equation with respect to the R-H wave was formulated into the eigenvalue-eigenvector problem consisting of the huge sparse matrix by expanding the variables with the spherical harmonic functions. It was shown that the definite threshold R-H wave amplitude for instability could be obtained by the normal mode method. It was revealed that some unstable modes were stationary, which tend to amplify without the time change of the spatial structure. The maximum growth rate of the most unstable mode turned out to be in almost linear proportion to the R-H wave amplitude. As a whole, the growth rate of the unstable mode was found to increase with the zonal- and total-wavenumber. The most unstable mode turned out to consist of more-than-one zonal wavenumber, and in some cases, the mode exhibited a discontinuity over the local domain of weak or vanishing flow. The normal mode method developed here could be readily extended to the basic state comprised of multiple zonalwavenumber components as far as the same total wavenumber is given.

• Journal of the Korean earth science society
• /
• v.41 no.4
• /
• pp.344-355
• /
• 2020

Tropical cyclone scale vortices and associated Rossby waves were investigated numerically using high-resolution barotropic models on the global domain. The equations of the barotropic model were discretized using the spectral transform method with the spherical harmonics function as orthogonal basis. The initial condition of the vortex was specified as an axisymmetric flow in the gradient wind balance, and four types of basic zonal states were employed. Vortex tracks showed similar patterns as those on the beta-plane but exhibited more eastward displacement as they moved northward. The zonal-mean flow appeared to control not only the west-east translation but also the meridional translation of the vortex. Such a meridional influence was revealed to be associated with the beta gyre and the Rossby wave, which are formed around the vortex due to the beta effect. In the case of the basic zonal state of climatological mean, the meridional translation speed reached the maximum value when the vortex underwent recurving.