Combined analyses of chemometrics and kriging for identifying groundwater contamination sources and origins at the Masan coastal area in Korea

Hydrogeochemical analyses including the basic statistics of chemical components, Piper??s trilinear diagram, and Mazor??s compositional bivariate diagram revealed that the main source and origin of groundwater contamination was seawater intrusion in the study area. However, the other sources and origins of groundwater contamination could be found by the combined analyses of chemometrics and kriging. Cluster analysis was helpful for the classification on the basis of the contamination characteristics of groundwater quality; however, it was not sufficient for the apportionment of groundwater contamination sources. Factor analysis (FA) determined three factors with 81.07% in total variance: Factor 1 for seawater contamination, Factor 2 for nitrate contamination, and Factor 3 for iron contamination. Factor analysis determined the sources of groundwater contamination; however, it could not discover the origins of contaminants except Factor 1. In backward stepwise mode, discriminant analysis decreased the number of parameters from 18 to 6 in discriminating the contaminant type with 96.2% correctness. TDS, Ca, NO₃, Mn, Fe, and Br were the most significant parameters for the discrimination of contaminants. Kriging analysis was very useful for the understanding of correlation and similarity between contaminants and factors of FA, and for the investigation of contaminant origins. It also showed that the similarity between factor scores and contaminant concentrations was proportional to the magnitudes of factor loadings for contaminants. This study represented that the combined analyses of chemometrics and kriging were very indispensable to the identification of groundwater contamination sources and origins, as well as for the spatial classification and assessment of groundwater quality.     

A semi-continuous measurement of gaseous ammonia and particulate ammonium concentrations in PM<Subscript>2.5</Subscript> in the ambient atmosphere

Ammonia has a short residence time in the atmosphere and rapidly neutralizes acid gases that occur near its source, requiring a rapid measurement system for ammonia and particulate ammonium concentrations to better understand their sources, temporal variation of ammonia emissions, and the formation of secondary ammonium aerosols. A semi-continuous measurement system, consisting of a diffusion scrubber, a particle growth chamber, an air-liquid separator, and a fluorescent detector, was developed to determine both gaseous ammonia (NH₃) and particulate ammonium (NH ₄ ⁺ ) in PM_2.5 in the ambient atmosphere of Gwangju, South Korea, during the months of March, April, July, and September of 2007. During the sampling periods, the average concentrations of ammonia and ammonium were found to be 2.33?±?1.29 μg/m³ and 1.89?±?0.99 μg/m³, respectively. Although the average gaseous ammonia concentration was highest in March, the particulate ammonium concentration was higher during the warmer season, reaching 2.08?±?1.07 μg/m³ and 2.32?±?0.94 μg/m³ in April and July, respectively, while only 1.68?±?0.61 μg/m³ in March and 1.24?±?0.99 μg/m³ in September. It is proposed that the higher availability of acid species during the warmer months produced a significant amount of particulate ammonium sulfate. Diurnal fluctuation of ammonia and ammonium during the warmer months showed that their peak time occurred at approximately 10:00 am. Both ammonia and ammonium concentrations were better correlated during the warmer months than during the cooler months. Further, the data suggest that the ammonia and ammonium were measured under well dispersed conditions, and multiple sources contributed to the ammonia at the sampling site.     

Seasonality in the relationship between El Nino and Indian Ocean dipole

The seasonal change in the relationship between El Nino and Indian Ocean dipole (IOD) is examined using the European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40), and the twentieth century simulations (20c3m) from the Geophysical Fluid Dynamics Laboratory Coupled Model, version 2.1. It is found that, both in ERA-40 and the model simulations, the correlation between El Nino (Nino3 index) and the eastern part of the IOD (90?C110°E; 10°S-equator) is predominantly positive from January to June, and then changes to negative from July to December. Correlation maps of atmospheric and oceanic variables with respect to the Nino3 index are constructed for each season in order to examine the spatial structure of their seasonal response to El Nino. The occurrence of El Nino conditions during January to March induces low-level anti-cyclonic circulation anomalies over the southeastern Indian Ocean, which counteracts the climatological cyclonic circulation in that region. As a result, evaporation decreases and the southeastern Indian Ocean warms up as the El Nino proceeds, and weaken the development of a positive phase of an IOD. This warming of the southeastern Indian Ocean associated with the El Nino does not exist past June because the climatological winds there develop into the monsoon-type flow, enhancing the anomalous circulation over the region. Furthermore, the development of El Nino from July to September induces upwelling in the southeastern Indian Ocean, thereby contributing to further cooling of the region during the summer season. This results in the enhancement of a positive phase of an IOD. Once the climatological circulation shifts from the boreal summer to winter mode, the negative correlation between El Nino and SST of the southeastern Indian Ocean changes back to a positive one.     

How predictable is the northern hemisphere summer upper-tropospheric circulation?

