Towards a Method for Quantitative LA‐ICP‐MS Imaging of Multi‐Phase Assemblages: Mineral Identification and Analysis Correction Procedures

Here, we present an approach to laser ablation ICP‐MS mapping of multi‐phase assemblages that permits the use of different internal standard elements, concentration values and reference materials for each mineral. In this way, we obtain not only broad pictures of elemental distributions within samples but can also extract high accuracy concentration data for any user‐selected region. This is accomplished by assigning regions of an image to corresponding mineral phases on a pixel‐by‐pixel basis. In this way, accurate trace element concentrations can be determined for each mineral phase, despite potential variations in their ablation characteristics. We present an example where elemental maps are constructed from ablation of a gabbroic sample that includes the phases apatite, amphibole and plagioclase. This work represents an important first step towards development of a method to produce highly accurate LA‐ICP‐MS elemental maps of multi‐phase samples.     

Data Analysis of Chang’E-1 Gamma-Ray Spectrometer and Global Distribution of U, K, and Th Elemental Abundances

Gamma-ray spectrometer (GRS) is one of the main payloads on the Chang’E-1 (CE-1) lunar probe, mainly aimed to detect the elemental abundances and distributions on the lunar surface. At 03:58 on 28 November 2007, it performed the first observation of the lunar gamma rays. As of 24 October 2008, 2105?h of effective gamma rays spectra had been acquired by CE-1 GRS, which covers the whole surface of the moon. This paper mainly describes the data processing procedures and methods of deriving the elemental abundances by using the CE-1 GRS time series corrected spectra: first, to bin data into pixels for mapping; then, to perform a background deduction of the cumulative spectra and obtain a peak area of the elements; and finally, to use the elemental abundances inversion model to produce the elemental abundances. Based on these processing methods, the global abundance maps of U, K, and Th at a 5°×5° equal-area pixel are acquired by CE-1 GRS data. The paper gives a preliminary analysis of the uncertainties of the elemental abundances     

基于Surfer的有限元后处理等值线图批处理程序开发

利用Surfer8.0的二次开发平台Scripter以及其与VisualBasic良好的兼容性,编制了等值线图批处理程序,实现了一次性大批等值线图的自动绘制。该程序可依据用户需要,进行指定样式的等值线图形自动拼合。为了操作方便还建立了友好的操作界面,有效地提高了有限元后处理的效率。     

Delineation of Integrated Anomaly with Generative Adversarial Networks and Deep Neural Networks in the Zhaojikou Pb-Zn Ore District,Southeast China

Geochemical maps are of great value in mineral exploration. Integrated geochemical anomaly maps provide comprehensive information about mapping assemblages of element concentrations to possible types of mineralization/ore, but vary depending on expert’s knowledge and experience. This paper aims to test the capability of deep neural networks to delineate integrated anomaly based on a case study of the Zhaojikou Pb-Zn deposit, Southeast China. Three hundred fifty two samples were collected, and ea...     

DEFLOGIC: a method for assessment of groundwater potential in karst terrains: Gurpi Anticline,southwest Iran

A new hydro-tectonic model, includes eight layers that affect karst hydrogeology was proposed for mapping of groundwater potential in karst areas of Gurpi Anticline, southwest Iran. To produce the groundwater potential map, remote sensing (RS) and GIS techniques were combined with fuzzy logic modeling. Criterion maps include the distances from discharge sites (D), the elevation difference from discharge sites (E), the distance from fractures (F), the fracture length density (L), the slope (O), the lithology (G), the distance from fractures intersections (I), and the fractures intersection density (C) were produced using GIS and RS techniques (DEFLOGIC layers). The approach of fuzzy sets was used to commensurate criterion maps, then fuzzy algebraic sum and gamma operators were applied to aggregate them. The weights of parameters of DEFLOGIC proposed in the range of 1 to 5, which standardized between 0 to 1, based on their importance in karst hydrogeology, professional judgments, and available exploration data. The final groundwater potential maps were verified by geoelectric and well-drilling data. The potential map prepared using fuzzy gamma operator with γ?=?0.92, which it is a flexible distinctive parameters of sum and product of fuzzy operator, depicts the best coincidence with exploration data. The final DEFLOGIC map shows the high groundwater potential in karst formations between Hati and Pebde valleys. The results support the efficiency of DEFLOGIC model to evaluating of groundwater potential in karst terrains, especially in Zagros ranges.     

