Atmospheric propagational phase variations are the dominant source of error for InSAR timeseries analysis, generally exceeding uncertainties from poor SNR or signal correlation. The spatial properties of these errors have been well studied, but their temporal dependence and correction have received much less attention to date. We present here an evaluation of the magnitude of tropospheric artifacts in derived time series after compensation using an algorithm that requires only the InSAR data themselves. The level of artifact reduction equals or exceeds that from many weather model based methods, while avoiding the need to access fine-scale atmosphere parameters globally at all times. Our method consists of identifying all points in an InSAR stack with consistently high correlation, and computing, then removing, a fit of the phase at each of these points with respect to elevation. Comparison with GPS truth yields a reduction of 3, from an rms misfit of 5-6 cm to ~2 cm over time. This algorithm can be readily incorporated into InSAR processing flows without need for outside information.

The paper presents an analysis of the method of recording the magnetic component of the Earth’s natural pulse electromagnetic field in an urban environment. This method of recording has already proved itself to be a method that allocates the stressed sections of rock mass at mining and, therefore, authors suppose its effectiveness for allocating active tectonic disturbances and forecasting accidents at underground utilities, what will help reduce the potential environmental hazard of these objects.

The paper is dedicated to the topical problem of examining permafrost state and the processes of its geocryological changes by means of geophysical methods. To monitor the cryolithozone, we propose and scientifically substantiate a new technique of pulsed electromagnetic cross-well sounding. Based on the vector finite-element method, we created a mathematical model of the cross-well sounding process with a pulsed source in a three-dimensional spatially heterogeneous medium. A high-performance parallel computing algorithm was developed and verified. Through realistic geoelectric models of permafrost with a talik under a highway, constructed following the results of electrotomography field data interpretation, we numerically simulated the pulsed sounding on the computing resources of the Siberian Supercomputer Center of SB RAS. The simulation results suggest the proposed system of pulsed electromagnetic cross-well monitoring to be characterized by a high sensitivity to the presence and dimensions of the talik. The devised approach can be oriented to addressing a wide range of issues related to monitoring permafrost rocks under civil and industrial facilities, buildings and constructions.

1 Abstract The grain-to-grain stress vertically in sediments is given by the overburden less the pore fluid pressure, σ, divided by the fraction of the horizontal area which is the supporting matrix , (1 − φ), φ being the porosity. It is proposed that the fractional reduction of this ratio, Λ, with time is given by the product of φ 4m/3) , (1 − φ) 4n/3 , and one or more Arrhenius functions A exp(−E/RT ) with m and n close to 1. This proposal is tested for shale sections in six wells from around the world for which porosity-depth data are available. Good agreement is obtained above 30-40 C. A single activation energy of 23+-5 kJ/mole, indicating pressure solution of quartz, 24 kJ/mol, was obtained. The average value of m is 1, indicating fractal pore-matrix spaces and water-wet interfaces. Grain-to -grain interfaces may be fractal with m close to 1, but can have lower values suggesting smooth surfaces and even grain-to-grain welding. Results are independent of over- or under-pressure of pore water. This model explains shale compaction quantitatively.

Allanite minerals are the principal host of REEs in the Sin Quyen, Iron Oxide Copper Gold (IOCG) type deposit. The geochemical characteristics of these minerals are discussed in this work. The studied allanites have an unstable concentration of all major elements, such as REE (14-27 wt%), Ca (9-16 wt%), Al (8-19 wt%), Si (26-34 wt%) and Fe (12-21 wt%). Two different varieties of these minerals are documented, the older with higher REE concentrations ranging from 20 to 27 wt%, and younger with lower total REE concentration ranging from 14 to 19 wt%, which occur as a rim surrounding the older. Differences between the two groups of allanites are documented by Raman spectra and optical properties. The WDS chemical composition indicate that the allanites belong to the Ce-La-ferriallanite family, with low ƩHREE with an average of 0.21 wt.%. This work also supports the estimated timing of the deposit development focusing on detailed petrological study, and documented chemical composition of allanites confirmed by simplified statistical analysis. Temperature 355ºC which was calculated using value of δ34S isotopes is interpreted as a temperature of the second crystallization stage of allanite group. The pressure of crystallization solution was calculated and is ranging from 0.98 to 5.88 MPa.

Three-dimensional (3D) facies and petrophysical models were generated from previously published data of carbonate strata in the Dam Formation (eastern Saudi Arabia) to quantitatively investigate, describe, understand, model, and predict the permeability anisotropy of tidal flat carbonate within a sequence stratigraphic framework. The resulting 3D models were used to conduct fluid flow simulations to demonstrate how permeability anisotropy influences the production of hydrocarbons and ultimately affects decisions concerning future drilling in the exploration and development of carbonate reservoirs that have tidal flat strata. The constructed 3D facies model consists of four lithofacies associations, two of which were grain-dominated associations and two of which were mud-dominated associations. These lithofacies associations varied spatially in four reservoir zones (zones 1 to 4), which represent two fourth-order sequences in the uppermost part of the Dam Formation. Zones 1 and 3 consist of transgressive parasequences, and zones 2 and 4 consist of the regressive parasequences of these sequences. The 3D porosity and permeability models have a coherent match with the distribution of the lithofacies and the stratigraphic framework of the Dam Formation. The results suggested that the permeability anisotropy in zones 1 and 3 is controlled by the occurrence of the grain-dominated lithofacies associated with tidal flat channels. This lithofacies association overlies the sequence boundaries of sequences 1 and 3, forms reservoir bodies with relatively high permeability values, and is elongated perpendicular to the shoreline of the depositional environments. In contrast, permeability anisotropy in zones 2 and 4 is thought to be controlled by the occurrence of the grain-dominated lithofacies associated with the oolitic shoal. This lithofacies association overlies the maximum flooding surface of sequences 2 and 4, forms reservoir bodies with relatively high permeability values, and is elongated parallel to the shoreline of the depositional environments. Fluid flow simulation results suggested that the trend in hydrocarbon production from the constructed 3D models depends on permeability anisotropy in each reservoir zone. Thus, recognizing trends in permeability anisotropy, which might be predicted using sequences stratigraphy, could help to identify potential areas for future drilling.

The Doig Phosphate Zone (DPZ) is a phosphate-bearing marine unit located at the base of the Doig Formation, in the Triassic section of the Western Canada Sedimentary Basin. The DPZ has a maximum thickness of 90 m and extends across northeastern British Columbia and western central Alberta. In this study we characterize the significance and interpret the origin of apatite in the DPZ, through mineralogical and geochemical analyses, thin section study and field emission-scanning electron microscopy. The occurrence of apatite in the DPZ is not evenly distributed, but restricted to discrete 10 to 20 cm thick beds, located near the base of the DPZ. Phosphorites are of two types: grainstones composed primarily of unconformity-bounded coated grains, and intraclastic phosphorites composed of detrital silt-sized grains and apatite coated grains in a cryptocrystalline phosphatic matrix. The phosphorite beds are records of stratigraphic condensation due to low detrital input during transgression. The erosionally-truncated phosphatic coated grains and intraclasts are interpreted to be a result of various phases of phosphatization, exhumation, erosion reworking, winnowing and redeposition in alternating quiescence and storms or bottom currents. The abundance of pyrite and chalcophile trace elements, as well as the low concentration of proxy elements for organic matter productivity and preservation, are further evidences of stratigraphic condensation, with sulfidic pore water development and extensive organic recycling promoted by biological activity during the long exposure times. The phosphorites were formed under oxygenated water conditions, as suggested by the depletion in Ce and presence of a diverse benthic fauna.

During the period January 2014 – October 2018, four strong earthquakes occurred in the Ionian Sea, Greece. A rich aftershock sequence followed each event of them. More analytically, according to the manual solutions of National Observatory of Athens, the first event (K1), occurred on 26 January 2014 in Kefallinia Island with magnitude M_L = 5.8, which was followed by another in the same region (K2) on 3 February 2014 with magnitude M_L = 5.7. The third event occurred on 17 November 2015, M_L = 6.0 in Lefkas Island (L1) and the last on 25 October 2018, M_L = 6.6 in Zande Island (Z1). The first three of these earthquakes caused moderate structural damages mainly in houses and produced particular unrest to the local population. This work presents first the calculation of the source parameters of the strong events as well as for all earthquakes with magnitude M_L > 4.0, using the methodology of the Moment tensor inversion and secondary data from permanent GPS stations were analyzed to confirm the findings from seismological data and to investigate the displacement due to the earthquakes.

The paper is a kind of a small guide book through the textures and structures of river deposits in valleys. It presents principles of the flument systems with morphological river terraces and fluments (terrace bodies), different stages of the morphological terraces, structure of the flument, texture of terrace flights, terrace stacks and terrace rows. Special topics deal with flument overlaps and insight into the valley ground.

The storm of November 12th-13th, 2019 provoked the displacements of boulders in a central Mediterranean rocky coast; with reference to a selected area, prone to the boulder production and geomorphologically monitored for years, a field-oriented study approach was applied for the phenomenon, by collating data concerning pre-storm locations and kinematics of these boulders. The number of displaced boulders is 11, that is, in terms of morphological imprint of a specific storm, one of the major study cases for the Mediterranean. In addition, based on widely used hydrodynamic equations, the minimum wave height required to displace the boulders is assessed. The values conform with the expected values for the wave climate dominating during the causative meteorological event and give a measure of the energy of the storm slamming the coast. Boulder dislodgements usually have plays a key role in determining the rate of the coastal recession, likely also in the investigated area. In view of an adverse climate evolution with a possible increase of energy and frequency of severe storms, the results deriving from the study of this morphodynamics should be considered for the hazard assessment and coastal management.

We identify the source of the Mw = 5.6 earthquake that hit west-central Epirus on March 21, 2020 00:49:52 UTC. We use synthetic aperture radar interferograms tied to one permanent Global Navigation Satellite System (GNSS) station (GARD). We model the source by inverting the INSAR displacement data. The inversion model suggests a shallow source on a low-angle fault (39°) dipping towards east with a centroid depth of 8.5 km. The seismic moment deduced from our model agrees with those of the published seismic moment tensors. This geometry is compatible with the Margariti thrust fault within the collision zone between Apulia and Eurasia. We also processed new GNSS data and estimate a total convergence rate between Apulia and Eurasia of 8.9 mm yr-1 , of which shortening of the crust between the Epirus coastal GNSS stations and station PAXO in the Ionian Sea is equivalent to ~ 50% of it or 4.6 mm yr−1. A 60-km wide deformation zone takes up nearly most of the convergence between Apulia-Eurasia, trending N318°E. Its central axis runs along the southwest coast of Corfu, along the northeast coast of Paxos, heading toward the northern extremity of the Lefkada island.

The deep Cretaceous tight sandstone in Kuqa Depression of Tarim foreland basin is an ultra-low porosity and ultra-deep gas-bearing reservoir, which is characterized by small pores, fine throats, and poor connectivity. The wireline logging responses are so complex, and especially, it is difficult to identify fluid types from resistivity logs. Based on acoustic, density, and neutron logs response differences in gas and water layers, effective fluid sensitivity factors are constructed for gas layer identification. From conventional logs, acoustic-neutron porosity difference, density-neutron porosity difference, and triple-porosity ratio are all sensitive parameters to the gas layer. From the NMR logging response mechanism, the density and NMR porosity difference, and T2 geometric mean of the movable fluid are also two sensitive parameters to the gas layer. Based on these parameters, a series of fluid typing charts are constructed and their adaptabilities are analyzed and compared. By contrast, NMR log interpretation is better, and triple-porosity comprehensive method from conventional logs is also effective when NMR logging is not available. Finally, the comprehensive fluid typing strategy by combining some methods for ultra-low porosity tight sandstone is summarized and optimized. This study is another alternative for fluid identification using non-electrical logs.

The Ijen volcano has Pleistocene age (294.00 ± 0.03 Ma), and this stratovolcano was very acidic crater water with a pH of 0-1. The acidic crater water seeps into the Banyupait river flow. Asembagus is a research area located on the northern slope of the Ijen volcano, and the Banyupait River drains this. The acidic river water flows from the Ijen Crater Lake, so the pH of the water river was very acidic. This study used several different analytic methodologies with some previous researchers, namely using the method of geological mapping, pH measurement, spectrophotometry, IRMS (Isotope Ratio Mass Spectrometer), and the technique of Induced Coupled Mass Spectrometry (ICP-MS). Besides, the petrographic analysis is used to determine the composition of rock minerals due to rocks interacting with acidic water. Banyupait River water in the Asembagus area has a pH of river water around 3-7.3, SO4 (220-683 ppm), and the type of water is meteoric water. Also, concentrations of Ca, K, Mg in the west Banyupait river irrigation water flow showed higher levels when compared to the eastern Banyupait River water flow. Likewise, REE elements from the Asembagus region showed lower concentrations compared to Ijen Crater water. This change in the level of chemical elements is caused by the acidity of the Banyupait River being diluted or mixed with water from other water. However, the spring was not affected by acidic water. The process of acidic water interaction with rocks can also be observed from rocks traversed by the Banyupait River flow. Chalcedony and hematite replace the primary minerals of basaltic rocks. This research is expected to improve the quality of water needed by the Asembagus community so that people can live healthily.

Clay minerals surfaces potentially played a role in prebiotic synthesis through adsorption of organic monomers that give rise to highly concentrated systems; facilitate condensation and polymerization reactions; protection of early biomolecules from hydrolysis and photolysis; and surface-templating for specific adsorption and synthesis of organic molecules. This review presents processes of clay formation using saponite as a model clay mineral since it is shown to catalyze organic reactions, easy to synthesize in large and pure form, and has tunable properties. In particular, a method involving urea is presented as a reasonable analog of natural processes. The method involves a two-step process – 1) formation of the precursor aluminosilicate gel and 2) hydrolysis of a divalent metal (Mg, Ni, Co, Zn) by the slow release of ammonia from urea decomposition. The aluminosilicate gels in the first step forms a 4-fold-coordinated Al³⁺ similar to what is found in nature such as in volcanic glass. The use of urea, a compound figuring in many prebiotic model reactions, circumvents the formation of undesirable brucite, Mg(OH)₂, in the final product by slowly releasing ammonia thereby controlling the hydrolysis of magnesium. In addition, the substitution of B and Ga for Si and Al in saponite is also described The saponite products from this urea-assisted synthesis were tested as catalysts for several organic reactions including Friedel-Crafts alkylation, cracking and isomerization reactions.

Western Anatolia comprises a vast amount of various volcanic successions spanning from Eocene to Upper Miocene periods. These units mainly display southward younging in broad sense and display large amounts of chemical variation that spanned from basalt to rhyolite. The southward younging of magmatism and chemical variations have been largely attributed to the retreat and roll-back of the Hellenic slab and the western escape of the Anatolian microplate. However, there is still a lack of high precision data to pinpoint the exact nature of the magmatism and lithospheric tectonics. In this contribution we investigated a poorly known region along the Western Anatolia along Manisa district called Karakılıçlı volcanic field. We investigated two different volcanic sections (Kalpakkaya and Çamlık hill) that display the best volcano-sedimentary features in terms of geochronology and geochemistry. Samples acquired from the bottom, middle and upper portions of these sections display Early-Middle Miocene ages of 17.64±0.20, 17.22±0.15, 16.16±0.17 and 16.36±0.13, 15.79±0.71 and 13.61±0.20 Ma respectively. The results indicate that the volcanism in the region generated by the melting of the mantle and/or lithospheric mantle by slab retreat and roll-back of the Hellenic slab and evolved in the shallow magma chambers/mushes by fractional crystallization, magma mixing and crustal assimilation.

