Search

Precipitation plays a crucial role in the global energy and water cycle and has important implications for food, water, and energy security. To enhance our understanding of the water cycle, it is invaluable to have a comprehensive historical record of precipitation. However, obtaining such records, especially for the period before the Industrial Revolution, can be challenging. During the Joseon Dynasty, Korea established a network for measuring rainfall and recorded this information in historical documents known as Seungjeongwon Ilgi and Ilseongnok. Recently, these documents have been digitized, providing us with daily precipitation data for Seoul spanning 130 years, from 1778 to 1907. By combining and analyzing these two documents, we were able to address inconsistencies found in previous studies and improve the quality of the data. Notably, this dataset is free of any missing values, making it the longest daily precipitation record in the world before the Industrial Revolution. Its availability to the public holds great potential for climate research in East Asia during the late Little Ice Age. [ABSTRACT FROM AUTHOR]

Results from eight regional climate models (RCMs) participating in the Impact of Initialized Land Temperature and Snowpack on Sub-seasonal to Seasonal Prediction (LS4P) initiative of the Global Energy and Water Exchanges (GEWEX) are examined and compared with observations over the Tibetan Plateau (TP). The RCM common domain covers most areas of East Asia with a horizontal resolution of 20–30 km. The model simulation covers a period from April to September in each year between 1991 and 2015. This study explores the RCMs' ability for seasonal climate simulation over the TP, focusing on the summer monsoon climate as part of the LS4P initiative. An intercomparison is made among eight RCMs for precipitation, surface air temperature, mid-troposphere atmospheric circulation, moisture conditions, and surface energy fluxes. It shows that the downscaling characteristics differ significantly between two major RCM types. The RegCM4 models show positive precipitation biases over the entire TP, especially over the south and southeast TP, while the WRF models mostly show both positive and negative precipitation biases over the TP with relatively high spatial correlation between simulated and observed precipitation. The multi-model ensemble mean produces overall smaller precipitation biases than most individual RCMs, with the largest biases over the southeastern TP, and smaller surface air temperature biases over most areas of the TP, especially over the central and southwestern TP. Moreover, the ensemble mean can better reproduce the inter-annual variation of precipitation and surface air temperature than most RCMs with proper magnitude. Sensitivity analyses using RegCM4 with different physics parameterizations show that varying land and cumulus schemes may induce large precipitation differences over the TP by affecting moisture and atmospheric circulation conditions in the lower and upper troposphere, respectively. Moreover, turbulent heat and radiation fluxes differences are associated with the temperature differences between different RegCM4 models. [ABSTRACT FROM AUTHOR]

Based on eddy-covariance measurements over a glacier, we investigate the scalar flux dissimilarity between sensible heat and water vapor transport in a stably stratified flow. The scalar flux correlation coefficient R F is used as a measure of variable levels of the flux similarity, which are often elevated due to a rising degree of the kinetic anisotropy of turbulence. Moreover, sensible heat is transported more efficiently than water vapor; and transport efficiencies of these two scalars are separated in terms of their variability with the velocity aspect ratio. Compared with air temperature fluctuations, turbulence characteristics of the water vapor concentration are subject to a more pronounced modification because of distinct magnitudes of R F . An innovative method is employed for connecting quadrant analysis and cospectral analysis, so that the hyperbolic quadrant-hole size can be coupled to the natural frequency underlying the fast Fourier transform. Then, we introduce a hypothetical octant hole whose size is invoked as a metric for the amplitude scale of fluctuating scalar fluxes. The contributions to R F are quantified for a variety of eddy structures that are associated with different ranges of the amplitude scale. Regarding larger-amplitude fluxes due to heated drier air parcels in descending motions, reductions in R F correspond to increasing flux fractions for water vapor, whereas the flux fractions for sensible heat are largely unchanged. Overall, a more substantial portion of the changes in R F can be ascribed to smaller-amplitude fluxes due to cooled moister air parcels and heated drier air parcels being involved, respectively, in ascending and descending motions. Reductions in R F relate to the flux fractions of a decreasing magnitude for sensible heat but of an increasing magnitude for water vapor. Highlights: In a stably stratified flow, dissimilar transport of sensible heat and water vapor is associated with anisotropy properties of turbulence. A new approach is developed for a scrutiny of scalar transport, whereby quadrant analysis and cospectral analysis can be interconnected. Extending octant analysis to scalar turbulence identifies eddy structures that exhibit distinct behaviors indicative of the flux dissimilarity. Comparatively small-amplitude fluxes are accountable for a substantial portion of the changes in the overall level of scalar flux correlation. [ABSTRACT FROM AUTHOR]

