Complexation of mercury by dissolved organic matter in surface waters of Galveston Bay, Texas

The chemical speciation of dissolved mercury in surface waters of Galveston Bay was determined using the concentrations of mercury-complexing ligands and conditional stability constants of mercury-ligand complexes. Two classes of natural ligands associated with dissolved organic matter were determined by a competitive ligand exchange-solvent solvent extraction (CLE-SSE) method: a strong class (L_s), ranging from 19 to 93 pM with an average conditional stability constant (K_HgLs) of 10²⁸, and a weak class (L_w) ranging from 1.4 to 9.8 nM with an average K_HgLs of 10²³. The range of conditional stability constants between mercury and natural ligands suggested that sulfides and thiolates are important binding sites for dissolved mercury in estuarine waters. A positive correlation between the estuarine distribution of dissolved glutathione and that of mercury-complexing ligands supported this suggestion. Thermodynamic equilibrium modeling using stability constants for HgL, HgCl_x, Hg(OH)_x, and HgCl(OH) and concentrations of each ligand demonstrated that almost all of the dissolved mercury (> 99%) in Galveston Bay was complexed by natural ligands associated with dissolved organic matter. The importance of low concentrations of high-affinity ligands that may originate in the biological system (i.e., glutathione and phytochelatin) suggests that the greater portion of bulk dissolved organic matter may not be important for mercury complexation in estuarine surface waters.     

New variogram modeling method using MGGP and SVR

A critical step for kriging in geostatistics is estimation of the variogram. Traditional variogram modeling comprise of the experimental variogram calculation, appropriate variogram model selection and model parameter determination. Selecting of the variogram model and fitting of model parameters is the most controversial aspect of geostatistics. Shapes of valid variogram models are finite, and sometimes, the optimal shape of the model can not be fitted, leading to reduced estimation accuracy. In this paper, a new method is presented to automatically construct a model shape and fit model parameters to experimental variograms using Support Vector Regression (SVR) and Multi-Gene Genetic Programming (MGGP). The proposed method does not require the selection of a variogram model and can directly provide the model shape and parameters of the optimal variogram. The validity of the proposed method is demonstrated in a number of cases.     

Application of S-SEBI model for crop evapotranspiration using Landsat-8 data over parts of North India

Estimation and monitoring of crop evapotranspiration (ET_c) or consumptive water use over large-area holds the key to irrigation management plans and regional drought preparedness. The objective of this study was to estimate ET_c by applying the simplified-surface energy balance index (S-SEBI) model to Landsat-8 data for the 2014–2015 period in parts of North India. An average ET_c was estimated 2.72 and 2.47 in mm day^?1 with 0.22, 0.18 standard deviation and 0.11, 0.07 standard error for Kharif and Rabi crops, respectively. On validation part, a close relationship was observed between S-SEBI derived and scintillometer observed evaporative fraction with 0.85 correlation coefficient and 0.86 agreement index. The statistical analysis also endorses the results accuracy and reliability with 0.026 and 0.602, relative root-mean square errors and model efficiency for wheat crop, respectively. The study showed that normalized difference vegetation index and LST are closely related and serve as a proxy for qualitative representation of ET_c.     

Mössbauer studies of iron silicate spinel at high pressure

We have measured in situ Mössbauer transmission spectra of iron silicate spinel (γ-Fe₂SiO₄) in a diamond anvil cell at room temperature and pressures up to 16 GPa. The observed spectra show a doublet characteristic of the paramagnetic state. The isomer shift and quadrupole splitting at atmospheric pressure are 1.10 and 2.63 mm/s, respectively, which are smaller than those of fayalite (α-Fe₂SiO₄). Both the isomer shift and quadrupole splitting decrease linearly with pressure with slope of ?0.003(1) and ?0.020(1) mm/sec · GPa, respectively. This simple linear trend suggests that no electronic or polymorphic transitions occur under 16 GPa except for those due to the small and continuous changes of volume and local symmetry under pressure. On the basis of a crystalline field calculation, the negative pressure derivative of the quadrupole splitting is associated with a trend towards an ideal cubic symmetry of the oxygen sublattice.     

A simulation of Asian dust events in March 2010 by using the ADAM2 model

Asian dust events occurred in Asia during March 2010 were simulated using the Asian Dust Aerosol Model 2 (ADAM2). The performance of the model for simulations of surface dust concentrations and dust event occurrences was tested at several monitoring sites located in the dust source region and the downstream region of Korea. The observed and modeled dust event occurrences at each monitoring site were defined with the hourly observed and modeled dust concentrations that were used to evaluate the performance of the model by constructing a contingency table for the dust event occurrence. It was found that the model simulated quite well the starting and ending times of dust events with their peak dust concentrations for most dust events occurred both in the dust source region and the downstream region of Korea. However, the model failed to simulate a few dust events observed in both regions mainly due to the inaccurate simulations of the meteorological fields. Inaccurate simulations of wind speeds have caused for the model to simulate dust events poorly in the dust source region whereas poor simulations of precipitation of the fifth-generation mesoscale model (MM5) model have led to miss dust events in the downstream region of Korea. The contingency table made with the hourly data for the dust event occurrence made it possible to evaluate the ADAM2 model for the simulation of the dust event occurrence. It was found that the model has the probabilistic simulation capability for dust events of about 78% with the hit rate of more than 83% and the false alarm rate of about 27% for the dust events occurred during March in 2010. The probabilistic capability of the model could be much improved by improving the meteorological model (MM5 model).     

The Asian Dust Aerosol Model 2 (ADAM2) with the use of Normalized Difference Vegetation Index (NDVI) obtained from the Spot4/vegetation data

