• Journal of the Korean earth science society
• /
• v.36 no.5
• /
• pp.476-483
• /
• 2015

The high-order Laplacian-type filter, which is capable of providing isotropic and sharp cut-off filtering on the spherical domain, is essential in processing geophysical data. In this study, a spherical high-order filter was designed by combining the Fourier method with finite difference-method in the longitude and latitude, respectively. The regular grid system was employed in the filter, which has uniform angular spacing including the poles. The singularity at poles was eliminated by incorporating variable transforms and continuity-matching boundary conditions across poles. The high-order filter was assessed using the Rossby-Haurwitz wave, the observed geopotential, and observed wind field. The performance of the filter was found comparable to the filter based on the Galerkin procedure. The filter, employing the finite difference method, can be designed to give any target order of accuracy, which is an important advantage being unavailable in other methods. The computational complexity is represented with 2n-1 diagonal matrices solver with n being the target order of accuracy. Along with the availability of arbitrary target-order, it is also advantageous that the filter can adopt the reduced grid to increase computational efficiency.

• Proceedings of the KSRS Conference
• /
• 2003.11a
• /
• pp.1236-1238
• /
• 2003

The Earth's total magnetic field was extracted from on board TAM (Three Axis Magnetometer) observations of KOMPSAT-1 satellite between June 19th and 21st, 2000. In the pre-processing, the TAM's telemetry data were transformed from ECI (Earth Centered Inertial frame) to ECEF (Earth Centered Earth Fixed frame) and then to spherical coordination, and self-induced magnetic field by satellite bus itself were removed by using an on-orbit magnetometer data correction method. The 2-D wavenumber correlation filtering and quadrant-swapping method were applied to the pre-processed data in order to eliminate dynamic components and track-line noise, respectively. Then, the spherical harmonic coefficients are calculated from KOMPSAT-1 data. To test the validity of the TAM's geomagnetic field, Danish/NASA/French ${\phi}$rsted satellite's magnetic model and IGRF2000 model were used for statistical comparison. The correlation coefficient between ${\phi}$rsted and TAM is 0.97 and IGRF and TAM is 0.96. It was found that the data from on board magnetometer observations for attitude control of Earth-observing satellites can be used to determinate the Earth's total magnetic field and that they can be efficiently used to upgrade the global geomagnetic field coefficients, such as IGRF by providing new information at various altitudes with better temporal and spatial coverage.

• The Bulletin of The Korean Astronomical Society
• /
• v.17
• /
• pp.18.1-18.1
• /
• 1992

• Proceedings of the KSRS Conference
• /
• 2002.10a
• /
• pp.242-242
• /
• 2002

The Earth's magnetic field acquired from KOMPSAT-1's TAM (Three-Axis Magnetometer) between June 19th and 21st 2000 was analyzed. The TAM, one of the KOMPSAT-1's Attitude and Orbit Control Subsystems, plays an important role in determining and controlling the satellite's attitude. This also can provide new insight on the Earth's magnetic field. By transforming the satellite coordinate from ECI to ECEF, spherical coordinate of total magnetic field was achieved. These data were grouped into dusk (ascending) and dawn (descending) data sets, based on their local magnetic times. This partitioning is essential for performing 1-D WCA (Wavenumber Correlation Analysis). Also, this enhances the perception of external fields in the Kompsat-1's TAM magnetic maps that were compiled according to different local. The dusk and dawn data are processed independently and then merged to produce a total field magnetic anomaly map. To extract static and dynamic components, the 1-D and 2-D WCAs were applied to the sub-parallel neighboring tracks and dawn-dusk data sets. The static components were compared with the IGRF, the global spherical harmonic magnetic field model. The static and dynamic components were analyzed in terms of corefield, external, and crustal signals based on their origins.

• Journal of the Korean earth science society
• /
• v.38 no.2
• /
• pp.105-114
• /
• 2017

Two dimensional finite element method with quadrilateral basis functions was applied to the spherical high order filter on the spherical surface limited area domain. The basis function consists of four shape functions which are defined on separate four grid boxes sharing the same gridpoint. With the basis functions, the first order derivative was expressed as an algebraic equation associated with nine point stencil. As the theory depicts, the convergence rate of the error for the spherical Laplacian operator was found to be fourth order, while it was the second order for the spherical Laplacian operator. The accuracy of the new high order filter was shown to be almost the same as those of Fourier finite element high order filter. The two-dimension finite element high order filter was incorporated in the weather research and forecasting (WRF) model as a hyper viscosity. The effect of the high order filter was compared with the built-in viscosity scheme of the WRF model. It was revealed that the high order filter performed better than the built in viscosity scheme did in providing a sharper cutoff of small scale disturbances without affecting the large scale field. Simulation of the tropical cyclone track and intensity with the high order filter showed a forecast performance comparable to the built in viscosity scheme. However, the predicted amount and spatial distribution of the rainfall for the simulation with the high order filter was closer to the observed values than the case of built in viscosity scheme.