Organic geochemistry of the kerogen and bitumen from the Ohio Shale by pyrolysis-gas chromatography

An exposure of the Ohio Shale at Bristol Ridge, Ohio in the Bellefontaine Outlier was sampled into 53 channel units. Proximate and elemental (CHN) analyses were performed on these samples to determine the concentration of organic matter. The organic content averaged 6.82% with the highest both in the lower half (11.08%) and extreme upper part (10.43%) of the section.Fragmentation of the kerogen of each of the whole rock samples was accomplished by pyrolyzing for 10 s at 750°C in the injection port of a gas chromatograph with a FID detector. The concentration and distribution of the diagnostic hydrocarbons fragmented at this temperature were used to infer the molecular construction of the high molecular weight kerogens. The suite of identified hydrocarbons for each of the samples throughout the section was similar, but a very large fluctuation in the amount of pyrolysis products was observed. This fluctuation is strongly related to organic richness as defined by the proximate and elemental analyses. The identification of diagnostic normal alkanes and aromatics and SCI indicates that the kerogen assemblage is composed of both marine and terrestial components — most likely a mixed Type II kerogen. The contribution of individual aromatic and normal alkanes is variable, which suggests that there were minor fluctuations in the proportions of terrestrial and marine organic matter during decomposition.Bitumen pyrolized from the whole rock and bitumen removed by Soxhlet extraction using a 4:1 benzene/methanol solvent mixture were analyzed. Comparison of these two techniques yield similar pristane/phytane ratios and Carbon Preference Index (CPI) values. On decomposition the whole rock bitumen exhibited a mean pristane/phytane ratio of 1.5 and a mean CPI value of 1.4, whereas the Soxhlet-extracted bitumen exhibited a mean pristane/phytane ratio of 1.6 and a mean value of 1.3. The low pristane/phytane ratios from this outlier indicates a possible marine environment, with reducing conditions existing immediately below the sediment-water interface. The CPI values also indicate the sedimentary sequence is thermally immature. Reflectance of the vitrinite (%R_o 0.05) and fluorescence (Q 0.93) of the algal bodies indicate a maturation at the beginning of the oil window or at boundary between thermal immaturity and maturity.     

全国航放数据库的建设

全国航放数据库系统成功地解决了不同类型航放测量数据源的录入、重复数据检索等关键技术问题并已录入30多个测区的航放(部分含航磁)数据；研究、应用概率论数理统计学原理，解决连片(重复)数据的标准化方法与实现；完成了检索的库原始数据、查询计算的标准化数据或比值计算数据等直接网格化的程序设计；探索了图形数据与标量数据库的接口技术，开发了屏幕绘图及与相关软件接口程序，做到了图并茂。使全国航放数据库具有：灵活、可靠的数据录入／导入功能，方便、可视性强的查询(统计)功能，实用、快捷的数据处理应用功能，适时、完备的输出功能，易用、可扩展的系统维护功能，严密、安全的数据保密功能。     

Application of LA ICP-MS analysis of arsenopyrite to gold metallogeny of the Meguma Terrane,Nova Scotia,Canada

Investigation of gold metallogeny in the Paleozoic Meguma terrane (Canada) is conducted through LA-ICP-MS analysis of arsenopyrite collected from eight slate-belt style vein gold deposits using a novel approach integrating elemental distribution maps and their derived elemental paragenesis with multi-element binary plots. The data reveal two distinct gold events: 1) an early event characterized by a Co-Ni-Mo-Sb-Se elemental association related to initial growth of arsenopyrite that reflects the presence of invisible gold (>10 ppm); and 2) a second event, spatially associated with late fracture sets, that is characterized by an Al-Ti-V-Mn element association and reflects either remobilization or upgrading of primary invisible gold and is manifest as visible gold.The results of this study indicate a complex and protracted history of gold mineralization which has important ramifications for the Meguma gold deposits, as well as other orogenic gold districts globally. In the case of the Meguma Terrane, it involves an initial gold event that is followed by element mobilization and, in the case of precious metals (Au, Ag), an upgrading through a zone refining process. In addition, the variable coupling and decoupling of elements is only revealed using in-situ derived LA-ICP-MS data.     