The retrospective forecast skill of three coupled climate models (NCEP CFS, GFDL CM2.1, and CAWCR POAMA 1.5) and their multi-model ensemble (MME) is evaluated, focusing on the Northern Hemisphere (NH) summer upper-tropospheric circulation along with surface temperature and precipitation for the 25-year period of 1981–2005. The seasonal prediction skill for the NH 200-hPa geopotential height basically comes from the coupled models’ ability in predicting the first two empirical orthogonal function (EOF) modes of interannual variability, because the models cannot replicate the residual higher modes. The first two leading EOF modes of the summer 200-hPa circulation account for about 84% (35.4%) of the total variability over the NH tropics (extratropics) and offer a hint of realizable potential predictability. The MME is able to predict both spatial and temporal characteristics of the first EOF mode (EOF1) even at a 5-month lead (January initial condition) with a pattern correlation coefficient (PCC) skill of 0.96 and a temporal correlation coefficient (TCC) skill of 0.62. This long-lead predictability of the EOF1 comes mainly from the prolonged impacts of El Niño-Southern Oscillation (ENSO) as the EOF1 tends to occur during the summer after the mature phase of ENSO. The second EOF mode (EOF2), on the other hand, is related to the developing ENSO and also the interdecadal variability of the sea surface temperature over the North Pacific and North Atlantic Ocean. The MME also captures the EOF2 at a 5-month lead with a PCC skill of 0.87 and a TCC skill of 0.67, but these skills are mainly obtained from the zonally symmetric component of the EOF2, not the prominent wavelike structure, the so-called circumglobal teleconnection (CGT) pattern. In both observation and the 1-month lead MME prediction, the first two leading modes are accompanied by significant rainfall and surface air temperature anomalies in the continental regions of the NH extratropics. The MME’s success in predicting the EOF1 (EOF2) is likely to lead to a better prediction of JJA precipitation anomalies over East Asia and the North Pacific (central and southern Europe and western North America).     

Polybrominated diphenyl ethers and polychlorinated biphenyls in sediments of southwest Taiwan: regional characteristics and potential sources

Very little information is available on the contamination of coastal sediments of Taiwan by PBDEs and PCBs. In this study, we determined the concentrations of 19 PBDE and 209 PCB congeners in 57 surface sediment samples to identify the possible sources of PBDEs and PCBs. The total PBDE and PCB concentrations ranged from below detection limit to 7.73 ng/g and 0.88-7.13 ng/g, respectively; these values are within the ranges observed for most coastal sediments worldwide. The PBDE congeners were dominated by BDE-209 (50.7-99.7%), with minor contributions from penta- and octa-BDEs. The signatures of PCB congeners suggested that PCB residues in Kaohsiung coast may be the legacy of past use or the result of ongoing inputs from the maintenance, repair and salvage of old ships. Principal component analysis of the congener-specific composition of PBDEs and PCBs revealed distinct regional patterns that are related to the use of commercial products.     

Long-term simulations of the sulfur concentrations over the China, Japan and Korea: A model comparison study

Three comprehensive acid deposition models were used to simulate the sulfur concentrations over northeast Asia over the period covering entire year of 2002, and discussed the aggregated uncertainties and discrepancies of the three models. The participating models are from the countries participating in the project of Longrange Transboundary Air Pollutants in Northeast Asia (LTP): China, Japan and Korea. The Eulerian Model-3/CMAQ (by China), Regional Air Quality Model (RAQM, by Japan), and Comprehensive Acid Deposition Model (CADM, by Korea) were employed by each country with common emissions data established by the administrative agencies of China, Japan and Korea. The episodic simulation results between 1 to 15, March 2002 are also presented, during which aircraft measurements were carried out over the Yellow sea. The episodic results show both a wide short-term variability in simulations against measurements, and maximum concentration differences of 3～5 times among the three models, requiring that further attention before confidence among the three models can be claimed for short-term simulations. However, the year-long cumulative simulations showed almost the same general features, with lower aggregated uncertainties between the three models, produced by the long term integration over northeast Asia.     

Representation of tropical storms in the northwestern pacific by the Modern-Era Retrospective analysis for research and applications

This study examines the tropical storms simulated in the Modern-Era Retrospective analysis for Research and Applications (MERRA) global atmospheric reanalysis for the recent 12 years (1998–2009), focusing on the tropical storm activity over the Northwestern Pacific. For validation, the International Best Track Archive for Climate Stewardship (IBTrACS) dataset is used as an observational counterpart. Climatological-mean features of the tropical storm genesis, tracks and their maximum intensity are the primary interests in this study. Regarding the genesis location of tropical storms, MERRA is reasonable in resolving major development regions over the South China Sea and the Northwestern Pacific close to the Philippines. The seasonal variation of the number of storms is also reproduced in a realistic way in MERRA, with peak values occurring from July to September. In addition, MERRA tends to reproduce the observed interannual variation of the number of tropical storms during the 12-years, though with a limited accuracy. The simulated paths toward higher latitudes are also reasonable in MERRA, where the reanalysis corresponds well with the observations in resolving frequent paths of westward moving storms and recurving storms toward the northeast. Regarding the intensity, MERRA captures the linear relationship between the minimum center pressure and the maximum wind speed near the surface at the maximum development. Some discrepancies from the observed features are found in the reanalysis, such as less frequent development of storms over the South China Sea and less frequent paths over this region. The reanalysis also does not attain the observed maximum intensity for the resolved tropical storms, particularly underestimating the center pressure. These deficiencies are likely related to limitations in the horizontal resolution and the parameterized physics of the data assimilation system.     

