Self-Organizing Maps for identification of zeolitic diagenesis patterns in closed hydrologic systems on the Earth and its implications for Mars

We have conducted a survey of zeolite occurrences in saline-alkaline paleolake deposits on Earth to identify the most prominent zeolite alteration patterns and to characterize the most common authigenic minerals and their paragenetic relationships. We collected the bulk mineral assemblages (from previous and our studies) as identified by X-ray diffraction from zeolitic tuff beds and associated sedimentary beds from thirteen paleolake deposits from the USA, Mexico, Greece, and Tanzania. We applied the Kohonen Self-Organizing Maps (SOM) to look for interesting patterns in the tuff bed mineral assemblages without prescribing any specific interpretation, and for information reduction and classification. Decision Tree (DT) method was applied to characterize these clusters. We were able to define clear class boundaries between fresh glass, non-analcime zeolites, analcime, and K feldspar. The non-analcime zeolites were further grouped into several classes based on mineral type. We also discuss the potential implications for Mars, showing that the mineral assemblages of diagenetic facies identified by SOM and DT can be used to test or validate the orbital, in situ, or modeling results, while the trained SOM provides a robust generalized ability to classify the new mineral assemblage data into the most common diagenetic facies identified in saline-alkaline paleoenvironments that contain zeolites. The study concludes that generalizing the complex geochemical behaviors using unsupervised statistical learning methods can help to identify the most prominent geochemical behaviors.     

Unexploded ordnance discrimination using time-domain electromagnetic induction and self-organizing maps

Self-organizing maps (SOM) are implemented for discrimination of geologic noise, buried metal objects and unexploded ordnance using the geophysical method of time-domain electromagnetic induction. The learning and misfit measures are based on a Euclidean metric. The U*-matrix method is shown to be a reliable tool for determining data clusters and cluster boundaries. The performance of SOM for data-type discrimination was tested using three synthetic, idealized geophysical datasets consisting of exponential, multi-exponential and stretched-exponential decaying transients. In addition, experimental data were acquired using a modified Geonics EM63 instrument. Results from the synthetic examples show that SOM clusters the data based on their functional origin, when represented using U*-matrices. The percentage of correct classification is 100%. Unsupervised learning using the field dataset obtained with the Geonics EM63 succeeded in producing a multi-clustered map in which the background transients cluster themselves and are separated from clusters associated with metal clutter objects and UXO. Even though in some cases the SOM did not produce a single cluster for each type of causative body, it was able to separate clutter data from target data by producing several small clusters. The results are encouraging in view of the heterogeneity and sparsity of the training dataset.     

甘肃省地震应急专题图的设计与实现

地震应急期间,标准清晰的专题图因其快速、直观、信息量大等优点,能及时有效地提供灾区基本情况和灾情信息,在地震灾害应急指挥及救援管理中发挥着重要作用。地震应急图件的制作流程涉及基础数据的处理、总体设计、符号设计、制图综合、打印出图等多个环节,甘肃省地震应急专题图存在应急图件制作不规范、产品产出时效性不高的问题。通过不同比例尺数据库的建立、地震应急专题图和产品目录的设计、专题图模板的制作、本地化快速出图软件的部署,使震后应急期间各类专题图件产出快速、标准化、产品化,使甘肃省地震应急快速制图的能力和水平得到提升。     

Across-arc geochemical variation of Quaternary lavas in West Java, Indonesia: Mass-balance elucidation using arc basalt simulator model

The Sunda Arc of Indonesia developed along the convergent margin between the Eurasian and the Australian Plates. More than 100 Quaternary volcanic centers occur along the arc. The West Java Arc is a segment of the Sunda Arc in which more than 10 volcanic centers are located, corresponding to the 120 to 200 km depth contours of the Wadati–Benioff zone. The geochemistry of 207 Quaternary lavas from six centers across the arc was investigated. The lavas range from basalt to dacite. Incompatible element abundances increase from the volcanic front to the rear‐arc in response to a change from low‐K to high‐K suites. Nd–Sr isotope compositions of the basalts scatter between mid‐ocean ridge basalt (MORB) source mantle and Indian Ocean sediment (SED) compositions, with volcanic front low‐K basalts having more radiogenic Nd than the rear‐arc basalts. It is suggested that mixing between slab‐derived fluids mainly from the SED and melt from MORB source mantle played a significant role in determining the geochemistry of the West Java basalts. Incompatible element patterns in primitive mantle normalized multi‐element plots are almost identical across the arc, except for greater inclination and weaker positive Sr spikes in the rear‐arc basalts. This suggests a lower degree of partial melting in the rear‐arc mantle, accompanied by change in SED fluid composition between the volcanic front and the rear‐arc. The latter is confirmed by fluid‐fluxed melting model calculations using multiple trace elements and Nd and Sr isotopes. All the West Java Arc lavas require deficit of Sr from the slab SED. This may occur due to selective breakdown of Sr‐rich hydrous silicate minerals, such as zoisite, at shallower depths before the SED component reaches the depth of dehydration effective for magma genesis. The rear‐arc basalts need further Sr deficits along with lesser fluid. These features are commonly observed in many arc basalts, and are likely attributable to the same mechanism.     

