Projections of high resolution climate changes for South Korea using multiple-regional climate models based on four RCP scenarios. Part 1: surface air temperature

We projected surface air temperature changes over South Korea during the mid (2026-2050) and late (2076-2100) 21st century against the current climate (1981-2005) using the simulation results from five regional climate models (RCMs) driven by Hadley Centre Global Environmental Model, version 2, coupled with the Atmosphere- Ocean (HadGEM2-AO), and two ensemble methods (equal weighted averaging, weighted averaging based on Taylor’s skill score) under four Representative Concentration Pathways (RCP) scenarios. In general, the five RCM ensembles captured the spatial and seasonal variations, and probability distribution of temperature over South Korea reasonably compared to observation. They particularly showed a good performance in simulating annual temperature range compared to HadGEM2-AO. In future simulation, the temperature over South Korea will increase significantly for all scenarios and seasons. Stronger warming trends are projected in the late 21st century than in the mid-21st century, in particular under RCP8.5. The five RCM ensembles projected that temperature changes for the mid/late 21st century relative to the current climate are +1.54°C/+1.92°C for RCP2.6, +1.68°C/+2.91°C for RCP4.5, +1.17°C/+3.11°C for RCP6.0, and +1.75°C/+4.73°C for RCP8.5. Compared to the temperature projection of HadGEM2-AO, the five RCM ensembles projected smaller increases in temperature for all RCP scenarios and seasons. The inter-RCM spread is proportional to the simulation period (i.e., larger in the late-21st than mid-21st century) and significantly greater (about four times) in winter than summer for all RCP scenarios. Therefore, the modeled predictions of temperature increases during the late 21st century, particularly for winter temperatures, should be used with caution.     

Direct and semi-direct radiative effects of anthropogenic aerosols in the Western United States: Seasonal and geographical variations according to regional climate characteristics

The direct and semi-direct radiative effects of anthropogenic aerosols on the radiative transfer and cloud fields in the Western United States (WUS) according to seasonal aerosol optical depth (AOD) and regional climate are examined using a regional climate model (RCM) in conjunction with the aerosol fields from a GEOS-Chem chemical-transport model (CTM) simulation. The two radiative effects cannot be separated within the experimental design in this study, thus the combined direct- and semi-direct effects are called radiative effects hereafter. The CTM shows that the AOD associated with the anthropogenic aerosols is chiefly due to sulfates with minor contributions from black carbon (BC) and that the AOD of the anthropogenic aerosol varies according to local emissions and the seasonal low-level winds. The RCM-simulated anthropogenic aerosol radiative effects vary according to the characteristics of regional climate, in addition to the AOD. The effects on the top of the atmosphere (TOA) outgoing shortwave radiation (OSRT) range from ?0.2?Wm^?2 to ?1?Wm^?2. In Northwestern US (NWUS), the maximum and minimum impact of anthropogenic aerosols on OSRT occurs in summer and winter, respectively, following the seasonal AOD. In Arizona-New Mexico (AZNM), the effect of anthropogenic sulfates on OSRT shows a bimodal distribution with winter/summer minima and spring/fall maxima, while the effect of anthropogenic BC shows a single peak in summer. The anthropogenic aerosols affect surface insolation range from ?0.6?Wm^?2 to ?2.4?Wm^?2, with similar variations found for the effects on OSRT except that the radiative effects of anthropogenic BC over AZNM show a bimodal distribution with spring/fall maxima and summer/winter minima. The radiative effects of anthropogenic sulfates on TOA outgoing longwave radiation (OLR) and the surface downward longwave radiation (DLRS) are notable only in summer and are characterized by strong geographical contrasts; the summer OLR in NWUS (AZNM) is reduced (enhanced) by 0.52?Wm^?2 (1.14?Wm^?2). The anthropogenic sulfates enhance (reduce) summer DLRS by 0.2?Wm^?2 (0.65?Wm^?2) in NWUS (AZNM). The anthropogenic BC affect DLRS noticeably only in AZNM during summer. The anthropogenic aerosols affect the cloud water path (CWP) and the radiative transfer noticeably only in summer when convective clouds are dominant. Primarily shortwave-reflecting anthropogenic sulfates decrease and increase CWP in AZNM and NWUS, respectively, however, the shortwave-absorbing anthropogenic BC reduces CWP in both regions. Due to strong feedback via convective clouds, the radiative effects of anthropogenic aerosols on the summer radiation field are more closely correlated with the changes in CWP than the AOD. The radiative effect of the total anthropogenic aerosols is dominated by the anthropogenic sulfates that contribute more than 80% of the total AOD associated with the anthropogenic aerosols.     

