U-Pb Ages of the Late Early Cretaceous Magmatic Rocks in the Zhalantun Area of Inner Mongolia

A great deal of early-to-mid Early Cretaceous magmatic activities have been recorded in the Zhalantun area of Inner Mongolia,while the late Early to Late Cretaceous magmatic rocks have been barely reported(Guo et al.,2018;Zhang Xiangxin et al.,2017).At present,only a few Late Cretaceous magmatic activities were reported in the Arongqi area,such as volcanic rocks of the Gushanzhen Formation.However,the Gushanzhen Formation lacks accurate isotopic age,and contemporaneous intrusive rocks has not been reported yet.In this work,we collected the volcanic rocks from the Gushanzhen Formation and contemporaneous intrusive rocks in the Zhanlantun and nearby,and aim to figure out the formation ages of volcanic rocks of the Gushanzhen Formation and accompanied intrusive rocks by analyzing zircon U-Pb isotopes(Fig.1).     

Self-Similar and Self-Affine Properties of Two-Dimensional Fracture Patterns in Rocks

In this paper, we investigate the fractal properties of binary maps of rock fractures at different scales and different geological types, as well as different families of fracture patterns obtained from a two-dimensional Laplacian growth model (LGM). From these analyses we figure out which families of the LGM patterns match the structural properties of the fracture binary maps. The LGM is defined in terms of a nonlinear map that depends on two parameters, λ and \mathfraka\mathfrak{a}, that respectively define the area and shape of the elements of the aggregate that conforms the patterns. The fractal dimension and roughness exponent of the LGM patterns are found to depend on \mathfraka\mathfrak{a}, with 0 < \mathfrak a < 10<\mathfrak {a}<1. From a detailed statistical analysis of these patterns we found that the fractal dimensions of capacity, correlation and information decrease monotonically as \mathfraka\mathfrak{a} increases. We also found that the values of these three fractal dimensions tend to collapse on top of each other as \mathfraka\lessapprox1\mathfrak{a}\lessapprox1. Remarkably, the fractal properties of rock fractures in the scales from millimeters up to a few meters appear to be well represented by the fractal structure of the LGM families of patterns with \mathfraka=0.15\mathfrak{a}=0.15 and 0.30, while the fractal properties of rock fractures in the scale of kilometers seems to be well represented by the LGM family with \mathfraka=0.90\mathfrak{a}=0.90. In addition, the three fractal dimension values of fracture binary maps in the scales from millimeters up to meters were found to be different between them. Nonetheless, for fractures in the scale of kilometers, the values of the three fractal dimensions are very close to each other as an indication of self-similar behavior. Analysis of the corrections to the scaling of the roughness exponent, ζ, suggests that they are negligible for the LGM family of fracture patterns with \mathfraka=0.9\mathfrak{a}=0.9. This finding points to a self-affine structure for this family of patterns. In fact, the calculated roughness exponent results are in the range of values characteristic of rock fractures.     

New U-Pb Geochronological Constraints on Formation and Evolution of the Susong Complex Zone in the Dabie Orogen

<正>Objective The Susong complex zone(SCZ)is a relatively lowgrade metamorphic unit mostly with an epidoteamphibolite facies,located in the southernmost part of the Dabie orogen.However,its rock compositions,ages,metamorphic processes and petrogenesis are still     

Application of the EnKF and Localization to Automatic History Matching of Facies Distribution and Production Data

