Geochemistry and environmental effect of cadmium in the super-large Jinding Pb-Zn deposit,Yunnan Province,China

The Jinding Pb-Zn deposit in Yunnan Province is the representative of a Cd-enriched area and mining activities lead to the release of Cd into the hypergenic ecosystem, resulting in Cd pollution. The concentrations of Cd vary greatly from one type to another type of rocks in the mining district. In the host rock, Cd concentrations range from 50×10^-6 to 650×10^-6 with an average of 310×10^-6. In primary ores, Cd concentrations range from 14×10^-6 to 2800×10^-6 with an average of 767×10^-6. However, in oxidized ores, Cd concentrations are highest, varying within the range of 110×10^-6 to 8200×10^-6 , averaging 1661×10^-6. It is shown that the oxidized ores are the main carder and environmental source of Cd. Leaching test showed that Pb/Zn ores are easy to oxidize and thereafter release Cd and other harmful elements. These leached elements in the leachate may be precipitated rapidly in the order of Zn〉Pb〉Cd. As for the concentration distribution of Cd in the Bijiang River, it is estimated to be 15.7 μg/L Cd in water, 49.3 mg/L in suspended substances, and 203.7 mg/L in sediments. The average value of Cd in soil from the polluted area is 83.0 mg/kg. Natural weathering of Cd-rich rocks and minerals imposes a potential environmental risk on the aquatic ecosystem of the Bijiang catchment.     

Changes in nutrient uptake of phytoplankton under the interaction between sunlight and phosphate in the Changjiang （Yangtze） River Estuary

We conducted ship-board incubation experiments to investigate changes in nutrient uptake of phytoplankton under different phosphate concentrations and irradiances in the Changjiang River Estuary and its adjacent waters in China. Under 100% natural irradiance the uptake rates of phosphate, silicate, and nitrate were accelerated at high phosphate levels （1.84 μM）, while under low irradiance （about 50% natural irradiance） their uptake rates were restrained at the high but stimulated greatly at the intermediate phosphate concentrations （1.26 μM）, as the growth of phytoplankton, changes in nitrite and ammonium uptake didn＇t follow an obvious pattern. Our results also showed that there were linear relationships between nitrate, silicate and phosphate uptake at different phosphate concentrations under low and high irradiances, and the growth period of phytoplankton was prolonged both at the high phosphate concentrations under high irradiance and at the intermediate concentrations under low irradiance, suggesting that the limitation of phytoplankton growth mainly reflected changes in its growth period, and because no such environment （low irradiance and low phosphate concentrations） actually existed in a high turbidity zone, phytoplankton blooms hardly occurred there. In the absence of irradiance, denitrification occurred readily and phytoplankton was kept decreasing, which resulted in phosphate regeneration.     

Simulation experiment on hydrocarbon generating potential of various source rocks on the northern margin of the Qaidam Basin,Qinghai, China

The northern margin of the Qaidam Basin is one of the main oil-gas exploration areas in western China, where source rocks are composed mostly of Middle and Lower Jurassic dark mud shale, carbargillite and coal. A large number of subsurface and outcrop samples differing in lithology with different types of organic matter were selected for resource evaluation, research and calculation. And among them, 13 samples were used for simulation experiment on hydrocarbon generating potential of various source rocks. At first, two kinds of heating modes were compared through simulation experiment, including single temperature-step heating and continual heating. Perhaps, the process of hydrocarbon generation and expulsion occurred naturally between a close system and an open system. In addition, the first heating mode was convenient, and all its reaction products were involved in the whole thermal evolution, and the final simulation experimental results were basically in consistency with the natural evolution trend Therefore, the first heating mode was adopted and the hydrocarbon yield of every sample was worked out. According to the type and lithology of organic matter and the hydrocarbon yield of samples for simulation experiment, hydrocarbon generation and expulsion mode with three kinds of lithology and five types of source rock has been established to provide the basis for hydrocarbon generation evaluation, research and resource calculation.     

<Superscript>40</Superscript>Ar-<Superscript>39</Superscript>Ar dating of lamprophyre dykes in the Baimazhai nickel deposit,Yunnan Province,China, and its geological significance

The Baimazhai nickel deposit, Yunnan Province, China, is located in the southern part of the Sanjiang （Tri-river） alkali-rich intrusive rock belt （Sanjiang ARIR）. In this paper was conducted ^40Ar-^39Ar dating of two phlogopites in lamprophyres which are, as dikes, widely distributed in the orefield, and two plateau ages were acquired, i.e., 32.46±0.62 Ma and 32.01±0.60 Ma, respectively （averaging 32.23±0.60 Ma）. The ages are obviously younger than those of the ore-hosted complex and mineralization of the Baimazhai nickel deposit. In combination with the characteristics, it is indicated that lamprophyres in the orefield and those in the Sanjiang ARIR are similar in tectonic setting, mineral assemblage and geochemistry. It is considered that lamprophyres in the orefield are the important component of the Sanjiang ARIR, and the lamprophyres and ore-hosted complex in the orefield represent the products of two times of different magmatic activity from different mantle sources. On the other hand, the age of lamprophyres in the orefield is older than that of the strike-slip shearing of the Ailaoshan-Honghe fault belt, suggesting that the strike-slip shearing of the Ailaoshan-Honghe fault belt is not the factor leading to magmatic activity of lamprophyres in the orefield, while it is more possible that magmatic activity of the Sanjiang ARIR promoted strike-slip shearing of the fault belt.     

