Effect of wave amplitude on turbulent flow in a wavy channel by direct numerical simulation

The present study numerically investigates the characteristics of three-dimensional turbulent flow in a wavy channel. For the purpose of a careful observation of the effect of the wave amplitude on the turbulent flow, numerical simulations are performed at a various range of the wave amplitude to wavelength ratio (0.01?α/λ?0.05), where the wavelength is fixed with the same value of the mean channel height (H). The immersed boundary method is used to handle the wavy surface in a rectangular grid system, using the finite volume method. The Reynolds number (Re=U_bH/ν) based on the bulk velocity (U_b) is fixed at 6760. The present computational results for a wavy surface are well compared with those of references. When α/λ=0.02, the small recirculating flow occurs near the trough at the instant, but the mean reverse flow is not observed. In the mean flow field, the reverse flow appears from α/λ=0.03 among the wave amplitude considered in this study. The domain of the mean reverse flow defined by the locations of separation and reattachment depends strongly on the wave amplitude. The pressure drag coefficient augments with increasing the wave amplitude. The friction drag coefficient shows the increase and decrease behavior according to the wave amplitude. The quantitative information about the flow variables such as the distribution of pressure and shear stress on the wavy surface is highlighted.     

Accumulation features of arsenic species in various fishes collected from coastal cities in Korea

In this study, 36 fish species were collected from three coastal cities in Korea to investigate levels and patterns of six arsenicals (arsenite: As (III), arsenate: As (V), arsenocholine: AsC, arsenobetaine: AsB, monomethylarsonic acid: MMA, and dimethylarsinic acid: DMA). The levels of ∑₆ As in the different fish species varied substantially, ranging from 0.02 μg As/g ww (Islaeli carp) to 9.65 μg As/g ww (Skate ray) with a median of 0.40 μg As/g ww. All the arsenicals in marine fishes showed higher levels than those in freshwater fishes due to fish feed living in saline water. Overall, marine carnivorous fishes seem to be more contaminated with arsenic. For all the fish samples, AsB (mean fraction: 90.6%) was dominant among the six arsenicals, indicating biomethylation of inorganic arsenic and accumulation of AsB. Fish species with high water contents showed elevated levels of As (III), but there was no further significant correlations between arsenicals and water/lipid contents. Concentrations of As (V) were significantly lower than those of As (III), which implies that As (V) is reduced during biomethylation of inorganic arsenic. Consequently, we hypothesize that the toxicity of arsenic (mainly derived from As (III)) can be increased by the reduction of As (V), especially for the fish species with higher water contents.     

Classification of precipitation types using fall velocity-diameter relationships from 2D-video distrometer measurements

Fall velocity-diameter relationships for four different snowflake types (dendrite, plate, needle, and graupel) were investigated in northeastern South Korea, and a new algorithm for classifying hydrometeors is proposed for distrometric measurements based on the new relationships. Falling ice crystals (approximately 40 000 particles) were measured with a two-dimensional video disdrometer (2DVD) during a winter experiment from 15 January to 9 April 2010. The fall velocity-diameter relationships were derived for the four types of snowflakes based on manual classification by experts using snow photos and 2DVD measurements: the coefficients (exponents) for different snowflake types were 0.82 (0.24) for dendrite, 0.74 (0.35) for plate, 1.03 (0.71) for needle, and 1.30 (0.94) for graupel, respectively. These new relationships established in the present study (PS) were compared with those from two previous studies. Hydrometeor types were classified with the derived fall velocity-diameter relationships, and the classification algorithm was evaluated using 3× 3 contingency tables for one rain-snow transition event and three snowfall events. The algorithm showed good performance for the transition event: the critical success indices (CSIs) were 0.89, 0.61 and 0.71 for snow, wet-snow and rain, respectively. For snow events, the algorithm performance for dendrite and plate (CSIs = 1.0 and 1.0, respectively) was better than for needle and graupel (CSIs = 0.67 and 0.50, respectively).     

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.     

