In this study, the environmental behavior of malodor pollutants (MPs:　H2S, CH3SH, DMS, and DMDS) was investigated around areas influenced by strong anthropogenic processes based on observations and modeling study (a CALPUFF dispersion model). The MP emission concentrations were measured from 8 industrial source regions (tire plants (S1-S3), waste water disposal plant (S4), and oil refinery (S5) in an urban center area and paper mill/incineration plant (S6) and livestock feedlots (S7-S8) in Ungsang area) in Yangsan city during a fall period in 2008 (21 October 2008). Overall, the most MPs emitted from the urban center area were found to affect the malodor pollution in their downwind areas during early morning (06:00 LST) and nighttime (18:00 and 21:00 LST), compared with those in the Ungsang area. For malodor intensity, the most MPs in the urban center area (especially S1 and S2) were found to be a significant contributor, whereas CH3SH and H2S in the Ungsnag area (especially S6) were the dominant contributor. The model study showed agreement in the spatial distributions of simulated MPs with those of the observations. The largest impact of MPs in the urban center area on the malodor pollution in its residential areas occurred at S1, S2, and S3 sites during nighttime, while that of MPs in the Ungsang area occurred at S6 and S8 sites. This may be caused mainly by the high MP emissions and in part by wind conditions (prevailing northeasterly winds with low wind speeds of 2-3 m/s).

Development of an artificial neural network model was presented to predict the daily maximum SO2 concentration in the urban-industrial area of Ulsan. The network model was trained during April through September for 2000-2005 using SO2 potential parameters estimated from meteorological and air quality data which are closely related to daily maximum SO2 concentrations. Meteorological data were obtained from regional modeling results, upper air soundings and surface field measurements and were then used to create the SO2 potential parameters such as synoptic conditions, mixing heights, atmospheric stabilities, and surface conditions. In particular, two-stage clustering techniques were used to identify potential index representing major synoptic conditions associated with high SO2 concentration. Two neural network models were developed and tested in different conditions for prediction: the first model was set up to predict daily maximum SO2 at 5 PM on the previous day, and the second was 10 AM for a given forecast day using an additional potential factors related with urban emissions in the early morning. The results showed that the developed models can predict the daily maximum SO2 concentrations with good simulation accuracy of 87% and 96% for the first and second model. respectively, but the limitation of predictive capability was found at a higher or lower concentrations. The increased accuracy for the second model demonstrates that improvements can be made by utilizing more recent air quality data for initialization of the model.

We focused on effects on data assimilation of simulated wind fields by using upper-air observations (wind profiler and sonde data). Local Analysis Prediction System (LAPS), a type of data assimilation system, was used for wind field modeling. Five cases of simulation experiments for sensitivity analysis were performed : which are EXP0) non data assimilation, EXP1) surface data, EXP2) surface data and sonde data, EXP3) surface data and wind profiler data, EXP4) surface data, sonde data and wind profiler data. These were compared with observation data. The result showed that the effects of data assimilation with wind profiler data were found to be greater than sonde data. The delicate wind fields in complex coastal area were simulated well in EXP3. EXP3 and EXP4 using wind profiler data with vertically high resolution represented well sophisticated differences of wind speed compared with EXP1 and EXP2, this is because the effects of wind profiler data assimilation were sensitively adjusted to first guess field than those of sonde observations.

We employed two data assimilation techniques including MM5 Four Dimensional Data Asssimilation (FDDA) and Local Analysis and Prediction System (LAPS) to find out the effects of the changed initial conditions on the wind fields simulation according to the objective analysis methods. We designed 5 different modeling cases. EXP B used no data assimilation system. Both EXP F1 using surface observations and EXP F2 with surface and upper-air observations employed MM5 FDDA. EXP L1 using surface observations and EXP L2 with surface and upper-air observations used LAPS. As results of, simulated wind fields using MM5 FDDA showed locally characterized wind features due to objective analysis techniques in FDDA which is forcefully interpolating simulated results into observations. EXP F1 represented a large difference in comparison of wind speed with EXP B. In case of LAPS, simulated horizontal distribution of wind fields showed a good agreement with the patterns of initial condition and EXP L1 showed comparably lesser effects of data assimilation of surface observations than EXP F1. When upper-air observations are applied to the simulations, while MM5 FDDA could hardly have important effects on the wind fields simulation and showed little differences with simulations with merely surface observations (EXP F1), LAPS played a key role in simulating wind fields accurately and it could contribute to alleviate the overestimated winds in EXP L1 simulations.

