Using Remote Sensing Data to Delineate the Lineaments for Hydrothermal Mineral Prediction in Heqing Area, Northwest Yunnan Province, China

Lineament extraction and analysis is one of the routine work in mapping medium and large areas using remote sensing data, most of which are satellite images. Landsat Enhanced Thematic Mapper （ETM） of 945×1 232 pixels subscene acquired on 21 March 2000 covering the northwestern part of Yunnan Province has been digitally processed using ER Mapper software. This article aims to produce lineament density map that predicts favorable zones for hydrothermal mineral occurrences and quantify spatial associations between the known hydrothermal mineral deposits. In the process of lineament extraction a number of image processing techniques were applied. The extracted lineaments were imported into MapGIS software and a suitable grid of 100 m×100 m was chosen. The Kriging method was used to create the lineament density map of the area. The results show that remote sensing data could be useful to extract the lineaments in the area. These lineaments are closely correlated with the faults obtained through other geological investigation methods. On comparing with field data the lineament-density map identifies two important high prospective zones, where large-scale deposits are already existing. In addition the map highlights unrecognized target areas that require follow up investigation.     

Digital Geological Mapping of Sinus Iridum Area of the Moon Based on the Chang’E-I Data

Based on the CCD images,IIM data and DEM data of China＇s lunar exploration project （Chang＇E-Ⅰ） and related processed and analytic results,an integrated study of regional geology of Sinus Iridum and its adjacent area was conducted,and a series of relevant researches and analyses were carried out,including analysis of impact craters and their extrusive and accumulative materials,division of stratigraphic and tectonic units and classification of rock types,integrated analysis of chronology and lunar evolution history.In consideration of crater＇s shape features,quantity and preserving status of filling materials,the lunar impact craters can be divided into 7 types and 11 subtypes,and the accumulative materials of craters are divided into 6 types and 9 accumulative groups.According to the content and distribution of TiO2 and image characteristics,the basalts are divided into high-TiO2,medium-TiO2 and low-TiO2 basalts.Discussion was made on division of tectonic units and evolution features in the study area.The geological map of the Sinus Iridum Quadrangle （LQ-4） at a scale of 1∶2.5 M was preliminarily compiled with the ArcGIS system,and the spatial database of the map was established.Related technical specification,procedure and method for lunar geological mapping have been worked out,so as to lay a foundation for the forthcoming geological mapping of the global Moon in China by using the data of Chang＇E-Ⅱ and also for comprehensive study and geological mapping of other celestial bodies in the future.     

Edge Effect Correction in the S-A Method for Geochemical Anomaly Separation

Anomaly separation using geochemical data often involves operations in the frequency domain, such as filtering and reducing noise/signal ratios. Unfortunately, the abrupt edge truncation of an image along edges and holes (with missing data) often causes frequency distribution distortion in the frequency domain. For example, bright strips are commonly seen in frequency distribution when using a Fourier transform. Such edge effect distortion may affect information extraction results; sometimes severely, depending on the edge abruptness of the image. Traditionally, edge effects are reduced by smoothing the image boundary prior to applying a Fourier transform. Zero-padding is one of the most commonly used smoothing methods. This simple method can reduce the edge effect to some degree but still distorts the image in some cases. Moreover, due to the complexity of geoscience images, which can include irregular shapes and holes with missing data, zero-padding does not always give satisfactory results. This paper proposes the use of decay functions to handle edge effects when extracting information from geoscience images. As an application, this method has been used in a newly developed multifractal method (S-A) for separating geochemical anomalies from background patterns. A geochemical dataset chosen from a mineral district in Nova Scotia, Canada was used to validate the method.     

Vertical Electrical Sounding Investigation in East River Nile Area （Khartoum State）, Sudan

Vertical electrical sounding (VES) was carried out in the east of the River Nile. The main objectives of the resistivity survey are to determine the types and thicknesses of sedimentary units in the area, to define the contact separating the sediments from the crystalline basement complex, and to determine the structural features of the subsurface formations. Several local depressions, whose maximum depth to the Imsement surface is about 160 m, are revealed as an outcome of the VES method, and suggested to have been infilled with undifferentiated units of the Nubian Group in particular Omdurman Formation. Thus, a depth to the basement complex is calculated and the associated structural map of the east of the River Nile is drawn. The map is useful for the groundwater drilling, as far as the presence or absence of the aquifer is concerned.     

Migration of Magnetotelluric Data in Two——Dimensional Model