Landslides impact on human activities and socio-economic development especially in mountainous areas. This study focuses on the comparison of the prediction capability of advanced machine learning techniques for rainfall-induced shallow landslide susceptibility of Deokjeokri catchment and Karisanri catchment in South Korea. The influencing factors for landslides i.e. topographic, hydrologic, soil, forest, and geologic factors are prepared from various sources based on availability and a multicollinearity test is also performed to select relevant causative factors. The landslide inventory maps of both catchments are obtained from historical information, aerial photographs and performing field survey. In this study, Deokjeokri catchment is considered as a training area and Karisanri catchment as a testing area. The landslide inventories content 748 landslide points in training and 219 points in testing areas. Three landslide susceptibility maps using machine learning models i.e. Random Forest (RF), Extreme Gradient Boosting (XGBoost) and Deep Neural Network (DNN) are prepared and compared. The outcomes of the analyses are validated using the landslide inventory data. A receiver operating characteristic curve (ROC) method is used to verify the results of the models. The results of this study show that the training accuracy of RF is 0.757 and the testing accuracy is 0.74. Similarly, training accuracy of XGBoost is 0.756 and testing accuracy is 0.703. The prediction of DNN revealed acceptable agreement between susceptibility map and the existing landslides with training and testing accuracy of 0.855 and 0.802, respectively. The results showed that, the DNN model achieved lower prediction error and higher accuracy results than other models for shallow landslide modeling in the study area

This communication reports an improvement of the quality of the electrical data obtained from the application of electrical resistivity tomography method on archaeological studies. The electrical contact between ground and electrode enhances significantly by using carbomer-based gel during the electrical resistivity tomography measurements. Not only does the gel promote the conservation of the building surface under investigation, but it also virtually eliminates the necessity of conventional spike electrodes, which in many archaeological studies are inadequate or not permitted. Results evidenced an enhancement in the quality of the electrical data obtained in the order of thousands of units compared with those without using the carbomer-based gel. The potential and capabilities of this affordable gel make it appropriate to be applied to other geoelectrical studies beyond archaeological investigations. Moreover, it might solve corrosion issues on conventional spike electrodes, and electrical multicore cables usually provoked for added saltwater attempting to improve the electrical contact.

The study of hydrothermal travertines contributes to the understanding of the interaction between physico-chemical processes and the role played by microbial mats and biofilms in influencing carbonate precipitation. Three active travertine sites were investigated in Central Italy to identify the types of carbonate precipitates and the associated microbial mats at varying physico-chemical parameters. Carbonate precipitated fabrics at the decimetre- to millimetre-scale and microbial mats vary with decreasing water temperature: a) at high temperature (55-44°C) calcite or aragonite crystals precipitate on microbial mats of sulphide oxidizing, sulphate reducing and anoxygenic phototrophic bacteria forming filamentous streamer fabrics, b) at intermediate temperature (44-40°C), rafts, coated gas bubbles and dendrites are associated with Spirulina cyanobacteria and other filamentous and rod-shaped cyanobacteria, c) low temperature (34-33°C) laminated crusts and oncoids in a terraced slope system are associated with diverse Oscillatoriales and Nostocales filamentous cyanobacteria, sparse Spirulina and diatoms. At the microscale, carbonate precipitates are similar in the three sites consisting of prismatic calcite (40-100 µm long, 20-40 µm wide) or acicular aragonite crystals organized in radial spherulites, overlying or embedded within biofilm EPS (Extracellular Polymeric Substances). Microsparite and sparite crystal size decreases with decreasing temperature and clotted peloidal micrite dominates at temperatures < 40°C, also encrusting filamentous microbes. Carbonates are associated with gypsum and Ca-phosphate crystals; EPS elemental composition is enriched in Si, Al, Mg, Ca, P, S and authigenic aluminium-silicates form aggregates on EPS. This study confirms that microbial communities in hydrothermal travertine settings vary as a function of temperature. Carbonate precipitate types at the microscale do not vary considerably, despite different microbial communities suggesting that travertine precipitation, driven by CO2 degassing, is influenced by biofilm EPS acting as template for crystal nucleation (EPS-mediated mineralization) and affecting the fabric types, independently from specific microbial metabolism.

Hydrogen gas is seeping from the sedimentary basin of São Franciso, Brazil. The seepages of H₂ are accompanied by helium whose isotopes reveal a strong crustal signature. Geophysical data indicates that this intra-cratonic basin is characterized by i) a relatively high geothermal gradient, ii) deep faults delineating a horst and graben structure and affecting the entire sedimentary sequence, iii) an archean to paleoproterozoïc basements enriched in radiogenic elements and displaying mafic and ultramafic units, and iv) a possible karstic reservoir located 400 m below the surface. The high geothermal gradient could be due to a thin lithosphere enriched in radiogenic elements, which can also contribute to a massive radiolysis process of water at depth, releasing an important amount of H₂. Alternatively, ultramafic rocks that may have generated H₂ during their serpentinization are also documented in the basement. The seismic profiles show that the faults seen at the surface are deeply rooted in the basement, and can drain deep fluids to shallow depths in a short time scale. The carbonate reservoirs within the Bambuí group which forms the main part of the sedimentary layers are crossed by the fault system and represent good candidates for temporary H₂ accumulation zones. The formation by chemical dissolution of sinkholes located at 400 m depth might explain the presence of sub-circular depressions seen at the surface. These sinkholes might control the migration of gas from temporary storage reservoirs in the upper layer of the Bambuí formation to the surface. The very high fluxes of H₂ escaping out of these structures which have been recently documented are, however, in disagreement with the newly developed H₂ production model in the Precambrian continental crust. They either question the validity of these models or the measurement methodology.

The geometry of mesoscopic inhomogeneities plays an important role in determining the macroscopic propagation behaviors of elastic waves in a heterogeneous medium. Non-equiaxed inhomogeneities can lead to anisotropic wave velocity and attenuation. Developing an accurate scattering theory to describe the quantitative relation between the microstructure features and wave propagation parameters is of fundamental importance for seismology and ultrasonic nondestructive characterization. This work presents a multiple scattering theory for strongly scattering elastic media with general tri-axial heterogeneities. A closed analytical expression of the shape-dependent singularity of the anisotropic Green’s tensor for the homogeneous reference medium is derived by introducing a proper non-orthogonal ellipsoidal coordinate. Renormalized Dyson’s equation for the coherent wave field is then derived with the help of Feynman’s diagram technique and the first-order-smoothing approximation. The exact dispersion curves and the inverse Q-factors of coherent waves in several representative medium models for the heterogeneous lithosphere are calculated numerically. Numerical results for small-scale heterogeneities with the aspect ratio varying from 1 to 7 show satisfactory agreement with those obtained from real earthquakes. The results for velocity dispersion give rise to a novel explanation to the formation mechanism of different seismic phases. The new model has potential applications in seismology and ultrasonic microstructure characterization.

Lake Urmia in northwestern Iran is the largest lake in Iran and the second largest saltwater lake in the world. The water level in Lake Urmia has decreased dramatically in recent years, due to drought, climate change, and overuse of water resources for irrigation. This shrinking of the lake may affect local climate conditions, assuming that the lake itself affects the local climate. In this study, we quantified the lake’s impact on the local climate by analyzing hourly time series of data on climate variables (temperature, vapor pressure, relative humidity, evaporation, and dewpoint temperature for all seasons, and local lake/land breezes in summer) for the period 1961-2016. For this, we compared high quality, long-term climate data obtained from Urmia and Saqez meteorological stations, located 30 km and 185 km from the lake center, respectively. We then investigated the effect of lake level decrease on the climate variables by dividing the data into 1961-1995 (normal lake level) and 1996-2016 (low lake level). The results showed that at Urmia station (close to the lake), climate parameters displayed fewer fluctuations and were evidently affected by Lake Urmia compared with those at Saqez station. The effects of the lake on the local climate increased with increasing temperature, with the most significant impact in summer and the least in winter. The results also indicated that, despite decreasing lake level, local climate conditions are still influenced by Lake Urmia, but to a lesser extent.

Subsurface drainage systems remove excess water from the soil profile thereby improving crop yields in poorly drained farmland. Knowledge of the position of the buried drain lines is important: 1) to improve understanding of leaching and offsite release of nutrients and pesticides, and 2) for the installation of a new set of drain lines between the old ones for enhanced soil water removal efficiency. Traditional methods of drainage mapping involve the use of tile probes and trenching equipment. While these can be effective, they are also time-consuming, labor-intensive, and invasive, thereby entailing an inherent risk of damaging the drainpipes. Non-invasive geophysical soil sensors provide a potential alternative solution. Previous research has focused on the use of time-domain ground penetrating radar (GPR), with variable success depending on local soil and hydrological conditions and the central frequency of the specific equipment employed. The objectives of this study were 1) to test the use of a stepped-frequency continuous wave (SFCW) 3D-GPR (GeoScope Mk IV 3D-Radar with DXG1820 antenna array) for subsurface drainage mapping, and 2) to evaluate the performance of a 3D-GPR with the use of a single-frequency multi-receiver electromagnetic induction (EMI) sensor (DUALEM) in-combination. The 3D-GPR system offers more flexibility for application to different (sub)surface conditions due to the coverage of wide frequency bandwidth. The EMI sensor simultaneously provides information about the apparent electrical conductivity (EC_a) for different soil volumes, corresponding to different depths. This sensor combination was evaluated on twelve different study sites with various soil types with textures ranging from sand to clay till. While the 3-D GPR showed a high success rate in finding the drainpipes at five sites (sandy, sandy loam, loamy sand, and organic topsoils), the results at the other seven sites were less successful due to limited penetration depth (PD) of the 3D-GPR signal. The results suggest that the electrical conductivity estimates produced by the inversion of EC_a data measured by the DUALEM sensor could be a useful proxy to explain the success achieved by the 3D-GPR in finding the drain lines. The high attenuation of electromagnetic waves in highly conductive media limiting the PD of the 3D-GPR can explain the findings obtained in this research.

Space-borne Synthetic Aperture Radar (SAR) Interferometry (InSAR) is now a key geophysical tool for surface deformation studies. The European Commission’s Sentinel-1 Constellation began acquiring data systematically in late 2014. The data, which are free and open access, have global coverage at moderate resolution with a 6 or 12-day revisit, enabling researchers to investigate large-scale surface deformation systematically through time. However, full exploitation of the potential of Sentinel-1 requires specific processing approaches as well as the efficient use of modern computing and data storage facilities. Here we present LiCSAR, an operational system built for large-scale interferometric processing of Sentinel-1 data. LiCSAR is designed to automatically produce geocoded wrapped and unwrapped interferograms and coherence estimates, for large regions, at 0.001° resolution (WGS-84 system). The products are continuously updated in a frequency depending on prioritised regions (monthly, weekly or live update strategy). The products are open and freely accessible and downloadable through an online portal. We describe the algorithms, processing, and storage solutions implemented in LiCSAR, and show several case studies that use LiCSAR products to measure tectonic and volcanic deformation. We aim to accelerate the uptake of InSAR data by researchers as well as non-expert users by mass producing interferograms and derived products.

The basic principles of geological risk calculation through Probability of Success (PoS) are mostly applied for numerical estimation of additional hydrocarbon existence in proven reservoirs or potential hydrocarbon discoveries in selected geological regional subsurface volume. It can be tailored and validated for a comprehensive input dataset collected in the selected petroleum province, adapted by dividing up geological events into several probability categories and classes. The most applied categories are (existence of) reservoir rocks, traps and isolators, source rocks, migration pathways and preservation conditions for hydrocarbons The methodology results in unique probability values as multiplication of independent statistical events, which can also be applied in the assessment of a potential hydrocarbon discovery of desired minimal volume and its value in any virtual currency like risk-neutral dollars. Such methodology has been extensively developed in the last decades in the Croatian subsurface, mostly in the Croatian part of the Pannonian Basin System (CPBS). Through the adaptation of geological categories, it was also applied in hybrid, i.e., stochastical, models developed in the CPBS (Drava Depression). Stochastically estimation of porosity was already applied. As the robustness of this methodology is very high, it was also modified to estimate the influence of water-flooding in increasing oil recovery in some proven Neogene sandstone reservoirs in the CPBS (Sava Depression). This new modification is being presented to be applied to geological risk calculation, intending to assess the safety of geological environment in deep wells, where depleted radioactive fuel would be disposed, a subject of great importance. The case study encompassed the magmatic and metamorphic rocks in the pre-Neogene basement of the CPBS. For disposal purpose, these are regionally lithologies considered as to be the safest ones considering petrophysical values, water saturation, recent weathering and tectonic activity.

A quite energetic seismic excitation consisting of one main and additional three distinctive earthquake clusters that occurred in the transition area between the Kefalonia Transform Fault Zone (KTFZ) and the continental collision between Adriatic and Aegean microplates, is thoroughly studied after high–precision aftershocks’ relocation. The activated fault segments are in an area where historical and instrumental data have never claimed the occurrence of a catastrophic (M>6.0) earthquake. The relocated seismicity initially defines an activated structure extending from the northern segment of the Lefkada branch of KTFZ with the same NNE–SSW orientation and dextral strike-slip faulting and then keeping the same sense of motion its strike becomes NE–SW and its dip direction NW. This provides unprecedented information on the link between the KTFZ and the Collision front and sheds more light on the regional geodynamics. The earthquake catalog, which is specially compiled for this study, starts one year before the occurrence of the Mw5.4 mainshock and adequately provides the proper data source for investigating the temporal variation of the b–value, which might be used for discriminating foreshock and aftershock behavior.

Geomathematics is extremely important in geosciences, particularly in the geology. The key for any geomathematical analysis is the definition of a typical model to be applied for further prognosis, either through deterministical or stochastical approaches. The selection of the appropriate procedure is presented in this paper. Two different geomathematical subfield datasets were used in subsurface geological mapping and palaeontology and different biostatistics applications, representing important geomathematical subfields in the Croatian geology. The different subsurface interpolation methods, tested, validated and recommended for application, were used to obtain the best possible outcome in reservoir modelling, in the cases with small datasets. Cross-validation may be selected as the main selection criteria, applied to the Croatian part of the Pannonian Basin System (abbr. CPBS). Recent advances in biostatistics applied in palaeontology and case studies from Croatia are also presented, where biometric studies are of significant importance in fossil biota. Data, methods and problems in geosciences is a vast subject, and address a wide spectrum of fundamental science. Because geology includes subsurface and surface geology, and very different datasets regarding variable and number of data, here are chosen two representative case study groups with original samples from Northern Croatia. Subsurface mapping has been presented on limited petrophysical datasets from the Northern Croatian, Miocene, hydrocarbon reservoirs. Biostatistics has been presented on very different samples, allowing to achieve paleoenvironmental reconstructions of size of relevant fossils, as dinosaurs or other species and their paleoenvironments. All examples highlight examples of the valuable application of geomathematical tools in geology. The results, cautiously validated and correlated with other, non-numerical (indicator, categorical) geological knowledge, are of enormous assistance in creating better geological models.