全球变化中的地球化学基线值研究及中国的研究现状

介绍了“全球地球化学制图”(IGM)项目及其在中国的进展情况。展望了我国地球化学基线值研究工作方向及研究重点。     

Tectonic terranes, metallogeny and regional geochemical surveys: an example from northern British Columbia

This study utilizes three major data sources: distribution of geological units; density, type, age and distribution of mineral deposits; and elemental analyses from regional geochemical stream sediment surveys to define parameters that ‘characterize’ tectonic terranes in northern British Columbia. A similar approach could be applied anywhere in the Canadian Cordillera.This area, NTS map sheets 104N, 104O and 104P along the British Columbia-Yukon border, forms a transect through allochthonous terranes into North American rocks. These are: the allochthonous island-arc Stikine, oceanic Cache Creek, cataclastic Yukon/Tanana, and island-arc Quensel terranes, the pericratonic Dorsey terrane; the parautochthonous oceanic Sylvester allochthon; and the autochthonous miogeoclinal North American Cassiar terrane. Plutonic rocks of Jurassic-Cretaceous to Tertiary age intrude all terranes.Data sources used in the study are geological base maps and reports, the Ministry of Energy Mines and Petroleum Resources' mineral deposit database (MINFILE) and analytical data from the National Regional Geochemical Survey stream sediment and water sampling program.Geological maps were compiled from various sources and plotted to act as bases for geochemical and mineral deposit overlays for analysis and interpretation.Geochemical samples were separated into background and anomalous populations and compared according to their source terranes. We found that mean concentrations from background sample populations for some elements are statistically distinctive for different terranes. Unfortunately, elemental correlation coefficients for the terranes are similar so cannot be used to characterize each terrane.Data on mineral deposits and occurrences were compiled from minfile and other sources. Particular attention was paid to deposits with histories of production or significant reported reserves. Deposits were sorted by type and commodity to produce synoptic metallogenic maps.The combined data from geological, geochemical and mineral deposit databases form a strong tool for interpreting and predicting patterns of mineralization.     

Seismic microzonation of Bangalore,India

In the present study, an attempt has been made to evaluate the seismic hazard considering local site effects by carrying out detailed geotechnical and geophysical site characterization in Bangalore, India to develop microzonation maps. An area of 220 km², encompassing Bangalore Mahanagara Palike (BMP) has been chosen as the study area. Seismic hazard analysis and microzonation of Bangalore are addressed in three parts: in the first part, estimation of seismic hazard is done using seismotectonic and geological information. Second part deals with site characterization using geotechnical and shallow geophysical techniques. In the last part, local site effects are assessed by carrying out one-dimensional (1-D) ground response analysis (using the program SHAKE2000) using both standard penetration test (SPT) data and shear wave velocity data from multichannel analysis of surface wave (MASW) survey. Further, field experiments using microtremor studies have also been carried out for evaluation of predominant frequency of the soil columns. The same has been assessed using 1-D ground response analysis and compared with microtremor results. Further, the Seed and Idriss simplified approach has been adopted to evaluate the soil liquefaction susceptibility and liquefaction resistance assessment. Microzonation maps have been prepared with a scale of 1:20,000. The detailed methodology, along with experimental details, collated data, results and maps are presented in this paper.     

Regional distribution of Al,B, Ba,Ca, K,La, Mg,Mn, Na,P, Rb,Si, Sr,Th, U and Y in terrestrial moss within a 188,000 km2 area of the central Barents region: influence of geology,seaspray and human activity