Since the wave equation of magnetotelluric(MT)field is similar to the one of seis-mic,the migration techniques used in seismic can be applied to MT data.In thispaper we make use of the principle of reflector mapping(i.e.U/D imaging principle)to image MT data.That is,the MT wavefield observed on the surface of the earth canbe resolved into upgoing and downgoing waves,the waves are extrapolated downward bythe phase-shift method or the phase-shift plus interpolation(PSPI)method.Conductivi-ty interfaces of the medium could be found by using the time coincidence of the upgoingand downgoing waves.Theoretical calculations show that the migration technique of MTdata presented here is very effective.It can not only enhance the lateral resolution of MTdata,but also obtain the visual image of subsurface interfaces.As compared with theconventional 2-D inversion,this procedure is more simple in calculation and can be easi-ly put into practice on a personal computer and is able to obtain the MT depth section,which is     

Solution of Multiple——Point Statistics to Extracting Information from Remotely Sensed Imagery

Two phenomena of similar objects with different spectra and different objects with similar spectrum often result in the difficulty of separation and identification of all types of geographical objects only using spectral information．Therefore,there is a need to incorporate spatial structural and spatial association properties of the surfaces of objects into image processing to improve the accuracy of classification of remotely sensed imagery．In the current article,a new method is proposed on the basis of the principle of multiple-point statistics for combining spectral information and spatial information for image classification．The method was validated by applying to a case study on road extraction based on Landsat TM taken over the Chinese YeHow River delta on August 8,1999． The classification results have shown that this new method provides overall better results than the traditional methods such as maximum likelihood classifier (MLC)     

Deep learning based classification of rock structure of tunnel face

The automated interpretation of rock structure can improve the efficiency,accuracy,and consistency of the geological risk assessment of tunnel face.Because of the high uncertainties in the geological images as a result of different regional rock types,as well as in-situ conditions(e.g.,temperature,humidity,and construction procedure),previous automated methods have limited performance in classification of rock structure of tunnel face during construction.This paper presents a framework for classifying multiple rock structures based on the geological images of tunnel face using convolutional neural networks(CNN),namely Inception-ResNet-V2(IRV2).A prototype recognition system is implemented to classify 5 types of rock structures including mosaic,granular,layered,block,and fragmentation structures.The proposed IRV2 network is trained by over 35,000 out of 42,400 images extracted from over 150 sections of tunnel faces and tested by the remaining 7400 images.Furthermore,different hyperparameters of the CNN model are introduced to optimize the most efficient algorithm parameter.Among all the discussed models,i.e.,ResNet-50,ResNet-101,and Inception-v4,Inception-ResNet-V2 exhibits the best performance in terms of various indicators,such as precision,recall,F-score,and testing time per image.Meanwhile,the model trained by a large database can obtain the object features more comprehensively,leading to higher accuracy.Compared with the original image classification method,the sub-image method is closer to the reality considering both the accuracy and the perspective of error divergence.The experimental results reveal that the proposed method is optimal and efficient for automated classification of rock structure using the geological images of the tunnel face.     

http://www.sciencedirect.com/science/article/pii/S1674987112000400

Monitoring of regional groundwater levels provides important information for quantifying groundwater depletion and assessing impacts on the environment Historically,groundwater level monitoring wells in Beijing Plain,China,were installed for assessing groundwater resources and for monitoring the cone of depression.Monitoring wells are clustered around well fields and urban areas.There is urgent need to upgrade the existing monitoring wells to a regional groundwater level monitoring network to acquire information for integrated water resources management.A new method was proposed for designing a regional groundwater level monitoring network.The method is based on groundwater regime zone mapping.Groundwater regime zone map delineates distinct areas of possible different groundwater level variations and is useful for locating groundwater monitoring wells.This method was applied to Beijing Plain to upgrade a regional groundwater level monitoring network.     

结合全球气候模型预测结果概率分析融化深度

A probabilistic approach may be adopted to predict freeze and thaw depths to account for the variability of (1) material properties, and (2) contemporary and future surface energy input parameters (e.g.air temperatures, cloud cover, snow cover) predicted with global climate models. To illustratc the probabilistic approach, an example of the prediction of thaw depths in Fairbanks, Alaska, is considered, More specifically, the Stefan equation is used together with the Monte Carlo simulation technique to make a probabilistic prediction of thaw penetration. The simulation results indicate that the variability in material properties, surface energy input parameters, and temperature data can lead to significant uncertainty in predicting thaw penetration. The Taylor series method was performed to determine the mean and standard deviation of thaw penetration and the results were compared to the Monte Carlo simulation results. The close comparison of the results suggests that the simpler Taylor series method may be applied to many cold regions problems to account for the variability of input parameters.     