Geochemical magma type discrimination: application to altered and metamorphosed basic igneous rocks

Five minor and trace elements, known to be chemically stable during alteration and metamorphism, have been combined in a set of binary diagrams that distinguish fresh tholeiites from alkali basalts. Of the five elements: Ti, P, Zr, Y, Nb, only P shows slight mobility during metamorphism, which is not sufficient to alter greatly the point distribution on the binary diagrams. Using these stable elements altered basaltic rocks: greenstones, spilites and amphibolites may be distinguished in the same way as fresh basalts, and their original magma may be identified as tholeiitic or alkaline basalt. All five elements are readily and rapidly determined, using XRF, thus this method may be applied as a rapid, easy way of discriminating the magma types of altered basaltic rocks. Using this method it can be demonstrated that alkali basalt magma was produced in minor quantities in the Precambrian.     

几种全岩样品40Ar-39Ar年龄谱的探讨

探讨了玄武岩、辉绿岩和断层泥中粘土矿物3种全岩样品在40Ar-39Ar分析中的某些问题。利用快中子照射过程中不同矿物发生不均匀Ar丟失的模式解释了玄武岩40Ar-39Ar全熔年龄与K-Ar年龄间的差异。对辉绿岩样品的40Ar-39Ar分析确定了辉绿岩脉的侵入时代,并说明该地区其它辉绿岩样品K-Ar年龄的离散是受后期变质作用的影响。对于断层泥中粘土矿物,快中子照射中39Ar的反冲丢失导致<1μm粒级部分的40Ar-39Ar全熔年龄偏高。<1μm和2～10μm粒级样品的40Ar-39Ar年龄谱中,低温区较低的阶段年龄可能是断层后期较强烈活动的结果。     

Application of integrated back‐propagation network and self organizing map for flood forecasting

Much of the nonlinearity and uncertainty regarding the flood process is because hydrologic data required for estimation are often tremendously difficult to obtain. This study employed a back‐propagation network (BPN) as the main structure in flood forecasting to learn and to demonstrate the sophisticated nonlinear mapping relationship. However, a deterministic BPN model implies high uncertainty and poor consistency for verification work even when the learning performance is satisfactory for flood forecasting. Therefore, a novel procedure was proposed in this investigation which integrates linear transfer function (LTF) and self‐organizing map (SOM) to efficiently determine the intervals of weights and biases of a flood forecasting neural network to avoid the above problems. A SOM network with classification ability was applied to the solutions and parameters of the BPN model in the learning stage, to classify the network parameter rules and to obtain the winning parameters. The outcomes from the previous stage were then used as the ranges of the parameters in the recall stage. Finally, a case study was carried out in Wu‐Shi basin to demonstrate the effectiveness of the proposal. Copyright © 2009 John Wiley & Sons, Ltd.     

Improvement of HVSR technique by self-organizing map (SOM) analysis

Data interpretation is one of the most important and thorny tasks in geosciences. Difficulties occur especially in non-invasive geophysical techniques and/or when the data that have to be analyzed are multidimensional, non-linear and highly noisy. Another important task is to ensure an efficient automatic data analysis, in order to allow a data interpretation as independent as possible from any a priori knowledge. This paper describes the post-processing application of a kind of neural network (self-organizing map, SOM) to the identification of the fundamental HVSR frequency of a given site. SOM results can be represented as two-dimensional maps, with a non-parametric mapping that projects the high dimensional original dataset in a fashion that provides both an unsupervised clustering and a highly visual representation of the data relationships. This innovative application of the SOM algorithm is presented with a case study related to the characterization of a mineral deposit.     

Basalt geochemistry and tectonic discrimination within continental flood basalt provinces