The performance of the ensemble Kalman filter (EnKF) for continuous updating of facies location and boundaries in a reservoir model based on production and facies data for a 3D synthetic problem is presented. The occurrence of the different facies types is treated as a random process and the initial distribution was obtained by truncating a bi-Gaussian random field. Because facies data are highly non-Gaussian, re-parameterization was necessary in order to use the EnKF algorithm for data assimilation; two Gaussian random fields are updated in lieu of the static facies parameters. The problem of history matching applied to facies is difficult due to (1) constraints to facies observations at wells are occasionally violated when productions data are assimilated; (2) excessive reduction of variance seems to be a bigger problem with facies than with Gaussian random permeability and porosity fields; and (3) the relationship between facies variables and data is so highly non-linear that the final facies field does not always honor early production data well. Consequently three issues are investigated in this work. Is it possible to iteratively enforce facies constraints when updates due to production data have caused them to be violated? Can localization of adjustments be used for facies to prevent collapse of the variance during the data-assimilation period? Is a forecast from the final state better than a forecast from time zero using the final parameter fields?To investigate these issues, a 3D reservoir simulation model is coupled with the EnKF technique for data assimilation. One approach to enforcing the facies constraint is continuous iteration on all available data, which may lead to inconsistent model states, incorrect weighting of the production data and incorrect adjustment of the state vector. A sequential EnKF where the dynamic and static data are assimilated sequentially is presented and this approach seems to have solved the highlighted problems above. When the ensemble size is small compared to the number of independent data, the localized adjustment of the state vector is a very important technique that may be used to mitigate loss of rank in the ensemble. Implementing a distance-based localization of the facies adjustment appears to mitigate the problem of variance deficiency in the ensembles by ensuring that sufficient variability in the ensemble is maintained throughout the data assimilation period. Finally, when data are assimilated without localization, the prediction results appear to be independent of the starting point. When localization is applied, it is better to predict from the start using the final parameter field rather than continue from the final state.     

Easy Path Wavelet Transform on Triangulations of the Sphere

The computation of sparse representations of data on the sphere (e.g. topographical data) is a crucial step for further processing such as multiple separation, migration, imaging and sparsity promoting data-recovery. The Easy Path Wavelet Transform (EPWT) is a new tool for sparse data representation that has recently been introduced for image compression. In this paper we consider the EPWT on spherical triangulations. It is a locally adaptive transform that works along pathways through the array of function values and exploits the local correlations of the data in a simple appropriate manner. In our approach the usual discrete one-dimensional orthogonal or biorthogonal wavelet transform can be applied. The EPWT can be used for defining a multiresolution analysis on the sphere in which the scaling spaces and the wavelet spaces depend adaptively on the data. Issues of implementation of the EPWT are also considered.     

Temperature Distribution in the Earth’s Mantle

Doklady Earth Sciences - The influence of the compressibility and sphericity of the mantle on the temperature distribution in various models of mantle convection is analyzed. A model providing a...     

The Use of Neural Networks for the Prediction of Zeta Potential of Kaolinite

The sign and the magnitude of the zeta potential must be known for many engineering applications. For clay soils, it is usually negative, but it is strongly dependent on the pore fluid chemistry. However, measurement of zeta potential time is time-consuming and requires special and expensive equipment. In this study, the prediction of zeta potential of kaolinite has been investigated by artificial neural networks (ANNs) and multiple regression analyses (MRAs). To achieve this, ANN and MRA models based on zeta potential measurements of kaolinite in the presence of salt and heavy metal cations at different pH values have been developed. The results of the models were compared with the experimental results. The performance indices, including coefficient of determination, root mean square error, mean absolute error, and variance, were used to assess the performance of the prediction capacity of the models developed in this study. The obtained indices make it clear that the constructed ANN models were able to predict zeta potential of kaolinite quite efficiently and outperformed the MRA models. Results showed that ANN models can be used satisfactorily to predict zeta potential of kaolinite as a rapid inexpensive substitute for laboratory techniques.     

Reaction rim growth in the system MgO-Al2O3-SiO2 under uniaxial stress

