Research Article  Open Access
SeungJae Lee, Juwon Lee, Steven J. Greybush, Minseok Kang, Joon Kim, "Spatial and Temporal Variation in PBL Height over the Korean Peninsula in the KMA Operational Regional Model", Advances in Meteorology, vol. 2013, Article ID 381630, 16 pages, 2013. https://doi.org/10.1155/2013/381630
Spatial and Temporal Variation in PBL Height over the Korean Peninsula in the KMA Operational Regional Model
Abstract
Spatial and temporal variations in planetary boundary layer height (PBLH) over the Korean Peninsula and its surrounding oceans are investigated using a regional grid model operated at the Korea Meteorological Administration (KMA). Special attention is placed on daily maximum mixing height for evaluation against two radiosonde observation datasets. In order to construct a new highresolution PBLH database with 3hour time and 10 km spatial resolution, shortterm integrations with the regional model are carried out for a oneyear period from June 2010 to May 2011. The resulting dataset is then utilized to explore the seasonal patterns of horizontal PBLH distribution over the peninsula for one year. Frequency distributions as well as monthly and diurnal variations of PBLH at two selected locations are examined. This study reveals specific spatiotemporal structure of boundary layer depth over the Korean Peninsula for the first time at a relatively highresolution scale. The results are expected to provide insights into the direction for operational tuning and future development in the model boundary layer schemes at KMA.
1. Introduction
The diurnally evolving structure of the planetary boundary layer (PBL) for a typical synoptic highpressure system is described by Stull [1]. The depth of the PBL provides important information for numerical weather prediction (NWP) and atmospheric dispersion models. It has been used as a key parameter in the trigger function for convection in cumulus parameterization schemes in many NWP models (e.g., [2, 3]). Accurate prediction of PBL’s vertical extent is crucial in determining whether harmful gases (e.g., those erupted from a volcano) would reach the ground or not [4]. It affects nearsurface atmospheric pollutant concentrations (e.g., [5]), lowlevel moisture availability, and updraft conditions prior to thundershowers (e.g., [6]). In particular, daytime mixedlayer (ML) height has been regarded as the location of a capping temperature inversion atop the convective boundary layer. Raupach et al. [7] and Denmead et al. [8] formularized the relationship between carbon dioxide concentration and ML depth through entrainment processes.
To identify the ML top, there have been several methods for which radiosondes, wind profilers, and ceilometer/light detection and ranging (LIDAR) were utilized (e.g., [9–14]). These kinds of techniques to extract PBLH information have been also applied over the Korean Peninsula to surface flux [15–17], radiosonde [18, 19], LIDAR [20], and wind profiler [21] data. Most of the researchers analyzed characteristics of PBLH at one location over the peninsula and obtained several important results. For example, about 60% of daytime maximum mixing depths estimated at Korean radiosonde sites have been found to vary in a range of 1000–2000 m altitude with higher (lower) values in spring and fall (summer and winter) [18]. Low values in the summertime are associated with the East Asia monsoon starting in June, and wintertime low values are related to the cold air masses of the Siberian high [22]. More recently, Lee et al. [19] developed an empirical algorithm for automated determination of mixing height from midday (06UTC or 15LST) radiosonde temperature, relative humidity, and watervapor mixing ratio profiles.
However, all these studies until now have used observations and estimations of mixing depths at local stations, and thereby the spatial distribution of PBL height and its temporal variation was not well known over the Korean Peninsula. The central goal of this paper is to use an operational regional grid model to obtain a typical PBLH fields over the Korean Peninsula with the least computational cost. A quantitative analysis is also performed for the frequency, horizontal distribution, and monthly and seasonal variation in PBL height over the Korean Peninsula and its surrounding ocean for the first time. Data and methods are presented in Section 2, model evaluation is described in Section 3, and the spatial distribution and temporal variation are shown in Sections 4 and 5, respectively. Summary and concluding remarks are given in Section 6.
2. Data and Methodology
2.1. Model PBL Height Data
This study follows Berbery and Collini [23] premise stating that “lacking regional analyses, shortterm forecasts with a regional model will serve to add information in regions and times that otherwise would not be available.” Based on this idea, we use the KWRF Model, a regional grid model operated at Korea Meteorological Administration (KMA), to fill in regions with no PBLH information. The KWRF initially started from the Weather Research and Forecasting (WRF) system [24]. The KMA implemented its 10 km resolution version in May 2007 and has continuously modified it for data assimilation, physical parameterization, higher resolution, and other diverse applications [25]. The current version of KWRF is initialized by the global analyses from the KMA UM model (N320/L50) and its boundary condition is updated every 6 h. Horizontal resolution is 10 km with 40 vertical layers extending from the surface to 10 hPa, and a 60 s time step is used. The inflow/outflowdependent relaxation method is used for the lateral boundary condition, the radiation method for the top boundary condition, and the daily NOAA SST data for the marine boundary condition. The model physics configuration is composed of the Dudhia shortwave radiation scheme [26], the RRTM longwave radiation scheme [27], the WRF Single Moment 6class microphysics [28], the New KainFritsch scheme [29], the YSU boundary layer scheme [30], the MoninObukhov Janjić surface layer scheme [31–33], and the Noah land surface model [34]. WRF model behavior for East Asia is reported in Lee et al. [35], Choi et al. [36], and Jung et al. [37].
In general, PBL height is not a primary atmospheric variable such as temperature, humidity, and wind, and an objective data analysis scheme does not assimilate it into a numerical model. This is one main reason why the KMA data assimilation and prediction systems do not include the PBL height as an analysis variable. This significantly hinders us from knowing specifically the spatiotemporal structure of PBL height over the Korean Peninsula. Therefore, a strategy is needed in order to construct a basic dataset for the PBL height as accurately as possible. To this end, this study uses shortperiod integrations of the KWRF. Basically, the dataset consists of a series of 24 h forecasts at 3 h intervals starting at 00UTC and 12UTC for each day. Then, a "proxy" analysis field of PBLH can be established by connecting the two 12–21 h forecasts, starting at 00UTC and 12UTC, respectively (see Figure 1). The 00–09 h forecasts are not used to avoid the model spinup effect.
The nonlocal YSU boundary layer scheme uses prescribed eddy diffusivity profiles and permits countergradient terms in the situation of convective boundary mixing. Also, as an upgraded version of the MRF PBL scheme [2, 38], it considers an explicit representation of entrainment at the top of the boundary layer. Stability functions for momentum and heat are determined according to atmospheric stability that has four regimes based on the critical BulkRichardson number (Rib_{cr}): free convection , forced convection , damped mechanical turbulence , and stable atmosphere . In the scheme, PBL height is diagnosed by where is acceleration due to gravity, and are virtual potential temperatures near the model surface and at the lowest sigma level, respectively, is the virtual potential temperature, and is the horizontal wind speed at height of the PBL. In the KWRF model, boundary layer depths at every time steps are always defined through (1).
Figure 2 shows the model domain and its topographic features. The model topography reproduces the mountain range and the sharp slopes along the east coast of the Korean Peninsula. Seoul, Daejeon, Daegu, Gwangju, and Busan are the largest five cities in the Republic of Korea. A high mountain among Daejeon, Daegu, and Gwangju is known as “Jiri Mountain.” Between the mountains, lowaltitude regions with terrain heights of no more than 100 m are spread out toward the west and south seacoasts.
2.2. Observed Data
The KWRFproduced PBLH can be evaluated using remotesensing equipment such as radiosonde and wind profiler. In this study, two types of radiosonde datasets at Osan site (See Figure 2) are employed for the model evaluation because the data are available for a longtime period (since the 1980s) over Republic of Korea. The first dataset, produced by Choi and Baek [18], is composed of average daily maximum mixing height (DMMH) data for 10 years (1983–1992). Their DMMH data were calculated using a parcel method used by Holzworth [9], where the DMMH is determined as a height where a dry adiabatic line passing through the day’s maximum 2 m temperature intersects a morning (00UTC or 09LST) temperature profile. This 10year average DMMH data (hereafter, CB98) play an important role as climatological boundarylayer depth data over Republic of Korea.
The second dataset is also composed of DMMH data, but it was created by more direct method to use only midday environmental profiles of virtual potential temperature, mixing ratio and relative humidity measured at 06UTC (15LST). The midday (15LST) radiosonde balloon launches have been conducted and archived since 2001 and were used in this study for creating a oneyear dataset (June 2010 to May 2011) of DMMH. Hereafter, this DMMH dataset is called RAOB and Lee et al. [19] describe it in more depth.
3. Comparison of Observed and Modeled PBLH at Osan
3.1. Osan Radiosonde Observations
While 00UTC and 12UTC profiles of thermodynamic variables over the Korean Peninsula are archived and their spatial and temporal variation are well known, the characteristics of 06UTC (15LST) profiles (a new dataset) have not yet been reported. Since the 06UTC profiles are heavily used for model evaluation, this subsection provides readers with the basic information on the fullrange and seasonal variability of the profiles. For this, we visualize them in a spaghetti plot that is useful to trace the individual patterns of vertical profiles at Osan and easily check their daily behavior by season. Figure 3 shows a spaghetti plot of the radiosonde temperature profiles at the Osan site in each of the four seasons. The thin lines correspond to daily 06UTC soundings, and the thick black line is the seasonal mean profile. Seasonal mean of all temperature data (over all dates and levels) ranged from −9.5°C (winter) to 14.6°C (summer). It is noticeable that daily differences in temperature profiles are relatively small in summer with standard deviation of 6.9°C, compared to 10.1°C in spring.
(a)
(b)
(c)
(d)
Empirical orthogonal function (EOF) analysis, also known as principal component analysis (PCA), can illuminate the patterns that explain the most variability within a dataset. It is important to be cognizant of these vertical structures in the dataset that may be masked by examining a single metric, such as DMMH, in isolation. The EOF expansion was truncated to 4, because these vectors together explain 98% of the variance of the dataset; further details on the methodology for this section are described in the appendix. Figure 4 shows the EOF pattern vectors for the temperature in each season. The first pattern vector in all four seasons slowly increases its magnitude with height and has uniform sign. This represents a shift (warmer or cooler) from the profile mean, which is more pronounced in the upper atmosphere than lower atmosphere. It is noteworthy that the first eigenvector explains considerably less of the variance during summertime (79% compared to 94, 95, and 91%), as shown in Figure 3 where the summer season exhibits the smallest spread. This indicates that summertime profiles may exhibit more complex and small scale variations in the vertical. The second eigenvector is nearly monotonically increasing (decreasing) from the surface to the upper atmosphere and represents changes in lapse rate (and hence stability) in the profiles. The third pattern has its maximum at around 1500 m, with a sharply opposite sign near the surface, and could be indicative of a temperature inversion at that level. The fourth pattern has two peaks and can explain finer scale vertical structure in the data. Many of the eigenvectors change sign or slope at around 1000 m, indicating that this level is critical in separating two vertical regimes (i.e., PBL and free atmosphere) in temperature profile shape.
(a)
(b)
(c)
(d)
3.2. Model Evaluation at Osan
To measure the performance of the PBLH produced by the KWRF Model, a pointtopoint comparison was conducted at the Osan site where the two types of observation datasets are available. The results are summarized in the form of a scatter diagram (Figure 5) where each point represents the monthly DMMH value. The linear regression equations on the plot provide information on the relationship between the model and the observed DMMHs and indicate how well the KWRF Model DMMHs are compared to the observed.
(a)
(b)
Figure 5(a) shows the scatter plot that compares the monthly KWRF DMMH with the 10year average data (CB98). One merit of the CB98 dataset is its very low data loss because it was computed on a thermodynamic diagram using the parcel method. It can be seen that while the model DMMH tended to be less than the observed, they have remarkable correlation (0.7). It shows that 70% of the variance in the KWRF DMMH values is explained by the 10year climatological values, and 30% of those differences remain unexplained. But the regression equation has a slope greater than 1 and a large offset (~830 m). This means that the model simulates the observed pattern quite well, but not for the magnitude. The large offset should not be considered critical because CB98 is much smoothed data by using a 10year averaging.
Figure 5(b) shows the scatter plot between the monthly KWRF DMMH and the RAOB DMMH data for one year except the summertime rainy season (July–August 2010), when RAOB estimates of DMMH are nearly unavailable and substantially inaccurate due to bad weather conditions like precipitation and strong winds. It shows that 79% of the variance in the KWRF DMMH values is explained by the 1year radiosondeestimated values. The line in the graph has a slope of ~1, and an offset (~353 m) smaller than CB98. This means that, compared to CB98, the KWRF model exhibits better agreement with RAOB in terms of magnitude and pattern.
In order to determine the probability that the computed correlation occurred by chance, we conducted a significance test using a twotailed test based on the mutually exclusive hypothesis with the common significance level (=0.05) and the degrees of freedom (=10) and found that the critical value is ~0.58. Since the two correlation values in Figure 5 are greater than 0.58, we can conclude that they are statistically significant and the probability of a chance occurrence is less than 5 out of 100.
It is noted from Figure 5 that overall the model DMMHs tend to be lower than the radiosonde estimates with negative offsets. The model’s underestimation of PBL height in the present study is firstly because the sounding data were mainly from clear days, while the modeled PBLH data were from all days. Second, the radiosondeestimated DMMH is a point value, whereas the model DMMH represents a grid box () average, so modelgridded values of DMMH and large locality of radiosonde sounding can cause the lower model DMMH than the observed. Third, the KWRF model itself can have systematic errors in simulating boundary layer depth. Lastly, there can be potential error sources including an incorrect landcover type and topographic effects in the model. This underestimating tendency is opposite to some studies (e.g., [39]) where the WRF with YSU PBL tend to overestimate the PBL depth. However, the overestimation of PBL height is somewhat unclear in other studies. For example, Shin and Hong [40] showed that YSU PBL had the second lowest depth during daytime among two nonlocal schemes and three local schemes. Such an overestimating or underestimating tendency seems to be dependent on observational system as well (e.g., radiosonde, wind profilers, LIDAR, etc.) used in the PBL height evaluation. When all these things are considered, the levels of correlation and offset shown in this section are acceptable.
4. Horizontal Distribution of the Model Boundary Layer Height
4.1. Diurnal Range and Morning Growth Rates
For simplicity, the diurnal range (daily maximum minus daily minimum) of PBLH is defined by 15LST PBLH minus 03LST PBLH. Figure 6 shows the spatial distribution of the diurnal range of PBLH by season. Overall, it shows that the PBLH field is influenced by seasonality and topographic characteristics. Compared to land PBLH, the maritime PBLH shows the lowest diurnal range values because of the ocean’s large heat capacity and the marine boundary layer it develops, and these two effects can even lead to negative values at some local points (e.g., southeastern sea of the peninsula in fall) due to the timelagged cycles of PBLH over land and ocean.
(a) Spring
(b) Summer
(c) Fall
(d) Winter
However, the diurnal range increases sharply inland and exceeds 1000 m over relatively low terrain including over the cities of Seoul, Daejeon, and Daegu in the warm season. Over land, springtime has the largest diurnal amplitude throughout the peninsula. Summertime also showed large diurnal amplitudes, but the spatial coverage was restricted within Republic of Korea. In summer, the spatial gradients and horizontal heterogeneity of diurnal range become large and this results in the largest diurnal range occurring around Seoul. A strong gradient was seen over the western shoreline near Seoul city where the diurnal amplitude varies from a low of ~25 m to values of 1000–1500 m. The difference in spatial heterogeneity of diurnal range is related to heavy rainfall in the rainy season during summertime and frequent clear skies in springtime. Wintertime has the smallest diurnal amplitude especially over mountainous areas. Fall has features intermediate between summer and winter in time and space.
Figure 7 shows the morning growth rate of PBLH, which is defined by 12LST PBLH minus 06LST PBLH. The growth rate of the PBLH during the morning hours is an important factor in the development of high ozone concentrations at ground level (Berman et al., [41]). Overall, the growth rate has similar patterns to the diurnal range in its horizontal distribution. Both Koreas show similar growth rates in spring, but in summer Republic of Korea showed more rapid growth rates than Democratic People’s Republic of Korea. The largest growth rates appeared inland over large cities and mountain valleys, where they exceeded 1100 m at some locations. Growth rates over high mountain ranges were lower (~500 m) and declined rapidly as the shoreline is approached, with values of a few tens of meters along the coastline. This pattern is consistent with the findings of Berman et al. [41] that the effect of the marine boundary layer is to depress surface temperature and thermal convection. Values of zero or negative growth rates appearing over much of the ocean indicate that the direct effect of incoming solar radiation is very small on the PBL development at those areas over the sea. Spatial distributions of nocturnal PBLH showed very little seasonal variation and they were not shown.
(a) Spring
(b) Summer
(c) Fall
(d) Winter
4.2. Daily Maximum Mixing Height (DMMH)
Since the PBLH data used in this study have a 3hour resolution, 06UTC (15LST) PBLH can be regarded as DMMH, which usually corresponds to maximum afternoon mixing depth in the case of noprecipitation days. Figure 8 displays the horizontal distribution of DMMH for the seasonal mean. In general, the daytime boundary layer is deeper over land as expected except winter. Large DMMHs can be found in nonmountainous regions such as the Seoul and Daegu metropolitan cities. This is related to the urban heat island effect as shown by Kim and Baik (2005) [42] and Lee and Kim [43]. In winter, however, the boundary depths are larger over the ocean, and this can be easily understood from the relatively high sea surface temperature due to the larger heat capacity of the sea.
(a) Spring
(b) Summer
(c) Fall
(d) Winter
Over the sea, the wintertime marine PBL depths are smaller in the West Sea than in the East Sea and tend to be smaller at high latitudes. This is affected mainly by the sea surface temperature distribution. Over land, the topographic effect is clearly seen regardless of season, especially along the high mountain ranges. This results in relatively low DMMHs in the mountainous areas in all seasons and reflects that the mountains have cooler daytime surface temperatures.
Figure 9 shows seasonal standard deviation of DMMH. It should be noticed that seacoast areas, especially the southeastern areas, show large variability in DMMH in spring. In the other seasons, lowaltitude regions show large variability compared to highaltitude regions. Over sea, the variability of DMMH is generally small compared to land, except for spring.
(a) Spring
(b) Summer
(c) Fall
(d) Winter
5. Frequency Distribution and Monthly Variation of the Model Boundary Layer Height at Two Local Points
5.1. Frequency Distribution of PBLH and DMMH
In order to see the distribution of PBLH data at all times of day (3 h interval) during the oneyear period, histograms of PBLHs were calculated in Figure 10 for Osan and Gwangju (see Figure 2). It shows the frequency of PBLH more quantitatively using the same bin width as Choi and Baek [18]. First of all, it can be seen that the occurrence frequency is high (low) for smaller (larger) PBLH bins. The frequency distribution of PBLH appears to follow approximately exponential distributions regardless of season. Both sites show the largest frequency in 100–500 m bin, occupying a relative frequency of ~44%.
(a) Spring
(b) Summer
(c) Fall
(d) Winter
(e) 1 year
Figure 11 shows histograms of DMMHs for the two sites and provides detailed information on occurrence frequency and relative frequency of DMMH for each site. They reveal that the frequency of DMMH at Osan (Gwangju) had the most occurrences in the 500–1000 m (1000–1500 m) bin, occupying relative frequency of ~62(56)%. Unlike Figure 10, the histograms are somewhat close to Gaussian distributions and this characteristic is more clear at Gwangju in the fall season. In summer, the shape of the distribution is considerably broader, implying larger variability. In winter, the distribution is skewed toward low DMMH values and deviates from the Gaussianlike shape.
(a) Spring
(b) Summer
(c) Fall
(d) Winter
(e) 1 year
The oneyear mean DMMH at the two sites shows again a Gaussianlike distribution (Figure 11(e)). Gwangju shows relatively lower (higher) frequency at small (large) PBLH bins, compared to Osan. Possible reasons include that Gwangju is geographically located south of Osan and thus climatologically warmer than Osan or that Gwangju is more unstable and convective than Osan during the daytime. Also, the Gwangju site showed a pattern more similar to a Gaussian distribution and had its mode near 1100 m. The Osan site shows larger skewness than Gwangju site and has its first mode at 900 m. The second peak at a height of 300 m at both sites is an interesting feature and remains for further analysis.
5.2. Monthly and Diurnal Variation of DMMH
In order to see the temporal trend in model DMMH, monthly and seasonal mean and standard deviation were calculated. Table 1 shows that the highest DMMH occurs in the spring season at both sites, and the spring mean value is ~1315 m (1406 m) for the Osan (Gwangju) site. The lowest DMMH appears in the winter season at both sites, and the winter mean value is ~630 m (691 m) for Osan (Gwangju) site. Summer and fall seasons show intermediate values of DMMH between the spring and winter values. Summer DMMH values are most spread out and the standard deviation is large, indicating that summer has the largest variability in the DMMH dynamics compared to the other seasons, whereas winter DMMHs are least spread out and show the smallest dispersion in DMMH. Compared with Choi and Baek [18], the most different aspect is that the KWRF DMMH is much smaller in winter. This indicates that the model PBL mixing is generally weak in the cold season and presumably nighttime stable regimes.

