Volume 9, Issue 1

We examined the petrography, petrology, and geochemistry of two types of gabbroic xenoliths (A- and B-type xenoliths) in olivine basalt and biotite rhyolite units among the dominantly rhyolitic rocks in Niijima volcano, northern Izu-Bonin volcanic arc, central Japan. A-type gabbroic xenoliths consisting of plagioclase, clinopyroxene, and orthopyroxene with an adcumulate texture were found in both olivine basalt and biotite rhyolite units, and B-type gabbroic xenoliths consisting of plagioclase and amphibole with an orthocumulate texture were found only in biotite rhyolite units. Geothermal- and barometricmodelling based on mineral chemistry indicated that the A-type gabbro formed at higher temperatures (899–955°C) and pressures (3.6–5.9 kbar) than the B-type gabbro (687–824°C and 0.8–3.6 kbar). These findings and whole-rock chemistry suggest different parental magmas for the two types of gabbro. The A-type gabbro was likely formed from basaltic magma, whereas the B-type gabbro was likely formed from an intermediate (andesitic) magma. The gabbroic xenoliths in erupted products at Niijima volcano indicate the presence of mafic to intermediate cumulate bodies of different origins at relatively shallower levels beneath the dominantly rhyolitic volcano.

This paper is devoted to the application of the Combined Resistivity Sounding and Profiling electrode configuration (CRSP) to detect underground cavities. Electrical resistivity surveying is among the most favorite geophysical methods due to its nondestructive and economical properties in a wide range of geosciences. Several types of the electrode arrays are applied to detect different certain objectives. In one hand, the electrode array plays an important role in determination of output resolution and depth of investigations in all resistivity surveys. On the other hand, they have their own merits and demerits in terms of depth of investigations, signal strength, and sensitivity to resistivity variations. In this article several synthetic models, simulating different conditions of cavity occurrence, were used to examine the responses of some conventional electrode arrays and also CRSP array. The results showed that CRSP electrode configuration can detect the desired objectives with a higher resolution rather than some other types of arrays. Also a field case study was discussed in which electrical resistivity approach was conducted in Abshenasan expressway (Tehran, Iran) U-turn bridge site for detecting potential cavities and/or filling loose materials. The results led to detect an aqueduct tunnel passing beneath the study area.

Vesuvius is a high-risk volcano and the 1631 Plinian eruption is a reference event for the next episode of explosive unrest. A complete stratigraphic and petrographic description of 1631 pyroclastics is given in this study. During the 1631 eruption a phonolite was firstly erupted followed by a tephritic phonolite and finally a phonolitic tephrite, indicating a layered magma chamber. We suggest that phonolitic basanite is a good candidate to be the primitive parental-melt of the 1631 eruption. Composition of apatite from the 1631 pyroclastics is different from those of CO 2 -rich melts indicating negligible CO 2 content during magma evolution. Cross checking calculations, using PETROGRAPH and PELE software, accounts for multistage evolution up to phonolite starting from a phonolitic basanite melt similar to the Vesuvius medieval lavas. The model implies crystal settling of clinopyroxene and olivine at 6 kbar and 1220°C, clinopyroxene plus leucite at a pressure ranging from 2.5 to 0.5 kbar and temperature ranging from 1140 to 940°C. Inside the phonolitic magma chamber K-feldspar and leucite would coexist at a temperature ranging from from 940 to 840°C and at a pressure ranging from 2.5 to0.5 kbar. Thus crystal fractionation is certainly a necessary and probably a sufficient condition to evolve the melt from phono tephritic to phonolitic in the 1631 magma chamber. We speculate that phonolitic tephrite magma refilling from deeper levels destabilised the chamber and triggered the eruption, as testified by the seismic precursor phenomena before 1631 unrest.

Earthquake activities along the Sumatra-Andaman Subduction Zone (SASZ) were clarified using the derived frequency-magnitude distribution in terms of the (i) most probable maximum magnitudes, (ii) return periods and (iii) probabilities of earthquake occurrences. The northern segment of SASZ, along the western coast of Myanmar to southern Nicobar, was found to be capable of generating an earthquake of magnitude 6.1–6.4 M w in the next 30–50 years, whilst the southern segment of offshore of the northwestern and western parts of Sumatra (defined as a high hazard region) had a short recurrence interval of 6-12 and 10-30 years for a 6.0 and 7.0 M w magnitude earthquake, respectively, compared to the other regions. Throughout the area along the SASZ, there are 70– almost 100% probabilities of the earthquake with M w up to 6.0 might be generated in the next 50 years whilst the northern segment had less than 50% chance of occurrence of a 7.0 M w earthquake in the next 50 year. Although Rangoon was defined as the lowest hazard among the major city in the vicinity of SASZ, there is 90% chance of a 6.0 M w earthquake in the next 50 years. Therefore, the effective mitigation plan of seismic hazard should be contributed.

Negative effects of tourism development in a destination are usually the consequence of the high concentration of tourists, accommodation facilities and the activities that are practiced in a relatively restricted area. One of the most important measures to protect the areas is to calculate the maximum number of tourists that can simultaneously reside in a region, i.e. the determination of the carrying capacity. This paper outlines a method for determining carrying capacity based on zoning of environmental resources and zoning within a region. The paper argues for a return to the idea of identifying maximum appropriate number of users. The main hypothesis of the paper is based on the statement that the development of tourism in Fruška Gora (Mountain) National Park in Northern Serbia must be in accordance with the basic principles of sustainability, including the determination of carrying capacity. The main research goal was to show the opinion of local residents about the uncontrolled development of tourism, and to determine the carrying capacity in four sports and recreational zones of the mountain. The carrying capacity of the area is calculated by Lavery and Stanev formulas.

Extraction of knickpoints or knickzones from a Digital Elevation Model (DEM) has gained immense significance owing to the increasing implications of knickzones on landform development. However, existing methods for knickzone extraction tend to be subjective or require time-intensive data processing. This paper describes the proposed Knickzone Extraction Tool (KET), a new raster-based Python script deployed in the form of an ArcGIS toolset that automates the process of knickzone extraction and is both fast and more user-friendly. The KET is based on multi-scale analysis of slope gradients along a river course, where any locally steep segment (knickzone) can be extracted as an anomalously high local gradient. We also conducted a comparative analysis of the KET and other contemporary knickzone identification techniques. The relationship between knickzone distribution and its morphometric characteristics are also examined through a case study of a mountainous watershed in Japan.

This study explores the influence of stereoscopic (real) 3D and monoscopic (pseudo) 3D visualization on the human ability to reckon altitude information in noninteractive and interactive 3D geovisualizations. A two phased experiment was carried out to compare the performance of two groups of participants, one of them using the real 3D and the other one pseudo 3D visualization of geographical data. A homogeneous group of 61 psychology students, inexperienced in processing of geographical data, were tested with respect to their efficiency at identifying altitudes of the displayed landscape. The first phase of the experiment was designed as non-interactive, where static 3D visual displayswere presented; the second phase was designed as interactive and the participants were allowed to explore the scene by adjusting the position of the virtual camera. The investigated variables included accuracy at altitude identification, time demands and the amount of the participant’s motor activity performed during interaction with geovisualization. The interface was created using a Motion Capture system, Wii Remote Controller, widescreen projection and the passive Dolby 3D technology (for real 3D vision). The real 3D visual display was shown to significantly increase the accuracy of the landscape altitude identification in non-interactive tasks. As expected, in the interactive phase there were differences in accuracy flattened out between groups due to the possibility of interaction, with no other statistically significant differences in completion times or motor activity. The increased number of omitted objects in real 3D condition was further subjected to an exploratory analysis.