Continental flood basalts are usually regarded as a single tectonomagmatic entity but frequently quoted examples exhibit a variety of tectonic settings. In one well-studied, classic, flood basalt province, the Mesozoic Karoo province of southern Africa, magmatism occurred in the following tectonic settings: (a) continental rifting leading to ocean-floor spreading in the South Atlantic Ocean (Etendeka suite of Namibia); (b) stretched continental lithosphere and rifting not leading directly to ocean-floor formation (Lebombo suite of southeastern Africa); and (c) an a-tectonic, within-plate, continental setting characterized by an absence of faulting or warping (Lesotho highlands and Karoo dolerites of South Africa). By means of spidergrams of the elements Rb, Ba, Th, Nb, K, La, Ce, Sr, Nd, P, Hf, Zr, Sm, Ti, Tb, Y, V, Ni and Cr, uncontaminated tholeiites from (c) above [i.e. the Lesotho-type continental flood basalts (LTCFB)] are compared with mid-ocean ridge basalts (MORB), ocean-island tholeiites (OIT), and tholeiites and calc-alkali basalts from subduction environments. The comparison reveals the LTCFBs are geochemically distinct. The differences are reflected in relative enrichments or depletions of the more incompatible elements (Rb-Ce) to less incompatible elements (Ce-Y), i.e. the overall slope of the spidergrams, and in anomalous enrichments or depletions of one or more of the elements Th, K, Nb, Sr, Ti, Hf, and Zr. The distinctive geochemical character of the Lesotho LTCFBs is interpreted in terms of a lithospheric mantle source for the basalts. This is supported by isotopic data. There are no major geochemical differences between Lesotho CFBs and basalts of the rift-related Etendeka and Lebombo suites, although the latter are somewhat enriched in Rb, Ba and K. However, unlike the Lesotho basalts, the Lebombo and Etendeka basalts are associated with voluminous silicic volcanics or intrusive centres and late-stage dolerites having MORB/OIT (i.e. asthenospheric) geochemical characteristics. The flood basalt/silicic magmatism/late-stage dyke swarm association is characteristic of several rift or thinned lithosphere environments (e.g., Ethiopia, Skye, eastern Greenland) but in many of these the flood basalts have ocean-island basalt (OIT) geochemical characteristics. The Lesotho-type CFB geochemistry is exhibited by the Grande Ronde Basalt of the Columbia River Group (a possible subduction-related flood basalt province) and the basic rocks associated with Mesozoic rifting in the North and South Atlantic. Basalt geochemistry alone is unhelpful in determining the tectonic setting of CFBs although the rift-related environments may be identified by the petrology and geochemistry of the whole igneous suite. A two-source model is proposed for the mantle-derived basic rocks in rift-related CFB provinces. Early enriched basalts are derived from the lithosphere and, following pronounced lithospheric attenuation or rifting, later MORB-like melts are emplaced from the rising asthenosphere. The presence of both Lesotho- and OIT-type geochemical patterns in rift-related CFBs suggests that the lithosphere exhibits different styles of enrichment.     

Identification and classification of multiple reflections with self-organizing maps

Artificial neural networks can be used effectively to identify and classify multiple events in a seismic data set. We use a specialized neural network, a self-organizing map (SOM), that tries to establish rules for the characterization of the physical problem. Selected seismic data attributes from CMP gathers are used as input patterns, such that the SOM arranges the data to form clusters in an abstract space. We show with synthetic and real data how the SOM can identify and classify primaries and multiples, and how it can classify the various types of multiple corresponding to a certain generating mechanism in the subsurface.     

Surface current patterns observed by HF radar: methodology and analysis of currents to the north of the Yaeyama Islands,East China Sea

A new method was developed to compare the classifications of different kind of data maps based on the self-organizing map (SOM) analysis. The surface current maps of the northern coast of Ishigaki and Iriomote Islands, East China Sea, observed by high-frequency (HF) radar were classified by SOM. Winds, sea surface temperatures (SST), and reanalysis data were also classified by SOM. The optimum area for classification was determined objectively by the new method, which relates these patterns to HF radar current patterns. We found two typical surface patterns: the first was that northeastward flows are dominant in the observation area and the second was that a clockwise eddy was also dominant in this area. The southwestward wind pattern was strongly related to the dominant clockwise eddy pattern of the HF radar current field.     

基于IGRF的地磁基准图技术

地磁导航需要获得目标区域的地磁基准图。研究了基于国际地磁参考场的地磁场要素计算和等值线绘制方法。以给定目标区域的经纬高为输入,计算高阶地磁球谐函数获取地磁矢量并计算地磁场要素,修正了因地球扁率带来的数据差异。采用具有虚节点的移动矩形网格算法来绘制地磁参数等值线,生成目标区域的地磁基准图。采用C#实现了地磁基准图软件,测试表明,获得的地磁等高线数据满足使用要求,等值线精度较高。     

Trace elements in ocean ridge basalts