Landslides are a frequent and recurrent problem in hilly regions of India and predicting them is always a challenging task. In this paper, an attempt was made to deal with this problem using advanced physical and numerical modeling methods. Detailed understanding of the initial slope failures is very interesting, and challenging at the same time, in the design and development of wireless sensor network based on early warning of landslide monitoring. A small scale physical model was developed to assess the instability through a sensor network with variable rain fall intensity. This was achieved by increasing the simulated rain water flow intensity in different time spans (dry condition, at t=0 to t= 30 min, 0.5 mm/min at t=30 to t= 60 min, 0.75 mm/min at t=60 to t=91 min and 1 mm/min at t=91 to t= 120 min). The water level and movement in the slope was recorded by rainfall sensor, vibration sensor, soil moisture sensor and a digital camera. The following changes were observed during the slope failure: a) movement of small particles at top of the slope; b) initial failure of medium size soil particle; c) scouring of soil mass; d) whole slope collapse. The obtained results clearly indicated the superiority and effectiveness of the proposed system in providing a factor of safety for the progressive slope.

This paper examines the influence of the petrographic characteristics of sandstones from Klepa Nafpaktias (Greece) on their suitability in construction (concrete) and energy storage applications. For this scope, ten sandstones were collected in order to study their petrographic characteristics using petrographic microscope and a GIS software as well as their basic physical, mechanical and physicochemical properties. Concrete specimens (C25/30) were produced with constant volume proportions, workability, mixing and curing conditions using different sizes of each aggregate type. Aggregates were mixed both in dry and water saturated states in concrete. Three different types of sandstone aggregates were examined and classified in three district groups according to their physicomechanical properties, petrographic characteristics and surface texture. The classification in groups after the concrete compressive strength test (UCS) verified the initial classification in the same three groups relative to their grain size from coarse to fine grained. As the grain size decreases their physicomechanical and physicochemical properties are getting better resulting in higher concrete strength values (25 to 32 MPa). Furthermore, the proposed ratio C/A (crystals/ mm2) seems to influence the aggregate properties which constitute critical factors for the final concrete strength, presenting the more fine grained sandstones as the most suitable for concrete aggregates. Concerning the use of Klepa Nafpaktias sandstones as potential energy reservoirs, the studied sandstones have the appropriate physicochemical properties for the implementation of a financially feasible CO₂ capture and storage scenario.

Legacy mining industry has left a large number of tailings ponds exposed to water and wind erosion that causes serious environmental and health problems. Prior to rehabilitation actions a deep sampling of the materials infilling the pond used to be necessary. Thus, the primary objective of this study is to demonstrate the usefulness of the Electrical Resistivity Tomography (ERT) method as a non-invasive tool to determine the physicochemical composition of mine tailings ponds, enabling more efficient and low-cost surveys. To achieve this objective, three ERT profiles and three boreholes in each profile were carried out, from each borehole three waste samples from differents depths were collected and a geochemical characterization of the samples was carried. In order to estimate the composition of the infilling wastes in tailing ponds from electrical resistivity measures, several regression models were calculated for different physicochemical properties and metal concentrations. As a result, a high resistivity area was depicted in profiles G2 and G3 while a non-resistive area (profile G1) was also found. Relationships among low resistivity values and high salinity, clay content and high metal concentrations and mobility were established. Specifically, calibrated models were obtained for electrical conductivity, particles sizes of 0.02-50 µm and 50-2000 µm, total Zn and Cd concentration, and bioavailable Ni, Cd and Fe. Therefore, the ERT technique could be considered as a useful tool for mine tailings ponds characterization, and it can be used to estmate some physicochemical properties and metal concentrations of this mine waste.

Morphometric analysis of Holocene pebbles from the Sava River gravels, in Zagreb alluvial aquifer system (NW Croatia), revealed distribution of their shapes along 30 km long observed watercourse. Limestones, dolomites and sandstones are determined as major (> 4%), and effusive magmatics, cherts and tuffs as minor lithotypes of the pebbles (up to 4%). Their distributions indicate mainly distant Alpine provenance for carbonate (limestones and dolomites) pebbles and local input for sandstones and minor lithotypes, laterally from the Samoborska gora and Medvednica Mts. Carbonates have predominately disc and sphere shapes, implying also their mainly distant sources. Scattered distributions of pebble shapes (sphere, disc, blade and rod) for sandstones and minor lithotypes indicate multiple sources, some of them probably local. Original sedimentary environments for main pebble lithotypes are tentatively interpreted from their flatness ratios, indicating predominant lake shore environments, followed by moraine and riverbed.

Tengchongite is a uranyl molybdate uranium mineral and it was found and named by Chen 1985.No more scholars studied on tengchongite after Chen’s work, The identification of this mineral has only been confirmed by single crystal X-ray diffraction. In the paper, micro laser Raman spectroscopy and infrared spectroscopy are used to identify the spectroscopy properties of tengchongite. The studies fill in the gaps of more than 30 years in terms of the molecular spectroscopy research of tengchongite. The mineral has an ideal model of Ca(UO₂)₆(MoO₄)₂O₅•12H₂O its bands attributed to the vibrating units are clearly identified in the Raman spectrum. Symmetric stretching modes at 812 cm^-1 and 839 cm^-1 are assigned to ν1 (UO₂)²⁺ ,The ν3 antisymmetric stretching modes of (UO₂)²⁺ are observed at 896 cm^-1, Symmetric stretching modes at 419^-1 and 440 cm^-1 are assigned to ν2 (UO₂)²⁺ .Symmetric stretching modes at 919cm-1 are assigned to ν1(MoO₄)^2-,The ν3 antisymmetric stretching modes of (MoO₄)^2- are observed at 752 cm^-1, the in-plane ν2(MoO₄)^2-and the out-of-plane ν4(MoO₄)^2- bending modes are at 169 cm^-1 and 254 cm^-1. IR spectrum of tengchongite shows the major uranyl band at 858 cm^-1 and 693 cm^-1, Mo-O bonds are observed at about 985 cm^-1 and 780 cm^-1, and H₂O groups are present with a wide range peaks from 3100 cm^-1 to 3500 cm^-1 and 1647 cm^-1, and the band at 1432.4 cm^-1 is probably due to the stretching vibration hydroxyl (OH^-1), therefore, tengchongite contains may include hydroxyl and its chemical formula needs to be modified .

Pore-fractures network play a key role in coalbed methane (CBM) accumulation and production, while the impacts of coal facies on the pore-fractures network performance are still poorly understood. In this work, the research on the pore-fracture occurrence of 38 collected coals from Sangjiang-Muling coal-bearing basins with multiple techniques including mercury intrusion porosimetry (MIP), micro-organic quantitative analysis, and optic microscopy, and its variation controlling of coal face were studied. The MIP curves of 38 selected coals indicating pore structures were subdivided into three typical types including type I of predominant micropores, type Ⅱ of predominant micropores and macropores with good connectivity and type Ⅲ of predominant micropores and macropores with poor connectivity. For coal facies, there are three various coal facies were distinguished, which include lake shore coastal wet forest swamp, the upper delta plain wet forest swamp, tidal flat wet forest swamp with Q-cluster analysis and tissue preservation index - gelification index (TPI-GI) and Wood index - groundwater influence index (WI -GWI). The results show there is positive relationship between tissue preservation index (TPI), wood index (WI) and mesopores (102nm-103nm), while a negative relationship between TPI, WI and macropores/fractures. In addition, groundwater level fluctuations can control the development of type C and D fractures, and the frequency of type C and D fractures shows an ascending trend with increasing GWI, which may be caused by the mineral hydration of the coal. Finally, from the perspective of the pore-fractures occurrence in CBM reservoirs, the wet forest swamp of upper delta plain is considered to be the optimization areas for Sanjiang-Mulinghe coal-bearing basins by a comparative study of various coal facies.

Land Surface Temperature (LST) and Land Use Land Cover (LULC) are the principal aspects of climate and environment studies. The object of the study is to assess spatial relationship between LST and remote sensing LULC indices at the global and continental scale. Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua daytime LST and eight LULC MODIS indices of 2018 prepared and processed using Earth Engine Code Editor. R squared and significance of the relationship values of randomly selected points computed in R program. The research observed the relationship between examined indices and LST is significant at the 0.001 level. Normalized Difference Water Index (NDWI) and Normalized Difference Snow Index (DSI) are the dominant drivers of LST in the world, Asia and North America. In Australia and Africa, Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) are the dominant drivers of LST. Albedo and Normalized Difference Soil Index (NDSI) have superior in Central America. In South America and Europe, the dominant driver of LST is NDWI. Relationship between albedo and LST is moderate inverse on a global scale. Observed relationship between LST and examined vegetation indices is positive in Europe and North America while inverse in Australia and Africa. All observed relationship between Normalized Difference Built-up Index (NDBI) and LST are positive. Association observed between NDSI and LST is positive in Australia, Africa and Central America.

This study presents an exploration into identifying the interactions between ocean waves and the continental margin in the origination of double-frequency (DF, 0.1-0.5 Hz) microseisms recorded at 33 stations across East Coast of USA (ECUSA) during a ten-day period of ordinary ocean wave climate. Daily primary vibration directions are calculated in three frequency bands and projected as great circles passing through each station. In each band, the great circles from all stations exhibit largest spatial density primarily near the continental slope in the western North Atlantic Ocean. Generation mechanisms of three DF microseism events are explored by comparing temporal and spatial variations of the DF microseisms with the migration patterns of ocean wave fronts in Wavewatch III hindcasts. Correlation analyses are conducted by comparing the frequency compositions of and calculating the correlation coefficients between the DF microseisms and the ocean waves recorded at selected buoys. The observations and analyses lead to a hypothesis that the continental slope causes wave reflection, generating low frequency DF energy and that the continental shelf is where high frequency DF energy is mainly generated in ECUSA. The hypothesis is supported by the primary vibration directions being mainly perpendicular to the strike of the continental slope.

Here is analysed the process of water re-injection in the two hydrocarbon reservoirs/fields in the Upper Miocene sandstone reservoirs, located in the western part of the Sava Depression (Croatia). Namely, those are the "A" field with "L" reservoir and the "B" field with "K" reservoir. Both currently produce hydrocarbons using a secondary recovery method, i.e. water injection (in fact re-injection of the field waters). Three regional reservoir variables had been analysed, namely porosity, permeability and injected water volumes. The number of data was small in all three cases. For porosity: reservoir “L” included 25 data, reservoir “K” 19 data; for permeability: reservoir “L” 10 data, reservoir “K” 18 data; for injected volumes of water: reservoir “L” 10 data; reservoir “K” 3 data. It defined selection of mapping algorithms mostly designed for small datasets (less than 20 points), i.e. Inverse Distance Weighting, Nearest and Natural Neighbourhood. Additionally, the Ordinary Kriging was used, but only with jack-knifed variograms, producing many “artificial points”. Results are extensively tested, using cross-validation and shape recognitions, and the Inverse Distance Weighting method is described as the most appropriate approach for mapping permeability and injected volumes in both reservoirs (“K” and “L”). The Kriging could be slightly outlined as the best approach for porosity. Obtained maps made possible application of the modified geological probability calculation as tools for prediction of successfulness of future injection (probability of 0.56). Consequently, results made possible to plan future injection more efficiently, with smaller injected volumes and same of higher hydrocarbon recovery. That could prevent useless injection, decrease number of injection wells, and save energy and funds invested in such processes.

Eight samples of limestones and marbles were studied by neutron diffraction to collect 2 Quantitative Texture (i.e., Crystallographic Preferred Orientations or CPO) of calcite deforming at 3 different depths in the crust. We studied the different Texture patterns developed in shear zones at 4 different depth and their influence on seismic anisotropies. Samples were collected in the French and 5 Italian Alps, Apennines, and Paleozoic Sardinian basement. They are characterized by isotropic to 6 highly anisotropic (e.g., mylonite shear zone) fabrics. Mylonite limestones occur as shear zone horizons 7 within the Cenozoic Southern Domain in Alpine thrust-and-fold belts (Italy), the Briançonnais domain 8 of the Western Alps (Italy-France border), the Sardinian Paleozoic back-thrusts or in the Austroalpine 9 intermediate units. The analyzed marbles were collected in the Carrara Marble, in the Austroalpine Units 10 in the Central (Mortirolo) and Western Alps (Valpelline). The temperature and depth of development of fabrics vary from < 100◦C, to 800◦C and depth from <10 km to about 30 km, corresponding from upper 12 to lower crust conditions. Quantitative Texture Analysis shows different types of patterns for calcite: 13 random to strongly textured. Textured types may be further separated in orthorhombic and monoclinic 14 (Types A and B), based on the angle defined with the mesoscopic fabrics. Seismic anisotropies were 15 calculated by homogenizing the single crystal elastic tensor, using the Orientation Distribution Function 16 calculated by the Quantitative Texture Analysis. The resulting P- and S-waves anisotropies show a wide 17 variability due to the textural types, temperature and pressure conditions, and dip of the shear planes.

The akimotoite--majorite--bridgmanite (Ak--Mj--Bm) triple point in MgSiO3 has been measured in large-volume press (LVP; COMPRES 8/3 assembly) and laser-heated diamond anvil cell (LHDAC). For the LVP data, we calculated pressures from the calibration by Leinenweber (2012). For the LHDAC data, we conducted \textit{in situ} determination of pressure at high temperature using the Pt scale by Dorogokupets and Dewaele (2007) at synchrotron. The measured temperatures of the triple point are in good agreement between LVP and LHDAC at 1990-2000 K. However, the pressure for the triple point determined from the LVP is 3.9+/-0.6~GPa lower than that from the LHDAC dataset. The triple point determined through these experiments will provide an important reference point in the pressure-temperature space for future high-pressure experiments and allow mineral physicists to compare the pressure--temperature conditions measured in these two different experimental methods.