Severe spatiotemporal heterogeneity of emissions sources and limited measurement networks have been hampering the monitoring and understanding of CO2 fluxes in large cities, a great concern in climate research as big cities are among the major sources of anthropogenic CO2 in the climate system. To understand the CO2 fluxes in Seoul, Korea, CO2 fluxes at eight surface energy balance sites, six urban (vegetation-area fraction < 15%) and two suburban (vegetation-area fraction > 60%), for 2017–2018 are analyzed and attributed to the local land-use and business types. The analyses show that the CO2 flux variations at the suburban sites are mainly driven by vegetation and that the CO2 flux differences between the urban and suburban sites originate from the differences in the vegetation-area fraction and anthropogenic CO2 emissions. For the CO2 fluxes at the urban sites; (1) vehicle traffic (traffic) and heating-fuel consumption (heating) contribute > 80% to the total, (2) vegetation effects are minimal, (3) the seasonal cycle is driven mainly by heating, (4) the contribution of heating is positively related to the building-area fraction, (5) the annual total is positively (negatively) correlated with the commercial-area (residential-area) fraction, and (6) the traffic at the commercial sites depend further on the main business types to induce distinct CO2 flux weekly cycles. This study shows that understanding and estimation of CO2 fluxes in large urban areas require careful site selections and analyses based on detailed consideration of the land-use and business types refined beyond the single representative land-use type widely-used in contemporary studies. [ABSTRACT FROM AUTHOR]

Based on observations over an alpine glacier, we investigate the turbulent flux dissimilarity between momentum and sensible heat transfer in a stably stratified katabatic flow. The flux correlation coefficient RF is employed as a measure of variable levels of the flux similarity, which are found to be modulated by the anisotropy of turbulence. In the katabatic wind regime over this glacier, heat is transported more efficiently than momentum. Based on quadrant analysis, the sweep–ejection differences in the flux fraction constitute turbulence characteristics in terms of the velocity aspect ratio, which unravel dissimilar transport of momentum and heat. Moreover, an innovative method is established for connecting quadrant analysis and cospectral analysis, whereby the hyperbolic quadrant hole size is coupled to the frequency underlying the fast Fourier transform. In extending applications of octant analysis, we introduce a hypothetical octant hole, whose size is solicited as a proxy for the amplitude scale of fluctuating fluxes. The contributions to R F can then be quantified for eddy structures that are associated with different amplitude scales. The katabatic flow structures identified from octant analysis differ in their behaviour so helping illuminate the outcome of flux dissimilarity. Exhibited as a statistical behaviour regardless of amplitude scale, along-wind rapid motions of heated air parcels can modify fractional contributions to the heat flux instead of the momentum flux, with reductions in R F related to decreasing heat-flux fractions. Besides, along-wind slow motions of cooled air parcels cannot modify the flux fractions for both momentum and heat. Thus, the flux dissimilarity due to low-speed low-temperature eddies cannot be explicable in terms of the flux fractions alone. These findings are an incipient step towards physical understanding of the turbulent flux dissimilarity for a stably stratified katabatic flow. [ABSTRACT FROM AUTHOR]

The dilution effect caused by boundary-layer evolution over land has strong influences on air quality. Accurate and continuous measurements of the boundary-layer height over urban areas are therefore needed for complete air-quality assessments. Commercial ceilometers, in combination with reliable and simple methodologies, can be used to retrieve the mixed-layer height, and represent a means of obtaining information on vertical mixing and atmospheric structure above cities. Here, we evaluate various retrieval algorithms based on the gradient method against high-temporal-resolution radiosonde observations. Based on the results, we propose a simple algorithm by using the gradient method, the correction of background noise and the moving averages, with the minimum number of parameters that need to be adjusted to the local properties and the instrument itself. The algorithm is adjusted for Seoul, Korea, and improves the retrieval performance by reducing high-frequency noise. The algorithm is used to investigate the relationship between the evolution of the daytime mixed-layer height and air pollution under a two-layer mixing model where changes in concentration depend only on the urban boundary-layer growth and air entrainment from the free atmosphere. Using 2 months of ceilometer retrievals of mixed-layer height and air-quality data from across the city, we find strong negative correlations for primary emitted pollutants such as NO2, CO, SO2, and particulate matter smaller than 10 µm, and a modest positive correlation for O3. The results provide insight into the significant influence of urban boundary-layer evolution on Seoul's air quality. [ABSTRACT FROM AUTHOR]