In this paper we analyse the dependence between gully networks, networks of valleys and faults in the area of Hronska pahorkatina Hill Land. The work is based on the analysis of directions and densities of these three networks in the study area and subunits of lower hierarchical level. The coincidence of all three networks is rare. This scenario occurs when networks of gullies are situated on the bottom of shallow valleys which are conditioned by the presence of faults. More often scenario is only the coincidence of network of valleys and gullies. The last scenario appears in areas with low fault density. The most specific scenario is perpendicularity between network of gullies and network of gullies. Gullies are situated on steep slopes of incised valleys in this case. The last scenario appears the most frequently and was also proven by findings from other studies. We propose three comprehensive explanations of the possible dependence between network of gullies and faults. We also suggest the draft of the possible dependence between network of valleys and gullies.

Nine geophagic material samples were analysed in order to characterise their mineralogical and elemental constituents, and infer on their health threats. Most abundant mineral in the samples was quartz, followed by smectite, kaolinite and muscovite in minor; and microcline, plagioclase, and goethite in trace quantities. Dominant major oxides were SiO 2 (43 - 74 wt%) and Al 2 O 3 (15 - 19 wt%). Chemical Indices of Alteration (79.37-99.34) and Weathering (94.38-99.92) values suggest moderate to extreme silicate weathering, and alkali and alkaline earth metals depletion. Based on molar proportions of Al 2 O 3 , CaO + Na 2 O, and K 2 O, weathering trend and mineralogical compositions of the soils showed more advanced argillic alteration with rare earth elements being more concentrated compared to Upper Continental Crust. Excessive amounts of quartz and heavy metals in the geophagic soils represent significant health threats to geophagic individuals, though heavy metals had low biaccessibility values. There is need for soils beneficiation before ingestion.

A neotectonic structural interpretation was conducted in the Fujian Province, west of the Taiwan Strait, by using computer image processing and 3D visualizations to enhance linear structural traces. The major faults in this area can be grouped into two conjugate shear fracture zones, with one trending to the northeast and the other trending to the northwest. PS-InSAR technology uses stable permanent target scatterer points to determine deformation rates and can effectively reduce the influence of spatiotemporal decorrelations and atmospheric anomalies that affect conventional D-InSAR techniques and prevent the formation of interference fringes. This study focuses on the fault zones located in the Quanzhou area of Fujian Province, where the 1604 M7.5–8.0 historic earthquake occurred. In total, 22 scenes of ERS SAR data from 1996 to 1999 were processed using PS-InSAR methods. The results show that the line of sight direction displacement rate of the main fault in the study area is 3–5 mm/yr, which indicates that the faults in this area are still active and subject to earthquake risk.

Organic shale is one of the most important unconventional resources all around the world. Total organic carbon (TOC) content is an important evaluation parameter of reservoir hydrocarbon source quality. The regular evaluation methods have higher requirements of well logs and core experiment data for statistical regression. Through analyzing the resistivity and gamma ray logging response characteristics of shale content and organic matters, combined with digital rock physics experiment simulation, we put forward and improve the dual-shale-content method for TOC content logging evaluation. The accuracy of this method is verified by actual data processing. The result shows the dual-shale-content method is simple to use and the error is small. And by comparing with the calculation results by using the Δ LogR method, it is revealed that the trend of the TOC content calculated by our method agrees with the core results better. This new method is suitable for the evaluation of TOC content in the area or interval where there is only conventional logs.

Dolerite dykes are widespread in the Mongo area within the granitic Guéra Massif (Chad, Central Africa). Dykes are several hundred metres to several kilometres long, a metre to decametre thick, and vertical, crosscutting the Pan-African granitic basement rocks. They are controlled by major Pan-African NNE-SSW, NE-SW and ENE-WSW faults. Rocks constituting the dykes exhibit typical doleritic textures ( i.e. intergranular, ophitic or subophitic). They are mainly composed of phenocrysts, microcrysts and microlites of clinopyroxene, amphibole, plagioclase, alkali feldspar and Fe-Ti oxides. Mongo dolerites are classified into two groups: Group 1 of basalt and trachybasalt and Group 2 of basaltic trachyandesite, which follow a trend with continental tholeiite affinities, and are differentiated through fractional crystallization. Trace element enrichment and REE contents show the high melting degree at shallow depths of enriched sub-continental lithospheric mantle, whose composition is intermediate between OIB and MORB, contaminated by small amounts of sediment during older subduction. Mongo dolerites are interpreted as imprints of either the final stage of stabilization of ancient continental crust, or the initiation of tectonic activity related to Pan-African mobile belts.

Non-Local means algorithm is a new and effective filtering method. It calculates weights of all similar neighborhoods’ center points relative to filtering point within searching range by Gaussian weighted Euclidean distance between neighborhoods, then gets filtering point’s value by weighted average to complete the filtering operation. In this paper, geometric distance of neighborhood’s center point is taken into account in the distance measure calculation, making the non-local means algorithm more reasonable. Furthermore, in order to better protect the geometry structure information of seismic data, we introduce structure tensor that can depict the local geometrical features of seismic data. The coherence measure, which reflects image local contrast, is extracted from the structure tensor, is integrated into the non-local means algorithm to participate in the weight calculation, the control factor of geometry structure similarity is added to form a non-local means filtering algorithm based on structure tensor. The experimental results prove that the algorithm can effectively restrain noise, with strong anti-noise and amplitude preservation effect, improving PSNR and protecting structure information of seismic image. The method has been successfully applied in seismic data processing, indicating that it is a new and effective technique to conduct the structure-preserved filtering of seismic data.

In the Es 3 formation (third section of the Shahejie) of the Eastern sag section of the Liaohe Depression, basalt and trachyte are predominant in the igneous rock. The reservoir consists of complex reservoir space types. Based on the porosity bins of nuclear magnetic logging and the porosity distribution of electric imaging logging, the pores’ sizes and distribution, as well as the mutual connectivity of the reservoir, were analyzed. Also, the characteristics of the different reservoirs were summarized. In regards to the oil reservoirs, large pores (PS>10) were found to account for the majority of the reservoir spaces, and the pore distribution was concentrated and well connected. However, for the poor oil reservoirs, the large and small pores were found to alternate, and the pore distribution was scattered and poorly connected. Within the dry layers, the smaller pores (PS<10) were predominant. The pore distributions were found to be influenced by lithology, facies, and tectonism. The reservoirs of the pyroclastic flow of the explosive facies had good connectivity, and the interlayer heterogeneity was relatively weak. This reservoir’s pore distributions were found to be mainly dominated by the larger pores (PS10-PS13), which displayed a concentrated distribution mainly in one porosity bin. Therefore, it was taken as a favorable facie belt in the eastern sag of the Liaohe Depression. The examination of the pore distribution characteristics of the igneous rock was the key to the evaluation of the properties and effectiveness of the igneous reservoirs in this study, which potentially has great significance to the future exploration and development of igneous rock.

The Qaidam basin, bounded by the Altyn Tagh fault in the north, is located in the northeast of the Tibet plateau, and it has important implications for understanding the history and mechanism of Tibetan plateau formation during the Cenozoic Indo-Eurasia collision. In this study, we constructed the main geological structures and surfaces in three dimensions through the interpolation of regularly spaced 2D seismic sections, constrained by wells data and surface geology of the Qaidam basin in northeast Tibet. Meanwhile the Cenozoic tectonic history of the Qaidam basin was reconstructed and the uplift mechanism of the Tibetan plateau was discussed. This study presents the subsurface data in conjunction with observations and analysis of the stratigraphic and sedimentary evolution. The Cenozoic deformation history of the Qaidam basin shows geologic synchroneity with uplifting history of the Tibet Plateau. It is therefore proposed that the deformation and uplifting in the south and north edges of the Tibet Plateau was almost synchronous. The total shortening and shortening rate during Cenozoic reached 25.5 km and 11.2% respectively across the Qaidam basin, indicating that the loss of the left-lateral strike slip rates of the Altyn Tagh fault has been structurally transformed into local crustal thickening across NW-trending folds and thrust faults. Meanwhile there is an about 11° vertical component along the strike-slip Altyn Tagh fault, the block oblique slip shows one more growth mechanism of the northeast Tibet.