Effect of snow cover on pan-Arctic permafrost thermal regimes

This study quantitatively evaluated how insulation by snow depth (SND) affected the soil thermal regime and permafrost degradation in the pan-Arctic area, and more generally defined the characteristics of soil temperature (T_SOIL) and SND from 1901 to 2009. This was achieved through experiments performed with the land surface model CHANGE to assess sensitivity to winter precipitation as well as air temperature. Simulated T_SOIL, active layer thickness (ALT), SND, and snow density were generally comparable with in situ or satellite observations at large scales and over long periods. Northernmost regions had snow that remained relatively stable and in a thicker state during the past four decades, generating greater increases in T_SOIL. Changes in snow cover have led to changes in the thermal state of the underlying soil, which is strongly dependent on both the magnitude and the timing of changes in snowfall. Simulations of the period 2001–2009 revealed significant differences in the extent of near-surface permafrost, reflecting differences in the model’s treatment of meteorology and the soil bottom boundary. Permafrost loss was greater when SND increased in autumn rather than in winter, due to insulation of the soil resulting from early cooling. Simulations revealed that T_SOIL tended to increase over most of the pan-Arctic from 1901 to 2009, and that this increase was significant in northern regions, especially in northeastern Siberia where SND is responsible for 50 % or more of the changes in T_SOIL at a depth of 3.6 m. In the same region, ALT also increased at a rate of approximately 2.3 cm per decade. The most sensitive response of ALT to changes in SND appeared in the southern boundary regions of permafrost, in contrast to permafrost temperatures within the 60°N–80°N region, which were more sensitive to changes in snow cover. Finally, our model suggests that snow cover contributes to the warming of permafrost in northern regions and could play a more important role under conditions of future Arctic warming.     

Radio Emission from Masuda Sources

We note that different models, providing comparably good interpretation of the hard X-ray properties of so-called Masuda sources, can make distinctly different predictions for the radio emission produced at the Masuda source by the same population of accelerated electrons. Accordingly, we calculate the radio emission within a few competing models, i.e., those involving magnetic, turbulent, and collisional trapping of the fast electrons in the coronal source. We show that even available incomplete radio observations of the classical Masuda event and a Masuda-like event on 31 December 2007, recently reported by Krucker et al. (Astrophys. J. 714, 1108, 2010) are highly valuable in restricting the physical model of the source. Furthermore, our study proposes that combination of more complete high-resolution X-ray and radio observations can allow unambiguous distinction between the competing Masuda source models.     

Effects of prescribed initial cloud droplet spectra on convective cloud and precipitation developments under different thermodynamic conditions: A modeling and observational study

A two-dimensional cloud model with bin microphysics was used to investigate the effects of cloud condensation nuclei (CCN) concentrations and thermodynamic conditions on convective cloud and precipitation developments. Two different initial cloud droplet spectra were prescribed based on the total CCN concentrations of maritime (300 cm^− 3) and continental (1000 cm^− 3) air masses, and the model was run on eight thermodynamic conditions obtained from observational soundings. Six-hourly sounding data and 1-hourly precipitation data from two nearby weather stations in Korea were analyzed for the year 2002 to provide some observational support for the model results.For one small Convective Available Potential Energy (CAPE) ( 300 J kg^− 1) sounding, the maritime and continental differences were incomparably large. The crucial difference was the production of ice phase hydrometeors in the maritime cloud and only water drops in the continental cloud. Ice phase hydrometeors and intrinsically large cloud drops of the maritime cloud eventually lead to significant precipitation. Meanwhile negligible precipitation developed from the continental cloud. For the three other small CAPE soundings, generally weak convective clouds developed but the maritime and continental clouds were of the same phases (both warm or both cold) and their differences were relatively small.Model runs with the four large CAPE ( 3000 J kg^− 1) soundings demonstrated that the depth between the freezing level (FL) and the lifting condensation level (LCL) was crucial to determine whether a cloud becomes a cold cloud or not, which in turn was found to be a crucial factor to enhance cloud invigoration with the additional supply of freezing latent heat. For two large CAPE soundings, FL–LCL was so deep that penetration of FL was prohibitive, and precipitation was only mild in the maritime clouds and negligible in the continental clouds. Two other soundings of similarly large CAPE had small FL–LCL, and both the maritime and continental clouds became cold clouds. Precipitation was strong for both but much more so in the maritime clouds, while the maximum updraft velocity and the cloud top were slightly higher in continental clouds. Although limited to small CAPE cases, more precipitation for smaller FL–LCL for a selected group of precipitation and thermodynamic sounding data from Korea was in support of these model results in its tendency.These results clearly demonstrated that the CCN effects on cloud and precipitation developments critically depended on the given thermodynamic conditions and not just the CAPE but the entire structure of the thermodynamic profiles had to be taken into account.     