Characteristics of canopy turbulence over a deciduous forest on complex terrain

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.     

Behavior of epidote at high pressure and high temperature: a powder diffraction study up to 10 GPa and 1,200 K

The thermo-elastic behavior of a natural epidote [Ca_1.925 Fe_0.745Al_2.265Ti_0.004Si_3.037O₁₂(OH)] has been investigated up to 1,200 K (at 0.0001 GPa) and 10 GPa (at 298 K) by means of in situ synchrotron powder diffraction. No phase transition has been observed within the temperature and pressure range investigated. P–V data fitted with a third-order Birch–Murnaghan equation of state (BM-EoS) give V ₀ = 458.8(1)ų, K _T0 = 111(3) GPa, and K′ = 7.6(7). The confidence ellipse from the variance–covariance matrix of K _T0 and K′ from the least-square procedure is strongly elongated with negative slope. The evolution of the “Eulerian finite strain” vs “normalized stress” yields Fe(0) = 114(1) GPa as intercept values, and the slope of the regression line gives K′ = 7.0(4). The evolution of the lattice parameters with pressure is slightly anisotropic. The elastic parameters calculated with a linearized BM-EoS are: a ₀ = 8.8877(7) Å, K _T0(a) = 117(2) GPa, and K′(a) = 3.7(4) for the a-axis; b ₀ = 5.6271(7) Å, K _T0(b) = 126(3) GPa, and K′(b) = 12(1) for the b-axis; and c ₀ = 10.1527(7) Å, K _T0(c) = 90(1) GPa, and K’(c) = 8.1(4) for the c-axis [K _T0(a):K _T0(b):K _T0(c) = 1.30:1.40:1]. The β angle decreases with pressure, β_P(°) = β_P0 −0.0286(9)P +0.00134(9)P ² (P in GPa). The evolution of axial and volume thermal expansion coefficient, α, with T was described by the polynomial function: α(T) = α₀ + α₁ T ^−1/2. The refined parameters for epidote are: α₀ = 5.1(2) × 10⁻⁵ K⁻¹ and α₁ = −5.1(6) × 10⁻⁴ K^1/2 for the unit-cell volume, α₀(a) = 1.21(7) × 10⁻⁵ K⁻¹ and α₁(a) = −1.2(2) × 10⁻⁴ K^1/2 for the a-axis, α₀(b) = 1.88(7) × 10⁻⁵ K⁻¹ and α₁(b) = −1.7(2) × 10⁻⁴ K^1/2 for the b-axis, and α₀(c) = 2.14(9) × 10⁻⁵ K⁻¹ and α₁(c) = −2.0(2) × 10⁻⁴ K^1/2 for the c-axis. The thermo-elastic anisotropy can be described, at a first approximation, by α₀(a): α₀(b): α₀(c) = 1 : 1.55 : 1.77. The β angle increases continuously with T, with β_T(°) = β_T0 + 2.5(1) × 10⁻⁴ T + 1.3(7) × 10⁻⁸ T ². A comparison between the thermo-elastic parameters of epidote and clinozoisite is carried out.     

Performances of two instrumented laboratory models for the study of rainfall infiltration into unsaturated soils

Rainfall infiltration and suction variation in unsaturated soils must be taken into consideration in the analysis of most slope stability problems, particularly in the tropical regions where the annual precipitation is high. The process of rainfall infiltration into unsaturated soils is an extremely complex problem attributed to the non-linearity of the hydraulic property functions of the unsaturated soils. This paper describes in detail two instrumented laboratory models, i.e. one-dimensional soil column, and two-dimensional slope model used to provide experimental evidences for the transient suction variations in the unsaturated soils under certain rainfall conditions. The performances of the laboratory models were tested on four typical types of residual soils, i.e. sand-gravel, silty gravel, sandy silt, and silt (kaolin), and a two-layered soil system, i.e. sandy silt underlain by silty gravel. The results showed that the suction distributions for the single-layered homogeneous soils obtained from the simpler one-dimensional soil column were almost identical to that of two-dimensional slope model. However, the two-dimensional slope model should be employed for the two-layered soil system because of the dominant effect of the lateral flow mechanism. The capillary barrier effect was observed when a less permeable soil layer was underlain by a more permeable soil layer. The minimum suction value in soil is governed by the rainfall intensity, rainfall duration, and the saturated permeability of soil. The infiltration rate of the fine-grained soils that subjected to shrink and crack was independent of the soil permeability, but was significantly governed by the preferential flows developed in the soils.