Hydrogeochemical processes of thallium in a rural area of Southwest Guizhou, China

The water system in a rural area of Lanmuchang in Southwest Guizhou is facing a risk of thallium (Tl) contamination due to Tl mineralization around the area. The major trace elements and Tl in the water system are studied to understand the hydrogeochemical processes of Tl constrained by Tl mineralization. The results showed that the dispersion pattern of Tl follows a descending order in concentration from mine groundwater (deep groundwater) →stream water→shallow groundwater→background water, reflecting the impact of Tl mineralization on the hydrogeochemical composition. Tl concentrations in stream water in both regimes are remarkably higher (2-30 fold) downstream than up- and mid-streams, probably caused by the unidentified discharge of deep groundwater. Low Tl levels are detected in the current drinking water, however, the highly elevated Tl in stream water and ground water may pose a potential environmental risk through daily washing and agricultural irrigation. This study suggests that human activities, such as agricultural irrigation, could intensify the environmental risk of Tl.     

Experiments on the primary migration of oil from source rocks

Simulation experiments on the primary migration of oil were carried out on massive samples. The results proved that oil generated from source rocks was expelled in the form of an independent oil phase. High oil-expulsion efficiency was observed. It follows that the primary migration of oil is not directly dependent on the quantity of oil generated from the source rocks. Therefore, the oil-expulsion proportion was high though some source rocks yielded only a limited amount of oil. A great deal of gas was produced at the same time of oil-generation. Thus, it can be concluded that the main expulsion energy for oil primary migration came from these gases.     

Refinement of a lunar FeO mapping method

Clementine UVVIS and NIR data from the lunar sampling sites (totaling 46 sampling sites) were processed and used to refine the iron determination method of Le Mouéic et al. (2000, 2002). We found that about 21 sampling sites are unsuitable to Le Mouéic et al's spectral parameters ("slope" and "depth l") because their 1500 nm filter could not be used as spectral parameters' right shoulders and to evaluate the depth of the 1-μm absorption feature accurately. We used the rest 25 sampling sites to refine the method developed by Le Mouéic et al. (2000, 2002 ) and obtained our own equation of FeO content determination. We tested our own equation, and the results are satifactory. In our work we also acquired some useful experiences in scientific applications of our own dataset of the Chang'E-1 mission.     

Geochemical characteristics of crude oils from Well Zheng-1 in the Junggar Basin, Xinjiang, China

Well Zheng-1 is located in the combined area of the central uplift and the north Tianshan piedmont depression in the Junggar Basin. Two oil-bearing beds are recognized at 4788–4797 m of the Lower Cretaceous Tugulu Formation (K1tg) and 4808.5–4812.5 m of the Lower Jurassic Sangonghe Formation (J1s). The geochemical characteristics of family composition, carbon isotopic composition, saturated hydrocarbons, sterane and terpane biomarkers and carotane of two crude oils are described in this paper. The results show that the geochemical characteristics of the two crude oils are basically similar to each other, indicating they were all derived mainly from the high mature, brine, algae-rich lake facies sediments. Oil-source correlation revealed that crude oils of the two beds were derived mainly from the source rocks of Permian and mixed by the oil derived from the source rocks of Jurassic and Triassic. This is consistent with the geological background with several sets of source rocks in the area studied.     

Fishing for diatoms: fish gut analysis reveals water quality changes over a 75-year period

Defining the reference state of some perturbed streams and rivers is challenging since their non-altered state is undocumented. Despite the near-absence of pristine sites in highly anthropogenic areas, information about aquatic communities that existed prior to human disturbance can be obtained from historic sources. Sediment coring is commonly used in paleolimnology to estimate past environmental conditions of lakes, but this technique is often not appropriate for studies of lotic systems. An alternative is to use diatom assemblages present in the guts of fish that were captured in the streams prior to significant human disturbance. The purpose of the present study was to evaluate the biological integrity of several streams in Ontario and Québec during the early twentieth century based on “paleo” diatom assemblages extracted from the guts of fish stored in museums. The Eastern Canadian diatom index (IDEC: Indice Diatomées de l’Est du Canada) was used to evaluate the biological status of “paleo” and “modern” diatom assemblages. The IDEC shows the position of diatom assemblages on a general pollution gradient. The comparison of IDEC values calculated for the 1925–1948 and the 2003–2007 periods showed that several streams were severely polluted in the early 1900s. In general, present water quality has declined compared to the early 1900s. The biological integrity of only three of the 22 sites has increased. IDEC values were not influenced by the species of fish studied.     

An Efficient Non-iterative Bulk Parametrization of Surface Fluxes for Stable Atmospheric Conditions Over Polar Sea-Ice