We synthesize reaction rims between thermodynamically incompatible phases in the system MgO-Al₂O₃-SiO₂ applying uniaxial load using a creep apparatus. Synthesis experiments are done in the MgO-SiO₂ and in the MgO-Al₂O₃ subsystems at temperatures ranging from 1150 to 1350 °C imposing vertical stresses of 1.2 to 29 MPa at ambient pressure and under a constant flow of dry argon. Single crystals of synthetic and natural quartz and forsterite, synthetic periclase and synthetic corundum polycrystals are used as starting materials. We produce enstatite rims at forsterite-quartz contacts, enstatite-forsterite double rims at periclase-quartz contacts and spinel rims at periclase-corundum contacts. We find that rim growth under the “dry” conditions of our experiments is sluggish compared to what has been found previously in nominally “dry” piston cylinder experiments. We further observe that the nature of starting material, synthetic or natural, has a major influence on rim growth rates, where natural samples are more reactive than synthetic ones. At a given temperature the effect of stress variation is larger than what is anticipated from the modification of the thermodynamic driving force for reaction due to the storage of elastic strain energy in the reactant phases. We speculate that this may be due to modification of the physical properties of the polycrystals that constitute the reaction rims or by deformation under the imposed load. In our experiments rim growth is very sluggish at forsterite-quartz interfaces. Rim growth is more rapid at periclase-quartz contacts. The spinel rims that are produced at periclase-corundum interfaces show parabolic growth indicating that reaction rim growth is essentially diffusion controlled. From the analysis of time series done in the MgO-Al₂O₃ subsystem we derive effective diffusivities for the Al₂O₃ and the MgO components in a spinel polycrystal as ${\rm D}_{MgO} = 1.4 \pm 0.2 \cdot 10^{-15}$  m²/s and ${\rm D}_{Al_2O_3} = 3.7 \pm 0.6 \cdot 10^{-16}$  m²/s for T?=?1350 °C and a vertical stress of 2.9 MPa.     

Two-dimensional Conditional Simulations Based on the Wavelet Decomposition of Training Images

Scale dependency is a critical topic when modeling spatial phenomena of complex geological patterns that interact at different spatial scales. A two-dimensional conditional simulation based on wavelet decomposition is proposed for simulating geological patterns at different scales. The method utilizes the wavelet transform of a training image to decompose it into wavelet coefficients at different scales, and then quantifies their spatial dependence. Joint simulation of the wavelet coefficients is used together with available hard and or soft conditioning data. The conditionally co-simulated wavelet coefficients are back-transformed generating a realization of the attribute under study. Realizations generated using the proposed method reproduce the conditioning data, the wavelet coefficients and their spatial dependence. Two examples using geological images as training images elucidate the different aspects of the method, including hard and soft conditioning, the ability to reproduce some non-linear features and scale dependencies of the training images.     

Reconstructing the event stratigraphy from the complex structural – stratigraphic architecture of an Archaean volcanic – intrusive – sedimentary succession: the Boorara Domain,Eastern Goldfields Superterrane,Western Australia

Two main deformational phases are recognised in the Archaean Boorara Domain of the Kalgoorlie Terrane, Eastern Goldfields Superterrane, Yilgarn Craton, Western Australia, primarily involving south-over-north thrust faulting that repeated and thickened the stratigraphy, followed by east-northeast – west-southwest shortening that resulted in macroscale folding of the greenstone lithologies. The domain preserves mid-greenschist facies metamorphic grade, with an increase to lower amphibolite metamorphic grade towards the north of the region. As a result of the deformation and metamorphism, individual stratigraphic horizons are difficult to trace continuously throughout the entire domain. Volcanological and sedimentological textures and structures, primary lithological contacts, petrography and geochemistry have been used to correlate lithofacies between fault-bounded structural blocks. The correlated stratigraphic sequence for the Boorara Domain comprises quartzo-feldspathic turbidite packages, overlain by high-Mg tholeiitic basalt (lower basalt), coherent and clastic dacite facies, intrusive and extrusive komatiite units, an overlying komatiitic basalt unit (upper basalt), and at the stratigraphic top of the sequence, volcaniclastic quartz-rich turbidites. Reconstruction of the stratigraphy and consideration of emplacement dynamics has allowed reconstruction of the emplacement history and setting of the preserved sequence. This involves a felsic, mafic and ultramafic magmatic system emplaced as high-level intrusions, with localised emergent volcanic centres, into a submarine basin in which active sedimentation was occurring.     

Porphyry Deposits of the South Urals: Rhenium Distribution

Please refer to the attachment(s) for more details     

On the Elimination of Bias Averaging-Errors in Proxy Records