A substantial over-prediction bias at low-to-moderate wind speeds in the Weather Research and Forecasting (WRF) model has been reported in the previous studies. Low-level wind fields play an important role in dispersion of air pollutants, including radionuclides, in a high-resolution WRF framework. By implementing two subgrid-scale orography parameterizations (Jimenez and Dudhia in J Appl Meteorol Climatol 51:300–316, 2012; Mass and Ovens in WRF model physics: problems, solutions and a new paradigm for progress. Preprints, 2010 WRF Users' Workshop, NCAR, Boulder, Colo. http://www.mmm.ucar.edu/wrf/users/workshops/WS2010/presentations/session%204/4-1%5fWRFworkshop2010Final.pdf, 2010), we tried to compare the performance of parameterizations and to enhance the forecast skill of low-level wind fields over the central western part of South Korea. Even though both subgrid-scale orography parameterizations significantly alleviated the positive bias at 10-m wind speed, the parameterization by Jimenez and Dudhia revealed a better forecast skill in wind speed under our modeling configuration. Implementation of the subgrid-scale orography parameterizations in the model did not affect the forecast skills in other meteorological fields including 10-m wind direction. Our study also brought up the problem of discrepancy in the definition of "10-m" wind between model physics parameterizations and observations, which can cause overestimated winds in model simulations. The overestimation was larger in stable conditions than in unstable conditions, indicating that the weak diurnal cycle in the model could be attributed to the representation error. [ABSTRACT FROM AUTHOR]

Local aerodynamic roughness parameters (zero-plane displacement, zd, and aerodynamic roughness length, z0) are determined for an urban park and a suburban neighbourhood with a new morphometric parameterisation that includes vegetation. Inter-seasonal analysis at the urban park demonstrates zd determined with two anemometric methods is responsive to vegetation state and is 1-4 m greater during leaf-on periods. The seasonal change and directional variability in the magnitude of zd is reproduced by the morphometric methods, which also indicate z0 can be more than halved during leaf-on periods. In the suburban neighbourhood during leaf-on, the anemometric and morphometric methods have similar directional variability for both zd and z0. Wind speeds at approximately 3 times the average roughness-element height are estimated most accurately when using a morphometric method which considers roughness-element height variability. Inclusion of vegetation in the morphometric parameterisation improves wind-speed estimation in all cases. Results indicate that the influence of both vegetation and roughness-element height variability are important for accurate determination of local aerodynamic parameters and the associated wind-speed estimates. [ABSTRACT FROM AUTHOR]

Aerosol effects on clouds and precipitation account for a large portion of uncertainties in the prediction of the future course of global hydrologic circulations and climate. As a process of a better understanding of interactions between aerosol, clouds and precipitation, simulations are performed for a mixed-phase convective multiple-cloud system over the tropics. Studies on single-cloud systems have shown that aerosol-induced increases in freezing, associated increases in parcel buoyancy and thus the intensity of clouds (or updrafts) are a main mechanism which controls aerosol-cloud-precipitation interactions in convective clouds. However, in the multiple-cloud system that plays much more important roles in global hydrologic circulations and thus climate than single-cloud systems, aerosol effects on condensation play the most important role in aerosol-induced changes in the intensity of clouds and the effects on freezing play a negligible role in those changes. Aerosol-induced enhancement in evaporation intensifies gust fronts and increases the number of subsequently developing clouds, which leads to the substantial increases in condensation and associated intensity of convection. Although aerosol-induced enhancement in freezing takes part in the increases in condensation by inducing stronger convergence around cloud bottom, the increases in condensation are ~one order of magnitude larger than those in freezing. It is found that while aerosol-induced increases in freezing create intermittent extremely heavy precipitation, aerosol-induced increases in evaporation enhance light and medium precipitation in the multiple-cloud system here. This increase in light and medium precipitation makes it possible that cumulative precipitation increases with increasing aerosol concentration, although the increase is small. It is interesting that the altitude of the maximum of the time- and domain-averaged hydrometeor mass densities is quite robust to increases in aerosol concentration. This is because locations of gust fronts and homogeneous freezing do not vary significantly with changing aerosol concentration and this outweighs aerosol effects on hydrometeor size. [ABSTRACT FROM AUTHOR]

On 21 September 2010, heavy rainfall with a local maximum of 259 mm d occurred near Seoul, South Korea. We examined the ability of the Weather Research and Forecasting (WRF) model in reproducing this disastrous rainfall event and identified the role of two physical processes: planetary boundary layer (PBL) and microphysics (MPS) processes. The WRF model was forced by 6-hourly National Centers for Environmental Prediction (NCEP) Final analysis (FNL) data for 36 hours form 1200 UTC 20 to 0000 UTC 22 September 2010. Twenty-five experiments were performed, consisting of five different PBL schemes-Yonsei University (YSU), Mellor-Yamada-Janjic (MYJ), Quasi Normal Scale Elimination (QNSE), Bougeault and Lacarrere (BouLac), and University of Washington (UW)-and five different MPS schemes-WRF Single-Moment 6-class (WSM6), Goddard, Thompson, Milbrandt 2-moments, and Morrison 2-moments. As expected, there was a specific combination of MPS and PBL schemes that showed good skill in forecasting the precipitation. However, there was no specific PBL or MPS scheme that outperformed the others in all aspects. The experiments with the UW PBL or Thompson MPS scheme showed a relatively small amount of precipitation. Analyses form the sensitivity experiments confirmed that the spatial distribution of the simulated precipitation was dominated by the PBL processes, whereas the MPS processes determined the amount of rainfall. It was also found that the temporal evolution of the precipitation was influenced more by the PBL processes than by the MPS processes. [ABSTRACT FROM AUTHOR]

This study presents an investigation of the spin-up behavior of soil moisture content (SMC) and evapotranspiration (ET) in an offline Noah land surface model (LSM) for East Asia, focusing on its interplay with the Asian monsoon. The set of 5-year recursive runs is conducted to properly assess the spin-up behavior of land surface processes and consists of simulations initialized with (1) a spatially uniform soil moisture, (2) NCEP GDAS soil moisture data, and (3) ECMWF ERA-Interim soil moisture data. Each run starts either after or before the summer monsoon. Initial SMCs from GDAS and ERA-Interim data significantly deviate from the equilibrium state (spin-up state) with the given input forcing even though the same equilibrium is reached within 3-year spin-up time, indicating that spin-up of land surface process is necessary. SMC reaches the equilibrium much quickly when (1) the consistent LSMs have been used in the prediction and analysis systems and (2) the spin-up simulation starts before the onset of heavy rainfall events during summer monsoon. For an area with heavy monsoon rainfall, the total column SMC and ET spin up quickly. The spin-up time over dry land is about 2-3 years, but for monsoon rainfall area decreases dramatically to about 3 months if the spin-up run starts just before the onset of monsoon. Further scrutiny shows that the spin-up time is well correlated with evaporative fraction given by the ratio between the latent heat flux and the available energy at the land surface. [ABSTRACT FROM AUTHOR]

Land Data Assimilation Systems have been developed to generate the surface initial conditions such as soil moisture and temperature for better prediction of weather and climate. We have constructed Korea Land Data Assimilation System (KLDAS) based on an uncoupled land surface modeling framework that integrates high-resolution in-situ observation, satellite data, land surface information from the WRF Preprocessing System (WPS) and the MODIS land products over the East Asia. To present better surface conditions, the KLDAS is driven by atmospheric forcing data from the in-situ rainfall gauges and satellite. In this study, we 1) briefly introduce the KLDAS, 2) evaluate the meteorological states near the surface and the surface fluxes reproduced by the KLDAS against the in-situ observation, and then 3) examine the performance of the mesoscale model initialized by the KLDAS. We have generated a 5-year, 10 km, hourly atmospheric forcing dataset for use in KLDAS operating across East Asia. The KLDAS has effectively reproduced the observed patterns of soil moisture, soil temperature, and surface fluxes. Further scrutiny reveals that the numerical simulations incorporating the KLDAS outputs show better agreement in both the simulated near-surface conditions and rainfall distribution over the Korean Peninsula, compared to those without the KLDAS. [ABSTRACT FROM AUTHOR]