Water-filled fractures continue to grow after the excavation of karst tunnels, and the hydraulic pressure in these fractures changes along with such growth. This paper simplifies the fractures in the surrounding rock as flat ellipses and then identifies the critical hydraulic pressure values required for the occurrence of tensile-shear and compression-shear failures in water-filled fractures in the case of plane stress. The occurrence of tensile-shear fracture requires a larger critical hydraulic pressure than compression-shear failure in the same fracture. This paper examines the effects of fracture strike and lateral pressure coefficient on critical hydraulic pressure, and identifies compression-shear failure as the main failure mode of water-filled fractures. This paper also analyses the hydraulic pressure distribution in fractures with different extensions, and reveals that hydraulic pressure decreases along with the continuous growth of fractures and cannot completely fill a newly formed fracture with water. Fracture growth may be interrupted under the effect of hydraulic tensile shear.

This study attempts to explore the establishment of a low-carbon tourism community by public cognition, intention, and behaviour change analysis in a case study of a heritage site, China. Low carbon tourism advocates a way of travel with low energy consumption, low pollution and low CO 2 emissions during personal activities. Behaviour change is not only influenced by internal individual aspects including a person’s awareness, attitudes, and capacity to change, but is also driven by external social aspects including the culture and environments in which a person lives. In this paper, questionnaire surveys and field interviews were used to obtain basic information, and with reference to TPB, studied and analyzed the characteristics of cognition, intention and behaviour change practice by low carbon tourism community economy participants. With the help of SPSS analysis software, we found that a person’s educational level or occupation might affect cognition of low carbon tourism, and motives for participating in low carbon tourism could reflect the public’s perception of its emotional value, cognitive value and functional value. Most respondents knew about low carbon tourism; however, when putting it into practice, habitual behaviour was the main barrier for tourists while the residents were passive and followed the needs and choices of the tourists. Therefore, a comprehensive low carbon tourism community system was proposed not only for addressing the aspects of awareness, intention, and practice from individual behavior, but also for covering policy, infrastructure, institution systems and mechanisms at the community level.

Mangrove forests provide a wide range of socioeconomic and ecological services for coastal communities. Extensive aquaculture development of mangrove waters in many developing countries has constantly ignored services of mangrove ecosystems, leading to unintended environmental consequences. Monitoring the current status and distribution of mangrove forests is deemed important for evaluating forest management strategies. This study aims to delineate the density distribution of mangrove forests in the Gulf of Fonseca, Central America with Rapideye data using the support vector machines (SVM). The data collected in 2012 for density classification of mangrove forests were processed based on four different band combination schemes: scheme-1 (bands 1–3, 5 excluding the red-edge band 4), scheme-2 (bands 1–5), scheme-3 (bands 1–3, 5 incorporating with the normalized difference vegetation index, NDVI), and scheme-4 (bands 1–3, 5 incorporating with the normalized difference red-edge index, NDRI). We also hypothesized if the obvious contribution of Rapideye red-edge band could improve the classification results. Three main steps of data processing were employed: (1) data pre-processing, (2) image classification, and (3) accuracy assessment to evaluate the contribution of red-edge band in terms of the accuracy of classification results across these four schemes. The classification maps compared with the ground reference data indicated the slightly higher accuracy level observed for schemes 2 and 4. The overall accuracies and Kappa coefficients were 97% and 0.95 for scheme-2 and 96.9% and 0.95 for scheme-4, respectively.

Paleoclimatic data are essential for fingerprinting the climate of the earth before the advent of modern recording instruments. They enable us to recognize past climatic events and predict future trends. Within this framework, a conceptual and logical model was drawn to physically implement a paleoclimatic database named WDB-Paleo that includes the paleoclimatic proxies data of marine sediment cores of the Mediterranean Basin. Twenty entities were defined to record four main categories of data: a) the features of oceanographic cruises and cores (metadata); b) the presence/absence of paleoclimatic proxies pulled from about 200 scientific papers; c) the quantitative analysis of planktonic and benthonic foraminifera, pollen, calcareous nannoplankton, magnetic susceptibility, stable isotopes, radionuclides values of about 14 cores recovered by Institute for Coastal Marine Environment (IAMC) of Italian National Research Council (CNR) in the framework of several past research projects; d) specific entities recording quantitative data on δ 18 O, AMS 14 C (Accelerator Mass Spectrometry) and tephra layers available in scientific papers. Published data concerning paleoclimatic proxies in the Mediterranean Basin are recorded only for 400 out of 6000 cores retrieved in the area and they show a very irregular geographical distribution. Moreover, the data availability decreases when a constrained time interval is investigated or more than one proxy is required. We present three applications of WDB-Paleo for the Younger Dryas (YD) paleoclimatic event at Mediterranean scale and point out the potentiality of this tool for integrated stratigraphy studies.

Low impact development (LID)/green infrastructure (GI) practices have been identified as the sustainable practices of managing the stormwater in urban areas. Due to the increasing population, most of the cities are more developing which results in the change of natural area into impervious areas (roads, buildings etc.). Moreover, urbanization and climate change are causing many water-related problems and making over cities unsafe and insecure. Under these circumstances, there is a need to introduce new stormwater management practices into developed cities to reduce the adverse impacts of urbanization. For this purpose, retrofitting low impact development practices demands more attention to reduce these water-related problems and trying to make our cities sustainable. In developed areas, there is a little space is available for the retrofitting of LID practices for the stormwater management. Therefore, the selection of an appropriate place to retrofitting LID practices needs more concern. This paper describes the successfully applied retrofitting LID practices around the globe. It also includes the process of applying retrofitting LID practices at the suitable place with the suitable combination. Optimal places for the retrofitting of different LID practices are also mentioned. This paper also highlights the barriers and potential solutions of retrofitting LID practices in urban areas.

Geostatistics is a discipline that deals with the statistical analysis of regionalized variables. In this case study, geostatistics is used to estimate geoid undulation in the rural area of Guayaquil town in Ecuador. The geostatistical approach was chosen because the estimation error of prediction map is getting. Open source statistical software R and mainly geoR, gstat and RGeostats libraries were used. Exploratory data analysis (EDA), trend and structural analysis were carried out. An automatic model fitting by Iterative Least Squares and other fitting procedures were employed to fit the variogram. Finally, Kriging using gravity anomaly of Bouguer as external drift and Universal Kriging were used to get a detailed map of geoid undulation. The estimation uncertainty was reached in the interval [-0.5; +0.5] m for errors and a maximum estimation standard deviation of 2 mm in relation with the method of interpolation applied. The error distribution of the geoid undulation map obtained in this study provides a better result than Earth gravitational models publicly available for the study area according the comparison with independent validation points. The main goal of this paper is to confirm the feasibility to use geoid undulations from Global Navigation Satellite Systems and leveling field measurements and geostatistical techniques methods in order to use them in high-accuracy engineering projects.

The Ordovician paleokarst reservoirs in the Tahe oilfield, with burial depths of over 5300 m, experienced multiple phases of geologic processes and exhibit strong heterogeneity. Core testing can be used to analyse the characteristics of typical points at the centimetre scale, and seismic datasets can reveal the macroscopic outlines of reservoirs at the >10-m scale. However, neither method can identify caves, cave fills and fractures at the meter scale. Guided by outcrop investigations and calibrations based on core sample observations, this paper describes the interpretation of high longitudinal resolution borehole images, the identification of the characteristics of caves, cave fills (sedimentary, breccia and chemical fills) and fractures in single wells, and the identification of structures and fill characteristics at the meter scale in the strongly heterogeneous paleokarst reservoirs. The paleogeomorphology, a major controlling factor in the distribution of paleokarst reservoirs, was also analysed. The results show that one well can penetrate multiple cave layers of various sizes and that the caves are filled with multiple types of fill. The paleogeomorphology can be divided into highlands, slopes and depressions, which controlled the structure and fill characteristics of the paleokarst reservoirs. The results of this study can provide fundamental meter-scale datasets for interpreting detailed geologic features of deeply buried paleocaves, can be used to connect core- and seismic-scale interpretations, and can provide support for the recognition and development of these strongly heterogeneous reservoirs.