The operational Asian Dust Aerosol Model (ADAM)1 in Korea Meteorological Administration has been modified to the ADAM2 model to be used as an operational forecasting model all year round not only in Korea but also in the whole Asian domain (70-160°E and 5-60°N) using the routinely available World Meteorological Organization (WMO) surface reporting data and the Spot/vegetation Normalized Difference Vegetation Index (NDVI) data for the period of 9 years from 1998 to 2006. The 3-hourly reporting WMO surface data in the Asian domain have been used to re-delineate the Asian dust source region and to determine the temporal variation of the threshold wind speed for the dust rise. The dust emission reduction factor due to vegetation in different surface soil-type regions (Gobi, sand, loess, and mixed soil) has been determined with the use of NDVI data. It is found that the threshold wind speed for the dust rise varies significantly with time (minimum in summer and maximum in winter) and surface soil types with the highest threshold wind speed of 8.0 m?s^?1 in the Gobi region and the lowest value of 6.0 m?s^?1 in the loess region. The statistical analysis of the spot/vegetation NDVI data enables to determine the emission reduction factor due to vegetation with the free NDVI value that is the NDVI value without the effect of vegetation and the upper limit value of NDVI for the dust rise in different soil-type regions. The modified ADAM2 model has been implemented to simulate two Asian dust events observed in Korea for the periods from 31 March to 2 April 2007 (a spring dust event) and from 29 to 31 December 2007 (a winter dust event) when the observed PM₁₀ concentration at some monitoring sites in the source region exceeds 9,000 μg m^?3. It is found that ADAM2 model successfully simulates the observed high dust concentrations of more than 8,000 μg m^?3 in the dust source region and 600 μg m^?3 in the downstream region of Korea. This suggests that ADAM2 has a great potential for the use of an operational Asian dust forecast model in the Asian domain.     

A Measure of Magnetic Field Discontinuity

In this paper we introduce a measure of magnetic field discontinuity, MAD, defined as Maximum Angular Difference between two adjacent magnetic field vectors. To examine the characteristics of the MAD, we have considered several active region models having a quadrupolar field configuration and computed MADs over these active regions by approximating the 3-D magnetic fields as an ensemble of charge potential fields or linear force-free fields. The computed MAD fields are studied in comparison with other flare activity indicators such as separators. It is found that (1) the region of high level MAD corresponds well with the separator, or practically the intersection of the separator with the plane of measurement, (2) it singles out local discontinuities of magnetic fields, and (3) the MAD can also be a measure of the evolutionary status of an active region.An observational test has been made for 2-D MADs, using the Yohkoh SXT observation of a flare in AR 6919 and the vector magnetogram taken at the Mees Solar Observatory during this flare activity. The high level contours of 2-D MAD are found to trace well the observed soft X-ray bright points, which indicates that the MAD could serve as a good flare activity indicator.     

Generation of rotational discontinuities by magnetic reconnection associated with microflares

Magnetic reconnection can take place between two plasma regions with antiparallel magnetic field components. In a time-dependent reconnection event, the plasma outflow region consists of a leading bulge region and a trailing reconnection layer. Magnetohydrodynamic (MHD) discontinuities, including rotational discontinuities, can be formed in both the bulge region and the trailing layer. In this paper, we suggest that the rotational discontinuities observed in the solar wind may be generated by magnetic reconnection associated with microflares in coronal holes. The structure of the reconnection layer is studied by solving the one-dimensional Riemann problem for the evolution of an initial current sheet after the onset of magnetic reconnection as well as carrying out two-dimensional MHD simulations. As the emerging magnetic flux reconnects with ambient open magnetic fields in the coronal hole, rotational discontinuities are generated in the region with open field lines. It is also found that in the solar corona with a low plasma beta ( 0.01), the magnetic energy is converted through magnetic reconnection mostly into the plasma bulk-flow energy. Since more microflares will generate more rotational discontinuities and also supply more energy to the solar wind, it is expected that the number of rotational discontinuities observed in the solar wind would be an increasing function of solar wind speed. The observation rate of rotational discontinuities generated by microflares is estimated to be dN _RD/dt - f/63 000 s (f > 1) at 1 AU. The present mechanism favors the generation of rotational discontinuities with a large shock normal angle.     

Precise calculation of the magnetosphere surface for a tilted dipole

The shape of the magnetosphere has been calculated self-consistently for inclinations of the Earth's magnetic dipole from perpendicular to the solar wind. Inclination angles of 0–35° have been chosen in steps of 5° and various smooth trends in the surface characteristics with increasing inclination angle noted. The surface points and the complete field at the surface points have been calculated for the entire surfaces. The neutral point region has been given precise study in one degree steps and is found to be tangent to the solar wind velocity and to have a smooth continuous curvature. The inclusion of the Earth's bow shock pressure and other sources of current have not been included at this stage of our program.     

Preliminary investigation of Late Pleistocene fauna from Ryonggok Cave No. 1, Sangwon County,North Hwanghae Province,Democratic People's Republic of Korea

The investigation of a limestone cave on the Ryonggok-ri of Sangwon County, North Hwanghae Province, DPRK, has yielded many faunal remains including Late Pleistocene deer, brown bears and horses. Uranium series dating of fossil teeth yielded dates of ~72ka at the base of Layer 2 – the bone-bearing unit – and ~44ka at the top. The rich and diverse mammal remains of Ryonggok Cave No. 1 are indicative of warm and humid temperate conditions during the Late Pleistocene, and reflect the presence of both dense forests and open grassland. In addition, pollen and spores from the site include two families and 24 genera of trees, six families and a genus of grasses and herbs, and two families and six genera of ferns, corroborating the environments inferred from the faunal remains.