A study on the polarimetric properties of various features using SIR-C data

Spaceborne Imaging Radar (SIR-C) data acquired over Gujarat, India in 1994 were processed and analysed using differnet techniques applicable to polarimetric SAR data such as polarization signatures, polarization index, decomposition of the signal and polarization phase difference and limited groundtruth data. It has been observed that multi-frequency polarimetric data enhances the potential of retrieving geo-physical parameters. The polarization signatures are found to vary with the nature of the target. Target decomposition of the returned signal will be useful for the classification of various features. Polarization Phase Difference (PPD) gives good information about the vegetation parameters.     

The GeoForschungsZentrum Potsdam/Groupe de Recherche de Gèodésie Spatiale satellite-only and combined gravity field models: EIGEN-GL04S1 and EIGEN-GL04C

The recent improvements in the Gravity Recovery And Climate Experiment (GRACE) tracking data processing at GeoForschungsZentrum Potsdam (GFZ) and Groupe de Recherche de Géodésie Spatiale (GRGS) Toulouse, the availability of newer surface gravity data sets in the Arctic, Antarctica and North-America, and the availability of a new mean sea surface height model from altimetry processing at GFZ gave rise to the generation of two new global gravity field models. The first, EIGEN-GL04S1, a satellite-only model complete to degree and order 150 in terms of spherical harmonics, was derived by combination of the latest GFZ Potsdam GRACE-only (EIGEN-GRACE04S) and GRGS Toulouse GRACE/LAGEOS (EIGEN-GL04S) mean field solutions. The second, EIGEN-GL04S1 was combined with surface gravity data from altimetry over the oceans and gravimetry over the continents to derive a new high-resolution global gravity field model called EIGEN-GL04C. This model is complete to degree and order 360 and thus resolves geoid and gravity anomalies at half- wavelengths of 55 km at the equator. A degree-dependent combination method has been applied in order to preserve the high accuracy from the GRACE satellite data in the lower frequency band of the geopotential and to form a smooth transition to the high-frequency information coming from the surface data. Compared to pre-CHAMP global high-resolution models, the accuracy was improved at a spatial resolution of 200 km (half-wavelength) by one order of magnitude to 3 cm in terms of geoid heights. The accuracy of this model (i.e. the commission error) at its full spatial resolution is estimated to be 15 cm. The model shows a reduced artificial meridional striping and an increased correlation of EIGEN-GL04C-derived geostrophic meridional currents with World Ocean Atlas 2001 (WOA01) data. These improvements have led to select EIGEN-GL04C for JASON-1 satellite altimeter data reprocessing. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The Bayesian detection of discontinuities in a polynomial regression and its application to the cycle-slip problem

This paper deals with the problem of detecting and correcting cycle-slips in Global Navigation Satellite System (GNSS) phase data by exploiting the Bayesian theory. The method is here applied to undifferenced observations, because repairing cycle-slips already at this stage could be a useful pre-processing tool, especially for a network of permanent GNSS stations. If a dual frequency receiver is available, the cycle-slips can be easily detected by combining two phase observations or phase and range observations from a single satellite to a single receiver. These combinations, expressed in a distance unit form, are completely free from the geometry and depend only on the ionospheric effect, on the electronic biases and on the initial integer ambiguities; since these terms are expected to be smooth in time, at least in a short period, a cycle-slip in one or both the two carriers can be modelled as a discontinuity in a polynomial regression. The proposed method consists in applying the Bayesian theory to compute the marginal posterior distribution of the discontinuity epoch and to detect it as a maximum a posteriori (MAP) in a very accurate way. Concerning the cycle-slip correction, a couple of simultaneous integer slips in the two carriers is chosen by maximazing the conditional posterior distribution of the discontinuity amplitude given the detected epoch. Numerical experiments on simulated and real data show that the discontinuities with an amplitude 2 or 3 times larger than the noise standard deviation are successfully identified. This means that the Bayesian approach is able to detect and correct cycle-slips using undifferenced GNSS observations even if the slip occurs by one cycle. A comparison with the scientific software BERNESE 5.0 confirms the good performance of the proposed method, especially when data sampled at high frequency (e.g. every 1 s or every 5 s) are available.     

A data-driven approach to local gravity field modelling using spherical radial basis functions

We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.