We focus on the improvement of accuracy of sea surface wind over complex coastal area during the warm season. Local Analysis Prediction System (LAPS) was used to improve the initial values in Mesoscale Meteorological model (MM5). During the clear summer days with weak wind speed, sea surface wind simulated with LAPS was compared with the case without LAPS. The results of modeling with LAPS has a good agreement mesoscale circulation such as mountain and valley winds on land and in case of modeling without LAPS, wind speed overestimated over the sea in the daytime. And the results of simulation with LAPS indicated similar wind speed values to observational data over the sea under influence of data assimilation using BUOY, QuikSCAT, and AMEDAS. The present study suggests that MM5 modelling with LAPS showed more improved results than that of without LAPS to simulate sea surface wind over the complex coastal area.

We focused on improvement in simulation of wind fields for the complex coastal area. Local Analysis and Prediction System(LAPS) was used as a data assimilation method to improve initial conditions. Case studies of different LAPS inputs were performed to compare improvement of wind fields. Five cases have been employed : Ⅰ) non data assimilation, Ⅱ) all available data, Ⅲ) AWS, buoy, QuikSCAT, Ⅳ) AWS, buoy, wind profiler, Ⅴ) AWS, buoy, AMEDAS. Data assimilation can supplement insufficiency of the mesoscale model which does not represent detailed terrain effect and small scale atmospheric flow fields. Result assimilated all available data showed a good agreement to the observations rather than other cases and estimated well the local meteorological characteristics including sea breeze and up-slope winds. Result using wind profiler data was the next best thing. This implies that data assimilation with many high-resolution sounding data could contribute to the improvements of good initial condition in the complex coastal area. As a result, these indicated that effective data assimilation process and application of the selective LAPS inputs played an important role in simulating wind fields accurately in a complex area.

In predicting oxidants concentration, the most important fact is to select a suitable photochemical reaction mechanism. Sensitivity analysis of O3 and other important photochemical oxidants concentrations was conducted by using CBM-IV model. The predicted oxidants concentration was considerably related with the initial concentration of formaldehyde, [NO2]/[NO],NOx, RH and RCHO. As the initial concentration of formaldehyde increased, concentration of NO2 increased. O3 concentration was proportional to the [NO2]/[NO] ratio. When the initial concentrations of RH and RCHO were high, photochemical reaction was more reactive, including more rapid conversion of NO to NO2 and increased oxidants. Also, the sensitivities of ozone formation to rate constants, Kl,K2andK3 in the NO2 photolysis were studied.

The southeastern coastal area of the Korean peninsula has a complex terrain including an irregular coastline and moderately high mountains. This implies that mesoscale circulations such as mountain-valley breeze and land-sea breeze can play an important role in wind field and ocean forcing. In this study, to improve the accuracy of complex coastal wind field(surface wind and sea surface wind), we carried out the sensitivity experiments based on PBL schemes in PSU/NCAR Mesoscale Model (MM5), which is being used in the operational system at Korea Meteorological Administration. Four widely used PBL parameterization schemes in sensitivity experiments were chosen: Medium-Range Forecast (MRF), High-resolution Blackadar, Eta, and Gayno-Seaman scheme. Thereafter, case(2004. 8. 26 - 8. 27) of weak-gradient flows was simulated, and the time series and the vertical profiles of the simulated wind speed and wind direction were compared with those of hourly surface observations (AWS, BUOY) and QuikSCAT data. In the simulated results, the strength of wind speed of all schemes was overestimated in complex coastal regions, while that of about four different schemes was underestimated in islands and over the sea. Sea surface wind using the Eta scheme showed the highest wind speed over the sea and its distribution was similar to the observational data. Horizontal distribution of the simulated wind direction was very similar to that of real observational data in case of all schemes. Simulated and observed vertical distribution of wind field was also similar under boundary layer(about 1 km), however the simulated wind speed was underestimated in upper layer.