Aerobic oxidation of mackinawite (FeS) and its environmental implication for arsenic mobilization

Oxidation of mackinawite (FeS) and concurrent mobilization of arsenic were investigated as a function of pH under oxidizing conditions. At acidic pH, FeS oxidation is mainly initiated by the proton-promoted dissolution, which results in the release of Fe(II) and sulfide in the solution. While most of dissolved sulfide is volatilized before being oxidized, dissolved Fe(II) is oxidized into green rust-like precipitates and goethite (α-FeOOH). At basic pH, the development of Fe(III) (oxyhydr)oxide coating on the FeS surface inhibits the solution-phase oxidation following FeS dissolution. Instead, FeS is mostly oxidized into lepidocrocite (γ-FeOOH) via the surface-mediated oxidation without dissolution. At neutral pH, FeS is oxidized via both the solution-phase oxidation following FeS dissolution and the surface-mediated oxidation mechanisms. The mobilization of arsenic during FeS oxidation is strongly affected by FeS oxidation mechanisms. At acidic pH (and to some extent at neutral pH), the rapid FeS dissolution and the slow precipitation of Fe (oxyhydr)oxides results in arsenic accumulation in water. In contrast, the surface-mediated oxidation of FeS at basic pH leads to the direct formation of Fe (oxyhydr)oxides, which provides effective adsorbents for As under oxic conditions. At acidic and neutral pH, the solution-phase oxidation of dissolved Fe(II) accelerates the oxidation of the less adsorbing As(III) to the more adsorbing As(V). This study reveals that the oxidative mobilization of As may be a significant pathway for arsenic enrichment of porewaters in sulfidic sediments.     

Geostatistical integration using 2-D electrical resistivity and 3-D gravity methods for detecting cavities in a karst area

Two-dimensional (2-D) electrical resistivity and three-dimensional (3-D) gravity explorations were undertaken to estimate the 3-D distribution of karst cavities at the area of Yongweol-ri in Korea where ground subsidence has occurred. Although the gravity method is a low-cost way of analyzing a 3-D continuous structure, its vertical resolving power is poor. In contrast, the electrical resistivity method can provide a 2-D and/or 3-D subsurface structure with a much higher lateral and vertical resolution than the gravity method. Accordingly, geostatistical methods and density information were used to enhance the 2-D resistivity structure revealed by the electrical resistivity method into a 3-D structure. The assumptions are, first, that each data set senses the same underlying geological structure in terms of different material properties and, secondly, that two different material properties are correlated locally or globally throughout the entire target area. As a result, the distribution of limestone cavities can be estimated under the assumption that they are mostly filled with groundwater and clayey soils and have abnormally low levels of resistivity and density. The estimated distribution corresponds with the grouting, borehole imaging and monitoring data. In this example, it can be seen that the integration analysis of 2-D electrical resistivity and 3-D gravity methods is a very powerful tool for 3-D subsurface imaging and that the method can provide enhanced imaging capabilities for 3-D cavities.