Abstract: Landslide is a frequent natural hazard in Chittagong Hilly Areas (CHA), Bangladesh, which causes the loss of lives and damage to the economy. Despite that, an official landslide inventory is still lacking in this area. In this paper, we present a landslide inventory of this area prepared using the visual interpretation of Google Earth images (Google Earth Mapping), field mapping, and literature search. We mapped 730 landslides that occurred from 2001 to 2017. Different landslide attributes, including type, size, distribution, state, water content, future risk, and causes, are presented in the dataset. In this area, slide and flow were the two dominant types of landslides. Among five districts (Bandarban, Chittagong, Cox's Bazar, Khagrachari and Rangamati), most (54.93%) of the landslides occurred in Chittagong and Rangamati districts. About 45.09% of the landslides were small (<100 m2) in size while the maximum size of the detected landslides was 85201.6 m2. This dataset will help to understand the characteristics of landslides in CHA and provide useful guidance for policy implementation.

We present in this work a method for estimation of plasma bubble mean zonal drift velocities using keograms generated from images of the OI 6300.0 nm nightglow emission collected from an equatorial station -- Cariri (7.4$^\circ$S, 36.5$^\circ$W), and a mid-latitude station -- Cachoeira Paulista (22.7$^\circ$S, 45$^\circ$W), both in the Brazilian sector. The mean zonal drift velocities were estimated for 239 events recorded from 2000 to 2003 in Cariri, and for 56 events recorded over Cachoeira Paulista from 1998 to 2000. It was found that plasma bubble zonal drift velocities are smaller ($\sim$60 ms$^{-1}$) for events occurring later in the night compared to those occurring earlier ($\sim$150 ms$^{-1}$). The decreasing rate of the zonal drift velocity is of $\sim$10 ms$^{-1}$/h. We have also found that, in general, bubble events appearing first in the west-most region of the keogram are faster than those appearing first in the east-most region of the keograms. Larger zonal drift velocities occur from 19 LT to 23 LT in a longitude range from 37$^\circ$ to 33$^\circ$. The method of velocity estimation using keograms compares favorably against the mosaic method developed by \cite{Arruda:2005}, but the standard deviation of the residuals for the zonal drift velocities from the two methods is $\sim$15 ms$^{-1}$

Knowledge of the spatio-temporal occurrence of avalanche activity is critical for avalanche forecasting and hazard mapping. We present a near-real time automatic avalanche monitoring system that outputs detected avalanche polygons within roughly 10 min after Sentinel- 1 SAR data download. Our avalanche detection algorithm has an average probability of detection of 67.2 % with a false alarm rate averaging 45.9, with maximum POD's over 85 % and minimum FAR's of 24.9 % compared to manual detection of avalanches. The high variability in performance stems from the dynamic nature of snow in the Sentinel-1 data. After tuning parameters of the detection algorithm, we processed five years of Sentinel-1 images acquired over a 150 x 100 km large area in Northern Norway, with the best setup. Compared to a dataset of field-observed avalanches, 77.3 % were manually detectable. Using these manual detections as benchmark, the avalanche detection algorithm achieved an accuracy of 79 % with high POD's in cases of medium to large wet snow avalanches. For the first time, we can present a dataset of spatiotemporal avalanche activity over several winters from a large region. This unique dataset allows for research into the relationship between avalanche activity and triggering meteorological factors, mapping of avalanche prone areas and near-real time avalanche activity monitoring to assist public avalanche forecasting. Currently, the Norwegian Avalanche Warning Service is using our processing system for pre-operational use in three regions in Norway.

High frequency knowledge of the spatio-temporal distribution of the Downwelling Surface Shortwave Flux (DSSF) and its diffuse fraction (fd) at the surface is nowadays essential for understanding climate processes at the surface-atmosphere interface, plant photosynthesis and carbon cycle, and for the solar energy sector. The EUMETSAT Satellite Application Facility for Land Surface Analysis operationally delivers estimation of the MDSSFTD (Downwelling Surface Short-wave radiation Fluxes – Total and Diffuse fraction) product with an operational status since the year 2019. The method for the retrieval was presented in the companion paper [40]. The part 2 now focuses on the evaluation of the MDSSFTD algorithm and presents the comparison of the corresponding outputs, i.e. total DSSF and diffuse fraction (fd) components, against in-situ measurements acquired at four BSRN stations over a seven-month period. The validation is performed on an instantaneous basis. We show that the satellite estimates of DSSF and fd meet the target requirements defined by the user community for all-sky (clear and cloudy) conditions. For DSSF, the requirements are 20Wm^-2 for DSSF<200Wm^-2, and 10% for DSSF>=200Wm^-2. The MBE and rMBE compared to the ground measurements are 3.618Wm^-2 and 0.252%, respectively. For fd, the requirements are 0.1 for fd<0.5, and 20% for fd>=0.5. The MBE and rMBE compared to the ground measurements are -0.044 and -17.699%, respectively. The study also provides a separate analysis of the product performances for clear sky and cloudy sky conditions. The importance of representing the cloud-aerosol radiative coupling in the MDSSFTD method is discussed. Finally, it is concluded that the quality of the Aerosol Optical Depth (AOD) forecasts currently available is enough accurate to obtain reliable diffuse solar flux estimates. This quality of AOD forecasts was still a limitation a few years ago.

Dayaoshan, as an important metal ore producing area in China, is faced with the dilemma of resource depletion due to long-term exploitation. In this paper, remote sensing method is used to circle the favorable metallogenic areas and find new ore points for Gulong. Firstly, vegetation interference bas been removed by using mixed pixel decomposition method with hyperplane and genetic algorithm (GA) optimization; then, altered mineral distribution information has been extracted based on principal component analysis (PCA) and support vector machine (SVM) method; Thirdly, the favorable areas of gold mining in Gulong has been delineated by using ant colony algorithm (ACA) optimization SVM model to remove false altered minerals; Lastly, field survey verified that the extracted alteration mineralization information is correct and effective. The results show that the mineral alteration extraction method proposed in this paper has certain guiding significance for metallogenic prediction by remote sensing.

The origin and significance of the tonalite–trondhjemite–granodiorite (TTG) units and the familiar metabasite xenoliths they host in the Yangtze Craton, China, remain controversial, and resolving these issues is important if we are to understand the evolution of the early Yangtze Craton. We focused on biotite–tremolite schist xenoliths in the Archean TTG units of the Kongling high-grade metamorphic terrane, and U–Pb dating of their zircons yielded ²⁰⁷Pb/²⁰⁶Pb ages of ca. 3.00 Ga, which provides a minimum age for the formation of the pre-metamorphic basic igneous rock. The host TTGs and late intrusive granitic dikes yield three groups of upper intercept ages at 2.87–2.88, 2.91–2.94, and 3.07 Ga, and a concordant age at 2.94 Ga, which suggest that the Yangtze continental nucleus underwent three important metamorphic–magmatic events in the Mesoarchean at ca. 3.00, 2.94, and 2.87 Ga. The biotite–tremolite schists have high ratios of K₂O/Na₂O and high contents of CaO, Cr, and Ni, thus showing the characteristics of high-K calc-alkaline island-arc volcanic rocks (basalt–andesite) that form by the partial melting of subducted oceanic crust. The data also provide further proof that a Mesoarchean metamorphic basement exists in the Yangtze Plate. Derivation of the magmatic protoliths of the biotite–tremolite schist enclaves from an oceanic crust during slab subduction, and the presence of these xenoliths within the TTG suite, indicate the existence of the initiation of plate tectonics during the Mesoarchean (≤2.94 Ga).

Several studies have shown that changes in incoming solar radiation and variations of the diffuse fraction can significantly modify the vegetation carbon uptake. Hence, monitoring the incoming solar radiation at large scale and with high temporal frequency is crucial for this reason along with many others. The EUMETSAT Satellite Application Facility for Land Surface Analysis (LSA SAF) operationally disseminates in near real time estimates of the downwelling shortwave radiation at the surface since 2005. This product is derived from observations provided by the SEVIRI instrument onboard the Meteosat Second Generation series of geostationary satellites, which covers Europe, Africa, the Middle East, and part of South America. However, near real time generation of the diffuse fraction at the surface level has only recently been initiated. The main difficulty towards achieving this goal was the general lack of accurate information on the aerosol particles in the atmosphere. This limitation is nowadays less important thanks to the improvements in atmospheric numerical models. This study presents an upgrade of the LSA-SAF operational retrieval method, which provides the simultaneous estimation of the incoming solar radiation and its diffuse fraction from satellite every 15 minutes. The upgrade includes a comprehensive representation of the influence of aerosols based on physical approximations of the radiative transfer within an atmosphere-surface associated medium. This article explains the retrieval method, discusses its limitations and differences with the previous method, and details the characteristics of the output products. A companion article will focus on the evaluation of the products against independent measurements of solar radiation. Finally, the access to the source code is provided through an open access platform in order to share with the community the expertise on the satellite retrieval of this variable.

Sea level rise is one of the main factor of risk for the preservation of cultural heritage sites located along the coasts of the Mediterranean basin. Coastal retreat, erosion and storm surges are yet posing serious threats to archaeological and historical structures built along the coastal zones of this region. In order to assess the coastal changes by the end of 2100 under an expected sea level rise of about 1 m, a detailed determination of the current coastline position and the availability of high resolution DSM, is needed. This paper focuses on the use of very high-resolution UAV imagery for the generation of ultra-high resolution mapping of the coastal archaeological area of Pyrgi, near Rome (Italy). The processing of the UAV imagery resulted in the generation of a DSM and an orthophoto, with an accuracy of 1.94 cm/pixel. The integration of topographic data with two sea level rise projections in the IPCC AR5 2.6 and 8.5 climatic scenarios for this area of the Mediterranean, were used to map sea level rise scenarios for 2050 and 2100. The effects of the Vertical Land Motion (VLM) as estimated from two nearby continuous GPS stations located as much as close to the coastline, were included in the analysis. Relative sea level rise projections provide values at 0.30±0.15 cm by 2050 and 0.56±0.22 by 2100, for the IPCC AR5 8.5 scenarios and at 0.13±0.05 cm by 2050 and 0.17±0.22 by 2100, for the IPCC AR5 2.6 scenario. These values of rise will correspond to a potential beach loss between 12.6% and 23.5% in 2100 for RCP 2.6 and 8.5 scenarios, respectively, while during the highest tides the beach will be reduced up to 46.4%. With these sea level rise scenarios, Pyrgi with its nearby Etruscan temples and the medieval castle of Santa Severa will be soon exposed to high risk of marine flooding, especially during storm surges, thus requiring suitable adaptation strategies.

There is an important complementarity between experimental methods for the study of high-temperature viscoelasticity in the time and frequency domains, that has not always been fully exploited. Here we show that parallel processing of forced-oscillation data and microcreep records, involving consistent use of either Andrade or extended Burgers creep-function models, yields a robust composite modulus-dissipation dataset spanning a broader range of periods than either technique alone. In fitting this dataset, the alternative Andrade and extended Burgers models differ in their partitioning of strain between the anelastic and viscous contributions. The extended Burgers model is preferred because it involves a finite range of anelastic relaxation times, and accordingly a well-defined anelastic relaxation strength. The new strategy offers the prospect of better constraining the transition between transient and steady-state creep, or equivalently, between anelastic and viscous behaviour.

There is crescent demand for knowledge improvement of the integrated water vapor (IWV) distribution in regions affected by heat islands that are associated with extreme rainfall events such as in the metropolitan area of Rio de Janeiro (MARJ). This work assessed the suitability and distribution of IWV in the MARJ using products from the Global Navigation Satellite Systems (GNSS), MODerate Resolution Imaging Spectroradiometer (MODIS), and radiosonde. GNSS data were collected by the tracking station named RDJN, from the cooperation of the International GNSS Monitoring and Assessment System (iGMAS) and the National Observatory of Brazil (Observatório Nacional - ON), and the tracking stations ONRJ, RIOD, and RJCG belonging to the Brazilian Network for Continuous Monitoring (RBMC) in the period of January 2015–August 2018. High variability of the near surface air temperature (T) and relative humidity (RH) were observed among eight meteorological sites considered. The mean T differences between sites, up to 4.4 °C, led to mean differences as high as 3.1 K for weighted mean temperature (Tm) and hence 0.83 mm for IWV differences. The performance of the MODIS MOD07 and MYD07 products provided a reasonably good representation of the mean spatial distribution of IWV, especially during the daylight passages of the satellites TERRA and AQUA. Local grid points of MODIS IWV estimates had relatively good agreement with the GNSS-derived IWV, with mean differences from -2.4–1.1 mm considering only daytime passages of the satellites TERRA and AQUA. During nighttime, MODIS underestimated IWV (from -9–-3 mm) with respect to GNSS, due to attenuation of IR radiation by clouds. A contrasting behavior was found in the radiosonde IWV estimates compared with the estimates from GNSS. There were dry biases of 1.4 mm (3.7% lower than expected) by radiosonde IWV during the daytime considering that all other estimates were unbiased and the differences between IWV GNSS and IWV RADS were consistent. Based on the IWV comparisons between radiosonde and GNSS at nighttime, the atmosphere over the radiosonde site is about 1.2 (2.3) mm wetter than over RIOD (RDJN) station. The long time series of the comparisons between IWVRDJN and IWVRIOD showed that the highest values of IWV occurred from the afternoon to nocturnal hours. Further, the atmosphere over the site RIOD was consistently about 1 mm wetter than over RDJN. These results showed the feasibility of the iGMAS RDJN station data compared with the RBMC, MODIS, and radiosonde data to investigate IWV in a region with occurrence of heat islands, and the peculiar physiographic and meteorological characteristics as in the MARJ. This work recommended the usage of complete meteorological station data collocated near every GNSS receiver aiming improvements of local GNSS IWV estimates and serving as additional support for operational numerical assimilation, weather forecast, and nowcast of extreme rainfall events.

In South Korea, landslides are caused by localized heavy rainfall and typhoons, which often occur in the summer season at natural slopes in mountainous areas and artificial slopes in urban surroundings. Flow-type landslides frequently occur in mountainous areas. To evaluate flow-type landslides, it is essential to identify the physical characteristics of soil, giving focus to the soil on the top layers of various types of slope. This study conducts a survey and an analysis of the characteristics of landslides that occurred in the study area with different geological conditions of granite and gneiss. The characteristics of soil in the area and its surroundings that have or have not undergone landslides for every geological condition is also evaluated. Based on these characteristics and a statistics method, it extracts the triggering factors, permeability coefficients (k), and shearing strength with cohesion (c) and internal friction angel (φ) of soils that are highly related to landslides around weathered soil layers. As a result, the permeability coefficients show significant relevance with void ratio (e), the effective size of grains (D10), and uniformity coefficient (cu), while the shearing strength with the proportion of fine-grained soil (Fines), uniformity coefficient (cu), degree of saturation (S), dry weight density (rd), and void ratio (e). By obtaining this result, the study uses the regression analysis to suggest models to estimate the permeability coefficients and shearing strength. For the gneiss area, the statistics-based estimation model (SEM) is proposed as k_gn = (1.488 × 10^-02 × e) + (1.076 × D₁₀) + (-1.629 × 10^-04 × c_u) - (1.893 × 10^-02) for permeability coefficients; cgn = (-0.712 × Fines) + (-0.131 × c_u) + 15.335 for cohesion; and φgn = (27.01 × r_d) + (-12.594 × e) + 6.018 for internal frictional angle of soils. For the granite area, the statistics-based estimation model (SEM) is proposed as kgr = (8.281 × 10^-03 × e) + (0.639 × D₁₀) + (-2.766 × 10^-05 × c_u) - (9.907 × 10^-03) for permeability coefficients; c_gr = (-0.689 × Fines) + (-0.0744 × S) + 18.59 for cohesion; and φ_gr = (33.640 × r_d) + (-0.875 × e) - 9.685 for internal frictional angle of soils. The use of statistics-based estimation models (SEMs) for landslide-triggering factors that trigger landslides will support the simple calculation of permeability coefficient and shearing strength (cohesion and internal frictional angle), only requiring information about the physical properties of soil at the natural slopes that have different geological features such as gneiss and granite areas.

Understanding the physical properties of ultramafic rocks is important for evaluating awide variety of petrologic models of the oceanic lithosphere, particularly upper mantle and lower crust. Hydration of oceanic peridotites results in increasing serpentine content, which affects lithospheric physical properties and the global bio/geochemical cycles of various elements. In understanding tectonic, magmatic and metamorphic history of the oceanic crust, interpreting seismic velocities, rock composition and elastic moduli are of fundamental importance. In this study we show that as serpentine content increases, density decreases linearly with a slope of 7.85. We also correlate increase in serpentine content with a linear decline in shear, bulk and Young’s moduli with slopes of 0.48, 0.77, 0.45 respectively. Our results show that increase in serpentine content of lower crust and forearc mantle could decrease elasticity of lithospehere and result in break-offs. Therefore tectonic processes at peridotite rich slow spreading ridges may be strongly affected by serpentine content, particularly serpentinization may be responsible for discontinuities in thin crust, and formation of weak fault zones.

Temperature-index modeling is used to determine the magnitude of temperature depression in the northern Sawatch Range required to maintain steady-state mass balances of six reconstructed glaciers at their extent during the local Last Glacial Maximum (LLGM), dated at ~21 ka. Assuming no significant differences in precipitation compared to modern values, mean annual temperatures in the region were on average 8.8 +0.5/–0.8 °C cooler than they are today. Allowing for modest increases or decreases in precipitation, required temperature depressions only differ by ± 0.2 °C. Temperature depression in the northern Sawatch Range are consistent, although slightly greater, with those determined in other ranges in Colorado using similar approaches. The estimates presented here are, however, substantially less than those suggested by several downscaled simulations of global LGM climate, that might be due to the need for improved calibration of such downscalings, or the models from which they are derived. Our estimates of LGM temperature depression are considerably greater than that previously determined in the study area and those in two other ranges in Colorado derived using different methodologies, the latter being most likely responsible for the discrepancies.

Retrogressive breach failures or coastal flow slides occur naturally in the shoreface in fine sands near dynamic tidal channels or rivers. They sometimes retrogress into beaches, shoal margins and river banks where they can threaten infrastructure and cause severe coastal erosion and flood risk. Ever since the first reports were published in the Netherlands over a century ago, attempts have been made to understand the geo-mechanical mechanism of flow slides. In this paper we have established that events, observed during the active phase, are characterized by a slow and steady retrogression into the shoreline, often continuing for many hours. This can be explained by the breaching mechanism, as elaborated in this paper. Recently, further evidence has become available in the form of video footage of active events in Australia and elsewhere, often publicly posted on the internet. All these observations justify the new term ‘retrogressive breach failure’ (RBF event). The mechanism has been confirmed in small-scale flume tests and in a large-scale field experiment. With a better understanding of the geo-mechanical mechanism, current protection methods can be better understood and new defense strategies can be envisaged. In writing this paper, we hope that the coastal science and engineering communities will better recognize and understand these intriguing natural events.

The Greenland Ice Sheet is now the largest land ice contributor to global sea level rise, largely driven by increased surface meltwater runoff from the ablation zone, i.e. areas of the ice sheet where annual mass losses exceed gains. This small but critically important area of the ice sheet has expanded in size by ~50% since the early 1960s, and satellite remote sensing is a powerful tool for monitoring the physical processes that influence its surface mass balance. This review synthesizes key remote sensing methods and scientific findings from satellite remote sensing of the Greenland Ice Sheet ablation zone, covering progress in 1) radar altimetry, 2) laser (lidar) altimetry, 3) gravimetry, 4) multispectral optical imagery and, 5) microwave and thermal imagery. Physical characteristics and quantities examined include surface elevation change, gravimetric mass balance, reflectance, albedo, and mapping of surface melt extent and glaciological facies and zones. The review concludes that future progress will benefit most from methods that combine multi-sensor, multi-wavelength, and cross-platform datasets designed to discriminate the widely varying surface processes in the ablation zone. Specific examples include fusing laser altimetry, radar altimetry, and optical stereophotogrammetry to enhance spatial measurement density, cross-validate surface elevation change, and diagnose radar elevation bias; fusing optical imagery, radar imagery, and microwave scatterometry to discriminate between snow, liquid water, refrozen meltwater, and bare ice near the equilibrium line altitude; combining optical reflectance with laser altimetry to map supraglacial lake, stream, and crevasse bathymetry; and monitoring the inland migration of snowlines, surface melt extent, and supraglacial hydrologic features.

