We have investigated the numerical methods to calculate model atmosphere for the analysis of spectral lines emitted from the sun and stars. Basic equations used in our calculations are radiative transfer, statistical equilibrium and charge-particle conservations. Transfer equation has been solved to get emitting spectral line profile as an initial value problem using Adams-Bashforth-Moulton method with accuracy as high as 12th order. And we have calculated above non linear differential equations simultaneously as a boundary value problem by finite difference method of 3 points approximation through Feautrier elimination scheme. It is found that all computing programs coded by above numerical methods work successfully for our model atmosphere.

We have investigated the solar activity variation with period shorter than 1000 days, through Fourier transformation of solar cycle 21 and 22 data. And real time predictions of the flare maximum intensity have been made by multilinear regression method to allow the use of multivariate vectors of sunspot groups or active region characteristics. In addition, we have examined the evolution of magnetic field and current density in active regions at times before and after flare occurrence, to check short term variability of solar activity. According to our results of calculation, solar activity changes with periods of 27.1, 28.0, 52.1, 156.3, 333.3 days for solar cycle 21 and of 26.5, 27.1, 28.9, 54.1, 154, 176.7, 384.6 days for solar cycle 22. Periodic components of about 27, 28, 53, 155 days are found simultaneously at all of two solar cycles. Finally, from our intensive analysis of solar activity data for three different terms of 1977~1982,1975~1998,and1978~1982, we find out that our predictions coincide with observations at hit rate of 76%,63% 76%,63% , 59 respectively.

We have made intensive calculations on the maximum relative sunspot number and the date of solar maximum of 23rd solar cycle, by using the statistical and precursor methods to predict solar activity cycle. According to our results of solar data processing by statistical method, solar maximum comes at between February and July of 2000 year and at that time, the smoothed sunspot number will reach to 114.3~122.8. while precursor method gives rather dispersed value of 118~171 maximum sunspot number. It is found that prediction by statistical method using smoothed relative sunspot number is more accurate than by any method to use any data of 10.7cm radio fluxes and geomagnetic aa, Ap indexes, from the full analysis of solar cycle pattern of these data. In fact, current ascending pattern of 23rd solar cycle supports positively our predicted values. Predicted results by precursor method for Apavg,aa31−36 Apavg,aa31−36 indexes show similar values to those by statistical method. Therefore, these indexes can be used as new precursors for the prediction of 23rd or next solar cycle.

We have carried out high precision orbit calculation, by using various numerical techniques with accuracy of higher than fourth order, in order for exact prediction on position and velocity of celestial bodies and artificial satellites. General second order ordinary differential equation has been solved numerically to test the performance for each of numerical methods. We have compared computed values with exact solution obtained by using universal variables for two body problem and discussed overall results of numerical methods used in our calculation. As a result, it is found that high order difference table method called as Gauss-Jackson method is best one with easiness and efficiency in the increase of accuracy by number of initial values.

We have developed solar observational system in the department of Astronomy & Space Sciences of KyungHee University, in order to monitor solar activities and construct solar database for space weather forecasting at maximum of 23rd solar cycle, as well as an solar education and exercise for undergraduate students. Our solar observational system consists of the full disk monitoring system and the regional observation system for H a fine structure. Full disk monitoring system is made of an energy rejection filter, 16cm refractor, video CCD camera and monitor. Monitored data are recorded to VHS video tape and analog output of video CCD can be captured as digital images by the computer with video graphic card. Another system for regional observation of the sun is made of energy rejection filter, 21cm Schmidt-Cassegrain reflector, H a filter with 1.6A pass band width and 375×242 375×242 CCD camera. We can observe H a fine structure in active regions of solar disk and solar limb, by using this system. We have carried out intense solar observations for a test of our system. It is found that Quality of our H a image is as good as that of solar images provided by Space Environmental Center. In this paper, we introduce the basic characteristics of the KyungHee Solar Observation System and result of our solar observations. We hope that our data should be used for space weather forecasting with domestic data of RRL(Radio Research Laboratory) and SOFT(SOlar Flare Telescope).

Theoretical calculations of the combined radiative transfer and statistical equilibrium equation including the charge-particle conservations have been earned out for a multilevel hydrogen atom in quiescent prominences. Cool and dense models show the steep changes of population and radiation field in the vicinity of the surface, while these physical quantities remain unchanged for models with temperature of 7,300K, regardless of total densities. Ionization rate of hydrogen atom related with metallic line formation varies in considerable amounts from the surface to the center of model prominences cooler than 6,300K. However, such cool models cannot release enough hydrogen line emissions to explain observed intensities. Prominence models with a temperature higher than 8,000K can yield the centrally reversed Lyman line profiles confirmed by satellite EUV observations. We find that queiscent prominence with a density between 2 × 10 11 and 10 12 c m − 3 should be in temperature range between 6,300K and 8,300K, in order to explain consistently observed H alpha, beta line emissions and n p / n l ratio.

Fine structures of a quiescent prominence are studied by analyzing high resolution H alpha filtergrams and H alpha line spectra observed at the Hida Observatory of Kyoto University. We have found two kinds of downward motions in the prominence. One of them is a movement with a constant acceleration below the solar gravity( ≃ 1 / 4 g s ) and the other with an uniform velocity( ≃ 16 K m / s ). The average life time and the size of prominence knots are estimated to be about 7 minutes and 4000Km, respectively. Spatial and brightness distribution of knots are also presented in this paper. With the analytical solutions derived from magnetostatic equilibrium in the prominence, we have examined the filamentary structure based on the Kippenhahn-Schluter model. Sag angles of the magnetic fields supporting the prominence matter are predicted from the observed density profile.