Canopy turbulence plays an important role in mass and energy exchanges at the canopy-atmosphere interface. Despite extensive studies on canopy turbulence over a flat terrain, less attention has been given to canopy turbulence in a complex terrain. The purpose of this study is to scrutinize characteristics of canopy turbulence in roughness sublayer over a hilly forest terrain. We investigated basic turbulence statistics, conditionally sampled statistics, and turbulence spectrum in terms of different atmospheric stabilities, wind direction and vertical structures of momentum fluxes. Similarly to canopy turbulence over a homogeneous terrain, turbulence statistics showed coherent structure. Both quadrant and spectrum analysis corroborated the role of intermittent and energetic eddies with length scale of the order of canopy height, regardless of wind direction except for shift of peak in vertical wind spectrum to relatively high frequency in the down-valley wind. However, the magnitude of the momentum correlation coefficient in a neutral condition was smaller than typical value over a flat terrain. Further scrutiny manifested that, in the up-valley flow, temperature skewness was larger and the contribution of ejection to both momentum and heat fluxes was larger compared to the downvalley flow, indicating that thermal instability and weaker wind shear in up-valley flow asymmetrically affect turbulent transport within the canopy. [ABSTRACT FROM AUTHOR]

n expanded planar-fit (PF) method over complex terrain is presented and applied to coordinate rotation of the eddy-covariance (EC) flux and vertical velocity estimation. Theoretical analysis indicates that PF coefficients depend on wind direction, and an expression of vertical velocity is deduced. We applied the theory using 1 year of observations from the KoFlux site in the Gwangneung Forest in Korea and investigated the influence of wind direction on the PF method. Then, we performed an expansion of the PF method to consider dependence of PF coefficients on wind direction and applied the PF method to every sector. The results show that the PF coefficients and tilt angles over complex terrain vary with wind direction. Two hundred 30-min data sets are sufficient to derive stable PF coefficients over hilly terrain for each sector. The relative difference in eddy-covariance flux between the general planar fit (GPF) and sector planar fit (SPF) is less than 10% for the scalar flux and about 18% for friction velocity. Vertical velocity and vertical advection (VA) terms were also calculated and compared using SPF and GPF methods, and a normal distribution and diurnal trend of real vertical velocity on clear days are presented. [ABSTRACT FROM AUTHOR]