We replace the current industry-standard empirical forecasts of oil production from hydrofractured horizontal wells in shales with a statistically and physically robust, accurate and precise approach, using the Bakken shale as an illustration. The proposed oil production forecasting method extends our previous work on predicting fieldwide gas production in the Barnett shale and merges it with our new scaling of oil production in shales. We first divide the existing 14,678 horizontal oil wells in the Bakken into 12 static samples in which depositional settings and completion technologies are similar. For each sample, we construct a purely data-driven P50 well prototype by merging the GEV distribution fits of oil production from appropriate well cohorts. We fit the parameters of our physics-based scaling curve to the statistical well prototypes, and obtain their smooth extrapolations to 30 years on production. By calculating the number of potential wells of each Bakken region, and scheduling future drilling programs, we stack up the extended well prototypes to achieve the most plausible forecast. We predict that Bakken will ultimately produce 5 billion barrels of oil from the existing wells, with the possible increments of 2 and 6 billion barrels from core and noncore areas.

Along with density and mass variations of the oceans driven by global warming, Glacial Isostatic Adjustment (GIA) in response to the last deglaciation still contributes significantly to present-day sea-level change. Indeed, in order to reveal the impacts of climate change, long term observations at tide gauges and recent absolute altimetry data need to be decontaminated from the effects of GIA. This is now realized by means of global models constrained by the observed evolution of the paleo-shorelines since the Last Glacial Maximum, which account for the complex interactions between the solid Earth, the cryosphere and the oceans. In the recent literature, past and present-day effects of GIA are often expressed in terms of fingerprints describing the spatial variations of several geodetic quantities like crustal deformation, the harmonic components of the Earth's gravity field, relative and absolute sea level. However, since it is driven by the sluggish readjustment occurring within the viscous mantle, GIA shall taint the pattern of sea-level variability also during the forthcoming centuries. The shapes of the GIA fingerprints reflect inextricable deformational, gravitational, and rotational interactions occurring within the Earth system. Using up-to-date numerical modeling tools, our purpose is to revisit and to explore some of the physical and geometrical features of the fingerprints, their symmetries and intercorrelations, also illustrating how they stem from the fundamental equation that governs GIA, i.e., the Sea Level Equation.

Scattering of elastic waves in heterogeneous media has become one of the most important problems in the field of wave propagation due to its broad applications in seismology, natural resource exploration, ultrasonic nondestructive evaluation and biomedical ultrasound. Nevertheless, it is one of the most challenging problems because of the complicated medium inhomogeneity and the complexity of the elastodynamic equations. A widely accepted model for the propagation and scattering of elastic waves, which properly incorporates the multiple scattering phenomenon and the statistical information of the inhomogeneities is still missing. In this work, the author developed a multiple scattering model for heterogeneous elastic continua with strong property fluctuation and obtained the exact solution to the dispersion equation under the first-order smoothing approximation. The model establishes an accurate quantitative relation between the microstructural properties and the coherent wave propagation parameters and can be used for characterization or inversion of microstructures. Starting from the elastodynamic differential equations, a system of integral equation for the Green functions of the heterogeneous medium was developed by using Green’s functions of a homogeneous reference medium. After properly eliminating the singularity of the Green tensor and introducing a new set of renormalized field variables, the original integral equation is reformulated into a system of renormalized integral equations. Dyson’s equation and its first-order smoothing approximation, describing the ensemble averaged response of the heterogeneous system, are then derived with the aid of Feynman’s diagram technique. The dispersion equations for the longitudinal and transverse coherent waves are then obtained by applying Fourier transform to the Dyson equation. The exact solution to the dispersion equations are obtained numerically. To validate the new model, the results for weak-property-fluctuation materials are compared to the predictions given by an improved weak-fluctuation multiple scattering theory. It is shown that the new model is capable of giving a more robust and accurate prediction of the dispersion behavior of weak-property-fluctuation materials. Numerical results further show that the new model is still able to provide accurate results for strong-property-fluctuation materials while the weak-fluctuation model is completely failed. As applications of the new model, dispersion and attenuation curves for coherent waves in the Earth’s lithosphere, the porous and two-phase alloys, and human cortical bone are calculated. Detailed analysis shows the model can capture the major dispersion and attenuation characteristics, such as the longitudinal and transverse wave Q-factors and their ratios, existence of two propagation modes, anomalous negative dispersion, nonlinear attenuation-frequency relation, and even the disappearance of coherent waves. Additionally, it helps gain new insights into a series of longstanding problems, such as the dominant mechanism of seismic attenuation and the existence of the Mohorovičić discontinuity. This work provides a general and accurate theoretical framework for quantitative characterization of microstructures in a broad spectrum of heterogeneous materials and it is anticipated to have vital applications in seismology, ultrasonic nondestructive evaluation and biomedical ultrasound.

Multimedia materials represent a promising approach for the Geoheritage promotion. Despite Geology is normally associated with natural environments, new tendencies are projected towards a better knowledge of the “geological reason” for the location and the development of urban settlements. The urban environment is, in fact, a perfect “laboratory” for opening the scientific topics to a broad audience. In this paper the experience of a geological exhibition organized in the town of Perugia city (Umbria, central Italy) is illustrated, highlighting the SECRET (SEe and CREaTe) for an effective dissemination activity. Panels, interactive tools, laboratories and trekking tours outside the museum represent the main activities, which counted more than eight thousand visitors in a few months. Moreover, the exhibition was the starting point for ongoing projects on geotourism in the city with important consequences in terms of visibility and financial return.

The shear failure of surface methane capture borehole (SMCB) is the main cause of shortening life cycle of SMCB but lack of lithological analysis. In order to improve the stability of SMCB and improve efficient drainage period, it is of great significance to investigate the lithology performances for shear failure of SMCB. Based on the direct shear tests and geological method, the results shows that the shear displacement increases as the grain size decreases. Mechanical jump occurs at the lithological boundaries, which is mainly determined by the composition of rock specimens. The cohesion is the mainly possible reason for the step change of shear strength. Lithology with high quartz and low clay may effectively improve shear strength and failure resistance. Boreholes drilled into the weaker siltstone and mudstone sections may potentially experience preferential damage due to the larger shear displacement and shear strength. Protective measures at these sections may improve the stability of the borehole casing. The probing data where it was found that boreholes closure validated the prediction.

Based on the geochemical analysis of the volcanic material from the sediment core AMK-340, central zone of the Reykjanes Ridge, we could detect two ash-bearing sediment units accumulated during the Termination I. They correlate to the Ash Zone I in the North Atlantic Late Quaternary sediments having an age of 12170-12840, within the Younger Dryas cold chronozone, and 13600-14540 years, within and Bølling-Allerød warm chronozone. The ash of the Younger Dryas unit is presented mostly by the mafic and persilicic material originated from the Icelandic volcanoes; Vedde Ash is presented in one sediment sample from this unit. The ash of the Bølling-Allerød unit is presented mostly by the mafic shards which are related to the basalts of the rift zone on the Reykjanes Ridge, having presumably the local origin. A detection of Vedde Ash helped to specify the timing of the previously reconstructed paleoceanographic changes for the Termination I in the point of study: a significant warming in the area could occur as early as 300 years prior to the end of the conventional Younger Dryas cold chronozone.

Currently, the use of numerical models for reproducing the landscape evolution of a river basin is part of the day-by-day research activities of fluvial engineers and geomorphologists. However, despite landscape modelling is based on a rather long tradition, and scientists and practitioners are trying to schematize the processes involved in the evolution of a landscape since decades, there is still the need for improving both the knowledge of the physical mechanisms and their numerical coding. The present review focuses on the first aspect, discussing the main components of a landscape evolution model and their more common schematizations, presenting possible open questions to be addressed towards an improvement of the reliability of such kind of models in describing the fluvial geomorphology.

Groundwater is the most extracted raw material in the world with global annual withdrawal rates of 800–1500 km3/year. In East and Southern Africa, 70 % of the population are reliant on shallow groundwater as their primary drinking water source. With increased population growth, intensification of agriculture and industrialization, conflicting demands on groundwater present a challenge to achieve the Sustainable Development Goals (6,3,11,12,15). Between 2015 and 2018, the Horn of Africa was affected by a series of climatic induced events, namely El Nino, La Nina and the Indian Ocean Diopole. These events modified the variability of rainfall patterns and resulted in long periods of low rainfall, low recharge and high evapotranspiration. As a result, shallow aquifers in alluvial deposits of Somali region have low yields and produce brackish and saline water. That situation prompted humanitarian water professionals to finance the transportation of water from selected locations with high groundwater potential through water trucks to areas facing groundwater depletion and drought. To address this challenge, UNICEF explored alternative, sustainable deeper groundwater sources that could be extracted using solar water pumping technology for multi water use. This paper describes a three-phase methodology of deep groundwater development of wells in the Ogaden Jesoma sandstone aquifers of the Somali region of the Horn of Africa to a depth of 600 meters below ground level. The results concluded that the deep sandstone aquifer of Jesoma can provide fresh water with yields of 15 l/s to the local population of Somali region. to the study provided insights into deep groundwater identification and development as well as adaptive deep boreholes drilling as a source for climate resilient water supplies.

Surface urban heat island (SUHI) depicts the deteriorating thermal environment in high-density cities and local climate zone (LCZ) classification provides a universal protocol for SUHI identification. In this study, taking the central urbanized area of Guangzhou in the humid subtropical region of China as the study area, the maps or images of LCZ, land surface temperature (LST), SUHI and urban design factors were achieved by using Landsat satellite data, GIS database and a series of retrieval and classification algorithms, and the urban design factors influencing SUHI were investigated based on 625 samples of LCZs. The results show that in the summer daytime under the clear sky condition, the LST varied greatly from 26 °C to 40 °C and the SUHI changed in a wide range of -6 °C to 8 °C in the LCZs of the study area. Seven and five urban design factors influencing the summer daytime SUHI were identified for the two dominant LCZ of LCZs 1-5 (LCZ 1 to LCZ 5) and the mixed LCZ (containing at least three types of LCZs), respectively. The summer daytime SUHI prediction models were obtained by using the step-wise multiple linear regression, with the performance of R2 of 0.697, RMSE of 1.21 °C, and the d value of 0.81 for the model of LCZs 1-5, and the values of 0.666, 1.66 °C, and 0.76 for the model of the mixed LCZ, indicating that the models can predict the changes of SUHI with LCZs to a large and satisfactory extent. This study presents a methodology to efficiently achieve a large sample of SUHI and urban design factors of LCZs in the largely urbanized cities, and provides information beneficial to the urban designs and regenerations in the humid subtropical region.