Discovering Crustal Deformation Bands by Processing Regional Gravity Field

Objectives: This article presents a new computational procedure to discover scratches buried in the earth's crust. We also validate this new interdisciplinary analysis method with regional gravity data located in a well-known Dabie orogenic zone for test. Methods: Based on the scratch analysis method evolved with mathematical morphology of surfaces, we present a procedure that extracts information of the crustal scratches from regional gravity data. Because the crustal scratches are positively and highly correlated to crustal deformation bands, it can be used for delineation of the crustal deformation belts. The scratches can be quantitatively characterized by calculation of the ridge coefficient function, whose high value traces delineate the deformation bands hidden in the regional gravity field. In addition, because the degree of crustal deformation is an important indicator of tectonic unit divisions, so the crust can be further divided according to the degree of crustal deformation into some tectonic units by using the ridge coefficient data, providing an objective base map for earth scientists to build tectonic models with quantitative evidence. Results: After the ridge coefficients are calculated, we can further enhance the boundary of high ridge-coefficient blocks, resulting in the so-called ridge-edge coefficient function. The high-value ridge-edge coefficients are well correlated with the edge faults of tectonic units underlay, providing accurate positioning of the base map for compilation of regional tectonic maps. In order to validate this new interdisciplinary analysis method, we select the Dabie orogenic zone as a pilot area for test, where rock outcrops are well exposed on the surface and detailed geological and geophysical surveys have been carried out. Tests show that the deformation bands and the tectonic units, which are conformed by tectonic scientists based on surface observations, are clearly displayed on the ridge and ridge-edge coefficient images obtained in this article. Moreover, these computer-generated images provide more accurate locations and geometric details. Conclusions: This work demonstrates that application of modern mathematical tools can promote the quantitative degree in research of modern geosciences, helping to open a door to develop a new branch of mathematical tectonics.     

A Study of Spatial Variability on Aquifer Hydraulic Conductivity K

A structural analysis of K of an aquifer system in the study area is presented, and the main direction and degree of the variability of K are found by using the unstationary regionalized variable theory of geostatistics. Optimal estimation of K has been made by universal kriging method (U K M ). Both spatial variability distribution map and division map of K are given.     

New Recognized Intellect for Prospecting Large‐superlarge Mineral Deposits

Based on the data base of 1285 mineral deposits of 22 commodities in 121 countries of 6 continents of the world, the authors use the linear trend analysis for their reserves to determine the cut-off limited order of reserves to select 36 exceptional superlarge (as peak mineral), 95 superlarge and 314 large deposits as new recognized intellect for their quantitative change. We have projected above 445 large-superlarge deposits on (1:5 M) global tectonic background map and divided 4 metallogenic domains, 21 metallogenic belts. Global metallogeny of large-superlarge deposits are: unity by endogenic, exogenic metamorphic and epigenetic in origin; speciality in different metallogenic domains and belts; preferentiality to ore-forming elements of Cu, Au, Fe, Ag, Cr, Mn, Zn, Pb, Sb, Hg, to continental margins or plate convergent belts, to Intra-continental tectono-magmatic complex belts and Large ductile shear zones; abnormality by the global oxyatmversion (excess oxygen atmospheric event) in Archean, redoxyatmversion (lack oxygen atmospheric event) in Proterozoic-Paleozoic, and tectonosphere thermal erosion (great amount of tectonic magmatic event) in Mesozoic-Cenozoic.     

A relational database of global U-Pb ages

Enhanced understanding of how sampling techniques affect estimates of the global U-Pb age-distribution have, in turn, constrained U-Pb database design. Recent studies indicate that each continent has a unique age-distribution, as determined by zircon ages dated by the U-Pb isotope method. Likewise, broad regions within a continent also exhibit diverse age-distributions. To achieve a reliable estimate of the global distribution, the heterogenous composition of the continental crust requires sampling as many regions as feasibly possible. To attain this goal, and to provide a method for calculating age histograms, the records from a recent global U-Pb compilation are supplemented with 281,631 new records. These additions increase the database size to 700,598 records. In addition, the data are now restructured and made available as a relational database. After filtering the records by the six age-models included with the database, the results reveal two problems that might generally be unrecognized. First, an abrupt switch in the best-age at any given point(such as 1000 Ma) from ~(206)Pb/~(238)U ages to ~(207)Pb/~(206)Pb ages artificially depresses the age-distribution at the cutoff point. Second, rejecting analyses based on either absolute discordance or the magnitude of 2σ precision errors artificially depresses the age-distribution between 900 Ma and 2000 Ma. The results indicate that, when estimating the global U-Pb age-distribution, the methods for determining best-age and for rejecting records both require some attention. Possible solutions include using either an Accuracy Model or a Precision Model for estimating best-age, and then including all U-Pb records in the estimate, rather than rejecting any of them.     