The article describes how the Structure-from-Motion (SfM) method can be used to calculate the volume of anthropogenic microtopography. In the proposed workflow, data is obtained using mass-market devices such as a compact camera (Canon G9) and a smartphone (iPhone5). The volume is computed using free open source software (VisualSFMv0.5.23, CMPMVSv0.6.0., MeshLab) on a PCclass computer. The input data is acquired from video frames. To verify the method laboratory tests on the embankment of a known volume has been carried out. Models of the test embankment were built using two independent measurements made with those two devices. No significant differences were found between the models in a comparative analysis. The volumes of the models differed from the actual volume just by 0.7‰ and 2‰. After a successful laboratory verification, field measurements were carried out in the same way. While building the model from the data acquired with a smartphone, it was observed that a series of frames, approximately 14% of all the frames, was rejected. The missing frames caused the point cloud to be less dense in the place where they had been rejected. This affected the model’s volume differed from the volume acquired with a camera by 7%. In order to improve the homogeneity, the frame extraction frequency was increased in the place where frames have been previously missing. A uniform model was thereby obtained with point cloud density evenly distributed. There was a 1.5% difference between the embankment’s volume and the volume calculated from the camera-recorded video. The presented method permits the number of input frames to be increased and the model’s accuracy to be enhanced without making an additional measurement, which may not be possible in the case of temporary features.

Coal and gas outbursts are the result of several geological factors related to coal seam gas (coal seam gas pressure P , coal seam sturdiness coefficient f and coal seam gas content W ), and these parameters can be used to classify the outburst hazard level of a coal seam. To classify the outburst hazard level of a coal seam by means of statistical methods, this study considered the geological parameters of coal seam gas and statistical data on the amount of material involved in coal outbursts. Through multivariate regression analysis, a multivariate regression equation between the outburst coal quantity and P , f and W was established. Using a significance evaluation of the aforementioned factors, the relative contributions of the gas-related geological parameters to the outburst hazard level of a coal seam were found to follow the order P > f > W . This work provides a scientific basis for evaluating the outburst hazard level of a coal seam and adopting feasible and economical outburst-prevention measures.

The Neogene strata of the onshore West Baram Province of NW Borneo contain organic rich rock formations particularly within the Sarawak basin. This basin is a proven prolific oil and gas province, thus has been a subject of great interest to characterise the nature of the organic source input and depositional environment conditions as well as thermal maturation. This study is performed on outcrop samples of Lambir, Miri and Tukau formations, which are of stratigraphic equivalence to the petroleum bearing cycles of the offshore West Baram delta province in Sarawak. The investigated mudstone samples are organic rich with a total organic carbon (TOC) content of more than 1.0 wt.%. The integration of elemental and molecular analyses indicates that there is no significant variation in the source input between these formations. The investigated biomarkers parameters achieved from acyclic isoprenoids, terpanes and steranes biomarkers of a saturated hydrocarbon biomarkers revealed that these sediments contain high contribution of land plants with minor marine organic matter input that was deposited and preserved under relatively oxic to suboxic conditions. This is further supported by low total sulphur (TS), high TOC/TN ratios, source and redox sensitive trace elements (V, Ni, Cr, Co and Mo) concentrations and their ratios, which suggest terrigenous source input deposited under oxic to suboxic conditions. Based on the analysed biomarker thermal maturity indicators, it may be deduced that the studied sediments are yet to enter the maturity stage for hydrocarbon generation, which is also supported by measured vitrinite reflectance values of 0.39-0.48% R o .

With the rapid economic development, the serious air pollution in Beijing attracts increasing attention in the last decade. Seen as one whole complex and grey system, the causal relationship between the social development and the air pollution in Beijing has been quantitatively analyzed in this paper. By using the grey relational model, the aim of this study is to explore how the socio-economic and human activities affect on the air pollution in the city of Beijing, China. Four air pollutants, as the particulate matter with size 2.5 micrometers or less (PM 2.5 ), particulate matter with size 10 micrometers or less (PM 10 ), sulfur dioxide (SO 2 ) and nitrogen dioxide (NO x ), are selected as the indicators of air pollution. Additionally, fifteen socio-economic indicators are selected to account for the regional socio-economic characteristics (economy variables, energy consumption variables, pollution emissions variables, environment and construction activity variables). The results highlight that all variables are associated with the concentrations of the four selected air pollutants, but with notable differences between the air pollutants. Most of the socio-economic indicators, such as industrial output, total energy consumption are highly correlated with PM 2.5 , while PM 10 , SO 2 , and NO x present in general moderate correlations with most of the socio-economic variables. Contrary to other studies and reports this study reveals that vehicles and life energy do not have the strongest effect on air pollution in Beijing. This study provides useful information to reduce air pollution and support decision-making for sustainable development.

Surfaces weathering of rocks in which mineral materials may be similar to or quite different from the minerals in the underlying parent rock completely control the reflectance spectra of the terrain. Our study of typical weathered and fresh rock samples from the Xiemisitai metallogenic belt, Western Junggar region, Xinjiang, found that weathering results in the formation of new materials that cause differences in the spectral features of fresh and weathered rock surfaces. Alterations induce variations in spectrum brightness, presence and intensity of characteristic absorption features, and spectral slope. Spectral differences between weathered and fresh rock surfaces are small for rhyolite, granite, and tuffaceous sandstone, but large for andesite, basalt, and diorite. Spectral changes in the 350–1000 nm wavelength region are attributed to alteration of iron oxides by atmospheric processes or secondary alteration of iron-rich minerals. Spectral features between 1000–2500 nm are caused by O–H vibrations, with features at 2200–2500 nm solely attributed to hydroxyl groups. The strongest Al–OH bands appear near 2200 nm, while Mg–OH bands are found near 2300 nm and 2350 nm. Results from this study can be used to better characterize and discriminate lithological units and potential mineral zones using hyperspectral and multispectral remote sensing techniques.

Geochemical compositions of twenty-four sandstone and shale samples from the Ecca Group were analysed to decipher their provenance, paleoweathering conditions and tectonic setting. The shales have high Fe 2 O 3 , K 2 O, TiO 2 , Ce, Cu, Ga, La, Nb, Nd, Rb, Sc, Sr, Th and Y content more than the sandstones, whereas, sandstones are higher in SiO 2 , Hf and Zr than the shales. The positive correlations of Al 2 O 3 with other elements as well as the abundance of Ba, Ce, Th, Rb, Zn and Zr suggest that these elements are primarily controlled by the dominant clay minerals. Tectonic discrimination diagrams revealed that the sandstones and shales are mostly of quartzose sedimentary provenance, suggesting that they were derived from a cratonic interior or recycled orogen. The binary plots of TiO 2 versus Ni, TiO 2 against Zr and La/Th versus Hf as well as the ternary diagrams of V-Ni-Th*10 indicate that the shales and sandstones were derived from felsic igneous rocks. A-CN-K (Al 2 O 3 -CaO-K 2 O) ternary diagram and indices of weathering (CIA, CIW and PIS) suggest that the granitic source rocks underwent moderate to high degree of chemical weathering. The CIA values range between 24.41% and 83.76%, indicating low to high weathering conditions. The CIW values for the studied sandstones and shales range from 25.90 to 96.25%, suggesting moderate to high intensive chemical weathering. ICV values for the sandstones and shales vary from 0.71 to 3.6 (averaging 1.20) and 0.41 to 1.05 (averaging 0.82), respectively. The k 2 O/Na 2 O ratios for the studied samples vary from 0.71 to 8.29, which reveal moderate to high maturity. The plot of CIA against ICV shows that most of the shales are geochemically mature and were derived from both weak and intensively weathered source rocks. The tectonic setting discrimination diagrams support passive-active continental margin setting of the provenance.

