• Geophysics and Geophysical Exploration
• /
• v.16 no.3
• /
• pp.131-138
• /
• 2013

We developed a method for inverting magnetic data to image 2D susceptibility models. The major difficulty in the inversion of the potential data is the nonuniqueness. Furthermore, generally the number of inversion blocks are greater than the number of the magnetic data available, and thus the magnetic inversion leads to under-determined problem, which aggravates the nonuniqueness. When the magnetic data were inverted by the general least-squares method, the anomalous susceptibility would be concentrated near the surface in the inverted section. To overcome this nonuniqueness problem, we propose a new resolution model constraint that is calculated from the parameter resolution. The model constraint imposes large penalty on the model parameter with good resolution, on the other hand small penalty on the model parameter with poor resolution. Thus, the deep-seated model parameter, generally having poor resolution, can be effectively resolved. The developed inversion algorithm is applied to the inversion of the synthetic data for typical models of magnetic anomalies and is tested on real airborne data obtained at the Okcheon belt of Korea.

• Geophysics and Geophysical Exploration
• /
• v.10 no.1
• /
• pp.60-68
• /
• 2007

Conventional analysis of spectral induced polarization (SIP) data consists of measuring impedances over a range of frequencies, followed by spectral analysis to estimate spectral parameters. For the quantitative and accurate estimation of subsurface SIP parameter distribution, however, a sophisticated and stable inversion technique is required. In this study, we have developed a two-step inversion approach to obtain the two-dimensional distribution of SIP parameters. In the first inversion step, all the SIP data measured over a range of frequencies are simultaneously inverted, adopting cross regularisation of model complex resistivities at each frequency. The cross regularisation makes it possible to enhance the noise characteristics of the inversion by imposing a strong assumption, that complex resistivities should show similar characteristics over a range of frequencies. In numerical experiments, we could verify that our inversion approach successfully reduced inversion artefacts. As a second step, we have also developed an inversion algorithm to obtain SIP parameters based on the Cole-Cole model, in which frequency-dependent complex resistivities from the first step are inverted to obtain a two-dimensional distribution of SIP parameters. In numerical tests, the SIP parameter images showed a fairly good match with the exact model, which suggests that SIP imaging can provide a very useful subsurface image to complement resistivity.

• Geophysics and Geophysical Exploration
• /
• v.9 no.4
• /
• pp.299-303
• /
• 2006

Geologic noise, especially located at shallow depth, can be a great obstacle in the interpretation of geophysical data. Thus, it is important to suppress geologic noise in order to accurately detect major anomalous bodies in the survey area. In the inversion of geophysical data, model parameters at shallow depth, which have small size and low contrast of physical property, can be regarded as one of geologic noise. The least-squares method with smoothness constraint has been widely used in the inversion of geophysical data. The method imposes a big penalty on the large model parameter, while a small penalty on the small model parameter. Therefore, it is not easy to suppress small anomalous boies. In this study, we developed a new inversion scheme which can effectively suppress geologic noise by imposing a big penalty on the slowly varying model parameter and a small penalty on the largely varying model parameter. We call the method MTE (main-target emphasizing) inversion. Applying the method to the inversion of 2.5D small loop EM data, we can ensure that it is effective in suppressing small anomalous boies and emphasizing major anomalous bodies in the survey area.

• Atmosphere
• /
• v.25 no.1
• /
• pp.1-18
• /
• 2015

The most objective way to overcome the limitation of numerical weather prediction model is to represent the uncertainty of prediction by introducing probabilistic forecast. The uncertainty of the numerical weather prediction system developed due to the parameterization of unresolved scale motions and the energy losses from the sub-scale physical processes. In this study, we focused on the growth of model errors. We performed ensemble forecast to represent model uncertainty. By employing the multi-physics scheme (PHYS) and the stochastic kinetic energy backscatter scheme (SKEBS) in simulating typhoon Rusa (2002), we assessed the performance level of the two schemes. The both schemes produced better results than the control run did in the ensemble mean forecast of the track. The results using PHYS improved by 28% and those based on SKEBS did by 7%. Both of the ensemble mean errors of the both schemes increased rapidly at the forecast time 84 hrs. The both ensemble spreads increased gradually during integration. The results based on SKEBS represented model errors very well during the forecast time of 96 hrs. After the period, it produced an under-dispersive pattern. The simulation based on PHYS overestimated the ensemble mean error during integration and represented the real situation well at the forecast time of 120 hrs. The displacement speed of the typhoon based on PHYS was closest to the best track, especially after landfall. In the sensitivity tests of the model uncertainty of SKEBS, ensemble mean forecast was sensitive to the physics parameterization. By adjusting the forcing parameter of SKEBS, the default experiment improved in the ensemble spread, ensemble mean errors, and moving speed.

• Journal of the Korean earth science society
• /
• v.30 no.7
• /
• pp.834-844
• /
• 2009