A correlary of sea floor spreading is that the production rate of ocean ridge basalts exceeds that of all other volcanic rocks on the earth combined. Basalts of the ocean ridges bring with them a continuous record in space and time of the chemical characteristics of the underlying mantle. The chemical record is once removed, due to chemical fractionation during partial melting. Chemical fractionations can be evaluated by assuming that peridotite melting has proceeded to an olivine-orthopyroxene stage, in which case the ratios of a number of magmaphile elements in the extracted melt closely match the ratios in the mantle. Comparison of ocean ridge basalts and chondritic meteorites reveals systematic patterns of element fractionation, and what is probably a double depletion in some elements. The first depletion is in volatile elements and is due to high accretion temperatures of a large percentage of the earth from the solar nebula. The second depletion is in the largest, most highly charged lithophile elements (“incompatible elements”), probably because the mantle source of the basalts was melted previously, and the melt, enriched in these elements, was removed. Migration of melt relative to solid under ocean ridges and oceanic plates, element fractionation at subduction zones, and fractional melting of amphibolite in the Precambrian are possible mechanisms for depleting the mantle in incompatible elements. Ratios of transition metals in the mantle source of ocean ridge basalts are close to chondritic, and contrast to the extreme depletion of refractory siderophile elements, the reason for which remains uncertain. Variation of ocean ridge basalt chemistry along the length of the ridge has been correlated with ridge elevation. Thus chemically anomalous ridge segments up to 1000 km long appear to broadly coincide with regions of high magma production (plumes, hot spots). Basalt heterogeneity at a single location indicates mantle heterogeneity on a smaller scale. Variation of ocean ridge basalt chemistry with time has not been established, in fact, criteria for recognizing old oceanic crust in ophiolite terrains are currently under debate. The similarity of rare earth element patterns in basalt from ocean ridges, back-arc basins, some young island arcs, and some continental flood basalts illustrates the dangers of tectonic labeling by rare earth element pattern.     

Olivine basalts at the Karymskii Volcanic Center: Mineralogy, petrogenesis, and magma sources

The olivine basalts of the Karymskii Volcanic Center (KVC) can be traced during the history of the area from the Lower Pleistocene until recently (the 1996 events); they are typical low-and moderate-potassium tholeiite basalts of the geochemical island-arc type. We have investigated the compositions of phenocryst minerals represented by plagioclase, olivine, clinopyroxene, as well as solid-phase inclusions of spinel in olivine, and more rarely in anorthite. The evolutionary trends of the rock-forming minerals provide evidence of the comagmaticity of these basalts, and thus of a long-lived intermediate magma chamber in the interior of the structure. The activity of this chamber is related to periodic transport of high temperature basalt melts to the surface. The geochemistry of the basalts is controlled by their origin at the same depleted magma source close to N-MORB, by successive crystallization of the primary melt, and by restricted mixing with magma components that are crystallizing at different depths. It is hypothesized that the solid-phase inclusions of high alumina spinel (hercynite?) found in olivine (and anorthite) of the basalts in the KVC north sector are of relict origin.     

Geochemistry of peridotite xenoliths in alkali basalts from Jeju Island, Korea

Abstract   Ultramafic xenoliths in alkali basalts from Jeju Island, Korea, are mostly spinel lherzolites with subordinate amounts of spinel harzburgites and pyroxenites. The compositions of major oxides and compatible to moderately incompatible elements of the Jeju peridotite xenoliths suggest that they are residues after various extents of melting. The estimated degrees of partial melting from compositionally homogeneous and unfractionated mantle to form the residual xenoliths reach 30%. However, their complex patterns of chondrite-normalized rare earth element, from light rare earth element (LREE)-depleted through spoon-shaped to LREE-enriched, reflect an additional process. Metasomatism by a small amount of melt/fluid enriched in LREE followed the former melt removal, which resulted in the enrichment of the incompatible trace elements. Sr and Nd isotopic ratios of the Jeju xenoliths display a wide scatter from depleted mid-oceanic ridge basalt (MORB)-like to near bulk-earth estimates along the MORB–oceanic island basalt (OIB) mantle array. The varieties in modal proportions of minerals, (La/Yb)_N ratio and Sr-Nd isotopes for the xenoliths demonstrate that the lithospheric mantle beneath Jeju Island is heterogeneous. The heterogeneity is a probable result of its long-term growth and enrichment history.     

Enforced self‐organizing map neural networks for river flood forecasting