Figure 12 shows variations in monthly mean DMMH at Osan and Gwangju with standard deviations denoted by the vertical error bars. Overall, the two sites show a similar behavior in the pattern and magnitude. The monthly mean DMMHs are high in March and June at both sites, forming two peaks in the graphs, due to fair weather conditions in spring and early summer. From December to January both sites exhibit the lowest values, and this is basically attributed to low surface temperature due to small insolation in the winter season. The lowest values are ~600 m and the highest values go over 1500 m, making the maximum seasonal differences as large as ~900 m.
It should be noted that there are two depressions of the mean DMMHs in May and July during the warm season. The one depression in July is clearly related to the increase of cloud cover and precipitation during the monsoon season (normally late June through early August in Korea), as shown in Figure 13(a). The large amount of precipitation at Osan in September 2010 was caused by the 21September heavy rainfall event which occurred in the central part of the Korean peninsula (see the satellite image in Figure 13(b)). It can be seen that the recordbreaking heavy rainfall contributed to lower the mean DMMH more at Osan than at Gwangju. This can be confirmed also by examining vertical radiosonde soundings measured at both sites on September 21, 2010. Figure 14 shows that Gwangju had a normal mixed layer with a depth of 1385 m on the day, whereas Osan had a highly disturbed boundary layer with large wind speed and fluctuating vertical profiles.
(a)
(b)
(a)
(b)
On the other hand, the other depression in May is not related to monsoon rainfall. Figure 12 displays a decrease in monthly mean DMMH from March to May. This is associated with typical springtime changes in turbulent heat fluxes and surface wind over the Korean Peninsula. Figure 15 shows changes in monthly mean values of several micrometeorological variables measured at the Haenam Farmland site (KoFlux/AsiaFlux site code: HFK; refer to the web site http://asiaflux.net and http://ncam.kr/page/koflux/database for detailed information) from March to May 2011. It shows that evapotranspiration increased by ~60% while sensible heat fluxes gradually decreased from March to May, resulting in decreasing in Bowen ratio. This indicates that plants started active growth and photosynthesis with transpiration depending on more insolation, allotting more energy toward latent heat rather than sensible heat as time goes by. The reduction in sensible heat fluxes leads to reduction in buoyancy. Wind speed also plays an important role in PBL depth through mechanical turbulence generation. Wind speed is one of the main variables that can identify the PBL depth under strongly convective conditions (e.g., Stull, [1]). Wind speed near the surface decreases from March, when the strongest nearsurface winds occur due to frequently passing highpressure systems over the Korean Peninsula (e.g., [44, 45]), to May. Thus, a reduction in both sensible heat fluxes and nearsurface wind speed accounts for the depression of mean DMMH in May.
6. Summary and Concluding Remarks
Planetary boundary layer height is an important variable in modeling weather phenomena and air quality. Current understanding of the spatial distribution and temporal change of PBLH over the Korean Peninsula is very partial and incomplete. This paper uses shortterm forecasts from the operational regional KWRF model to obtain information on the spatial distribution and temporal variations in PBLH that are unavailable with the current observations and analyses over the region. To validate the model, we make use of previously unpublished 15LST Korean radiosonde profiles. Their seasonal and vertical modes of variability are presented using EOF analysis.
Model integrations for June 2010–May 2011 suggest that PBLHs exhibited the largest diurnal range, morning growth rate and daytime maximum values over lower terrain inland in spring, with the smallest values over higher terrain in winter. Strong horizontal gradients of PBLH were observed along the coastline especially in summer. The surrounding ocean showed very low (even slightly negative) values in the diurnal range, morning growth rate, and daytime maximum fields of PBLHs due to its low surface temperature and thermally depressive marine boundary layers. Large DMMHs were located over nonmountainous regions (especially over the cities of Seoul, Daejeon, and Daegu as a result of the urban heat island) in spring, while the highest DMMHs were located over the East Sea and the South Sea in winter. Regarding the variability of DMMH, it was large in the southeastern seacoast areas in spring and was small over all seas surrounding the peninsula in summer.
Local analyses of model PBLHs and DMMHs at the Osan and Gwangju points showed that histograms of model PBLH tended to have exponential distributions of frequency at both locations, while those of model DMMH exhibited almost Gaussian distributions especially at Gwangju. Monthly variation of DMMH had its main peaks in March and June and its lowest values in winter, with the largest variability in summer. There were two secondary minimums of DMMHs in May and July as well. The minimum DMMHs in winter and July are mainly attributed to the lowest surface temperature and the rainy season, respectively. It was suggested that the secondary minimum of mean DMMH in May can be accounted for by the decreasing trend in both Bowen ratio and nearsurface wind speed during spring over the peninsula.
The current study provides meso and regional scale analyses of PBLH fields over the Korean Peninsula for the first time. The complex structure of PBLH variation in space and time is reasonably shown using 3hourly model forecasts for one year, but additional work must include a comparison with observation data over a longerterm period. It is likely that the use of higherresolution model PBLH estimates with a more frequent output interval (e.g., hourly) may be necessary; in this regard, data from a new operational regional model (with a horizontal resolution of 1.5 km) being developed at KMA may play an important role toward better estimates of PBLH in the region. The revealed spatiotemporal characteristics of PBLH from the KWRF Model will be used as an important reference for future operational tuning and development of a new boundary layer scheme in the model. The constructed PBLH dataset can be applied to many related research including boundary layer CO_{2} dynamics and air pollutants dispersion.
Appendix
EOF Analysis of Radiosonde Temperature Profiles
Empirical orthogonal function (EOF) analysis, also known as principal component analysis (PCA), can illuminate the patterns that explain the most variability within a dataset. In this paper, EOF analysis is applied to 06UTC (15LST) Osan radiosonde temperature profiles. The data was first divided by season: June to August 2010 (summer), September to November 2010 (autumn), December 2010 to February 2011 (winter), and March to May 2011 (spring). For EOF analysis, data must be on a common vertical grid. Thus, the data were interpolated linearly to every 100 feet from 100 to 4000 feet, and records with missing data were ignored. Next, the mean profile was subtracted from the dataset so that the EOFs explain variability about the mean. Then, EOF analysis was performed separately on each variable for each season by computing the covariance matrix in space and its eigenvalue decomposition, and ordering the eigenvalues and eigenvectors by explained variance. For brevity, only temperature results are reported here. EOF pattern vectors can be rotated to relax the orthogonality constraint, which aids the physical interpretation of results [46, 47]. Without rotation, the orthogonality constraint of EOFs leads to the emergence of Buell patterns, which are characteristic of domain shape and can appear regardless of the field under consideration [46]; here we use the varimax technique. It should also be noted that the sign of the pattern vector is irrelevant for physical interpretation, only the shape.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This study was supported by the Korea Meteorological Administration Research and Development Program under Grant CATER 20133030. The HFK data were provided by KoFlux from the projects funded by Korea Meteorological Administration and the A3 program of National Research Foundation of Korea. The first author thanks Juyeol Yun and Aastha (Bindu Malla Thakuri) for their valuable KoFlux site management.
References
 R. B. Stull, An Introduction to Boundary Layer Meteorology, Kluwer Academic Publishers, Norwell, Mass, USA, 1988.
 S.Y. Hong and H.L. Pan, “Nonlocal boundary layer vertical diffusion in a mediumrange forecast model,” Monthly Weather Review, vol. 124, no. 10, pp. 2322–2339, 1996. View at: Google Scholar
 S.H. Shin and K.J. Ha, “Effects of spatial and temporal variations in PBL depth on a GCM,” Journal of Climate, vol. 20, no. 18, pp. 4717–4732, 2007. View at: Publisher Site  Google Scholar
 H. Nagai, H. Terada, and M. Chino, “Relation between nasty smell at Kanto area on 28 Aug 2000 and eruption at Miyake Island: examination by numerical simulation,” Tenki, vol. 48, pp. 227–230, 2001 (Japanese). View at: Google Scholar
 D. J. Jacob, Introduction to Atmospheric Chemistry, Princeton University Press, 1999.
 R. Latha and B. S. Murthy, “Boundary layer signatures of consecutive thunderstorms as observed by Doppler sodar over western India,” Atmospheric Research, vol. 99, no. 2, pp. 230–240, 2011. View at: Publisher Site  Google Scholar
 M. R. Raupach, O. T. Denmead, and F. X. Dunin, “Challenges in linking atmospheric CO_{2} concentrations to fluxes at local and regional scales,” Australian Journal of Botany, vol. 40, no. 45, pp. 697–716, 1992. View at: Google Scholar
 O. T. Denmead, M. R. Raupach, F. X. Dunin, H. A. Cleugh, and R. Leuning, “Boundary layer budgets for regional estimates of scalar fluxes,” Global Change Biology, vol. 2, no. 3, pp. 255–264, 1996. View at: Google Scholar
 G. C. Holzworth, “Estimates of mean maximum mixing depths in the contiguous United States,” Monthly Weather Review, vol. 92, pp. 235–242, 1963. View at: Google Scholar
 J. L. Heffter, “Transport layer depth calculations,” in Proceedings of the 2nd Joint Conference on Applications of Air Pollution Meteorology, New Orleans, La, USA, 1980. View at: Google Scholar
 W. M. Angevine, A. B. White, and S. K. Avery, “Boundarylayer depth and entrainment zone characterization with a boundarylayer profiler,” BoundaryLayer Meteorology, vol. 68, no. 4, pp. 375–385, 1994. View at: Publisher Site  Google Scholar
 P. Seibert, F. Beyrich, S.E. Gryning, S. Joffre, A. Rasmussen, and P. Tercier, “Review and intercomparison of operational methods for the determination of the mixing height,” Atmospheric Environment, vol. 34, no. 7, pp. 1001–1027, 2000. View at: Publisher Site  Google Scholar
 S. Liu and X. Z. Liang, “Observed diurnal cycle climatology of planetary boundary layer height,” Journal of Climate, vol. 23, no. 21, pp. 5790–5809, 2010. View at: Publisher Site  Google Scholar
 M. Haeffelin, F. Angelini, Y. Morille et al., “Evaluation of mixingheight retrievals from automatic profiling lidars and ceilometers in view of future integrated networks in Europe,” BoundaryLayer Meteorology, vol. 143, no. 1, pp. 49–75, 2012. View at: Publisher Site  Google Scholar
 S.U. Park and I.H. Yoon, “Estimation of atmospheric boundary layer parameters using routinely available meteorological data,” Journal of the Korean Meteorological Society, vol. 27, pp. 32–54, 1991. View at: Google Scholar
 J.B. Lee, “A study on the air pollution potential in the central part of Korea,” Journal of Korea Air Pollution Research Association, vol. 7, pp. 41–47, 1991 (Korean). View at: Google Scholar
 Y.G. Kim, “Calculation of the convective mixed layer by estimation of sensible heat flux,” Journal of the Korean Environmental Sciences Society, vol. 7, pp. 639–645, 1998 (Korean). View at: Google Scholar
 J.S. Choi and S.O. Baek, “An approach to estimate daily maximum mixing height (DMMH) in Pohang, Osan, and Kwangju areas—analysis of 10 years data from 1983 to 1992,” Journal of Korea Air Pollution Research Association, vol. 14, pp. 379–385, 1998 (Korean). View at: Google Scholar
 S.J. Lee, J. Kim, and C.H. Cho, “An automated monitoring of atmospheric mixing height from routine radiosonde profiles over South Korea using a webbased data transfer method,” Environmental Monitoring and Assessment. In press. View at: Google Scholar
 S.C. Yoon, J.G. Won, S.W. Kim, and G.M. Lim, “Measurement of mixed layer height using Lidar,” in Proceedings of Korean Society for Atmospheric Environment, pp. 436–437, 1999. View at: Google Scholar
 S.J. Lee and H. Kawai, “Mixing depth estimation from operational JMA and KMA windprofiler data and its preliminary applications: examples from four selected sites,” Journal of the Meteorological Society of Japan, vol. 89, no. 1, pp. 15–28, 2011. View at: Publisher Site  Google Scholar
 K.D. Min, S.H. Kim, K.E. Kim, and B.H. Kwon, “Seasonal and local characteristics of atmospheric mixed layer over Kyungpook Province,” Journal of Korean Meteorological Society, vol. 35, pp. 539–548, 1999. View at: Google Scholar
 E. H. Berbery and E. A. Collini, “Springtime precipitation and water vapor flux over southeastern South America,” Monthly Weather Review, vol. 128, no. 5, pp. 1328–1346, 2000. View at: Google Scholar
 W. C. Skamarock, J. B. Klemp, J. Dudhia et al., “A Description of the Advanced Research WRF Version 3,” NCAR Technical Note NCAR/TN475+STR, National Center for Atmospheric Research, Boulder, Colo, USA, 2008, http://www.mmm.ucar.edu/wrf/users/docs/arw_v3.pdf. View at: Google Scholar
 I.H. Cho, H.D. Yoo, and J.O. Lim, “Implementation and test run of the next generation regional forecast model (KWRF),” in Proceedings of Korean Meteorological Society, pp. 368–369, 2005 (Korean). View at: Google Scholar
 J. Dudhia, “Numerical study of convection observed during the winter monsoon experiment using a mesoscale twodimensional model,” Journal of the Atmospheric Sciences, vol. 46, no. 20, pp. 3077–3107, 1989. View at: Google Scholar
 E. J. Mlawer, S. J. Taubman, P. D. Brown, M. J. Iacono, and S. A. Clough, “Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlatedk model for the longwave,” Journal of Geophysical Research D, vol. 102, no. 14, pp. 16663–16682, 1997. View at: Google Scholar
 S.Y. Hong and J.O. Lim, “The WRF singlemoment 6class microphysics scheme (WSM6),” Journal of the Korean Meteorological Society, vol. 42, pp. 129–151, 2006. View at: Google Scholar
 J. S. Kain and J. Kain, “The KainFritsch convective parameterization: an update,” Journal of Applied Meteorology, vol. 43, no. 1, pp. 170–181, 2004. View at: Google Scholar
 S.Y. Hong, Y. Noh, and J. Dudhia, “A new vertical diffusion package with an explicit treatment of entrainment processes,” Monthly Weather Review, vol. 134, no. 9, pp. 2318–2341, 2006. View at: Publisher Site  Google Scholar
 A. S. Monin and A. M. Obukhov, “Osnovnye zakonomernosti turbulentnogo peremeshivanija v prizemnom sloe atmosfery,” Trudy Geofizicheskogo Instituta, Akademiya Nauk SSSR, vol. 24, no. 151, pp. 163–187, 1954. View at: Google Scholar
 Z.I. Janjić, “The surface layer in the NCEP eta model,” in Proceedings of the 11th Conference on Numerical Weather Prediction, Norfolk, Va, USA, August 1996. View at: Google Scholar
 Z.I. Janjić, “Nonsingular implementation of the MellorYamada level 2.5 scheme in the NCEP Meso model,” NCEP Office Note 437, 2002. View at: Google Scholar
 F. Chen and J. Dudhia, “Coupling and advanced land surfacehydrology model with the Penn StateNCAR MM5 modeling system—part I: model implementation and sensitivity,” Monthly Weather Review, vol. 129, no. 4, pp. 569–585, 2001. View at: Google Scholar
 D.K. Lee, D.Y. Eom, J.W. Kim, and J.B. Lee, “Highresolution summer rainfall prediction in the JHWC realtime WRF system,” AsiaPacific Journal of Atmospheric Sciences, vol. 46, no. 3, pp. 341–353, 2010. View at: Publisher Site  Google Scholar
 H.Y. Choi, J.H. Ha, D.K. Lee, and Y.H. Kuo, “Analysis and simulation of mesoscale convective systems accompanying heavy rainfall: the goyang case,” AsiaPacific Journal of Atmospheric Sciences, vol. 47, no. 3, pp. 265–279, 2011. View at: Publisher Site  Google Scholar
 S.H. Jung, E.S. Im, and S.O. Han, “The effect of topography and sea surface temperature on heavy snowfall in the Yeongdong region: a case study with high resolution WRF simulation,” AsiaPacific Journal of Atmospheric Sciences, vol. 48, pp. 259–273, 2012. View at: Google Scholar
 S.J. Lee, D. F. Parrish, S.Y. Park et al., “Effects of 2m air temperature assimilation and a new nearsurface observation operator on the NCEP gridpoint statisticalinterpolation system,” AsiaPacific Journal of Atmospheric Sciences, vol. 47, no. 4, pp. 353–376, 2011. View at: Publisher Site  Google Scholar
 G. Svensson, A. A. M. Holtslag, V. Kumar et al., “Evaluation of the diurnal cycle in the atmospheric boundary layer over land as represented by a variety of singlecolumn models: the second GABLS experiment,” BoundaryLayer Meteorology, vol. 140, no. 2, pp. 177–206, 2011. View at: Publisher Site  Google Scholar
 H. H. Shin and S.Y. Hong, “Intercomparison of planetary boundarylayer parametrizations in the WRF model for a single day from CASES99,” BoundaryLayer Meteorology, vol. 139, no. 2, pp. 261–281, 2011. View at: Publisher Site  Google Scholar
 S. Berman, J.Y. Ku, and S. T. Rao, “Spatial and temporal variation in the mixing depth over the Northeastern United States during the summer of 1995,” Journal of Applied Meteorology, vol. 38, no. 12, pp. 1661–1673, 1999. View at: Google Scholar
 Y.H. Kim and J.J. Baik, “Maximum urban heat island intensity in Seoul,” Journal of Applied Meteorology, vol. 41, pp. 651–659, 2002. View at: Google Scholar
 S.H. Lee and H.D. Kim, “Effects of regional warming due to urbanization on daytime local circulations in a complex basin of the Daegu metropolitan area, Korea,” Journal of Applied Meteorology and Climatology, vol. 47, no. 5, pp. 1427–1441, 2008. View at: Publisher Site  Google Scholar
 J.H. Kim and I.U. Chung, “Study on mechanisms and orographic effect for the springtime downslope windstorm over the Yeongdong region,” Atmosphere, vol. 16, pp. 67–83, 2006 (Korean). View at: Google Scholar
 S.J. Lee and Y.C. Kim, “A numerical forecast and verification of the aircraft turbulence observed over South Korea,” AsiaPacific Journal of Atmospheric Sciences, vol. 38, pp. 493–507, 2002 (Korean). View at: Google Scholar
 M. B. Richman, “Rotation of principal components,” Journal of Climatology, vol. 6, no. 3, pp. 293–335, 1986. View at: Google Scholar
 S. J. Greybush, S. E. Haupt, and G. S. Young, “The regime dependence of optimally weighted ensemble model consensus forecasts of surface temperature,” Weather and Forecasting, vol. 23, no. 6, pp. 1146–1161, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 SeungJae Lee et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.