Technology and Its Application： MIS with Geographical Information

Despite the constant expansion of the application of GIS, in recent years, the data process-ing ability of the application system developed with GIS_Tool is quite limited. Management information system with the geographical information can not only be used to deal with data commonly, but also deal with digitized maps as GIS_Tool does. Thus it can help the user develop an application system high-qual-ified both in maps and traditional data. This paper, by taking the method of integrated MapInfo in the application of Visual Basic as an example, mainly illustrates two technological problenm-integrated map technology and the multiuserization of Maplnfo geographical information. And also the concrete solu-tions to the problems are provided.     

Application of Spectral Angle Mapper Classification to Discriminate Hydrothermal Alteration in Southwest Birjand,Iran, Using Advanced Spaceborne Thermal Emission and Reflection Radiometer Image Processing

The purpose of this study is to evaluate the Spectral Angle Mapper (SAM) classification method for determining the optimum threshold (maximum spectral angle) to unveil the hydrothermal mineral assemblages related to mineral deposits. The study area indicates good potential for Cu-Au porphyry, epithermal gold deposits and hydrothermal alteration well developed in arid and semiarid climates, which makes this region significant for Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) image processing analysis. Given that achieving an acceptable mineral mapping requires knowing the alteration patterns, petrochemistry and petrogenesis of the igneous rocks while considering the effect of weathering, overprinting of supergene alteration, overprinting of hypogene alteration and host rock spectral mixing, SAM classification was implemented for argillic, sericitic, propylitic, alunitization, silicification and iron oxide zones of six previously known mineral deposits: Maherabad, a Cu-Au porphyry system; Sheikhabad, an upper part of Cu-Au porphyry system; Khoonik, an Intrusion related Au system; Barmazid, a low sulfidation epithermal system; Khopik, a Cu-Au porphyry system; and Hanish, an epithermal Au system. Thus, the investigation showed that although the whole alteration zones are affected by mixing, it is also possible to produce a favorable hydrothermal mineral map by such complementary data as petrology, petrochemistry and alteration patterns.     

The impact of crustal rheology on natural seismicity: Campi Flegrei caldera case study

We analyze the crustal rheology beneath the active resurgent Campi Flegrei caldera(CFc) in Southern Italy by modelling the 3 D brittle-ductile(B/D) transition, based on available thermal, geological and geophysical data. Firstly, the thermal field in the conductive physical regime is modeled using a finite element method; based on an optimization tool, this method is applied to evaluate the location and dimensions of the deep thermal source beneath the caldera. A horizontally-extended thermal anomaly located at about 5000 m depth below sea level is identified beneath Pozzuoli Bay, a part of the CFc. The same isotherm is located at a depth of 20,000 m beyond the caldera. This indicates a higher horizontal temperature gradient in the caldera with respect to the surrounding area. Next, we utilize this thermal model to image the 3D rheological stratification of the shallow crust below the caldera with two different values of strain rates. Within the caldera, the B/D transitions with ε equal to 10^-12 s^-1 and 10^-8 s^-1 are located at 3000 m and 5000 m depths, respectively. Outside the caldera, the transition is very deep(15,000-20,000 m), seemingly uninfluenced by the thermal state of the CFc volcanism. Finally, we compare these results with the spatial distribution of earthquake hypocenters, Benioff strain release and b-value distribution to investigate the relationship between crustal rheology and seismicity characteristics. Our analysis reveals that the image of the B/D transition is in agreement with the distribution of earthquake hypocenters, constraining the potential seismogenic volume of the region. Our study demonstrates that knowledge of the rheological state of a volcanic system is an important element to interpret its dynamic, forecast future activity and improve evaluation of the associated seismic hazard.     

Illite Crystallinity Mapping of Very Low Grade Metamorphism of Triassic Metapelites in the Zoigê Area, Western China