Five hundred and ninety-eight samples of terrestrial moss (Hylocomium splendens and Pleurozium schreberi) collected from a 188,000 km² area of the central Barents region (NE Norway, N Finland, NW Russia) were analysed by ICP-AES and ICP-MS. Analytical results for Al, B, Ba, Ca, K, La, Mg, Mn, Na, P, Rb, Si, Sr, Th, U and Y concentrations are reported here. Graphical methods of data analysis, such as geochemical maps, cumulative frequency diagrams, boxplots and scatterplots, are used to interpret the origin of the patterns for these elements. None of the elements reported here are emitted in significant amounts from the smelting industry on the Kola Peninsula. Despite the conventional view that moss chemistry reflects atmospheric element input, the nature of the underlying mineral substrate (regolith or bedrock) is found to have a considerable influence on moss composition for several elements. This influence of the chemistry of the mineral substrate can take place in a variety of ways. (1) It can be completely natural, reflecting the ability of higher plants to take up elements from deep soil horizons and shed them with litterfall onto the surface. (2) It can result from naturally increased soil dust input where vegetation is scarce due to harsh climatic conditions for instance. Alternatively, substrate influence can be enhanced by human activity, such as open-cast mining, creation of ‘technogenic deserts’, or handling, transport and storage of ore and ore products, all of which magnify the natural elemental flux from bedrock to ground vegetation. Seaspray is another natural process affecting moss composition in the area (Mg, Na), and this is most visible in the Norwegian part of the study area. Presence or absence of some plant species, e.g., lichens, seems to influence moss chemistry. This is shown by the low concentrations of B or K in moss on the Finnish and Norwegian side of the (fenced) border with Russia, contrasting with high concentrations on the other side (intensive reindeer husbandry west of the border has selectively depleted the lichen population).     

The application of geochemical pattern recognition to regional prospecting: A case study of the Sanandaj–Sirjan metallogenic zone,Iran

In regional exploration programs, the distribution of elements in known mineral deposits can be used as a guide for the classification of deposits, search for new prospects and modeling ore deposit patterns. The Sanandaj–Sirjan Zone (SSZ) is a major metallogenic zone in Iran, containing lead and zinc, iron, gold, copper deposits. In the central part of the SSZ, lead and zinc mineralization is widespread and hitherto exploration has been based on geological criteria. In this study, we used clustering techniques applied to element distribution for classification lead and zinc deposits in the central part of the SSZ. The hierarchical clustering technique was used to characterize the elemental pattern. Elements associated with lead and zinc deposits were separated into four clusters, encompassing both ore elements and their host rock-forming elements. It is shown that lead and zinc deposits in the central SSZ belong to two genetic groups: a MVT type hosted by limestone and dolomites and a SEDEX type hosted by shale, volcanic rocks and sandstone. The results of elemental clustering were used for pattern recognition by the K-means method and the respective deposits were classified into four distinct categories. K-means clustering also reveals that the elemental associations and spatial distribution of the lead and zinc deposits exhibit zoning in the central part of the SSZ. The ratios of ore-forming elements (Sb, Cd, and Zn) vs. (Pb and Ag) show zoning along an E–W trend, while host rock-forming elements (Mn, Ca, and Mg) vs. (Ba and Sr) show a zoning along a SE–NW trend. Large and medium deposits occur mainly in the center of the studied area, which justify further exploration around occurrences and abandoned mines in this area. The application of a pattern recognition method based on geochemical data from known mineralization in the central SSZ, and the classification derived from it, uncover elemental zoning, identify key elemental associations for further geochemical exploration and the potential to discover possible target areas for large to medium size ore deposits. This methodology can be applied in a similar way to search for new ore deposits in a wide range of known metallogenic zones.     

Lateral lithospheric conductivity patterns,from European synchronous observatory data

Much information on the regional lithospheric structure may come from MTS data acquired by synchronous 2D arrays and processed with regard to the nonlocal response of a laterally inhomogeneous subsurface. We suggest to invert the nonlocal MT responses applying correlation of all surface horizontal and vertical components of the geomagnetic field recorded simultaneously at all stations. The inversion algorithm has been applied to 2004–2005 European observatory data of diurnal S_q variations for first five harmonics and yielded lateral conductivity patterns for different periods. The maps show spatial correlation between conductivity maxima and lithospheric thickness minima and, specifically, highlight the contours of the Pannonian basin, where lithosphere is as thin as ～50 km, from seismic data.     