Modelling Korean extreme rainfall using a Kappa distribution and maximum likelihood estimate

Summary Attempts to use the 4-parameter Kappa distribution (K4D) with the maximum likelihood estimates (MLE) on the summer extreme daily rainfall data at 61 gauging stations over South Korea have been made to obtain reliable quantile estimates for several return periods. A numerical algorithm for searching MLE of K4D by minimizing the negative log-likelihood function with penalty method has been described. The isopluvial maps of estimated design values corresponding to selected return periods have been presented. The highest return values are centered at sites in the south-western part of the Korean peninsula. The distribution of return values for annual maxima of 2-day precipitation (AMP2) is more similar to the climatological features of annual total precipitation of Korea than that of annual maxima of daily precipitation (AMP1). Our results of return values delineate well the horizontal patterns of the heavy precipitation over the Korean peninsula. Received January 15, 2001 Revised October 8, 2001     

A dissection of the surface temperature biases in the Community Earth System Model

Based upon the climate feedback-responses analysis method, a quantitative attribution analysis is conducted for the annual-mean surface temperature biases in the Community Earth System Model version 1 (CESM1). Surface temperature biases are decomposed into partial temperature biases associated with model biases in albedo, water vapor, cloud, sensible/latent heat flux, surface dynamics, and atmospheric dynamics. A globally-averaged cold bias of ?1.22 K in CESM1 is largely attributable to albedo bias that accounts for approximately ?0.80 K. Over land, albedo bias contributes ?1.20 K to the averaged cold bias of ?1.45 K. The cold bias over ocean, on the other hand, results from multiple factors including albedo, cloud, oceanic dynamics, and atmospheric dynamics. Bias in the model representation of oceanic dynamics is the primary cause of cold (warm) biases in the Northern (Southern) Hemisphere oceans while surface latent heat flux over oceans always acts to compensate for the overall temperature biases. Albedo bias resulted from the model’s simulation of snow cover and sea ice is the main contributor to temperature biases over high-latitude lands and the Arctic and Antarctic region. Longwave effect of water vapor is responsible for an overall warm (cold) bias in the subtropics (tropics) due to an overestimate (underestimate) of specific humidity in the region. Cloud forcing of temperature biases exhibits large regional variations and the model bias in the simulated ocean mixed layer depth is a key contributor to the partial sea surface temperature biases associated with oceanic dynamics. On a global scale, biases in the model representation of radiative processes account more for surface temperature biases compared to non-radiative, dynamical processes.     

Retrieval of outgoing longwave radiation from COMS narrowband infrared imagery

Hourly outgoing longwave radiation(OLR) from the geostationary satellite Communication Oceanography Meteorological Satellite(COMS) has been retrieved since June 2010. The COMS OLR retrieval algorithms are based on regression analyses of radiative transfer simulations for spectral functions of COMS infrared channels. This study documents the accuracies of OLRs for future climate applications by making an intercomparison of four OLRs from one single-channel algorithm(OLR12.0using the 12.0 μm channel) and three multiple-channel algorithms(OLR10.8+12.0using the 10.8 and 12.0 μm channels; OLR6.7+10.8using the 6.7 and 10.8 μm channels; and OLR All using the 6.7, 10.8, and 12.0 μm channels). The COMS OLRs from these algorithms were validated with direct measurements of OLR from a broadband radiometer of the Clouds and Earth's Radiant Energy System(CERES) over the full COMS field of view [roughly(50°S–50°N, 70°–170°E)] during April 2011.Validation results show that the root-mean-square errors of COMS OLRs are 5–7 W m-2, which indicates good agreement with CERES OLR over the vast domain. OLR6.7+10.8and OLR All have much smaller errors(～ 6 W m-2) than OLR12.0and OLR10.8+12.0(～ 8 W m-2). Moreover, the small errors of OLR6.7+10.8and OLR All are systematic and can be readily reduced through additional mean bias correction and/or radiance calibration. These results indicate a noteworthy role of the6.7 μm water vapor absorption channel in improving the accuracy of the OLRs. The dependence of the accuracy of COMS OLRs on various surface, atmospheric, and observational conditions is also discussed.     