The moon always use the same side to face toward the earth, but there is a dead angle in the mainstream theory of explaining this phenomenon. That is, it cannot explain why the moon doesn't rotate around the axis which is a straight line to connect the mass centers of the earth and moon. Because the numerous meteorite impact craters on the lunar surface indicate that the moon is completely possible to obtain external momentums and rotate around this axis. This paper proposes a plain explanation, that is, the universal gravitation between the earth and moon as well as the earth's magnetic field have formed a trinity restraint mechanism on the moon. According to this explanation, the moon's rotation can be locked, and the mechanism of lunar libration has been revealed out, which can also confirm mutually with the natural phenomenon that the moon has sought a balance in the swing. In addition, with the help of all kinds of detection data from the Apollo moon landings and other circumlunar spacecraft, as well as the studies and analysis of lunar soil samples, the conclusion is that as far as a whole for the moon, it belongs to paramagnetic substances, and its relative permeability is between 1.008 and 1.03. Although the magnetic flux density of the earth on the lunar orbit has been dropped below 0.0008125 nT or lower due to the impact of the solar wind, but it can be used as a reason to lock the moon without rotating around the axis which is a straight line to connect the mass centers of the earth and moon. If another main reason to cause the existence of this fact cannot be found, even if the magnetic flux density of the geomagnetism in lunar orbit is very small, it also should not be artificially ignored. In this regard, we can artificially change the intensity of the earth's magnetic field, and carefully observe the lunar libration and in the distance between the earth and the moon, to verify the arguments in this paper.

The UNESCO World Heritage Hadrian’s Villa lies over the Colli Albani volcanic district near Rome. Magnetic, paleomagnetic, radar, and electric resistivity surveys were performed in the Plutonium–Inferi sector to detect buried buildings and outline a segment of the underground system of tunnels that link different zones of the villa. In particular, a paleomagnetic analysis of the bedrock unit allowed to accomplish an accurate geomagnetic field modelling and characterize the archaeological sources of the magnetic field anomalies. We used a computer-assisted forward modelling procedure to generate a structural model of the sources of the observed anomalies. The intrinsic ambiguity of the magnetic field modelling was reduced with the support of ground penetrating radar amplitude slices and an analysis of radar and electric resistivity profiles. The bedrock lithology in this area is an ignimbrite tuff characterized by abundant iron oxides. The high-amplitude magnetic anomalies observed in the Plutonium–Inferi area are due to strong bedrock remnant magnetization and susceptibility contrasts between topsoil infill of cavities and the surrounding tuff. The resulting magnetization model of the Plutonium–Inferi complex shows that the observed anomalies are mostly due to the presence of tunnels, skylights and a system of ditches excavated in the tuff.

Portugal Mainland has hundreds of thousands of people living in the Atlantic coastal zone, with numerous high economic value activities and a high number of infrastructures that must be protected from natural coastal hazard, namely extreme storms and sea level rise (SLR). In the context of climate change adaptation strategies, a reliable and accurate assessment of the physical vulnerability to SLR is crucial. This study is a contribution to the implementation of flooding standards imposed by the European Directive 2007/60/EC, which requires each member state to assess the risk associated to SLR and floods caused by extreme events. Therefore, coastal hazard in the Continental Atlantic coast of Portugal Mainland was evaluated for 2025, 2050 and 2100 in the whole coastal extension with different sea level scenarios for different extreme event return periods and due to SLR. A coastal flooding probabilistic map was produced based on the developed methodology using Geographic Information Systems (GIS) technology. The Extreme Flood Hazard Index (EFHI) was determined on flood probabilistic bases through five probability intervals of 20% of amplitude. For a given SLR scenario, the EFHI is expressed, on the probabilistic flooding maps for an extreme tidal maximum level, by five hazard classes ranging from 1 (Very Low) to 5 (Extreme).

The controls of depositional environments on reservoir quality have been evaluated in terms of porosity and permeability of the Gabo Field, Niger Delta, Nigeria. Data used in this research include Well logs, Core data and photos, and grain size analysis for Wells 51 and 52 in the study area. Standard methods as applicable in petrophysical and sedimentological analysis has been adopted. Thirteen reservoir units have been identified in wells 51 and 52 which had 5 reservoirs cored each. The lithofacies units of the identified reservoirs across the study area, comprise pebbly sands, coarse -, medium -, fine- and very fine-grained sands, sandy mud, silty sands and heteroliths. The heteroliths – very fine-grained silty muds are highly bioturbated. Ophiomorpha and skolithos are the major trace fossils with sedimentary structures (ripple lamination, wavy lenticular and planar beds, cross bedded sands, coarsening and fining upward). The facies associations interpreted for the study area are Channel and Coastal barrier systems and the environment of deposition as distributary channel, upper and lower shoreface. The sedimentary processes that deposited facies ranged from high energy regimes, reworking by waves to low energy with periodic influx of silts and muds. The average porosity and permeability for reservoirs in Well 51 is 16.7% and 1317 Md, reservoirs in Well 52 is 28.2% and 2330Md whereas porosity range for the study area is 2% - 32% and permeability is 1.2 – 10600 Md. The reservoir quality reservoir of the sand units in Well 51 (7, 9 and 13) and Well 52 (5, 7, 9, 11 and 13) is excellent - good, this is because of the dynamics environments of deposition (upper shoreface and distributary channel) as well as the mechanisms that play out during deposition such as bioturbation, sorting, sedimentary structures formed. Whereas the poor quality across the reservoirs especially the lower shoreface and prodelta facies is as result of lack bioturbation, connectivity, multiplicity of burrows that may have been plugged by clay and intercalation of shale and sand (heteroliths). This research has shown that environments of deposition have direct influence the reservoir quality in terms of porosity and permeability.

To estimate The potential inventory of natural gas hydrates in the Levant Basin we correlated the gas hydrate stability zone (GHSZ), modeled with locally estimated thermodynamic parameters, with seismic indicators of gas. Compilation of oceanographic measurements define the deep-water temperature and salinity to 13.8°C and 38.8‰ respectively, predicting the top GHSZ at a water depth of 1250±5 m. Assuming beneath the seafloor a hydrostatic pore-pressure, the water body salinity, and geothermal gradients ranging between 20 to 28.5°C/km, yields a useful first-order base-GHSZ approximation. Our model predicts that the entire northwestern half of the Levant Basin lies within the GHSZ, with a median thickness of ~150 m. High amplitude seismic reflectivity (HASR) imaged on an extensive 3D seismic dataset, consistently correlates with verified active seafloor gas seepage and is pervasively distributed across the deep-sea fan of the Nile within the Levant. Two main trends observed for the distribution of HASR are suggested to represent: (1) shallow gas and possibly hydrates, within buried channel-lobe systems 25 to 100 m beneath the seafloor; and (2) a regionally discontinuous bottom simulating reflection (BSR) broadly matching the modeled base GHSZ. We therefore estimate the potential methane hydrates reserve within the Levant Basin at ~4 Tcf.

The signal-to-noise ratio (SNR) data is part of the global navigation satellite systems (GNSS) observables. In a marine environment, the oscillation of the SNR data can be used to derive reflector heights. Since the attenuation of the SNR oscillation is related to the roughness of the sea surface, the significant wave height (SWH) of the water surface can be calculated from the analysis of the attenuation. The attenuation depends additionally on the relation between the coherent and the incoherent part of the scattered power. The latter is a function of the correlation length of the surface waves. Because the correlation length changes with respect to the direction of the line of sight relative to the wave direction, the attenuation must show an anisotropic characteristic. In this work, we present a method to derive the wave direction from the anisotropy of the attenuation of the SNR data. The method is investigated based on simulated data as well by the analysis of experimental data from a GNSS station in the North Sea.

Gas hydrate provinces are present in at least in two areas along Brazil’s continental margin: (1) the Rio Grande Cone in the southeast, and (2) the Amazon deep-sea fan in the equatorial region. The occurrence of gas hydrates in these depocentres was first detected geophysically and has recently been proven by seafloor sampling of gas vents, detected as water column acoustic anomalies rising from seafloor depressions (pockmarks) and/or mounds, many associated with seafloor faults. The gas vents include typical features of cold seep systems, including shallow sulphate reduction depths (<4 m), authigenic carbonate pavements and chemosynthetic ecosystems. In both areas, gas sampled in hydrate and in sediments is dominantly formed by biogenic methane. Calculation of the methane hydrate stability zone for water temperatures in the two areas shows that gas vents occur along its feather edge (water depths between 510-760 m in the Rio Grande Cone and 500-670 m in the Amazon deep-sea fan) but also in deeper waters within the stability zone. Gas venting along the feather edge of the stability zone could reflect gas hydrate dissociation and release to the oceans, as inferred on other continental margins, or upward fluid flow through the stability zone facilitated by tectonic structures recording the gravitational collapse of both depocentres. The potential quantity of venting gas on the Brazilian margin under different scenarios of natural or anthropogenic change require further investigation. The studied areas provide a natural laboratory where these critical processes can be analysed and quantified.

Knowledge about the liquefaction vulnerability in Pariaman city which is prone to an earthquake is very much needed in disaster mitigation based spatial planning. The liquefaction is an event of loss of the strength of the sandy soil layer caused by the vibration of the earthquake, where the liquefaction occurs in the sandy soil layer which has loose material in the form of sand that is not compact or not solid. This research was conducted by analyzing the potential of liquefaction vulnerability based on the Conus penetration to produce Microzonation of the susceptibility of subsidence due to liquefaction at 4 locations in Pariaman city, i.e., Marunggi village, Taluak village, Pauh Timur village and Padang Birik-Birik village. The Conus penetration testing is carried out using the piezocone (CPTU) method and mechanical Cunos penetration, and approach using Geographic Information System (GIS). The results showed that the potential of liquefaction was found at the sandy soil layer of sand and a mixture of sand and silt, which is characterized by the value of Cunos resistance and local resistance each smaller than 15 MPa and 150 kPa at varying depths. Based on the results of the analysis using this method, the critical conditions of liquefaction found in the medium sandy soil to solid. The fine sand layer which has the potential for liquefaction is in sand units formed from coastal deposits, coastal ridges and riverbanks. This liquefaction vulnerability zones analysis is limited to a depth of 6.00 m due to the limitations of the equipment used. The results of the analysis show that the fine sand layer which has the potential for liquefaction occurs at a depth of> 1.00-6.00 m with the division of zones, i.e., 1) High liquefaction in the sandy soil layer which has a critical acceleration (a) <0.10 g with shallow groundwater surface; 2) Medium liquefaction in the sandy soil layer which has a critical acceleration (a) between 0.10–0.20 g with shallow groundwater surface; and 3) Low and very low liquefaction in the sandy soil layer which has a critical acceleration (a) between 0.20–0.30 g with an average groundwater deep enough surface.

From 1964 to the early 1970s, Orson Anderson led a research program at the Lamont Geological Observatory in the newly-emerging field of “mineral physics”. In collaboration with colleagues Edward Schreiber and Naohiro Soga, Orson exploited the techniques of physical acoustics to study the behavior of the sound velocities of minerals at elevated pressures and temperatures. This research program also included semi-empirical studies of relationships between the bulk modulus and the molar volume of solids and the use of lattice dynamics to calculate the elastic moduli of cubic structures as a function of pressure to predict instabilities, as well as theoretical investigations of the Lagrangian vs. Eulerian formulations of finite strain equations of state.

Here are analysed data taken in two hydrocarbon fields ("A" and "B"), located in the western part of Sava Depression (North Croatia). They are in the secondary phase of production. The selected reservoirs "L" (in the “A” Field) and "K" (“B”) are of the Lower Pontian (Upper Miocene) age and belong to Kloštar-Ivanić Formation. Due to strong tectonics, there are numerous tectonic block, relatively rarely sampled with well and laboratory tests. Here are selected two variables for interpolation - reservoirs permeabilities and the injected volumes of field water. The following interpolation methods are described, compared and applied: Nearest Neighbourhood, Natural Neighbour (the first time in the Sava Depression) and Inverse Distance Weighting. The last one has been proven as the most appropriate for datasets with size lower than 20 points.

The Middle East spans the transition between temperate Mediterranean climate in the Levant to hyper-arid sub-tropical deserts in the southern part of the Arabian Peninsula, with the complex alpine topography in the northeast feeding the Euphrates and Tigris rivers which support life in the southeastern Fertile Crescent (FC). Climate projections predict severe drying in major parts much of the ME in response to global warming, making it important to understand the controls of hydro-climate perturbations in the region. Here we discuss 23 ME speleothem stable oxygen isotope (δ¹⁸Occ) records from 16 sites from the SISAL_v1 database, which provide a record of past hydro-climatic variability. Sub-millennial changes in ME speleothem δ¹⁸Occ values primarily indicate changes in past precipitation amounts superimposed on variations of the main synoptic pattern in the region, specifically Mediterranean cyclones. The coherency (or lack thereof) between regional records is reviewed from Pleistocene to present, covering the Last Glacial Maximum (LGM), prominent events during deglaciation, and transition into the Holocene. The available speleothem δ¹⁸Occ time-series are investigated by binning and normalizing at 25-year and 200-year time windows over the Holocene. Important Holocene climatic oscillations are discussed, such as the 8.2 ka, 4.2 ka and 0.7 ka (the Little Ice Age) Before Present events. Common trends in the standardized anomalies are tested against different climate archives. Finally, recommendations for future speleothem-based research in the region are given along with comments on the utility and completeness of the SISAL database.

The mantle mass flow interconnected with the process of Glacial Isostatic Adjustment (GIA) and the reformation of the Earth’s crust constantly perturbs the observed gravity field towards a hypothetic isostatic state. We analyse the temporal changes of the gravity field from the GRACE data, using different mathematical and/or statistical methods to detect the GIA amidst other gravity signals. A number of gravimetric post-glacial land uplift rate (LUR) modelling methods are investigated and compared with the data from a total number of 515 GPS stations and preferred GIA forward models in Fennoscandia and North America. We investigate three mathematical methods, namely regression, principal component, and independent component analysis (ICA) to extracting the GIA signal from the GRACE monthly geoid heights. We use some regularization techniques to exploit the GRACE monthly data to their maximum spatial resolution and to increase the Signal to Noise Ratio of their short wavelengths. Near the centres of the study areas the gravimetric LUR model using the fast-ICA algorithm of Hyvärinen and Oja (2000) is shown to be in a complete agreement with the GPS data and the predictions of the GIA forward models, and for the whole areas, subject to epeirogeny movement of the two regions, their discrepancies reach to the extrema at -1.8 and +3.3, and -4.5 and +7.5 mm/a, respectively. We show that the largest discrepancies between the gravimetric model using the ICA method and the GIA forward model, occur for the sub-regions likely collocated with strong ice mass change signals.

The Koka gold deposit is located in the Elababu shear zone between the Nakfa terrane and the Adobha Abiy terrane, NW Eritrea. Based on the paragenetic study two main stages of gold mineralization were identified in the Koka gold deposit: 1) an early stage of pyrite-chalcopyrite-sphalerite-galena-gold-quartz vein; and 2) a second stage of pyrite-quartz veins. NaCl-aqueous inclusions, CO2-rich inclusions, and three-phase CO₂-H₂O inclusions occur in the quartz veins at Koka. The ore-bearing quartz veins formed at 268℃, from NaCl-CO₂-H₂O(-CH₄) fluids averaging 5 wt% NaCl eq. The ore-forming mechanisms include fluid immiscibility during stage I, and mixing with meteoric water during stage II. Oxygen, hydrogen and carbon isotopes suggest that the ore-forming fluids originated as mixtures of metamorphic water, meteoric water and magmatic water, whereas sulfur isotope suggest an igneous origin. Features of geology and ore-forming fluid at Koka deposit are similar to those of orogenic gold deposits, suggesting the Koka deposit might be an orogenic gold deposit related to granite.