• Journal of Astronomy and Space Sciences
• /
• v.16 no.2
• /
• pp.159-166
• /
• 1999

As the spatial resolution of remote sensing satellites becomes higher, very accurate determination of the position of a LEO (Low Earth Orbit) satellite is demanding more than ever. Non-symmetric Earth gravity is the major perturbation force to LEO satellites. Since the orbit propagation is performed in the celestial frame while Earth gravity is defined in the terrestrial frame, it is required to convert the coordinates of the satellite from one to the other accurately. Unless the coordinate conversion between the two frames is performed accurately the orbit propagation calculates incorrect Earth gravitational force at a specific time instant, and hence, causes errors in orbit prediction. The coordinate conversion between the two frames involves precession, nutation, Earth rotation and polar motion. Among these factors, unpredictability and uncertainty of Earth rotation, called UTI-UTC, is the largest error source. In this paper, the effect of UTI-UTC on the accuracy of the LEO propagation is introduced, tested and analzed. Considering the maximum unpredictability of UTI-UTC, 0.9 seconds, the meaningful order of non-spherical Earth harmonic functions is derived.

• Korean Journal of Remote Sensing
• /
• v.18 no.3
• /
• pp.171-179
• /
• 2002

Sampling rates become inconsistent when spatial data in the spherical coordinate are resampled with respect to latitudinal or longitudinal degree for mathematical processes such as Fourier Transform, and this results in distortions of the processed data in the wavenumber domain. These 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 ERS-1 satellite radar altimeter data from the Barents Sea in the Russian Arctic, and a method is presented to minimize the distortion using the Lambert Conformal Conic map projection. This approach was found to enhance the free-air gravity anomalies in both data and wavenumber domains.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.26 no.10
• /
• pp.924-931
• /
• 2015

When a target locates at low altitude over a curved earth surface and far away from a radar, we examine the accuracy of the conventional formulations to compute the radar echo signal. The 4-ray model is used to calculate the scattering by the target to consider the ground plane effect. In this paper, the diffracted wave is not included. Based on the required parameters able to be calculated by known equations for the estimation of the wave propagation phenomena, the radar echo signal is computed and verified by comparing with measured data sets.

• International Journal of Aeronautical and Space Sciences
• /
• v.13 no.4
• /
• pp.484-490
• /
• 2012

This paper is devoted to investigate the feasibility of using a medium power ground-based laser to produce a torque on LEO satellites of various shapes. The laser intensity delivered to a satellite is calculated using a simple model of laser propagation in which a standard atmospheric condition and linear atmospheric interaction mechanism is assumed. The laser force is formulated using a geocentric equatorial system in which the Earth is an oblate spheroid. The torque is formulated for a cylindrical satellite, spherical satellites and for satellites of complex shape. The torque algorithm is implemented for some sun synchronous low Earth orbit cubesats. Based on satellites perigee height, the results demonstrate that laser torque affecting on a cubesat has a maximum value in the order of $10^{-9}$ which is comparable with that of solar radiation. However, it has a minimum value in the order of $10^{-10}$ which is comparable with that of gravity gradient. Moreover, the results clarify the dependency of the laser torque on the orbital eccentricity. As the orbit becomes more circular it will experience less torque. So, we can conclude that the ground based laser torque has a significant contribution on the low Earth orbit cubesats. It can be adjusted to obtain the required control torque and it can be used as an active attitude control system for cubesats.

• Economic and Environmental Geology
• /
• v.49 no.4
• /
• pp.269-280
• /
• 2016

I presented three vector crustal magnetic anomaly components and six gradients by using spherical Slepian functions over the cap area of $20^{\circ}$ of radius centered on the South Pole. The Swarm mission, launched by European Space Agency(ESA) in November of 2013, was planned to put three satellites into the low-Earth orbits, two in parallel in East-West direction and one in cross-over of the higher altitude. This orbit configuration will make the gradient measurements possible in North-South direction, vertical direction, as well as E-W direction. The gravity satellites, such as GRACE and GOCE, have already implemented their gradient measurements for recovering the accurate gravity of the Earth and its temporal variation due to mass changes on the subsurface. However, the magnetic gradients have little been applied since Swarm launched. A localized magnetic modeling method is useful in taking an account for a region where data availability was limited or of interest was special. In particular, computation to get the localized solutions is much more efficient and it has an advantage of presenting high frequency anomaly features with numbers of solutions fewer than the global ones. Besides, these localized basis functions that were done by a linear transformation of the spherical harmonic functions, are orthogonal so that they can be used for power spectrum analysis by transforming the global spherical harmonic coefficients. I anticipate in scientific and technical progress in the localized modeling with the gradient measurements from Swarm and here will do discussion on the results of the localized solution to represent the three vector and six gradient anomalies over the Antarctic area from the synthetic data derived from a global solution of the spherical harmonics for the crustal magnetic anomalies of Swarm measurements.