In climate and weather prediction models the near-surface turbulent fluxes of heat and momentum and related transfer coefficients are usually parametrized on the basis of Monin–Obukhov similarity theory (MOST). To avoid iteration, required for the numerical solution of the MOST equations, many models apply parametrizations of the transfer coefficients based on an approach relating these coefficients to the bulk Richardson number \(Ri_{b}\). However, the parametrizations that are presently used in most climate models are valid only for weaker stability and larger surface roughnesses than those documented during the Surface Heat Budget of the Arctic Ocean campaign (SHEBA). The latter delivered a well-accepted set of turbulence data in the stable surface layer over polar sea-ice. Using stability functions based on the SHEBA data, we solve the MOST equations applying a new semi-analytic approach that results in transfer coefficients as a function of \(Ri_{b}\) and roughness lengths for momentum and heat. It is shown that the new coefficients reproduce the coefficients obtained by the numerical iterative method with a good accuracy in the most relevant range of stability and roughness lengths. For small \(Ri_{b}\), the new bulk transfer coefficients are similar to the traditional coefficients, but for large \(Ri_{b}\) they are much smaller than currently used coefficients. Finally, a possible adjustment of the latter and the implementation of the new proposed parametrizations in models are discussed.     

Fault plane solutions of the january 26th, 2001 bhuj earthquake sequence

A 12-station temporary microearthquake network was established by the Geological Survey of India for aftershock monitoring of the January 26th, 2001 Bhuj earthquake (M _w 7.6) in the Kutch district of Gujarat state, western India. The epicentres of the aftershocks show two major trends: one in the NE direction and the other in the NW direction. Fault-plane solutions of the best-located and selected cluster of events that occurred along the NE trend, at a depth of 15–38 km, show reverse faulting with a large left-lateral strike-slip motion, which are comparable with the main-shock solution. The NW trending upper crustal aftershocks at depth <10 km, on the other hand, show reverse faulting with right-lateral strike-slip motion, and the mid crustal and lower crustal aftershocks, at a depth of 15–38 km, show pure reverse faulting as well as reverse faulting with right-lateral and left-lateral strike-slip motions; these solutions are not comparable with the main-shock solution. It is inferred that the intersection of two faults has been the source area for stress concentration to generate the main shock and the aftershocks.     

An efficient solution of real-time data processing for multi-GNSS network

Global navigation satellite systems (GNSS) are acting as an indispensable tool for geodetic research and global monitoring of the Earth, and they have been rapidly developed over the past few years with abundant GNSS networks, modern constellations, and significant improvement in mathematic models of data processing. However, due to the increasing number of satellites and stations, the computational efficiency becomes a key issue and it could hamper the further development of GNSS applications. In this contribution, this problem is overcome from the aspects of both dense linear algebra algorithms and GNSS processing strategy. First, in order to fully explore the power of modern microprocessors, the square root information filter solution based on the blocked QR factorization employing as many matrix–matrix operations as possible is introduced. In addition, the algorithm complexity of GNSS data processing is further decreased by centralizing the carrier-phase observations and ambiguity parameters, as well as performing the real-time ambiguity resolution and elimination. Based on the QR factorization of the simulated matrix, we can conclude that compared to unblocked QR factorization, the blocked QR factorization can greatly improve processing efficiency with a magnitude of nearly two orders on a personal computer with four 3.30 GHz cores. Then, with 82 globally distributed stations, the processing efficiency is further validated in multi-GNSS (GPS/BDS/Galileo) satellite clock estimation. The results suggest that it will take about 31.38 s per epoch for the unblocked method. While, without any loss of accuracy, it only takes 0.50 and 0.31 s for our new algorithm per epoch for float and fixed clock solutions, respectively.     

Impact of Indian ocean dipole on the coastal upwelling features off the southwest coast of India

A three-dimensional regional ocean model is used to examine the impact of positive Indian ocean dipole (pIOD) events on the coastal upwelling features at the southwest coast of India (SWCI). Two model experiments are carried out with different surface boundary conditions that prevailed in the normal and pIOD years from 1982 to 2010. Model experiments demonstrate the weakening of coastal upwelling at the SWCI in the pIOD years. The reduced southward meridional wind stress off the SWCI leads to comparatively lower offshore Ekman transport during August–October in the pIOD years to that in normal years. The suppressed coastal upwelling results in warmer sea surface temperature and deeper thermocline in the pIOD years during June–September. The offshore spatial extent of upwelled colder (<?22 °C) water was up to 75.5° E in August–September in normal years that was limited up to 76.2° E in pIOD years. The heat budget analysis reveals the decreased contribution of vertical entrainment process to the mixed layer cooling in pIOD years which is almost half of that of normal years in October. The net heat flux term shows warming tendency during May–November with a higher magnitude (+?0.4 °C day^?1) in normal years than pIOD years (+?0.28 °C day^?1). The biological productivity is found to reduce during the pIOD years as the concentration of phytoplankton and zooplankton decreases over the region of coastal upwelling at SWCI. Nitrate concentration in the pIOD years dropped by half during August–September and dropped by an order of magnitude in October as compared to its ambient concentration of 13 μmol L^?1 in normal years.     

Multiple Random Walk Simulation: A Fast Method to Map Grade Uncertainty with Large Datasets

Multiple Random Walk Simulation consists of a methodology adapted to run fast simulations if close-spaced data are abundant (e.g., short-term mining models). Combining kriging with the simulation of random walks attempts to approximate traditional simulation algorithm results but at a computationally faster way when there is a large amount of conditioning samples. This paper presents this new algorithm illustrating the situations where the method can be used properly. A synthetic study case is presented in order to illustrate the Multiple Random Walk Simulation and to analyze the speed and goodness of its results against the ones from using Turning Bands Simulation and Sequential Gaussian Simulation.