X-ray diffraction methods for estimating the metamorphic grade of diagenetic, anchizone and epizone in metapelites are reviewed and applied to samples from a 7000?m+ borehole in western China and surface samples from the surrounding Zoigê area. Kübler’s illite crystallinity (IC) measurements provide more consistent results than calculated values of percentage of illite in the I/S mixed layers and percentage of I/S mixed layers. Down-borehole IC values display a typical burial metamorphic relationship between stratigraphic level and IC. A method for preparing very low grade metamorphic maps is described, and isograds plotted on a regional geological map at selected values of IC, delineating a high temperature diagenetic zone, an anchizone, and an epizone. The map shows that IC values are controlled by stratigraphic level in the north of the study area (i.e. burial metamorphism), and proximity to an igneous intrusive body in the south (i.e. contact metamorphism).     

An alternative interpretation for the map expression of “abrupt” changes in lateral stratigraphic level near transverse zones in fold-thrust belts

The map expression of "abrupt" changes in lateral stratigraphic level of a thrust fault has been traditionally interpreted to be a result of the presence of(1) a lateral(or oblique) thrust-ramp,or(2) a frontal ramp with displacement gradient、and/or(3) a combination of these geometries.These geometries have been used to interpret the structures near transverse /ones in fold-thrust belts(FTB).This contribution outlines an alternative explanation that can result in the same map pattern by lateral variations in stratigraphy along the strike of a low angle thrust fault.We describe the natural example of the Leamington transverse zone.which marks the southern margin of the Pennsy lvanian—Permian Oquirrh basin with genetically related lateral stratigraphic variations in the North American Sevier(TB.Thus.the observed map pattern at this zone is closely related to lateral stratigraphic variations along the strike of a horizontal fault.Even though the present-day erosional level shows the map pattern that could be interpreted as a lateral ramp.the observed structures along the Leamington zone most likely share the effects of the presence of a lateral(or oblique) ramp,lateral stratigraphic variations along the fault trace.and the displacement gradient.     

Decoding earth's plate tectonic history using sparse geochemical data

Accurately mapping plate boundary types and locations through time is essential for understanding the evolution of the plate-mantle system and the exchange of material between the solid Earth and surface environments.However,the complexity of the Earth system and the cryptic nature of the geological record make it difficult to discriminate tectonic environments through deep time.Here we present a new method for identifying tectonic paleo-environments on Earth through a data mining approach using global geochemical data.We first fingerprint a variety of present-day tectonic environments utilising up to 136 geochemical data attributes in any available combination.A total of 38301 geochemical analyses from basalts aged from 5-0 Ma together with a well-established plate reconstruction model are used to construct a suite of discriminatory models for the first order tectonic environments of subduction and mid-ocean ridge as distinct from intraplate hotspot oceanic environments,identifying 41,35,and 39 key discriminatory geochemical attributes,respectively.After training and validation,our model is applied to a global geochemical database of 1547 basalt samples of unknown tectonic origin aged between 1000-410 Ma,a relatively ill-constrained period of Earth’s evolution following the breakup of the Rodinia supercontinent,producing 56 unique global tectonic environment predictions throughout the Neoproterozoic and Early Paleozoic.Predictions are used to discriminate between three alternative published Rodinia configuration models,identifying the model demonstrating the closest spatio-temporal consistency with the basalt record,and emphasizing the importance of integrating geochemical data into plate reconstructions.Our approach offers an extensible framework for constructing full-plate,deeptime reconstructions capable of assimilating a broad range of geochemical and geological observations,enabling next generation Earth system models.     

Integrated model for earthquake risk assessment using neural network and analytic hierarchy process:Aceh province,Indonesia

Catastrophic natural hazards,such as earthquake,pose serious threats to properties and human lives in urban areas.Therefore,earthquake risk assessment(ERA)is indispensable in disaster management.ERA is an integration of the extent of probability and vulnerability of assets.This study develops an integrated model by using the artificial neural network–analytic hierarchy process(ANN–AHP)model for constructing the ERA map.The aim of the study is to quantify urban population risk that may be caused by impending earthquakes.The model is applied to the city of Banda Aceh in Indonesia,a seismically active zone of Aceh province frequently affected by devastating earthquakes.ANN is used for probability mapping,whereas AHP is used to assess urban vulnerability after the hazard map is created with the aid of earthquake intensity variation thematic layering.The risk map is subsequently created by combining the probability,hazard,and vulnerability maps.Then,the risk levels of various zones are obtained.The validation process reveals that the proposed model can map the earthquake probability based on historical events with an accuracy of 84%.Furthermore,results show that the central and southeastern regions of the city have moderate to very high risk classifications,whereas the other parts of the city fall under low to very low earthquake risk classifications.The findings of this research are useful for government agencies and decision makers,particularly in estimating risk dimensions in urban areas and for the future studies to project the preparedness strategies for Banda Aceh.