Combining in situ isotopic, trace element and textural analyses of quartz from four magmatic-hydrothermal ore deposits

This study couples in situ ¹⁶O, ¹⁷O and ¹⁸O isotope and in situ trace element analyses to investigate and characterize the geochemical and textural complexity of magmatic-hydrothermal quartz crystals. Euhedral quartz crystals contemporaneous with mineralization were obtained from four magmatic-hydrothermal ore deposits: El Indio Au–Ag–Cu deposit; Summitville Au–Ag–Cu deposit; North Parkes Cu–Au deposit and Kingsgate quartz-Mo–Bi–W deposit. The internal features of the crystals were imaged using cathodoluminescence and qualitative electron microprobe maps. Quantitative isotopic data were collected in situ using 157 nm laser ablation inductively coupled plasma mass spectrometry (for 40 trace elements in quartz) and sensitive high-resolution ion microprobe (for 3 isotopes in quartz). Imaging revealed fine oscillatory zoning, sector zoning, complex “macromosaic” textures and hidden xenocrystic cores. In situ oxygen isotope analyses revealed a δ¹⁸O range of up to 12.4 ± 0.3 ‰ in a single crystal—the largest isotopic range ever ascribed to oscillatory zonation in quartz. Some of these crystals contain a heavier δ¹⁸O signature than expected by existing models. While sector-zoned crystals exhibited strong trace element variations between faces, no evidence for anisotropic isotope fractionation was found. We found: (1) isotopic heterogeneity in hydrothermal quartz crystals is common and precludes provenance analysis (e.g., δD–δ¹⁸O) using bulk analytical techniques, (2) the trace element signature of quartz is not an effective pathfinder toward noble metal mineralization and (3) in three of the four samples, both textural and isotopic data indicate non-equilibrium deposition of quartz.     

Aeromagnetic data analysis of the Southern Aravalli Fold Belt: Its implications in understanding the inter-relationship among the migmatites and Gneissic Rocks,Aravalli supracrustals and Godhra granite

The Southern Aravalli Fold Belt represents a geologically complex terrain comprising migmatites, gneisses, low-grade metasediments and intrusive granites. Processing and interpretation of the existing aeromagnetic data using modern techniques has helped in understanding the geological complexity of the area. High quality aeromagnetic maps — total intensity, vertical derivative, analytic signal and 3D Euler depth — have been produced from the digitally processed data to augment visualisation and interpretation of the data. Strong east-west to northeast-southwest long strike-length anomalies in the southern part within the gneissic/granitic terrain indicate presence of narrow bands of high-susceptibility (magnetic) materials probably along regional lineaments. The thick sequence (> 1 km) of magnetically transparent sedimentary cover in the Rajgad basin is better visualised in the analytic signal and Euler Depth maps. The sediments in the northern part of the basin are deposited in a grabben formed by a northwest-southeast trending marginal fault in the north and an intrabasinal east-west fault in the south as demonstrated in the Euler 3D map. The Lunavada and Champaner Groups are separated by a shear zone comprising en-echelon shears clearly indicated in the shaded relief and vertical derivative maps. Interpretation of magnetic relief of the aeromagnetic data reveals that the area comprises of five distinct litho-tectonic units. This study demonstrates that the old aeromagnetic data acquired few decades back can provide valuable geological information, if processed using modern techniques.     

Actinide crystal-liquid partitioning for clinopyroxene and Ca3(PO4)2

Using fission and alpha track radiography techniques, we have measured partition coefficients (D) for the actinide elements Th, U and Pu between diopsidic clinopyroxene, whitlockite [β-Ca₃ (PO₄)₂] and silicate liquid at 20kbar. Equilibrium partitioning at the crystal-liquid interface is assumed, and corrections for actinide zoning have been applied to the measured D values. Reproducibility for both actinide and minor element D values is carefully examined as a criterion for crystal-liquid interface equilibrium. The data are mostly compatible with interface equilibrium except for experiments at high cooling rates ( ? 30 deg/hr). Partition coefficients for Th/U/Pu of about 0.002/0.002/0.06 are measured for clinopyroxene and 1.2/0.5/3.4 for whitlockite. At an oxygen fugacity of 10^?8.5, Pu is much more readily incorporated into the crystalline phases than is U or Th because of the importance of trivalent Pu. The D_Pu(cpx) is similar to D(cpx) of the light rare earths supporting the concept of Pu/(rare earth) dating.     