In this study, we focused on the improvements in the simulation of sea surface wind over the complex coastal area. MM5 model being currently used to predict sea surface wind at Korea Meteorological Administration, was used to verify the accuracy to estimate the local wind field. A case study was performed on clear days with weak wind speed(4 m/s), chosen by the analysis of observations. The model simulations were conducted in the southeastern area of Korea during the selected periods, and observational data such as AWS, buoy and QuikSCAT were used to compare with the calculated wind components to investigate if simulated wind field could follow the tendency of the real atmospheric wind field. Results showed that current operational model, MM5, does not estimate accurately sea surface wind and the wind over the coastal area. The calculated wind speed was overestimated along the complex coastal regions but it was underestimated in islands and over the sea. The calculated diurnal changes of wind direction could not follow well the tendency of the observed wind, especially at nighttime. In order to exceed the limitations, data assimilation with high resolution data and more specificated geographical information is expected as a next best policy to estimate accurately the environment of local marine wind field.

This study presents the characteristics of nocturnal inversion layer and their effect on the concentration variations of surface air pollutants using tethersonde and automatic weather station (AWS, 2 layer tower) system in Ulsan during 2003. The method for the distinction of inversion intensity was decided based on the sum of nocturnal temperature gradient. As the results, there was a close correlation (correlation coefficient of 0.76) between the maximum inversion height obtained from tethersonde and the sum of nocturnal temperature gradient. The air pollutant concentration was also directly proportional to the inversion intensity. When the inversion intensity was strong in the nighttime, ozone (O3) concentration was lower, while nitrogen dioxide (NO2) concentration was higher. The carbon monoxide (CO) concentration was gradually higher according to the nocturnal inversion intensity, whereas sulfur dioxide (SO2) concentration was relatively constant. In addition, we found that there was no correlation between the inversion intensity and TSP concentration.

It has been researched the relationship between deposition velocity and factors which could affect the deposition phenomena and deposition velocity also has been estimated for several land-use types. The typical deposition velocities are complex functions of surface types, atmospheric stabilities, friction velocities, air pollutants and so on. The canopy resistance is major contribution to the model's total resistance for O3. Canopy wetness is also an important factor to calculate deposition velocity. We considered the canopy wetness as canopy water content(CWC) in our Model. But, it is not easy to observe CWC over each land-use types. In this study, we use CWC observed by EMEFS(CANADA Environment Service, 1988) to examine the influence of CWC in estimation of O3 dry deposition velocity(Vd) in summertime. The value of O3 Vd range 0.2～0.7 cms-1 on dry surface and 0.01～0.35 cms-1 on wet surface in daytime.

We calculated dose rate using radiative transfer equations to consider radiative processes distinctly. The dose rate at Pohang(36°02'N, 129°23'E) was calculated using measured ozone and meteorological data and two-stream approximations(quadrature, Eddington, delta Eddington, PIFM(practical improved flux method), discrete ordinate, delta discrete ordinate) are used in solving equation. The purpose of this study is to determine the most compatible radiative transfer approximation for simulating the radiative and photochemical processes of atmosphere through comparision between calculated and measured values. Dose rate of the biologically effective irradiance in the region 0.28-0.32μm showed the highest value when quadrature and Eddington was used and lower value on condition that delta scaling was applied. Correlation coefficient between dose rate at surface using radiation transfer equation and measured UV-B at Pohang was 0.78, 0.79 and 0.81 when detla Eddington, PIFM and delta discrete ordinate were used. Also, in case of above approximations were used, MBE(Mean Bias Error) was within -0.3MED/30min and RMBE(Relative Mean Bias Error) was below 10% between 1200 LST and 1400 LST. Approximations which are compatible in estimating radiative process are detla Eddington, PIFM and delta discrete ordinate. Especially, in case that radiative process is considered more detail, delta discrete ordinate increased the number of stream is proper.