Presented is semivariogram and the Ordinary Kriging analyses of porosity data from the Sava Depression (Northern Croatia), as part of the Croatian part of the Pannonian Basin System. Data are taken from hydrocarbon reservoirs of the Lower Pontian (Upper Miocene) age, which belongs to the Kloštar-Ivanić Formation. Original datasets had been jack-knifed with purpose to “artificially” increased data and calculate the more reliable semivariograms. The results showed that such improvements can assist in the interpolation of more reliable maps. The both sets, made by original and jack-knifed data, need to be compared using geological recognition of non-allowed shapes (“bull-eyes”, “butterfly effects”) as well as cross-validation results. That comparison made possible to select the most appropriate porosity interpolation.

The detection of seismic events at regional and teleseismic distances is critical to Nuclear Treaty Monitoring. Traditionally, detecting regional and teleseismic events has required the use of an expensive multi-instrument seismic array; however in this work, we present DeepPick, a novel seismic detection algorithm capable of array-like performance from a single trace. We achieve this directly, by training our single-trace detector against labeled events from an array catalog, and by utilizing a deep temporal convolutional neural network. The training data consists of all arrivals in the International Seismological Centre Catalog for seven seismic arrays over a five year window from 1 Jan 2010 to 1 Jan 2015, yielding a total training set of 608,362 detections. The test set consists of the same seven arrays over a one year window from 1 Jan 2015 to 1 Jan 2016. We report our results by training the algorithm on six of the arrays and testing it on the seventh, so as to demonstrate the transportability and generalization of the technique to new stations. Detection performance against this test set is outstanding. Fixing a type-I error rate of 1%, the algorithm achieves an overall recall rate of 73% on the 141,095 array beam picks in the test set, yielding 102,394 correct detections. This is more than 4 times the 23,259 detections found in the analyst-reviewed single-trace catalogs over the same period, and represents an 8dB improvement in detector sensitivity over current methods. These results demonstrate the potential of our algorithm to significantly enhance the effectiveness of the global treaty monitoring network.

This paper introduces an application of R programming language for geostatistical data processing with a case study of the Mariana Trench, Pacific Ocean. The formation of the Mariana Trench, the deepest among all hadal oceanic depth trenches, is caused by complex and diverse geomorphic factors affecting its development. Mariana Trench crosses four tectonic plates: Mariana, Caroline, Pacific and Philippine. The impact of the geographic location and geological factors on its geomorphology has been studied by methods of statistical analysis and data visualization using R libraries. The methodology includes following steps. Firstly, vector thematic data were processed in QGIS: tectonics, bathymetry, geomorphology and geology. Secondly, 25 cross-section profiles were drawn across the trench. The length of each profile is 1000-km. The attribute information has been derived from each profile and stored in a table containing coordinates, depths and thematic information. Finally, this table was processed by methods of the statistical analysis on R. The programming codes and graphical results are presented. The results include geospatial comparative analysis and estimated effects of the data distribution by tectonic plates: slope angle, igneous volcanic areas and depths. The innovativeness of this paper consists in a cross-disciplinary approach combining GIS, statistical analysis and R programming.

The Malyy (Little) Murun massif of the Aldan Shield of the Siberian Craton has long been a kind of the geologists’ Siberian Mecca. It attracted thousands of geologists, prospectors and mineral collectors despite of its remote location. It is famous for a dozen of new and rare minerals, including gemstones charoite and dianite (the latter is the market name for strontian potassicrichrerite), as well as for specific alkaline igneous rocks. Despite of this, the age of the Malyy Murun igneous complex and associated metasomatic and hydrothermal mineral associations remained poorly constrained. In this paper, we provide extensive 40Ar/39Ar geochronological data to reveal its age and temporal history. It appeared that being unique in terms of rocks and constituent minerals, the Malyy Murun is one of multiple alkaline massifs and lavas emplaced in Early Cretaceous (~137-128 Ma) within a framework of extensional setting of the Aldan Shield and nearby Transbaikalian region. The extension has occurred 40-60 million years after the supposed closure of the Mongolia-Okhotsk Ocean and peak of orogeny in Early-Middle Jurassic.

Computation of Common Middle Point seismic sections and their subsequent time migration and diffraction imaging provides very important knowledge about the internal structure of 3D heterogeneous geological media and are key elements for successive geological interpretation. Full-scale numerical simulation, that computes all single shot seismograms, provides a full understanding of how the features of the image reflect the properties of the subsurface prototype. Unfortunately, this kind of simulations of 3D seismic surveys for realistic geological media needs huge computer resources, especially for simulation of seismic waves’ propagation through multiscale media like cavernous fractured reservoirs. In order to significantly reduce the query of computer resources we propose to model these 3D seismic cubes directly rather than the shot-by-shot simulation with subsequent CMP stacking. To do that we modify the well known "exploding reflectors principle" for 3D heterogeneous multiscale media by use of the finite-difference technique on the base of grids locally refined in time and space. To be able to simulate realistic models and acquisition we develop scalable parallel software, which needs reasonable computational costs. Numerical results for simulation of Common Middle Points sections and their time migration are presented and discussed.

Geologists employ high-spatial-resolution (HR) remote sensing (RS) data for many diverse applications as they effectively reflect detailed geological information, enabling high-quality and efficient geological surveys. Applications of HR RS data to geological and related fields have grown recently. By analyzing these applications, we can better understand the results of previous studies and more effectively use the latest data and methods to efficiently extract key geological information. HR optical remote sensing data are widely used in geological hazard assessment, seismic monitoring, mineral exploitation, glacier monitoring, and mineral information extraction due to high accuracy and clear object features. Compared with optical satellite images, synthetic-aperture radar (SAR) images are stereoscopic and exhibit clear relief, strong performance, and good detection of terrain, landforms, and other information. SAR images have been applied to seismic mechanism research, volcanic monitoring, topographic deformation, and fault analysis. Furthermore, a multi-standard maturity analysis of the geological applications of HR images using literature from the Science Citation Index reveals that optical remote sensing data are superior to radar data for mining, geological disaster, lithologic, and volcanic applications, but inferior for earthquake, glacial, and fault applications. Therefore, geological remote sensing research needs to be truly multidisciplinary or interdisciplinary, ensuring more detailed and efficient surveys through cross-linking with other disciplines. Moreover, the recent application of deep learning technology to remote sensing data extraction has improved automatic processing and data analysis capabilities.

The acquisition of high resolution morpho-bathymetric data on the Calabro Tyrrhenian continental margin (Southern Italy) enabled us to identify several mass-wasting processes, including shallow gullies, shelf-indenting canyons and landslides. In particular, we focus our attention on submarine landslides occurring from the coast down to -1700 m, with mobilized volumes ranging from some hundreds up to tens of millions of cubic meters. These landslides also show a large variability of geomorphic features in the headwall, translational and toe domain. Based on their morphology and distribution, four main classes of coastal/submarine landslides have been recognized: a) rocky coastal/shallow-water failures characterized by large hummocky deposits offshore; b) large-size and isolated scars with associated landslide deposits, mostly occurring on open slope environment and lower part of tectonically-controlled escarpments; c) a linear array of coalescent and nested landslide scars occurring in the upper part of tectonically-controlled escarpments and canyon flanks; d) a cauliflower array of small and coalescent scars occurring in canyon headwall. The latter two classes of landslides are also characterized by a marked retrogressive evolution and their landslide deposits are generally not recognizable on the morpho-bathymetric data. By integrating the morpho-bathymetric dataset with the results of previous studies, we also discuss the main factors controlling the variability in size and morphology of these submarine landslides to provide insights on their failure and post-failure behavior.

Carbonate precipitation induced by microorganism has become a hot spot in the field of carbonate sedimentology, while the effect of different magnesium on biominerals has rarely been studied. Therefore, magnesium sulfate and magnesium chloride were used to investigate the significant role played on carbonate minerals. In this study, Staphylococcus epidermidis Y2 was isolated and identified by 16S rDNA homology comparison. The ammonia, pH, carbonic anhydrase, carbonate and bicarbonate ions were investigated. The mineral phase, morphology and elemental composition were analyzed by XRD and SEM-EDS. The ultrathin slices of bacteria were analyzed by HRTEM-SAED and STEM. The result showed that this bacterium could release ammonia and carbonic anhydrase to increase pH, and elevate the supersaturation via a large number of carbonate and bicarbonate ions released through carbon dioxide hydration catalyzed by carbonic anhydrase. The crystal cell density of monohydrocalcite was lower in magnesium chloride medium than that in magnesium sulfate medium. The crystal grew in a mode of spiral staircas in magnesium sulfate medium, while in a concentric circular pattern in magnesium chloride medium. There was no obvious intracellular biomineralization. This study may be helpful to further understand the biomineralization mechanism, may also provide some references for the reconstruction of paleogeological environment.

A finite-difference approach with non-uniform meshes was presented for simulating magnetotelluric responses in 2D structures. We presented the formulation of this scheme and gave some sights into its successful implementation, and compared finite-difference solution with known numerical results and simple analytical solutions. First, a homogeneous half-space model was tested and the finite-difference approach can provide very good accuracy for 2D magnetotelluric modeling. Then we compared to the analytical solutions for the two-layered model, the relative errors of the apparent resistivity and the impedance phase were both increased when the frequency was increased. In the end, we compare our finite-difference simulation results with COMMEMI 2D-0 model with the finite-element solutions. Both results are in close agreement to each other. These comparisons confirm the validity and reliability of our finite-difference algorithm.

The active debate about the processes governing the organic-rich sediment deposition generally involves the relative roles of elevated primary productivity and enhanced preservation related to anoxia. However, other less spotlighted factors could have a strong impact on such deposits: e.g. residence time into the water column (bathymetry), sedimentation rate, transport behavior of organo-mineral floccules on the sea floor. They are all strongly interrelated and may be obscured in the current conceptual models inspired from most representative modern analogues (i.e. upwelling zones and stratified basins). To improve our comprehension of organic matter distribution and heterogeneities, we conducted a sensitivity analysis on the processes involved in organic matter production and preservation which have been simulated within a 3D stratigraphic forward model. The Lower-Middle Toarcian of the Paris Basin was chosen as a case study as it represents one of the best documented example of marine organic matter accumulation. The relative influence of the critical parameters (bathymetry, diffusive transport, oxygen mixing rate and primary production) on the output parameters (Total Organic Carbon, and oxygen level), determined performing a Global Sensitivity Analysis, shows that, in the context of a shallow epicontinental basin, a moderate primary productivity (> 175 gC.m^-².yr^-1) can led to local anoxia and organic matter accumulation. We argue that, regarding all the processes involved, the presence and distribution of organic-rich intervals is linked as a first-order parameter to the morphology of the basin (e.g. ramp slope, bottom topography). These interpretations are supported by very specific ranges of critical parameters which allowed to obtain output parameter values in accordance with the data. This quantitative approach and its conclusions open new perspectives about the understanding of global distribution and preservation of organic-rich sediments.

Large amounts of gas hydrate are present in marine sediments offshore Taitao Peninsula, near the Chile Triple Junction. Here, marine sediments on the forearc contain carbon that is converted to methane in a zone of very high heat flow and intense rock deformation above the downgoing oceanic spreading ridge separating the Nazca and Antarctic plates. This regime enables vigorous fluid migration. Here we present an analysis of the spatial distribution, concentration, estimate of gas phases (gas hydrate and free gas) and geothermal gradients in the accretionary prism and forearc sediments offshore Taitao (45.5° - 47° S). Velocity analysis of Seismic Profile RC2901-751 indicates gas hydrate concentration values <10% of the total rock volume, and extremely high geothermal gradients (<190 °Ckm^-1). Gas hydrates are located in shallow sediments (90-280 meters below the seafloor). The large amount of hydrate and free gas estimated (7.21x10¹¹ m³ and 4.1x10¹⁰ m³, respectively), the high seismicity, the mechanically unstable nature of the sediments, and the anomalous geothermal conditions, set the stage for potential massive releases of methane to the ocean mainly through hydrate dissociation and/or migration directly to the seabed through faults. We conclude that the Chile Triple Junction is an important methane seepage area and should be the focus of novel geological and ecological research.

Earth's global surface temperature shows variability on an extended range of temporal scales and satisfies an emergent scaling symmetry. Recent studies indicate that scale invariance is not only a feature of the observed temperature fluctuations, but an inherent property of the temperature response to radiative forcing, and a principle that links the fast and slow climate responses. It provides a bridge between the decadal- and centennial-scale fluctuations in the instrumental temperature record, and the millennial-scale equilibration following perturbations in the radiative balance. In particular, the emergent scale invariance makes it possible to infer equilibrium climate sensitivity (ECS) from the observed relation between radiative forcing and global temperature in the instrumental era. This is verified in ensembles of Earth system models (ESMs), where the inferred values of ECS correlate strongly to estimates from idealized model runs. For the range of forcing data explored in this paper, the method gives best estimates of ECS between 2.3 and 3.4 K.

Selenium-rich Au-Ag mineralization has been discovered in the Kremnica ore district, central Slovakia. The mineralization is hosted by a single quartz-dolomite vein hosted by Neogene propyllitized andesites of the Kremnica stratovolcano. Ore mineralogy and crystal chemistry of individual ore minerals have been studied here. The early base-metal ore mineralization composed of pyrite, sphalerite and chalcopyrite lacks selenium, whereas the superimposed Au-Ag paragenesis is Se-enriched. The Au-Ag alloys, uytenbogaardtite, minerals of the galena-clausthalite series, acanthite-naumannite series, diaphorite, miargyrite, pyrargyrite-proustite, polybasite group, minerals of the tetrahedrite group and andorite branch (andorite IV, andorite VI, Ag-excess fizélyite), freislebenite, and rare Pb-Sb sulphosalts (scaiinite, robinsonite, plagionite) have been identified here. Besides selenides, the most Se-enriched phases are miargyrite, proustite-pyrargyrite, and polybasite-pearceite, whose Se contents are among the highest reported worldwide. In addition, one new phase has been found, corresponding to a Se-analogue of pearceite containing 2.08 - 3.54 apfu Se. The style of mineralization, paragenetic situation, and chemical trends observed in individual minerals are comparable to those of Au-Ag low-sulphidation epithermal Au-Ag mineralizations of the Kremnica and neighboring Štiavnica and Hodruša-Hámre Ore Districts, however, the pronounced enrichment in selenium is a specific feature of the studied vein only.

Extraction of high-resolution surface waves is essential in surface-wave survey. Because reflections usually interfere with surface waves on X component in a multicomponent seismic exploration, it is difficult to extract dispersion curves of surface waves. The situation goes more serious when the frequencies and velocities of higher-mode surface waves are close to those of PS-waves. A method for surface-wave extraction is proposed based on the morphological differences between reflections and surface waves. Frequency-domain high-resolution linear Radon transform (LRT) and time-domain high-resolution hyperbolic Radon transform (HRT) are used to represent surface waves and reflections respectively. Then, the sparse representation problem based on the morphological component analysis (MCA) is built and optimally solved to obtain high-fidelity surface waves. An advantage of our method is its ability to extract surface waves when their frequencies and velocities are close to those of reflections. Furthermore, results of synthetic and field examples confirm that the proposed method can attenuate the distortion of surface-wave dispersive energy caused by reflections, which contributes to extracting accurate dispersion curves.