Landslide susceptibility maps comparing frequency ratio and artificial neural networks: a case study from the Nepal Himalaya

This study considers landslide susceptibility mapping by means of frequency ratio and artificial neural network approaches using geographic information system (GIS) techniques as a basic analysis tool. The selected study area was that of the Panchthar district, Nepal. GIS was used for the management and manipulation of spatial data. Landslide locations were identified from field survey and aerial photographic interpretation was used for location of lineaments. Ten factors in total are related to the occurrence of landslides. Based on the same set of factors, landslide susceptibility maps were produced from frequency ratio and neural network models, and were then compared and evaluated. The weights of each factor were determined using the back-propagation training method. Landslide susceptibility maps were produced from frequency ratio and neural network models, and they were then compared by means of their checking. The landslide location data were used for checking the results with the landslide susceptibility maps. The accuracy of the landslide susceptibility maps produced by the frequency ratio and neural networks is 82.21 and 78.25%, respectively.     

Applications of a GIS-based geotechnical tool to assess spatial earthquake hazards in an urban area

A geotechnical information system (GTIS) was constructed within a spatial geographic information system (GIS) framework to reliably predict geotechnical information and accurately estimate site effects at Gyeongju, an urban area in South Korea. The system was built based on both collected and performed site investigation data in addition to acquired geo-knowledge data. Seismic zoning maps were constructed using the site period (T _G) and mean shear wave velocity to a depth of 30 m (V _S30), and these maps were presented as a regional strategy to mitigate earthquake-induced risks in the study area. In particular, the T _G distribution map indicated the susceptibility to ground motion resonance in periods ranging from 0.2 to 0.5 s and the corresponding seismic vulnerability of buildings with two to five stories. Seismic zonation of site classification according to V _S30 values was also performed to determine the site amplification coefficients for seismic design and seismic performance evaluation at any site and administrative subunit in the study area. In addition, we investigated the site effects according to subsurface and surface ground irregularities at Gyeongju by seismic response analyses in time domains based on both two- and three-dimensional spatial finite element models, which were generated using spatial interface coordinates between geotechnical subsurface layers predicted by the GTIS. This practical study verified that spatial GIS-based geotechnical information can be a very useful resource in determining how to best mitigate seismic hazards, particularly in urban areas.     

A 2200-yr record of hydrologic variability from Foy Lake, Montana, USA, inferred from diatom and geochemical data

A 2200-yr long, high-resolution (∼5 yr) record of drought variability in northwest Montana is inferred from diatoms and δ¹⁸O values of bio-induced carbonate preserved in a varved lacustrine core from Foy Lake. A previously developed model of the diatom response to lake-level fluctuations is used to constrain estimates of paleolake levels derived from the diatom data. High-frequency (decadal) fluctuations in the de-trended δ¹⁸O record mirror variations in wet/dry cycles inferred from Banff tree-rings, demonstrating the sensitivity of the oxygen-isotope values to changes in regional moisture balance. Low frequency (multi-centennial) isotopic changes may be associated with shifts in the seasonal distribution of precipitation. From 200 B.C. to A.D. 800, both diatom and isotope records indicate that climate was dry and lake level low, with poor diatom preservation and high organic carbon: nitrogen ratios. Subsequently, lake level rose slightly, although the climate was drier and more stable than modern conditions. At A.D. 1200, lake level increased to approximately 6 m below present elevation, after which the lake fluctuated between this elevation and full stage, with particularly cool and/or wetter conditions after 1700. The hydrologic balance of the lake shifted abruptly at 1894 because of the establishment of a lumber mill at the lake's outlet. Spectral analysis of the δ¹⁸O data indicates that severe droughts occurred with multi-decadal (50 to 70 yr) frequency.