The Granjeno Schist is a meta-volcanosedimentary upper Paleozoic complex in northeastern Mexico. We suggest different tectonic settings for metamorphism of its serpentinite and talc-bearing rocks based on petrographic and geochemical compositions. According to the REE ratios (La N /Yb N = 0.51 –20.0 and La N /Sm N = 0.72–9.1) and the enrichment in the highly incompatible elements Cs (0.1 ppm), U (2.8 ppm), and Zr (60 ppm) as well as depletion in Ba (1 – 15 ppm), Sr (1 –184 ppm), Pb (0.1 –14 ppm), and Ce (0.1 –1.9 ppm) the rocks have mid-ocean ridge and subduction zones characteristics. The serpentinite contains Al-chromite, ferrian chromite and magnetite. The Al-chromite is characterized by Cr# of 0.48 to 0.55 suggesting a MORB origin, and Cr# of 0.93 to 1.00 for the ferrian chromite indicates a prograde metamorphism. We propose at least two serpentinization stages of lithospheric mantle for the ultramafic rock of the Granjeno Schist, (1) a first in an ocean-floor environment at sub-greenschist to greenschist facies conditions and (2) later a serpentinization phase related to the progressive replacement of spinel by ferrian chromite and magnetite at greenschist to low amphibolite facies conditions during regional metamorphism. The second serpentinization phase took place in an active continental margin during the Pennsylvanian. We propose that the origin of the ultramafic rocks is related to an obduction and accretional event at the western margin of Pangea.

A part of the Upper Oligocene sand-rich turbidite systems of the Central Carpathian Basin is represented by the Zuberec Formation. Sand/mud-mixed deposits of this formation are well exposed in the northern part of the basin, allowing us to interpret the turbidite succession as terminal lobe deposits of a submarine fan. This interpretation is based on the discrimination of three facies associations that are comparable to different components of distributive lobe deposits in deep-water fan systems. They correspond to the lobe off-axis, lobe fringe and lobe distal fringe depositional subenvironments, respectively. The inferences about the depositional paleoenvironment based on sedimentological observations are verified by statistical analyses. The bed-thickness frequency distributions and vertical organization of the facies associations show cyclic trends at different hierarchical levels that enable us to reconstruct architectural elements of a turbidite fan. First, small-scale trends correspond with shift in the lobe element centroid between successive elements. Differences in the distribution and frequency of sandstone bed thicknesses as well as differences in the shape of bed-thickness frequency distributions between individual facies associations reflect a gradual fining and thinning in a down-dip direction. Second, meso-scale trends are identified within lobes and they generally correspond to the significant periodicity identified by the time series analysis of the bed thicknesses. The meso-scale trends demonstrate shifts in the position of the lobe centroid within the lobe system. Both types of trends have a character of a compensational stacking pattern and could be linked to autogenic processes. Third, a largescale trend documented by generally thickening-upward stacking pattern of beds, accompanied by a general increase of the sandstones/mudstones ratio and by a gradual change of percentage of individual facies, could be comparable to lobe-system scale. This trend probably indicates a gradual basinward progradation of lobe system controlled by allogenic processes related to tectonic activity of sources and sea-level fluctuations.

The structural order of Cretaceous-Tertiary kaolins of the Douala Sub-Basin in Cameroon was determined in this study. This was achieved using Fourier Transform Infrared spectroscopy (FTIR) with attenuated total reflectance (ATR) on the Bomkoul, Dibamba, Ediki, Logbaba, Missole and Yatchika kaolins. Transmittance spectra of 20 samples were recorded in the mid-infrared regions (MIR). Results show that most of the kaolins had the four distinguishable bands in the hydroxyl (OH) stretching region, though the second transmittance band (3670 cm −1 ) had low intensities. The hydroxyl stretching of water bands (3457 and 3442 cm −1 for KGa-1b and KGa-2, respectively) were mainly observed in all Logbaba samples at 3443 cm −1 and 3445 cm −1 . The two bands reflecting OH deformation of inner-surface hydroxyl and inner hydroxyl groups (937 and 915 cm −1 ) were quite visible in the Bomkoul, Dibamba, Logbaba and Missole II samples; slightly visible in all Yatchika and two of the Ediki; and not visible in Missole I and one of the Ediki samples. Therefore, based on the empirical classification of the degree of order of kaolinite, 10 of the studied samples had ordered structures. Three displayed partially ordered structures; four had poorly ordered structures, and three exhibited disordered structures. This study revealed that kaolins from Dibamba, Logbaba and Missole II in the Douala Sub-Basin had the best crystallinity, which is one of the important characteristics needed for industrial applications.

Satellite imaging systems have known limitations regarding their spatial and temporal resolution. The approaches based on subpixel mapping of the Earth’s environment, which rely on combining the data retrieved from sensors of higher temporal and lower spatial resolution with the data characterized by lower temporal but higher spatial resolution, are of considerable interest. The paper presents the downscaling process of the land surface temperature (LST) derived from low resolution imagery acquired by the Advanced Very High Resolution Radiometer (AVHRR), using the inverse technique. The effective emissivity derived from another data source is used as a quantity describing thermal properties of the terrain in higher resolution, and allows the downsampling of low spatial resolution LST images. The authors propose an optimized downscaling method formulated as the inverse problem and show that the proposed approach yields better results than the use of other downsampling methods. The proposed method aims to find estimation of high spatial resolution LST data by minimizing the global error of the downscaling. In particular, for the investigated region of the Gulf of Gdansk, the RMSE between the AVHRR image downscaled by the proposed method and the Landsat 8 LST reference image was 2.255°C with correlation coefficient R equal to 0.828 and Bias = 0.557°C. For comparison, using the PBIM method, it was obtained RMSE = 2.832°C, R = 0.775 and Bias = 0.997°C for the same satellite scene. It also has been shown that the obtained results are also good in local scale and can be used for areas much smaller than the entire satellite imagery scene, depicting diverse biophysical conditions. Specifically, for the analyzed set of small sub-datasets of the whole scene, the obtained RSME between the downscaled and reference image was smaller, by approx. 0.53°C on average, in the case of applying the proposed method than in the case of using the PBIM method.

Fractional differentiation has adequate use for investigating real world scenarios related to geological formations associated with elasticity, heterogeneity, viscoelasticity, and the memory effect. Since groundwater systems exist in these geological formations, modelling groundwater recharge as a real world scenario is a challenging task to do because existing recharge estimation methods are governed by linear equations which make use of constant field parameters. This is inadequate because in reality these parameters are a function of both space and time. This study therefore concentrates on modifying the recharge equation governing the EARTH model, by application of the Eton approach. Accordingly, this paper presents a modified equation which is non-linear, and accounts for parameters in a way that it is a function of both space and time. To be more specific, herein, recharge and drainage resistance which are parameters within the equation, became a function of both space and time. Additionally, the study entailed solving the non-linear equation using an iterative method as well as numerical solutions by means of the Crank-Nicolson scheme. The numerical solutions were used alongside the Riemann-Liouville, Caputo-Fabrizio, and Atangana-Baleanu derivatives, so that account was taken for elasticity, heterogeneity, viscoelasticity, and the memory effect. In essence, this paper presents a more adequate model for recharge estimation.