Weather elements were observed by the AWS (Automatic Weather System) and dustfall particles were collected by the modified American dust jar (wide inlet bottle type) at 4 sampling sites in Busan area from March, 1999 to February, 2000. Thirteen chemical species (Al, Ca, Cd, Cr, Cu, Fe, K, Mg, Mn, Ni, Pb, Si, and Zn) were analyzed by AAS and ICP. The purposes of this study were to estimate qualitatively various bulk deposition flux of dustfall and insoluble components by applying regional and seasonal wind intensity. Frequency of wind speed were found in order of low(1-3㎧), very low(<1㎧), medium(3-8㎧) and high(>8㎧), and annual mean had higher range at low(1-3㎧) for 56.3%. Strong negative linear correlation were observed between dustfall and wind direction (northeastern and eastern), but strong positive linear correlation were observed between dustfall and wind direction (western and northwestern) at industrial, commercial and coastal zone(p<0.05). While a negative correlation were observed between wind speed frequency of very low(<1㎧) and dustfall, and positive correlation were observed between wind speed frequency of low(1-3㎧) and dustfall in coastal zone(p<0.05). The correlation coefficient was observed 0.556 between wind speed frequency of low(1-3㎧) and Ni by commercial zone(p<0.05). The correlation coefficient show well-defined insoluble trace metals (Al, Ca, Cr, Cu, Fe, Pb, and Zn) and wind speed frequency of low(1-3㎧) at coastal zone, which was found significant difference(p<0.01).

A case of air pollution study of estimation of emission rate for source to railroad rolling stock and emissive character about is unusual. Recent emission rate of railroad rolling stock was estimated with emission factor of EPA in three region(Seoul, Incheon, and Gyeonggido). But this EPA factor could be incorrect because Korea and America have a different railroad environment in the variety of fuel and character of railroad rolling stock. In this study, emission rate of each line, car and region(district) with emission factor of Korean case(National Railroad Administration, 1997) was estimated. In Busan Metropolitan City, railroad rolling stocks were divided accoding to, each line, car and service. Particularly, the idle fuel rate omitted in the preceding study was included in calculation fuel rate. Total emission rate of Busan Metropolitan City was 887.41t/year. Each emission rates of Kyeongbu line, Donghaenambu line, Gaya line, Uam line, Bujeon line, and idle was 489.15t, 196.46t, 33.94t, 12.66t, 6.47t, and 148.72t, respectively.

Dustfall particles were collected by the modified American dust jar (wide inlet bottle type) at 6 sampling sites in Pusan area from March, 1999 to February, 2000. Thirteen chemical species (Al, Ca, Cd, Cr, Cu, Fe, K, Mg, Mn, Ni, Pb, Si, and Zn) were analyzed by AAS and ICP. The purposes of this study were to estimate qualitatively various bulk deposition flux of dustfall and insoluble components by applying regional and seasonal distribution. Dustfall amount of regional variations were found in order of coastal zone, industrial zone, commercial zone, agricultural zone and residential zone, and seasonal total dustfall had higher concentrations during spring for 6.741 ton/㎢/season, lower concentrations during summer for 1.989 ton/㎢/season, and annual total concentration was 17.742 ton/㎢/year. The regional distributions of enrichment factor show well-defined anthropogenic metals (Cd, Cu, Pb, and Zn) at industrial and agricultural zone, and contribution rate of soil particles were found in order of summer, fall, winter and spring. Factor loading effects of chemical composition of dustfall were found in order of road traffic emission source and combustion processed source, industrial activity source, soil source and marine source.