The Greenland ice sheet is a major contributor to sea level rise, adding an estimated 0.47 +/− 0.23 mm/yr to global mean sea level between 1991 and 2015 (van den Broeke et al., 2016). Making sea level rise projections for the future and understanding the processes controlling current observed rates of sea level rise are crucially dependent on understanding the present-day state of the ice sheet. Here, we provide an overview of the current state of the mass budget of Greenland based on satellite gravimetry and remote sensing observations of surface elevation change, ice sheet velocity and calving front positions. We also combine these essential climate variables with a regional climate model (RCM) output from an ice sheet model (ISM) to gain insight into poorly understood ice sheet dynamical and surface mass processes. On average from 1992 to 2017 the ice sheet in some locations has lost up −2.65 m/yr in elevation based on ESA Radar altimetry analysis. Calving fronts have retreated all around Greenland since the 1990s and in only two out of 28 study locations have they remained stable. The locations of grounding lines at 5 key glaciers with floating ice tongues have remained stable over the observation period. However a detailed case study at Petermann glacier with an ice fracture model shows the sensitivity of these floating ice shelves to future climate change. GRACE gravimetric mass balance (GMB) data allows us to tie together disparate lines of evidence showing that Greenland has lost about 265 +/− 25 Gt/yr of ice over the period 2002 to 2015. RCM and ISM simulations show that surface mass processes dominate the overall Greenland ice sheet mass budget except for areas of fast ice sheet flow but marked differences between models and between models and observations indicate that not all processes are captured accurately, indicating areas of greater uncertainty and directions of future research for future sea level rise projections.

Porphyry copper-molybdenum deposits (PCDs) are the world’s most important source of copper, molybdenum and rhenium. Previous studies have reported that some PCDs can have sub-economic to economic grades of critical metals, i.e., those elements that are both essential for modern societies and subject to the risk of supply restriction (e.g., platinum group elements (PGE), rare earth elements (REE), In, Co, Te, Ge, Ga, among others). Even though some studies have reported measured concentrations of Pd and Pt in PCDs, their occurrence and mineralogical form remain poorly constrained. Furthermore, these reconnaissance studies have focused predominantly on porphyry Cu-Au deposits, but very limited information is available for porphyry Cu-Mo systems. In this contribution, we report the occurrence of critical metal (Pd, Pt, Au, Ag, and Te) inclusions in copper sulfides from the world’s largest PCD, the supergiant Río Blanco-Los Bronces deposit in central Chile. Field emission scanning electron microscope (FESEM) observations of chalcopyrite and bornite from the potassic alteration zone reveal the presence of micro- to nano-sized particles (<1-10 μm) of noble metals, most notably Pd, Au, and Ag. The high-resolution data show that these inclusions are mostly tellurides, such as merenskyite [PdTe2], Pd-rich hessite [Ag2Te], sylvanite [(Ag, Au)Te2] and petzite [Ag3AuTe2]. The data point to Pd (and probably Pt) partitioning in copper sulfides during the high-temperature potassic alteration stage, opening new avenues of research aimed at investigating not only the mobility of PGE during mineralization and partitioning into sulfides, but also at evaluating the potential of porphyry Cu-Mo deposits as a source for noble metals.

This work describes and interprets the presence of heavy minerals in the western Portuguese continental margin using a set of 78 bottom samples collected from 3 distinct areas of this margin: Porto, Aveiro and Nazaré canyon head areas. The main transparent heavy mineral suite (minerals with frequencies >10%), is composed by amphiboles (hornblende), mica (biotite), andalusite, tourmaline and garnet. A secondary suite (minerals with frequencies between 1 and 10%), is composed by pyroxene (enstatite, diopside and augite), staurolite, zircon and apatite. With very low frequency representing less than 1% we found rutile, olivine, kyanite, monazite, epidote, sphene, anatase, sillimanite and brookite. The main primary sources (igneous and metamorphic rocks) explain the presence of these minerals. However, the application of the principal component analysis, with a previous application of the centered log ratio transformation of the heavy mineral data, also stresses for the importance of the grain sorting as a process in controlling the heavy mineral occurrence. The importance of this process is mostly sustained by the distribution pattern of mica and of the most flattened amphibole grains in a way that these particles tend to have a hydraulic affinity to finer grained sediments.

This paper examines the effect of the aggregate type on concrete strength and more specifically how the petrographic characteristics of various aggregate rocks as well as their physico-mechanical properties influences the durability of C 25/30 strength class concrete. The studied aggregate rocks are derived from Veria-Naousa and Edessa ophiolitic complexes as well as granodiorite and albitite rocks from their surrounding areas in central Macedonia (Greece). Concretes are produced with constant volume proportions, workability, mixing and curing conditions using different sizes of each aggregate type. Aggregates were mixed both in dry and water saturated states in concretes. Six different types of aggregates were examined and classified in three district groups according to their physicomechanical properties, petrographic characteristics and surface texture. The classification in groups after the concrete compressive strength test verified the initial classification in the same three groups. Group I (ultramafic rocks) presents the lowest concrete strengths, depending on their high alteration degree and the low mechanical properties of ultramafic aggregates. Group II (mafic rocks and granodiorite) presents a wide range of concrete strengths, depending on different petrographic characteristics and mechanical properties. Group III (albite rocks) presents the highest concrete strengths, depending on their lowest alteration degree and their highest mechanical properties. Therefore mineralogy and microstructure of the coarse aggregates affects the final strength of the concrete specimens.

We propose a probabilistic graphical model for discriminative substrate characterization, to support geological and biological habitat mapping in aquatic environments. The model, called a fully connected conditional random field (CRF), is demonstrated using multispectral and monospectral acoustic backscatter from heterogeneous seafloors in Patricia Bay, British Columbia, and Bedford Basin, Nova Scotia. Unlike previously proposed discriminative machine learning algorithms, the CRF model considers both the relative backscatter magnitudes of different substrates and their relative proximities. The model therefore combines the statistical flexibility of a machine learning algorithm with an inherently spatial treatment of the substrate. The CRF model predicts substrates such that nearby locations with similar backscattering characteristics are likely to be in the same substrate class. The degree of proximity and allowable backscatter similarity are controlled by parameters that are learned from the data. CRF model results were evaluated against a popular generative model known as a Gaussian Mixture model that doesn't include spatial dependencies, only covariance between substrate backscattering response over different frequencies. Both models are used in conjunction with sparse bed observations/samples in a supervised classification. A detailed accuracy assessment, including a leave-one-out cross-validation analysis, was performed using both models. Using multispectral backscatter, the GMM model trained on 50% of the bed observations resulted in a 75% and 89% average accuracies in Patricia Bay and Bedford Basin, respectively. The same metrics for the CRF model were 78% and 95%. Further, the CRF model resulted in a 91% mean cross-validation accuracy across four substrate classes at Patricia Bay, and a 99.5% mean accuracy across three substrate classes at Bedford Basin, which suggest that the CRF model generalizes extremely well to new data. This analysis also showed that the CRF model was much less sensitive to the specific number and locations of bed observations than the generative model, owing to its ability to incorporate spatial autocorrelation in substrates. The CRF approach therefore may prove to be a powerful `spatially aware' alternative to other discriminative classifiers.

A waterbody detection technique is an essential part of digital elevation model (DEM) generation to delineate land-water boundaries and set flattened elevations. This paper describes the technical methodology for improving the initial tile-based waterbody data that are created during production of the ASTER GDEM, because without improvement such tile-based waterbody data are not suitable for incorporating into the new ASTER GDEM Version 3. Waterbodies are classified into three categories: sea, lake, and river. For sea-waterbodies, the effect of sea ice is removed to better delineate sea shorelines in high latitude areas, because sea ice prevents accurate delineation of sea shorelines. For lake-waterbodies, the major part of the processing is to set the unique elevation value for each lake using a mosaic image that covers the entire lake area. Rivers present a unique challenge, because their elevations gradually step down from upstream to downstream. Initially, visual inspection is required to separate rivers from lakes. A stepwise elevation assignment, with a step of one meter, is carried out by manual or automated methods, depending on the situation. The ASTER GWBD product consists of a global set of 1º latitude-by-1º longitude tiles containing water body attribute and elevation data files in geographic latitude and longitude coordinates and with one arc second posting. Each tile contains 3601-by-3601 data points. All improved waterbody elevation data are incorporated into the ASTER GDEM to reflect the improved results.

We examine the mass balance of the glaciers in the Novaya Zemlya Archipelago, located in the Russian High Arctic using time series of time-variable gravity from the NASA/DLR Gravity Recovery and Climate Experiment (GRACE) mission, laser altimetry data from the NASA Ice Cloud and land Elevation Satellite (ICESat) mission, and radar altimetry data from the ESA CryoSat-2 mission. We present a new algorithm for detecting changes in glacier elevation from these satellite altimetry data and evaluate its performance in the case Novaya Zemlya by comparing the results with GRACE. We find that the mass loss of Novaya Zemlya increased from 10±5 Gt/yr over 2003-2009 to 14±4 Gt/yr over 2010-2016, with a brief period of near mass balance between 2009 and 2011. The results are consistent across the gravimetric and altimetric methods. Furthermore, the analysis of elevation change from CryoSat-2 indicates that 60\% of the mass loss occurs at low elevation, where thinning rates are highest. We also find that marine-terminating glaciers in Novaya Zemlya are thinning significantly faster than land-terminating glaciers, which indicates an important role of ice dynamics of marine-terminating glaciers. We posit that the glacier changes have been caused by changes in atmospheric and ocean temperatures. We find that the increase in mass loss after 2010 is associated with a warming in air temperatures, which increased the surface melt rates. There is no enough information on the ocean temperature at the front of the glaciers to conclude on the role of the ocean, but we posit that the temperature of subsurface ocean waters must have increased during the observation period.

Field and numerical investigations at Happisburgh, East coast of England, UK, sought to characterize beach thickness and determine geologic framework controls on coastal change. After a major failure of coastal protection infrastructure, removal of about 1 km of coastal defense along the otherwise protected cliffed coastline of Happisburgh triggered a period of rapid erosion over 20 years of ca. 140 m. Previous sensitivity studies suggested that that beach thickness plays a major role in coastal recession. These studies were limited, however, by a lack of beach volume data. In this study, we have integrated the insights gained from our understanding of the Quaternary geology of the area, a novel non-intrusive passive seismic survey method, and a 3D novel representation of the subsurface source and transportable material into a coastal modelling environment to explore the role of beach thickness on the back wearing and down wearing of the cliffs and consolidated platform, respectively. Results show that beach thickness is non-homogeneous along the study site: we estimate that the contribution to nearshore sediment budget via platform down wearing is of a similar order of magnitude as sediment lost from the beach and therefore non-negligible. We have provided a range of evidence to support the idea that the Happisburgh beach is a relatively thin layer perched on a sediment rich platform of sand and gravel. This conceptualization differs from previous publications, which assumes that the platform was mostly till and fine material. This has direct implication on the regional sediment management of along this coastline. The present study contributes to our understanding of a poorly known aspect of coastal sediment budgeting and outlines a quantitative approach that allows for simple integration of geological understanding for coastline evolution assessments worldwide.

Subsurface temperature is the key parameter in geothermal exploration. An accurate estimation of the reservoir temperature is of high importance and usually done either by interpolation of borehole temperature measurement data or numerical modeling. However, temperature measurements at depths which are of interest for deep geothermal applications (usually deeper than 2 km) are generally sparse. A pure interpolation of such sparse data always involves large uncertainties and usually neglects knowledge of the 3D reservoir geometry or the rock and reservoir properties governing the heat transport. Classical numerical modeling approaches at regional scale usually only include conductive heat transport and do not reflect thermal anomalies along faults created by convective transport. These thermal anomalies however are usually the target of geothermal exploitation. Kriging with trend does allow including secondary data to improve the interpolation of the primary one. Using this approach temperature measurements of depths larger than 1,000 m of the federal state of Hessen/Germany have been interpolated in 3D. A 3D numerical conductive temperature model was used as secondary information. This way the interpolation result reflects thermal anomalies detected by direct temperature measurements as well as the geological structure. This results in a considerable quality increase of the subsurface temperature estimation.

Western Europe is the region with the highest density of published speleothem δ18O (δ18Ospel) records worldwide. Here we review these records in light of the recent publication of the Speleothem Isotopes Synthesis and Analysis (SISAL) database [1]. We investigate how representative the spatial and temporal distribution of the available records is for climate in Western Europe, and review potential sites and strategies for future studies. We show that spatial trends in precipitation δ18O are mirrored in the speleothems, providing means to better constrain the factors influencing δ18Ospel at a specific location. Coherent regional δ18Ospel trends are found over stadial-interstadial transitions of the last glacial, especially in high altitude Alpine records. Over the Holocene, regional trends are less clearly expressed, due to lower signal-to-noise ratios in δ18Ospel, but can potentially be extracted with the use of statistical methods. Overall, this first assessment highlights the potential of the European region for speleothem palaeoclimate reconstruction, while underpinning the importance of knowledge of local factors for a correct interpretation of δ18Ospel.

We discuss the influence played by different statistical models in the prediction of porosity and litho-fluid facies from logged and post-stack inverted acoustic impedance (Ip) values. We compare the inversion and classification results obtained under three different a-priori statistical assumptions: an analytical Gaussian distribution, an analytical Gaussian-mixture model and a non-parametric mixture distribution. The first model assumes Gaussian distributed porosity and Ip values, thus neglecting their facies-dependent behaviour caused by different lithologic and saturation conditions. Differently, the other two statistical models relate each component of the mixture to a specific litho-fluid facies, so that the facies-dependency of porosity and Ip values is taken into account. Blind well tests are used to validate the final predictions, whereas the analysis of the maximum-a-posteriori (MAP) solutions, the coverage ratio and the contingency analysis tools are used to quantitatively compare the inversion outcomes. This work points out that the correct choice of the statistical petrophysical model could be crucial in reservoir characterization studies. Indeed, for the investigated zone it turns out that the simple Gaussian model constitutes an oversimplified assumption, while the two mixture models provide more accurate results, although the non-parametric one yields slightly superior predictions with respect to the Gaussian-mixture assumption.

The application of geophysical methods for the characterization of geomorphologic subsurface features has become very popular in recent years. Convential geophysical applications on undifferentiated and weathered slopes have been less effective. Multichannel analysis of surface waves (MASW) is a new hope for the solution of the problem. MASW methods was used to the first step in order to reconstruct the ancient eruption of Mount Samalas. MASW was applied to the remnants of ancient volcanic sediment on Lombok Island to determine the physical and physical characteristics of seismic subsurface features. This paper is also discussion the application's capabilities of the MASW method on undifferentiated and weathered volcanic deposits on the slope area and classified them. The results of the research show that MASW technique is powerful in characterizing and delineating the stratum of the layers of undifferentiated and weathered volcanic precipitation areas. MASW is able to predict subsurface structures, especially for shallow and thin layers and presenting site class information for internal structure (vertical resolution). Based on S wave velocity structure, undifferentiated and weathered volcanic layer, genes are classified into site class C, D, and E. In the future, this method will completed with combination of other geophysical method and parameters variation (frequency and spaces) for isopach mapping of volcanic sediment and then analysis of ancient volcanic eruption dynamics.