Knowledge of and insight into past environmental conditions can be obtained by processing and analyzing proxies. The proxies need to be processed as precisely and accurately as possible, otherwise the conclusion of the analysis will be biased. A calibration method which reduces bias errors in the proxy measurements due to averaging is presented. Sampling with nonzero sample sizes causes an averaging of the true proxy signal over the volume of the sample. The method is applied on a linear synthetic record which results in an optimal correction for frequency components ranging from the dc-frequency (DC) to one half of the sample frequency (f _s /2). Next, the method is tested on non-linear synthetic data where the signal is reconstructed reasonably well. Finally, the method is applied to a real vessel density record of R. mucronata from Makongeni, Kenya, and to a real delta deuterium record of ice core EDC from dome C, Antarctica. The method discussed in this paper is a valuable tool for the calibration of proxy measurements; it can be applied as a correction for low resolution measurements and expanded to other types of samples and proxies. Working with small sample sizes (high resolution) amounts to working near the detection limit, where the signal-to-noise-ratio is low. This correction method provides an alternative in which low resolution measurements can be upgraded to minimize the loss of information due to larger sample sizes.     

Probabilistic Marcinkiewicz–Zygmund Inequalities on the Rotation Group

Approximation problems on the rotation group SO(3) naturally arise in various fields, like crystallography, chemistry, and biology. For example, in crystallographic texture analysis one is confronted with the problem of evaluating so-called orientation density functions (ODFs). In many situations one only has a finite number of measurements at scattered sampling nodes. In order to reconstruct ODFs over all rotations, so-called Marcinkiewicz–Zygmund inequalities on the rotation group are an important tool. These inequalities provide norm equivalences between polynomials on SO(3) and their sample values. Recently shown equivalences depend on a density parameter of the sampling set and the proven inequalities hold true for polynomials on SO(3) whose degree does not exceed an upper bound which is determined by this density parameter. In this paper, we show that we can enlarge this upper bound for the polynomial degree significantly if we are satisfied by such norm equivalences that hold with a given probability only. Moreover, we show that there are fixed sampling sets for which we get probabilistic Marcinkiewicz–Zygmund inequalities that hold for polynomials on SO(3) of all degrees.     

Another Look at the Chemical Relationships in the Dissolved Phase of Complex River Systems

Large rivers are a major pathway for the erosion products of continents to reach the oceans. The riverine transport of dissolved and particulate materials is generally related to a large number of interactions involving climate, hydrological, physico-chemical and biological aspects. Consequently, the investigation of large rivers allows the erosion processes at a global scale to be addressed, with information about biogeochemical cycles of the elements, weathering rates, physical erosion rates and CO₂ consumption by the acid degradation of continental rocks. Today, good databases exist for the major dissolved ions in the world’s largest rivers. Since concentration of ions in river waters has to be considered in a compositional context, it is necessary to study the implications of considering the simplex, with its proper geometry, as the natural sample space. Using the additive (alr) or the isometric (ilr) log-ratio transformations, a composition can be represented as a real vector; but only in the second case can these coordinates be mapped onto orthogonal axes. Using data related to the dissolved load of some of the most important rivers in the world, the relationships among the major ions frequently used in molar ratio mixing diagrams have been investigated with alternative tools. Following the balances approach, an investigation of the properties of aqueous solutions of electrolytes that are often treated in terms of equilibrium constants is undertaken. The role played by the source—rain water, weathering of silic, carbonatic and evaporitic rocks, pollution—from which elements and chemical species can potentially be derived, has been checked through an investigation of a probabilistic model able to describe the relationships among the different components that contribute to the chemical composition of a river water sample.     

Time-Series Reconstruction from Natural Archive Data with the Averaging Effect Taken into Account