For the purpose of predicting air pollutants concentration in Pusan coastal urban, we used an Eulerian model of flow and dispersion/chemistry/deposition process considering SST effects which estimate through POM. The results of air quality model including emission from various sources show that the seasonal variation pattern of respective pollutants was affected by the seasonal SST fields and local circulation. Horizontal deviation of diurnal SST was 2.5∼4K, especially large gradients in coastal region. Through numerical simulation of wind fields we predicted that local circulation prevailed during daytime in summer and nighttime in winter. So high concentration distribution showed toward inland in spring and summer seasons, while high concentration distribution showed at inland near coast in autumn and winter.

The atmospheric conditions and the transport mechanism of long-range transport of air pollutants from coastal area to inland area were investigated using regular meteorological data and air pollution data obtatined from the southeastern area of Korea. Daytime temperature over the inland area(Taegu) was higher than that over the coastal area(Pusan) and the temperature difference of about 5∼6℃ when the thermal low was most fully developed and the sea level pressure over Taegu was lower than that over Pusan by about 4∼5 hPa at that time. Therefore this low pressure appeared to the thermally induced low. Air mass polluted from the coastal area during the morning period was transported inland area, at first by the sea breeze and by the large scale wind system toward the thermal low generated in the mountainous inland region. This was explained by the fact that the concentration of air pollutants over Taegu increased throughtout the late afternoon.

In order to investigate horizontal wind field in the boundary layer around Pusan area, wind speed and wind direction measured at 14 AWS(Automatic Weather Station), 1997, was used. The wind direction at PRMO(Pusan Regional Meterological Office) was showed that southwest and northeast wind dominated for spring and summer, northeast wind for fall and northwest for winter. Anticline flow was showed at 'Gaekumm' which is located between Mt. Backyang(641 m) and Mt. Yumkwang(503 m) and affected on wind field at 'Pusanjin'. The low wind speed and various wind direction was represented at the basin topography, 'Buckgu', 'Jeasong', 'Ilkwang' and 'Kijang'. The annual mean wind speed at 14 sites, 2.5 ms-1, was lower than that measured at PRMO, 3.9 ms-1. The wind direction analysis showed that the case of same direction in compare with that measured at PRMO is about 54% and case of opposite direction is about 12%. Annual and seasonal mean windrose showed wind direction is affected by not only synoptic weather state but also topography.

Visibility reduction is a barometer of air pollution, which people can notice easily. First of all, we need to measure quantified visibility continuously in order to examine visibility reduction. Prevailing visibility is not practical to measure visibility depending on observer's expertness. Scattering visibility using Forward Scattering Meter (Belfort Visibility Sensor 6230) has been measured at Kwangan-Dong in Pusan and analysed since July, 1998. According to the analysis, the correlation coefficient(R) between prevailing visibility and scattering visibility was 0.7235. The visibility appeared that each frequency of poor visibility (under 6㎞) and good visibility(over 25㎞) was 10.6%, 9.7% on summertime in Pusan and the visibility range from 10km to 20km ranked high frequency as a half of whole ranges. The order of correlation coefficients between visibility and air pollutants are ranking CO, PM10 and NO2, that values are 0.5878, 0.5369, 0.5284 respectively. In meteorological factor, the case of poor visibility presented more weakly wind speed and higher relative humidity than the case of good visibility. The correlation coefficient between calculated visibility of multiple linear regression model and observed visibility was 0.7215. But the trend of calculated and observed visibility variation was similar with the exception of several good visibility cases.

Estimating dimensions of attractors are the most basic tools to analyze properties of chaotical dynamic systems. In this paper, we estimate correlation dimensions of meteorological variables, such as wind speed (v) and temperature (T) observed in Kimhae International Airport when the land-sea breeze circulation is appeared and find low non-integer values that reflect the deterministic chaos characterizing the dynamics. We compare the results with the correlation dimensions of 2-dimensional model that is calculated by finite element method. Though the correlation dimensions of the calculated wind speed (v) are less than those of the observed wind speed (v), we can suggest that the land-sea breeze circulation has not a unique mechanism. The land-sea breeze phenomenon is a complicated dynamics, which is constructed with various scale motions of atmosphere. In further research, we hope to find more accurate dynamics of land-sea breeze through wide observations and using of more sophisticated prediction models.