The research aimed at testing the combination of GAM/M-GAM models (Geosite Assessment Model/Modified Geosite Assessment Model) on selected geoheritage sites (geosites) of great scientific significance and geotourism potential. Testing was done on eight sites in the City of Belgrade (Serbian capital), an area which has significant potential for geotourism development in typical urban conditions. For this purpose, an assessment scale was used to highlight differences and similarities between main and additional values of the observed geosites. The modification of the original GAM model is based on the inclusion of visitors’ opinions regarding the importance of indicators in the assessment process. The assessment was done by using both GAM/M-GAM and its results were analyzed and compared afterwards. The analysis has successfully identified locations and features of geosites that require action for maintaining or increasing their overall value and function. Moreover, the principal aim of the paper was to analyze the relevance of each sub-indicator for the assessment process by introducing the importance factor in the modified model. The authors were able to point out those values of principal importance for geosite visitors, as well as to attach a different relevance to sub-indicators, which can influence the position of the geosites in the GAM/M-GAM matrices.

This work is based on the assumption that a resistivity meter can effectively monitor water saturation in earth levees and can be used as a warning system when saturation exceeds the expected seasonal maxima. We performed time-lapse ERT measurements to assess the capability of this method to detect areas where seepage is critical. These measurements were also very useful to design a prototype monitoring system with remarkable savings by customizing the specifications according to field observations. The prototype consists of a remotely controlled low-power resistivity meter with a spread of 48 stainless steel 20 × 20 cm plate electrodes buried at half-meter depth. We deployed the newly-designed permanent monitoring system on a critical levee segment. A weather station and an ultrasonic water level sensor were also installed in order to analyse the correlation of resistivity with temperature, rainfalls and water level seasonal variations. The preliminary analysis of the monitoring data shows that the resistivity maps follow a very reasonable trend related with the saturation/drying cycle of the levee caused by the seasonal variations of the water level in the irrigation channel. Sharp water level changes cause delayed and smooth resistivity variations. Rainfalls and, to a lesser extent, temperature seem to have an influence on the collected data but effects are apparently negligible beyond 1 m depth. The system is currently operating and results are continuously monitored.

Utilization of satellite data by means of remote sensing practices creates a wider window of opportunities to conduct robust outcomes Landsat OLI-8 data acquired over the peninsula of Halkidiki in Greece for the summer period of 2013, were analyzed to determine their utility to classify natural resources categories based on the determination of spectral bands combination. Consequently, spectral bands combination can be used to classify various categories based on their higher overall accuracy assessment. Spectral information contained in each of the utilized channels of Landsat OLI-8, statistics defining 12 cover type classes of interest were calculated and used as a basis for classification of the designated study area. Spectral bands combinations (4-5-6/ 3-4-5/ 3-4-5-6/ 4-5-6-7) fulfill the required accuracy for forest purpose as they conducted overall accuracies ranging between 7.85% and 96.29%. The spectral band combination 3-4-5-6-7 conducted the highest overall accuracy (98.15%). The poorest results were obtained by the analysis of spectral bands combination 4-6 (68.52%). Single spectral band 6 conducted the best overall results; however, spectral bands 5 and 6 seemed to be the most useful spectral bands combination for each category.

North Anatolian Fault Zone, on which large earthquakes have occurred in the past, migrates regularly from east to west, and it is one of the most active faults in the world. The purpose of this study is to estimate the coda wave quality factor ( Qc ) for each of the five sub regionsthat were determined according to the fault rupture of these large earthquakes and along the fault. 978 records have been analyzed for 1.5, 3, 6, 9, 12 and 18 Hz frequencies by Single Backscattering Method. Along the fault, the variations in the Qc with lapse time are determined via, Qc = (136±25)f (0.96±0.027) , Qc = (208±22)f (0.85±0.02) Qc = (307±28)f (0.72±0.025) at 20, 30, 40 sec lapse times, respectively. The estimated average frequency-dependence quality factor for all lapse time are; Qc (f) = (189±26)f (0.86±0.02) for Karliova-Tokat region; Qc (f) = (216±19)f (0.76±0.018) for Tokat-Çorum region; Qc (f) = (232±18)f (0.76±0.019) for Çorum-Adapazari region; Qc (f) = (280±28)f (0.79±0.021) for Adapazari-Yalova region; Qc (f) = (252±26)f (0.81±0.022) for Yalova-Gulf of Saros region. The coda wave quality factor at all the lapse times and frequencies is Qc(f) = (206±15)f (0.85±0.012) in the study area. The most change of Qc with lapse time is determined at Yalova-Saros region. The result may be related to heterogeneity degree of rapidly decreases towards the deep crust like compared to the other sub region. Moreover, the highest Qc is calculated between Adapazari – Yalova. It was interpreted as a result of seismic energy released by 1999 Kocaeli Earthquake. Besides, it couldn’t be established a causal relationship between the regional variation of Qc with frequency and lapse time associated to migration of the big earthquakes. These results have been interpreted as the attenuation mechanism is affected by both regional heterogeneity and consist of a single or multi strands of the fault structure.

Petrography of the sandstones of Ecca Group, Karoo Supergroup in the Eastern Cape Province of South Africa have been investigated on composition, provenance and influence of weathering conditions. Petrographic studies based on quantitative analysis of the detrital minerals revealed that the sandstones are composed mostly of quartz, feldspar and lithic fragments of metamorphic and sedimentary rocks. The sandstones have an average framework composition of 24.3% quartz, 19.3% feldspar, 26.1% rock fragments, and 81.33% of the quartz grains are monocrystalline. These sandstones are generally very fine to fine grained, moderate to well sorted, and subangular to subrounded in shape. In addition, they are compositionally immature and can be classified as feldspathic wacke and lithic wacke. The provenance characteristics suggest the influence of plutonic and metamorphic terrains (meta-magmatic arc) as the main source rock with minor debris derived from recycled sedimentary rocks. The latter revealed that the compositional immaturity of the sandstones is a result of weathering or recycling and short transport distance. The weathering diagrams and semi-quantitative weathering index indicate that the Ecca sandstones are mostly from a plutonic source area, with climatic conditions ranging from arid to humid. The detrital modal compositions of these sandstones are related to back arc to island and continental margin arc. These results, therefore, support previous studies that infer foreland basin setting for the Karoo Basin.

The aim of this study was to introduce a new morphometric index named Relief Index (RI). RI is the ratio of the total length of the contour lines and the surface area at which they occur. This easily calculated index provides an objective quantitative measure of relief variability as an important feature in geomorphological studies. To achieve this goal, a highly detailed morphometric analysis was carried out using a high-resolution (1m×1m) DEM. Twenty one sample areas in southern Poland were examined. These analyses showed RI, as a good tool for rapidly evaluating topography heterogeneity in division into relief classes. I distinguished 4 classes of the Relief Index that classify earth surface considering the variability of the relief. Results of the calculations demonstrated that there is a significant correlation between RI and the local relief and slopes, but there is no correlation between RI and planar curvatures and TWI. The relief of the sample areas were analysed using geomorphometric parameters (slopes, local relief, planar curvatures). Moreover the influence of the DEM resolution on Relief Index values was examined.

The main aim of this paper is to determine the depositional environments of an Upper-Eocene coarse-grained clastic succession in the Buda Hills, Hungary. First of all, we measured some commonly used parameters of samples (size, amount, roundness and sphericity) in a much more objective overall and faster way than with traditional measurement approaches, using the newly developed Rock Analyst application. For the multivariate data obtained, we applied Combined Cluster and Discriminant Analysis (CCDA) in order to determine homogeneous groups of the sampling locations based on the quantitative composition of the conglomerate as well as the shape parameters (roundness and sphericity). The result is the spatial pattern of these groups, which assists with the interpretation of the depositional processes. According to our concept, those sampling sites which belong to the same homogeneous groups were likely formed under similar geological circumstances and by similar geological processes. In the Buda Hills, we were able to distinguish various sedimentological environments within the area based on the results: fan, intermittent stream or marine.