Environmental information of the past can be obtained by processing and analyzing proxies recorded by environmental archives. Natural archives are sampled at a distance grid along their accretion axis. Starting from these distance series, a time series needs to be constructed, because comparison of different data records is only meaningful on a time grid. However, distance–time relationships are nonlinear as the accretion rate of natural archives is dependent on environmental and physiological factors. Furthermore, in environmental archives, samples are taken over a volume in distance, rather than over a point in distance. This implies that the sample-values will be averaged over the volume of the sample. In this paper a method is proposed, which establishes the nonlinear distance–time relationship and corrects for the averaging effects. The method is built upon the assumption that the proxy record on a time axis is harmonic. If this is not the case, then a harmonic approximation is made. As a consequence of the nonlinear distance–time relationship, this harmonic proxy signal is nonlinearly distorted on a distance axis. As such, a harmonic signal model with a nonlinear phase distortion and an averaging effect is fitted on the data. Since environmental records are short data records, the statistical performance of the estimator on noisy data is verified by means of Monte Carlo simulations. The applicability of the method is demonstrated on the measurement of the vessel density, in a mangrove tree, Rhizophora mucronata, which is an indicator of the rainfall in tropical coastal ecosystems.     

The Hilbert Transform on the Two-Sphere: A Spectral Characterization

The local analysis of signals arising on the sphere is a common task in earth sciences. On the real line the analytic signal turned out to be an important representation in local one-dimensional signal processing. Its generalization to two dimensions is the monogenic signal, and the properties of the analytic and the monogenic signal in the Fourier domain are well known. A generalization to the sphere is given by the Hilbert transform on the sphere known from Clifford analysis. To obtain a spectral characterization, the transform has to be decomposed into spherical harmonic functions. In this paper, we derive the spherical harmonic coefficients of the Hilbert transform on the sphere and give a series expansion. This will show that it acts as a differential operator on the spherical harmonic basis functions of the Laplace equation solution, analogously to the Riesz transform in two dimensions. This allows an interpretation of the Hilbert transform suitable for signal processing of signals naturally arising on the two-sphere. We show that the scale space naturally arising is a Poisson scale space in the unit ball. In addition, the obtained interpretation of the Hilbert transform is used for orientation analysis of plane waves. This representation is justified as a novel signal model on the sphere which can be used to construct intensity and rotation-invariant operators for local signal analysis in a scale-space concept.     

Basic Characteristics of the Earth''s Temperature Distribution in Southern China

Sourthern China referred to in this paper embraces all the provinces to the south of latitude 32° N within the territory of China, exclusive of Tibet and Taiwan.The paper is based on the data obtained from temperature measurements of 680 boreholes drilled for petroleum, coal and hydrogeology. The depth of these boreholes ranges between 1,000 and 2,000m in most cases, being generally not less than 400m and reaching 7,000m in maximum.Through the elaboration and analysis of the temperature data, isotherm maps showing the temperatures at depths of 1,000m, 2,000m and 3,000m respectively and a geothermal gradient map of southern China have been compiled. These maps indicate that the earth's temperature distribution in southern China is characterized by a elatively low value in areas of its central part and those to the west of the Anning River valley of the western part(30℃, 40℃ and 70℃ at depths of 1,000m, 2,000m and 3,000m respectively), and a comparatively high value in the southeastern coastal areas and western Yunnan(60℃, 80℃ and over 120℃ respectively), and that the geothermal gradient also increases from less than 1.5℃/100 m in the central part towards both the east and the west to 2.0-3.0℃/100m or even over 4.0℃/100m.The earth's temperature distribution is strictly controlled by the deep structure of the earth's crust and is very closely related to the character of the regional geological structure, The stable central areas remarkably differ in the earth's temperature distribution from the tectonicalty active areas such as the southeastern coastal areas and the Jinggu-Tengchong area in western Yunnan, the former being lower than the latter. Meanwhile, due to the well-developed karst features in southern China, the strong groundwater flow may exert certain cooling effect on the temperature down to a depth of about 1,500m.According to the regularity of the earth's temperature distribution, regional deep structure of the earth's crust, regional geological structure and hydrogeological conditions, the following four models are proposed to explain the mechanism of formation and distribution of the earth's temperature fields in southern China: 1) onduction type; 2) convection type; 3) deep heat-source conduction-convection type; and 4) precipitation infiltration cooling type.     

Structure of the Khatanga–Lomonosov Fracture Zone according to Seismic Data

Doklady Earth Sciences - Based on the results of geophysical data interpretation, the structural features of the Khatanga–Lomonosov fracture zone (KhLZ), adapted to the junction zone between...