Self‐organizing maps (SOMs) have been successfully accepted widely in science and engineering problems; not only are their results unbiased, but they can also be visualized. In this study, we propose an enforced SOM (ESOM) coupled with a linear regression output layer for flood forecasting. The ESOM re‐executes a few extra training patterns, e.g. the peak flow, as recycling input data increases the mapping space of peak flow in the topological structure of SOM, and the weighted sum of the extended output layer of the network improves the accuracy of forecasting peak flow. We have investigated an ESOM neural network by using the flood data of the Da‐Chia River, Taiwan, and evaluated its performance based on the results obtained from a commonly used back‐propagation neural network. The results demonstrate that the ESOM neural network has great efficiency for clustering, especially for the peak flow, and super capability of modelling the flood forecast. The topology maps created from the ESOM are interesting and informative. Copyright © 2007 John Wiley & Sons, Ltd.     

A field sampling strategy for semivariogram inference of fractures in rock outcrops

The stochastic continuum (SC) representation is one common approach for simulating the effects of fracture heterogeneity in groundwater flow and transport models. These SC reservoir models are generally developed using geostatistical methods (e.g., kriging or sequential simulation) that rely on the model semivariogram to describe the spatial variability of each continuum. Although a number of strategies for sampling spatial distributions have been published in the literature, little attention has been paid to the optimization of sampling in resource- or access-limited environments. Here we present a strategy for estimating the minimum sample spacing needed to define the spatial distribution of fractures on a vertical outcrop of basalt, located in the Box Canyon, east Snake River Plain, Idaho. We used fracture maps of similar basalts from the published literature to test experimentally the effects of different sample spacings on the resulting semivariogram model. Our final field sampling strategy was based on the lowest sample density that reproduced the semivariogram of the exhaustively sampled fracture map. Application of the derived sampling strategy to an outcrop in our field area gave excellent results, and illustrates the utility of this type of sample optimization. The method will work for developing a sampling plan for any intensive property, provided prior information for a similar domain is available; for example, fracture maps or ortho-rectified photographs from analogous rock types could be used to plan for sampling of a fractured rock outcrop.     

Quantifying uncertainty in remotely sensed land cover maps

Remotely sensed land cover maps are increasingly used as inputs into environmental simulation models whose outputs inform decisions and policy-making. Risks associated with these decisions are dependent on model output uncertainty, which is in turn affected by the uncertainty of land cover inputs. This article presents a method of quantifying the uncertainty that results from potential mis-classification in remotely sensed land cover maps. In addition to quantifying uncertainty in the classification of individual pixels in the map, we also address the important case where land cover maps have been upscaled to a coarser grid to suit the users’ needs and are reported as proportions of land cover type. The approach is Bayesian and incorporates several layers of modelling but is straightforward to implement. First, we incorporate data in the confusion matrix derived from an independent field survey, and discuss the appropriate way to model such data. Second, we account for spatial correlation in the true land cover map, using the remotely sensed map as a prior. Third, spatial correlation in the mis-classification characteristics is induced by modelling their variance. The result is that we are able to simulate posterior means and variances for individual sites and the entire map using a simple Monte Carlo algorithm. The method is applied to the Land Cover Map 2000 for the region of England and Wales, a map used as an input into a current dynamic carbon flux model.     

一种综合地学数据的矿产资源潜力制图方法

研究适应信息化时代特征的矿产资源潜力制图新技术、新方法对推动矿产资源评价理论与技术的发展具有重要的意义.笔者把GIS技术、图像分类算法和空间统计学理论进行有机集成,在空间统计学的空间结构分析技术和遥感图像纹理分类算法的基础上,提出了一种以综合地学数据（地质、地球物理、地球化学和遥感图像数据等）为基本数据源的矿产资源潜力自动制图方法.该方法的技术流程为：①数据准备,即对地球物理和地球化学勘探数据进行预处理,生成一个物化遥综合图像文件;②图像空间结构性分析和纹理图像生成,以综合地学图像为研究对象,用空间统计学的结构分析技术研究地学数据综合图像的空间结构性,生成纹理图像;③纹理图像多元分类,用实验变差函数纹理分类方法对研究区进行多元分类,生成分类专题图;④分类后处理,用叠置分析修正空间分类结果,生成区域矿产资源潜力分布图.     

Uranium in marine basalts: Concentration,distribution and implications

The uranium content of glass from chilled margins of oceanic tholeiitic basalt flows is generally ≤0.1ppm, even for old samples with highly altered crystalline interiors. Such low values represent the original whole rock concentrations, although subsequent to eruption low-temperature weathering has added uranium, and other elements, to the crystalline portions of these basalts. Consideration of the K/U ratios of altered samples suggests that basalt weathering may provide the major oceanic sink for these two elements.