Large-scale conglomerate fan-delta aprons were typical deposits on the slope of Mahu Depression during the Early Triassic. Without outcrops, it is difficult to study the lithofacies only by examining the limited cores from the main oil-bearing interval of the Baikouquan Formation. Borehole electrical imaging log provides abundant high-resolution geologic information that is obtainable only from real rocks previously. Referring to the lithology and sedimentary structure of cores, a case study of fan-deltas in the Lower Triassic Baikouquan Formation of the Mahu Depression presents a methodology for interpreting the complicated lithofacies utilizing borehole electrical images. Eleven types of lithologies and five types of sedimentary structures are summarized in borehole electrical images. The sediments are fining upward from gravel to silt and clay in the Baikouquan Formation. Fine-pebbles and granules are the main deposits in T 1 b 1 and T 1 b 2 , but sandstones, siltstones and mudstones are more developed in T 1 b 3 . The main sedimentary textures are massive beddings, cross beddings and scour-and-fill structures. Parallel and horizontal beddings are more developed in T 1 b 3 relatively. On integrated analysis of the lithology and sedimentary structure, eight lithofacies from electrical images, referred to as image lithofacies, is established for the fan-deltas. Granules to coarse-pebbles within massive beddings, granules to coarse-pebbles within cross and parallel beddings, siltstones within horizontal and massive beddings are the most developed lithofacies respectively in T 1 b 1 , T 1 b 2 and T 1 b 3 . It indicates a gradual rise of the lake level of Mahu depression during the Early Triassic, with the fan-delta aprons retrograding towards to the margin of the basin. Therefore, the borehole electrical imaging log compensate for the limitation of cores of the Baikouquan Formation, providing an effective new approach to interpret the lithofacies of fan-delta.

Grain size analysis is a vital sedimentological tool used to unravel the hydrodynamic conditions, mode of transportation and deposition of detrital sediments. In this study, detailed grain-size analysis was carried out on thirty-five sandstone samples from the Ecca Group in the Eastern Cape Province of South Africa. Grain-size statistical parameters, bivariate analysis, linear discriminate functions, Passega diagrams and log-probability curves were used to reveal the depositional processes, sedimentation mechanisms, hydrodynamic energy conditions and to discriminate different depositional environments. The grain-size parameters show that most of the sandstones are very fine to fine grained, moderately well sorted, mostly near-symmetrical and mesokurtic in nature. The abundance of very fine to fine grained sandstones indicate the dominance of low energy environment. The bivariate plots show that the samples are mostly grouped, except for the Prince Albert samples that show scattered trend, which is due to the either mixture of two modes in equal proportion in bimodal sediments or good sorting in unimodal sediments. The linear discriminant function analysis is dominantly indicative of turbidity current deposits under shallow marine environments for samples from the Prince Albert, Collingham and Ripon Formations, while those samples from the Fort Brown Formation are lacustrine or deltaic deposits. The C-M plots indicated that the sediments were deposited mainly by suspension and saltation, and graded suspension. Visher diagrams show that saltation is the major process of transportation, followed by suspension.

We presented new C and O isotope data of rockforming calcite in terrestrial carbonate deposits from Neogene basins of Attica (Greece), coupled with standard mineralogical and bulk geochemical results. Whereas both isotope datasets [ δ 18 O from −8.99 to −3.20‰(VPDB); δ 13 C from −8.17 to +1.40‰(VPDB)] could be interpreted in principle as indicative of a meteoric origin, the clear lack of a statistical correlation between them suggests diverse sources for the isotopic variation of the two elements. On the basis of broad correlations between lower carbon isotope data with increasing Fe and bulk organic carbon, we interpreted the light carbon isotope signatures and As enrichments as both derived mainly from a depositional process involving increased supply of metals and organic carbon to the original basins. Periodically augmented biological production and aerobic cycling of organic matter in the ambient lake waters, would have led to the precipitation of isotopically light calcite in concert with elevated fluxes of As-bearing iron oxy-hydroxide and organic matter to the initial terrestrial carbonate sediment. The terrestrial carbonate deposits of Attica therefore represented effective secondary storage reservoirs of elevated As from the adjacent mineralized hinterland; hence these and similar deposits in the region ought to be regarded as key geological candidates for anomalous supply of As to local soils, groundwater and related human activities.

Geolocated photos from google Panoramio are used as a proxy for evaluation of aesthetic values appreciation of different landscape types in Slovakia. We collected the photo’s metadata from years 2005 - 2014 and calculated the density of photos uploaded by unique user per square kilometre. Then we compared the photos density in different landscape types. The most appreciated are subalpine and alpine landscape types. The high photo density was also found in urban landscapes where most of the population live. Outside the urban area, we found that less intensive type of landscapes are visually more attractive. From the abiotic landscape categories the most aesthetically valuable are landscapes in giant highlands and glacial giant highlands. The lowland landscape used intensively for agricultural production is less attractive.

The geochemical characteristics of selected kaolins from Cameroon and Nigeria are presented, with an attempt to elucidate on their possible industrial applications by comparing them to world-known kaolin deposits. Major oxides concentrations were subjected to factor analyses in interpreting their relationships. Geochemical indices, including chemical index of alteration (CIA), chemical index of weathering (CIW) and the index of compositional variability (ICV) were computed and plotted on binary and ternary diagrams to determine the intensity of weathering of the kaolins and discriminate their different source rock types. Kaolinite was the major phase, followed by quartz, illite and goethite as minor phases. Minerals in trace phases included smectite, anatase, muscovite, gibbsite, microcline, palygorskite and calcite. Mean abundances of major oxides in wt% were: SiO 2 (56.96)>Al 2 O 3 (24.09)>Fe 2 O 3 (3.78)>TiO 2 (1.53)> K 2 O (1.26)> MgO (0.27)>CaO (0.20)>Na 2 O (0.17)>P 2 O 5 (0.05)>MnO (0.04). The CIW versus CIA and ICV versus CIA plots showed that most of the kaolins clearly depicted extreme silicate weathering. The current applications of kaolins from Cameroon and Nigeria include ceramics and manufacturing of bricks and tiles. Low MgO, CaO, Na 2 O, K 2 O and TiO 2 further position the kaolins for pharmaceutics, cosmetics, rubber and plastic applications. Thus, the studied kaolins have the potential to contribute to improved economic development of these countries.

The article presents the spatial pattern analysis of the ASG-EUPOS permanent GNSS stations in Poland. Using different methods and tools (nearest neighbour, Riplay’s K-function, morphology of Thiessen polygons) we proved that the station distribution model changes within scales. At short distances up to 65 km, which are typical lengths in the network, stations are irregularly dispersed. Increasing this distance to 130 km and over could result in a clustered pattern. The Thiessen polygon area in 72% depends on the level of urbanization, especially coverage of forested and built-up areas as well as the density of the transportation network. The smallest density of the ASG-EUPOS sites is one station over 10,000 sq. km, which is two times more than is stated in the national regulations. The mean distance from ASG-EUPOS location to the nearest station is about 41.5 km.

Surface area estimation is a widely used tool for resource evaluation in the physical world. When processing large scale spatial data, the input/output (I/O) can easily become the bottleneck in parallelizing the algorithm due to the limited physical memory resources and the very slow disk transfer rate. In this paper, we proposed a stream tilling approach to surface area estimation that first decomposed a spatial data set into tiles with topological expansions. With these tiles, the one-to-one mapping relationship between the input and the computing process was broken. Then, we realized a streaming framework towards the scheduling of the I/O processes and computing units. Herein, each computing unit encapsulated a same copy of the estimation algorithm, and multiple asynchronous computing units could work individually in parallel. Finally, the performed experiment demonstrated that our stream tilling estimation can efficiently alleviate the heavy pressures from the I/O-bound work, and the measured speedup after being optimized have greatly outperformed the directly parallel versions in shared memory systems with multi-core processors.