In boundary layer meteorology, surface layer similarity theory plays a critical role in measuring and modeling biospheric fluxes. In stable boundary layer, surface layer similarity called z-less stratification has been one of main research topics for over than two decades and the issue has yet to be settled in micrometeorology. In this scientific discussion on z-less turbulence, different turbulence statistics were used inconsistently and it was argued that z-less turbulence was valid if only any turbulence statistics were constant with different atmospheric stabilities. Consequently, such inconsistently tested turbulence statistics and misconception on z-less turbulence hinder us from correctly understanding turbulence structure in the stable boundary layer. This note revisits z-less turbulence and emphasizes that different dimensionless turbulence statistics generally do not exhibit a common behavior in the limit of z-less stratification. [ABSTRACT FROM AUTHOR]

A simple aerodynamic-variance method is proposed to fill gaps in continuous CO2 flux measurements in rainy conditions, when open-path analysers do not function. The method requires turbulent conditions (friction velocity greater than 0.1ms−1), and uses measurements of mean wind speed, and standard deviations of temperature and CO2 concentration fluctuations to complement, and at times replace, eddy-covariance measurements of friction velocity, sensible heat flux and CO2 flux. Friction velocity is estimated from the mean wind speed with a flux-gradient relationship modified for the roughness sublayer. Since normalised standard deviations do not follow Monin-Obukhov similarity theory in the roughness sublayer, a simple classification scheme according to the scalar turbulence scale was used. This scheme is shown to produce sensible heat and CO2 flux estimates that are well correlated with the measured values. [ABSTRACT FROM AUTHOR]

Background: It is important to quantify changes in CO2 sources and sinks with land use and land cover change. In the last several decades, carbon sources and sinks in East Asia have been altered by intensive land cover changes due to rapid economic growth and related urbanization. To understand impact of urbanization on carbon cycle in the monsoon Asia, we analyze net CO2 exchanges for various land cover types across an urbanization gradient in Korea covering high-rise high-density residential, suburban, cropland, and subtropical forest areas. Results: Our analysis demonstrates that the urban residential and suburban areas are constant CO2 sources throughout the year (2.75 and 1.02 kg C m−2 year−1 at the urban and suburban sites), and the net CO2 emission indicate impacts of urban vegetation that responds to the seasonal progression of the monsoon. However, the total random uncertainties of measurement are much larger in the urban and suburban areas than at the nonurban sites, which can make it challenging to obtain accurate urban flux measurements. The cropland and forest sites are strong carbon sinks because of a double-cropping system and favorable climate conditions during the study period, respectively (− 0.73 and − 0.60 kg C m−2 year−1 at the cropland and forest sites, respectively). The urban area of high population density (15,000 persons km−2) shows a relatively weak CO2 emission rate per capita (0.7 t CO2 year−1 person−1), especially in winter because of a district heating system and smaller traffic volume. The suburban area shows larger net CO2 emissions per capita (4.9 t CO2 year−1 person−1) because of a high traffic volume, despite a smaller building fraction and population density (770 persons km−2). Conclusions: We show that in situ flux observation is challenging because of its larger random uncertainty and this larger uncertainty should be carefully considered in urban studies. Our findings indicate the important role of urban vegetation in the carbon balance and its interaction with the monsoon activity in East Asia. Urban planning in the monsoon Asia must consider interaction on change in the monsoon activity and urban structure and function for sustainable city in a changing climate. [ABSTRACT FROM AUTHOR]