Winter precipitation variability over Korean Peninsula associated with ENSO

In this study, winter precipitation variability associated with the El Niño-Southern Oscillation (ENSO) over the Korean Peninsula was investigated using a 5-pentad running mean data because significant correlation pattern cannot be revealed using seasonal-mean data. It was found a considerably significant positive correlation between Niño3 sea-surface temperature and precipitation during early winter (from Mid-November to early-December), when the correlation coefficient is close to 0.8 in early-December; the correlation is distinctively weakened during late winter. It is demonstrated that such sudden intraseasonal change in relation to ENSO is associated with the presence of anticyclonic flow over the Kuroshio extension region (Kuroshio anticyclone). In early winter, there is strong southerly wind over the Korean Peninsula, which is induced by the Philippine Sea anticyclone and Kuroshio anticyclone. However, in January, although the Philippine Sea anticyclone develops further, the Kuroshio anticyclone suddenly disappears; as a result, the impact of ENSO is considerably weakened over the Korean Peninsula. These results indicate that the Kuroshio anticyclone during El Niño peak phase plays a critical role by strongly affecting Northeast Asia climate, including the Korean Peninsula. In addition, it is also found that there are distinctive interdecadal changes of the relationship between ENSO and precipitation over the Korean Peninsula. In particular, the strong correlation in early winter is clearer in the recent 30 years than that in the previous period of 1950–1979.     

Impact of climate change on the persistent turbidity issue of a large dam reservoir in the temperate monsoon region

Long-term discharge of turbid water from reservoirs after flood events is a major socioenvironmental problem in many countries, including Korea. This study used a suite of mathematical models to simulate the fate of turbidity flows in the Soyanggang Reservoir in Korea, an important source of drinking water for the Seoul Capital Area, in response to extreme floods based on the Representative Concentration Pathway 4.5 climate scenario. It evaluated the effectiveness of the selective withdrawal facility (SWF), installed recently in the Soyanggang Reservoir to control persistent turbidity. Extreme floods with a maximum daily inflow rate greater than the historical maximum observed in 2006 were projected to occur four times in this century. The fate and transport of turbidity flows were highly influenced by both the thermal stability of the reservoir and the season in which the flood event occurred. Thus, SWF operations should consider the timing of extreme events (i.e., the imminence of the autumn turnover) to mitigate the impact of high turbidity on the water supply and downstream ecosystem. It was found to be ineffective under extreme events if these occurred in two consecutive years. Current reservoir operations, which rely heavily on the SWF, are likely to be inadequate to overcome the negative effects of extreme-turbidity events on reliably providing safe water supplies. Coping with the worst event expected to occur in the future would require additional countermeasures such as bypassing high-turbidity water.     

The Korean Integrated Model (KIM) System for Global Weather Forecasting

The Korea Institute of Atmospheric Prediction Systems (KIAPS) began a national project to develop a new global atmospheric model system in 2011. The ultimate goal of this 9-year project is to replace the current operational model at the Korea Meteorological Administration (KMA), which was adopted from the United Kingdom’s Meteorological Office’s unified model (UM) in 2010. The 12-km Korean Integrated Model (KIM) system, consisting of a spectral-element non-hydrostatic dynamical core on a cubed sphere grid and a state-of-the-art physics parameterization package, has been launched in a real-time forecast framework, with initial conditions obtained via the advanced hybrid four-dimensional ensemble variational data assimilation (4DEnVar) over its native grid. A development strategy for KIM and the evolution of its performance in medium-range forecasts toward a world-class global forecast system are described. Outstanding issues in KIM 3.1 as of February 2018 are discussed, along with a future plan for operational deployment in 2020.     