In this paper, a location-based multiple point statistics method is developed to model a non-stationary reservoir. The proposed method characterizes the relationship between the sedimentary pattern and the deposit location using the relative central position distance function, which alleviates the requirement that the training image and the simulated grids have the same dimension. The weights in every direction of the distance function can be changed to characterize the reservoir heterogeneity in various directions. The local integral replacements of data events, structured random path, distance tolerance and multi-grid strategy are applied to reproduce the sedimentary patterns and obtain a more realistic result. This method is compared with the traditional Snesim method using a synthesized 3-D training image of Poyang Lake and a reservoir model of Shengli Oilfield in China. The results indicate that the new method can reproduce the non-stationary characteristics better than the traditional method and is more suitable for simulation of delta-front deposits. These results show that the new method is a powerful tool for modelling a reservoir with non-stationary characteristics.

The construction of tunnels and underground engineering in China has developed rapidly in recent years in both the number and the length of tunnels. However, with the development of tunnel construction technology, risk assessment of the tunnels has become increasingly important. Water inrush is one of the most important causes of engineering accidents worldwide, resulting in considerable economic and environmental losses. Accordingly, water inrush prediction is important for ensuring the safety of tunnel construction. Therefore, in this study, we constructed a three-dimensional discrete network fracture model using the Monte Carlo method first with the basic data from the engineering geological map of the Longmen Mountain area, the location of the Longmen Mountain tunnel. Subsequently, we transformed the discrete fracture networks into a pipe network model. Next, the DEM of the study area was analysed and a submerged analysis was conducted to determine the water storage area. Finally, we attempted to predict the water inrush along the Longmen Mountain tunnel based on the Darcy flow equation. Based on the contrast of water inrush between the proposed approach, groundwater dynamics and precipitation infiltration method, we conclude the following: the water inflow determined using the groundwater dynamics simulation results are basically consistent with those in the D2K91+020 to D2K110+150 mileage. Specifically, in the D2K91+020 to D2K94+060, D2K96+440 to D2K98+100 and other sections of the tunnel, the simulated and measured results are in close agreement and show that this method is effective. In general, we can predict the water inflow in the area of the Longmen Mountain tunnel based on the existing fracture joint parameters and the hydrogeological data of the Longmen Mountain area, providing a water inrush simulation and guiding the tunnel excavation and construction stages.

In this study, firstly, a practical and educational geostatistical program (JeoStat) was developed, and then example analysis of porosity parameter distribution, using oilfield data, was presented. With this program, two or three-dimensional variogram analysis can be performed by using normal, log-normal or indicator transformed data. In these analyses, JeoStat offers seven commonly used theoretical variogram models (Spherical, Gaussian, Exponential, Linear, Generalized Linear, Hole Effect and Paddington Mix) to the users. These theoretical models can be easily and quickly fitted to experimental models using a mouse. JeoStat uses ordinary kriging interpolation technique for computation of point or block estimate, and also uses cross-validation test techniques for validation of the fitted theoretical model. All the results obtained by the analysis as well as all the graphics such as histogram, variogram and kriging estimation maps can be saved to the hard drive, including digitised graphics and maps. As such, the numerical values of any point in the map can be monitored using a mouse and text boxes. This program is available to students, researchers, consultants and corporations of any size free of charge. The JeoStat software package and source codes available at: http://www.jeostat.com/JeoStat_2017.0.rar

The Khorat Plateau from northeast Thailand, the upstream part of the Mun River flows through clastic sedimentary rocks. A massive amount of sand was transported. We aimed to understand the evolution of fluvial system and to discuss the advantages of two shallow geophysical methods for describing subsurface morphology of modern and paleo-channels. We applied Electrical Resistivity Tomography (ERT) and Ground Penetrating Radar (GPR) to characterize the lateral, vertical morphological and sedimentary structures of paleo-channels, floodplain and recent point bars. Both methods were interpreted together with on-sites boreholes to describe the physical properties of subsurface sediments. As a result, we concluded that four radar reflection patterns including reflection free, shingled, inclined and hummocky reflections were appropriated to apply as criteria to characterize lateral accretion, the meandering rivers with channel-filled sequence and floodplain were detected from ERT profiles. The changes in resistivity correspond well with differences in particle size and show relationship with ERT lithological classes. Clay, silt, sand, loam and bedrock were classified by the resistivity data. Geometry of paleo-channel embayment and lithological differences can be detected by ERT, whereas GPR provides detail subsurface facies for describing point bar sand deposit better than ERT.

The Sidoarjo mudflow in East Java, Indonesia, has been erupting since May 29 th , 2006. The eruption has been known as the Lusi (lumpur Sidoarjo), which was previously considered as a remote seismic event consequence, but current geyser-like activities show an association with a geothermal phenomenon. A method of characterizing rare earth elements (REE) is commonly an effective tool for recognizing a geothermal system, and here it is adapted to particularly indicate the environmental origin of the Lusi mud. Results show that the Lusi hot mud is made of a porous smectite structure of a shale rock type, which becomes an ideal tank for trapping the REE, especially the light REE. Volcanic activities seem to be an important influence in the eruption; however, since there is a lack of significant isotopic evidences in the mobilization of the REE during the eruption, the chloride neutral pH water of the Lusi may hardly contain the REE. The moderate Ce and Eu anomalies found in the REE patterns of the mud strongly indicate a sea-floor basin as the most probable environment for the REE fractionation during the sedimentary rock formation, in which the weathering processes of volcanic rock origin enriched the Lusi shale with the REE.

Onitsha is one of the largest commercial cities in Africa with its population growth rate increasing arithmetically for the past two decades. This situation has direct and indirect effects on the natural resources including vegetation and water. The study aimed at assessing land use-land cover (LULC) change and its effects on the vegetation and landscape from 1987 to 2015 using geoinformatics. Supervised and unsupervised classifications including maximum likelihood algorithm were performed using ENVI 4.7 and ArcGIS 10.1 versions. The LULC was classified into 7 classes: built-up areas (settlement), waterbody, thick vegetation, light vegetation, riparian vegetation, sand deposit (bare soil) and floodplain. The result revealed that all the three vegetation types decreased in areas throughout the study period while, settlement, sand deposit and floodplain areas have remarkable increase of about 100% in 2015 when compared with the total in 1987. Number of dominant plant species decreased continuously during the study. The overall classification accuracies in 1987, 2002 and 2015 was 90.7%, 92.9% and 95.5% respectively. The overall kappa coefficient of the image classification for 1987, 2002 and 2015 was 0.98, 0.93 and 0.96 respectively. In general, the average classification was above 90%, a proof that the classification was reliable and acceptable.

We used 33 stations belonging of the Ecuador Continuous Monitoring GNSS Network (REGME) during the period 2008-2014, with aim to contribute with a methodological approach for the estimation of a new velocity model for Continental Ecuador. We used daily solutions to perform the analysis of GNSS time series, to obtain models of the series that best fit, taking into count the trend, seasonal variations and the type of noise. The sum of all these components represent the real-time series, and thus we can have a better estimation of the velocity parameter and its uncertainty. The velocities were calculated introducing the trend, seasonality and noise that were presented in each series into the overall model, which improved uncertainty by 12% and changed in magnitude up to 1.7 mm/yr and 2.5 mm/yr in the horizontal and vertical components, respectively, with respect to the initial velocities. The velocity field describes the crustal movement of the REGME stations in mainland Ecuador with uncertainty of 1 mm/yr and 2 mm/yr for the horizontal and vertical components, respectively. Finally, a velocity model has been developed using the kriging technique whose geostatistical approach has been based on the data to identify the spatial characteristics by examining the observations by peers. The mean square error (rms) of prediction obtained in this method is 1.78 mm/yr and 1.95 mm/yr in the east and north components, respectivaly. The vertical component could not be modeled due to its chaotic behavior.