Characteristics of tephra in Holocene lake sediments on King George Island, West Antarctica: implications for deglaciation and paleoenvironment

Several reworked tephra layers in gravity-flow deposits are present in lacustrine core sediments collected from Hotel and Rudy Lakes on King George Island, South Shetland Islands, maritime sub-Antarctica. This study tests the values of tephra for establishing regional tephrochronologies for lakes in ice-covered landscapes in the vicinity of volcanoes. The tephra record is more abundant in a long Hotel Lake core (515 cm long). This study uses volcanic glass samples from five tephra layers of Hotel Lake and from one tephra layer of Rudy Lake. Morphologically, tephras are mixtures of basaltic and pumice shards, both having various degrees of vesicularity. Major element analyses of glass shards reveal that the majority of the glass fragments belong to basic glass (<60 wt% SiO₂), compositionally ranging from basalt to low-silica andesite and subalkaline series medium-K tholeiites, probably sourced from Deception Island located 130 km southwest of King George Island. Less than 20% of tephra belongs to silicic glass and occurs in three tephra horizons E of Hotel Lake. However, source volcano(es) for about 10% of basic tephra and silicic tephra are not readily identified from nearby volcanic centers. Except for the studied tephra in Rudy Lake, all tephra samples in Hotel Lake are not ashfall deposits but reworked and redeposited pyroclasts derived from retreating ice sheet, resulting in the occurrence of geochemically equivalent tephra samples in different tephra horizons. The dating of the studied tephra horizons represents the timing of deglaciation rather than that of volcanic eruptions. The result of this study implies that combined with sedimentological information more chemical criterion is necessary to study tephrochronology and regional correlation and to understand paleoenvironmental changes using tephra.     

Discrimination of small earthquakes and artificial explosions in the Korean Peninsula using Pg/Lg ratios

The purpose of this study is to develop a technique to discriminate artificial explosions from local small earthquakes ( M ≤ 4.0) in the time–frequency domain. In order to obtain spectral features of artificial explosions and earthquakes, 3-D spectrograms (frequency, time and amplitude) have been used. They represent a useful tool for studying the frequency content of entire seismic waveforms observed at local and regional distances (Kim, Simpson & Richards 1994). P and S(L g ) waves from quarry blasts show that the frequency content associated with the dominant amplitude appears above 10  Hz and Rg phases are observed at close distances. P and S(L g ) waves from the Tongosan earthquake have strong amplitudes below 10  Hz. For the Munkyong earthquake, however, a broader frequency content up to 20  Hz is found.
  For discrimination between small earthquakes and explosions, Pg/L g spectral ratios are used below 10  Hz, and through spectrogram analysis we can see different frequency contents of explosions and earthquakes. Unfortunately, because explosion data recorded at KSRS array are digitized at 20  sps, we cannot avoid analysing below 10  Hz because of the Nyquist frequency. In order to select time windows, the group velocity was computed using multiple-filter analysis (MFA), and free-surface effects have been removed from all three-component data in order to improve data quality. Using FFT, a log-average spectral amplitude is calculated over seven frequency bands: 0.5 to 3, 2 to 4, 3 to 5, 4 to 6, 5 to 7, 6 to 8 and 8 to 10  Hz. The best separation between explosions and earthquakes is observed from 6 to 8  Hz. In this frequency band we can separate explosions with log ( Pg/L g ) above −0.5, except EXP1 recorded at SIHY1-1, and earthquakes below −0.5, except the Munkyong earthquake record at station KMH.     

Effect of spatial variability of cross‐correlated soil properties on bearing capacity of strip footing

Geotechnical engineering problems are characterized by many sources of uncertainty. Some of these sources are connected to the uncertainties of soil properties involved in the analysis. In this paper, a numerical procedure for a probabilistic analysis that considers the spatial variability of cross‐correlated soil properties is presented and applied to study the bearing capacity of spatially random soil with different autocorrelation distances in the vertical and horizontal directions. The approach integrates a commercial finite difference method and random field theory into the framework of a probabilistic analysis. Two‐dimensional cross‐correlated non‐Gaussian random fields are generated based on a Karhunen–Loève expansion in a manner consistent with a specified marginal distribution function, an autocorrelation function, and cross‐correlation coefficients. A Monte Carlo simulation is then used to determine the statistical response based on the random fields. A series of analyses was performed to study the effects of uncertainty due to the spatial heterogeneity on the bearing capacity of a rough strip footing. The simulations provide insight into the application of uncertainty treatment to geotechnical problems and show the importance of the spatial variability of soil properties with regard to the outcome of a probabilistic assessment. Copyright © 2009 John Wiley & Sons, Ltd.     

Estimating evapotranspiration from small watersheds using a water and energy balance approach

Meteorological and environmental data measured in semiarid watersheds during the summer monsoon and winter periods were used to study the interrelationships among flux, meteorological and soil water variables, and to evaluate the effects of these variables on the daily estimation of actual evapotranspiration (AET). The relationship between AET and potential evapotranspiration (PET) as a function of soil water content, as suggested by Thornthwaite–Mather and by Morton, was studied to determine its applicability to the study area. Furthermore, multiple linear regression (MLR) analysis was employed to evaluate the order of importance of the meteorological and soil water factors involved. The results of MLR analysis showed that the combined effects of available energy, soil water content and wind speed were responsible for more than 70% of the observed variations in AET during the summer monsoon period. The analyses also indicate that the combined effects of available energy, vapour pressure deficit and wind speed were responsible for more than 70% of the observed variations in AET during the winter period. However, the test results of two different approaches, using the relationships between AET and PET as a function of soil water content, indicated some inadequacy. Copyright © 2007 John Wiley & Sons, Ltd.     

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).     

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.