Two-dimensional magnetotelluric (MT) inversion algorithm using minimum support (MS) stabilizer functional was implemented in this study to enhance the contrast of inverted images. For this implementation, this study derived a formula in discrete form for creeping model updates in the least-squares linearized inversion. A spatially varying regularization parameter determination algorithm, which is known as ACB (Active Constraint Balancing), was also adopted to stabilize the inversion process when using MS stabilizer as a model constraint. Inversion experiments for a simple isolated body model show well the feature of MS stabilizer in concentrating the anomalous body compared with the second-order derivative model constraint. This study also compared MS stabilizer and the second-order derivative model constraints for a model having multiple anomalous bodies to show the applicability of the algorithm into field data.

• Journal of the Korean earth science society
• /
• v.25 no.6
• /
• pp.495-501
• /
• 2004

Using numerous ground-based or space-based telescopes by many astronomers, AGN spectroscopic monitoring campaigns have been carried out over many years to study the variability of continua and emission lines. We investigate the SWP IUE spectra of a large aperture configuration for the Seyfert 1 galaxies NGC 4151 and NGC 5548. We estimate the BLR electron number densities and their variation from the line ratios of C III] 1909 to Si III] 1892. With the ratios of C IV 1550 to C III] 1909 which give us the information on the ionization parameter of BLR, we try to find the physical conditions of the BLR and activities of he super massive black hole surroundings. The BLR density variations scale as 4 and 8 for NGC 4151 and NGC 5548, respectively. Based on the BLR size and C III] line profiles, we found both black hole masses as about $10^7$ $M_{\odot}$.

• Journal of the Korean earth science society
• /
• v.31 no.1
• /
• pp.88-94
• /
• 2010

A simple laboratory device and methodology are developed for deep understanding of the western boundary current (WBC). The apparatus consists of a rotating (count-clockwise) rectangular tank with a sloping bottom in order to simulate the beta effect, the variation of the Coriolis parameter with latitude. We also used a clockwise rotating disk at the surface water to mimic the wind stress forcing in mid-latitude oceans. Four experiments were carried out using some combination of a bottom type and a rotation of the tank. Experiment with the beta effect clearly demonstrated the WBClike flow as well as the Sverdrup interior. However, the water in a tank gave rise to an inertial motion under the influence of a constant Coriolis force alone. We also discussed a stiffening of the rotating fluid on the basis of the Taylor-Proudman effect. We believe that the apparatus and methods developed in this study help to understand the WBC due to the beta effect.

• Korean Journal of Remote Sensing
• /
• v.29 no.1
• /
• pp.137-149
• /
• 2013

Solar insolation is essential to understand the interaction between the earth and solar system, and it is a significant parameter that is utilized in various research fields including earth science, agriculture, and energy engineering. Although solar insolation is broadly measured in the ground-based observation station, it is difficult to identify the spatial distribution of solar insolation accurately. The remote sensing approach is known to have several benefits because it can provide continuous data sets for large area. In this study, we conducted the validation of solar insolation from COMS in the South Korea by comparing with flux tower observation. The results showed that the correlations between COMS and observation were high in both 30 minutes interval data and daily average data. Thus, we can identify that COMS can provide a reasonable estimate of solar insolation.

• Journal of the Korean earth science society
• /
• v.31 no.2
• /
• pp.139-150
• /
• 2010

This study investigated the depositional environment of the area of Anmyeondo in the mid-western coast of Korea. A total of 142 surface sediments were analyzed, and leveling survey was conducted as well in 12 transects at the sandy tidal flat in the area. The surface sediments are classified into several different types (slightly gravelly Sand, Sand, slightly gravelly sandy Mud, etc.) and they show fining-distribution toward the offshore. Textural parameters are characterized by being coarse, well sorted, positively skewed in the supratidal and intertidal flat, whereas being finer, poorly sorted, negatively skewed in the subtidal flat. According to the C-M diagram, the sediment transport modes are characterized by the mixed mode of bedload and graded suspension in supratidal and intertidal flat, and by graded suspension or uniform suspension in the subtidal flat. Histogram shows a bimodal distribution in the lower subtidal flat which is indicative of two sediment transport modes.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.86-92
• /
• 2006

The full-waveform inversion algorithm using normalised seismic wavefields can avoid potential inversion errors due to source estimation required in conventional full-waveform inversion methods. In this paper, we have modified the inversion scheme to install a weighted smoothness constraint for better resolution, and to implement a staged approach using normalised wavefields in order of increasing frequency instead of inverting all frequency components simultaneously. The newly developed scheme is verified by using a simple two-dimensional fault model. One of the most significant improvements is based on introducing weights in model parameters, which can be derived from integrated sensitivities. The model-parameter weighting matrix is effective in selectively relaxing the smoothness constraint and in reducing artefacts in the reconstructed image. Simultaneous multiple-frequency inversion can almost be replicated by multiple single-frequency inversions. In particular, consecutively ordered single-frequency inversion, in which lower frequencies are used first, is useful for computation efficiency.