Adsorption of Rb and Sr at the orthoclase (0 0 1)-solution interface

Adsorption of Rb⁺ and Sr²⁺ at the orthoclase (0 0 1)-solution interface is probed with high-resolution X-ray reflectivity and resonant anomalous X-ray reflectivity. Specular X-ray reflectivity data for orthoclase in contact with 0.01 m RbCl solution at pH 5.5 reveal a systematic increase in electron density adjacent to the mineral surface with respect to that observed in contact with de-ionized water (DIW). Quantitative analysis indicates that Rb⁺ adsorbs at a height of 0.83 ± 0.03 Å with respect to the bulk K⁺ site with a nominal coverage of 0.72 ± 0.10 ions per surface unit mesh (55.7 Å²). These results are consistent with an ion-exchange reaction in which Rb⁺ occupies an inner-sphere adsorption (IS) site. In contrast, X-ray reflectivity data for orthoclase in contact with 0.01 m Sr(NO₃)₂ solution at pH 5.3 reveal few significant changes with respect to DIW. Resonant anomalous X-ray reflectivity was used to probe Sr²⁺ adsorption and to image its vertical distribution. This element-specific measurement reveals that Sr²⁺ adsorbs with a total coverage of 0.37 ± 0.02 ions per surface unit mesh, at a substantially larger height (3.28 ± 0.05 Å) than found for Rb⁺, and with a relatively broad density distribution (having a root-mean-square width of 1.88 ± 0.08 Å for a single-peak model), implying that Sr²⁺ adsorbs primarily as a fully-hydrated outer-sphere (OS), species. Comparison to a two-height model suggests that 13 ± 5% of the adsorbed Sr²⁺ may be present as an IS species. This partitioning implies a ∼5 kJ/mol difference in free energy between the IS and OS Sr²⁺ on orthoclase. Differences in the partitioning of Sr²⁺ between IS and OS species for orthoclase (0 0 1) and muscovite (0 0 1) suggest control by the geometry of the IS adsorption site. Results for the OS distribution are compared to predictions of the Poisson-Boltzmann equation in the strong coupling regime, which predicts an intrinsically narrow vertical diffuse ion distribution; the OS distribution might thus be thought of as the diffuse ion profile in the limit of high surface charge.     

Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions

Precise and efficient numerical simulation of transport processes in subsurface systems is a prerequisite for many site investigation or remediation studies. Random walk particle tracking (RWPT) methods have been introduced in the past to overcome numerical difficulties when simulating propagation processes in porous media such as advection-dominated mass transport. Crucial for the precision of RWPT methods is the accuracy of the numerically calculated ground water velocity field. In this paper, a global node-based method for velocity calculation is used, which was originally proposed by Yeh (Water Resour Res 7:1216–1225, 1981). This method is improved in three ways: (1) extension to unstructured grids, (2) significant enhancement of computational efficiency, and (3) extension to saturated (groundwater) as well as unsaturated systems (soil water). The novel RWPT method is tested with numerical benchmark examples from the literature and used in two field scale applications of contaminant transport in saturated and unsaturated ground water. To evaluate advective transport of the model, the accuracy of the velocity field is demonstrated by comparing several published results of particle pathlines or streamlines. Given the chosen test problem, the global node-based velocity estimation is found to be as accurate as the CK method (Cordes and Kinzelbach in Water Resour Res 28(11):2903–2911, 1992) but less accurate than the mixed or mixed-hybrid finite element methods for flow in highly heterogeneous media. To evaluate advective–diffusive transport, a transport problem studied by Hassan and Mohamed (J Hydrol 275(3–4):242–260, 2003) is investigated here and evaluated using different numbers of particles. The results indicate that the number of particles required for the given problem is decreased using the proposed method by about two orders of magnitude without losing accuracy of the concentration contours as compared to the published numbers.