The origin of the Pedra Furada sandstone tubular structures (South of Lisbon, Portugal)

The Pedra Furada is a 12 m wide, 18 m high outcrop feature showing hundreds of ferruginised sandy tubes and looking in part like a giant organ. In this paper the origin of the tubes is explained on the basis of geochemical, petrographic and microscopic (optical and electronic) analytical data. The tubes are considered to represent vertical escape channels for overpressured water, exhibiting inward decreasing grain size due to water velocity gradients inside the escape channels. The ferruginisation is due to iron oxides associated with colloidal/clayey fine sediments and to goethite formed from solution. The overpressure of water may be due to seismically fluidised beds below the Pedra Furada outcrop or to artesian water ascent. In both cases, fault rupturing may have played a major role in the focussing of the ascending flow.     

High-dimensional geostatistical history matching

Hydrocarbon reservoir modelling and characterisation is a challenging subject within the oil and gas industry due to the lack of well data and the natural heterogeneities of the Earth’s subsurface. Integrating historical production data into the geo-modelling workflow, commonly designated by history matching, allows better reservoir characterisation and the possibility of predicting the reservoir behaviour. We present herein a geostatistical-based multi-objective history matching methodology. It starts with the generation of an initial ensemble of the subsurface petrophysical property of interest through stochastic sequential simulation. Each model is then ranked according the match between its dynamic response, after fluid flow simulation, and the observed available historical production data. This enables building regionalised Pareto fronts and the definition of a large ensemble of optimal subsurface Earth models that fit all the observed production data without compromising the exploration of the uncertainty space. The proposed geostatistical multi-objective history matching technique is successfully implemented in a benchmark synthetic reservoir dataset, the PUNQ-S3, where 12 objectives are targeted.     

Gravimetric and magnetic fabric study of the Sintra Igneous complex: laccolith-plug emplacement in the Western Iberian passive margin

International Journal of Earth Sciences - The geometry and emplacement of the ~ 96 km2, Late Cretaceous Sintra Igneous complex (SIC, ca. 80 Ma) into the West Iberian passive...     

Permian volcanic rocks from the Apuseni Mountains (Romania): Geochemistry and tectonic constraints

An Early Permian volcanic assemblage is well exposed in the central-western part of the Apuseni Mountains (Romania). The rocks are represented by rhyolites, basalts and subordinate andesites suggesting a bimodal volcanic activity that is intimately associated with a post-orogenic (Variscan) syn-sedimentary intra-basinal continental molasse sequences. The mafic and mafic-intermediate rocks belong to sub-alkaline tholeiitic series were separated in three groups (I–III) showing a high Th and Pb abundances, depletion in Nb, Ta and Sr, and slightly enriched in LREE patterns (La_N/Yb_N = 1.4–4.4). Isotopically, the rocks of Group I have the initial ratios ⁸⁷Sr/⁸⁶Sr_(i) = 0.709351–0.707112, ¹⁴³Nd/¹⁴⁴Nd_(i) = 0.512490–0.512588 and high positive ?_Nd270 values from 3.9 to 5.80; the rocks of Group II present for the initial ratios values ⁸⁷Sr/⁸⁶Sr_(i) = 0.709434–0.710092, ¹⁴³Nd/¹⁴⁴Nd_(i) = 0.512231–0.512210 and for ?_Nd270 the negative values from −1.17 to −1.56; the rocks of Group III display for the initial ratios the values ⁸⁷Sr/⁸⁶Sr_(i) = 0.710751–0.709448, ¹⁴³Nd/¹⁴⁴Nd_(i) = 0.512347–0.512411 and for ?_Nd270 the positive values from 1.64 to 2.35. The rocks resembling continental tholeiites, suggest a mantle origin and were further affected by fractionation and crustal contamination. In addition, the REE geochemistry (1 > Sm_N/Yb_N < 2.5; 0.9 > La_N/Sm_N < 2.5) suggests that these rocks were generated by high percentage partial melting of a metasomatized mantle in the garnet peridotite facies. The felsic rocks are enriched in Cs, Rb Th and U and depleted in Nb, Ta, Sr, Eu, and Ti. The REE fractionation patterns show a strong negative Eu anomaly (Eu/Eu* = 0.23–0.40). The felsic rocks show the initial ratios the values: ⁸⁷Sr/⁸⁶Sr_(i) = 0.704096–0.707805, ¹⁴³Nd/¹⁴⁴Nd_(i) = 0.512012–0.512021 and for ?_Nd270 the negative values from −5.27 to −5.44. They suggest to be generated within the lower crust during the emplacement of mantle-derived magmas that provided necessary heat to crustal partial melting.     

Generative Adversarial Network Applied to Ore Type Modeling in Complex Geological Environments

Mathematical Geosciences - In complex geological environments, the spatial domain characterization of different ore types has been one of the most important challenges in the assessment of mineral...     

Urban forest loss using a GIS-based approach and instruments for integrated urban planning:A case study of Jo(a)o Pessoa,Brazil

Urban forests play an important role in the thermal comfort and overall life of local populations in large-and medium-sized cities.This study analyzes urban forest loss and maps land use and land cover (LULC) changes between 1991 and 2018 by evaluating the use of urban planning instruments for the mitigation of urban forest loss in Jo～o Pessoa,Brazil.For this purpose,satellite-derived LULC images from 1991,2006,2010 and 2018 and data on urban forest loss areas obtained using the Google Earth Engine were used.In addi-tion,this paper also discusses the instruments used for integrated urban planning,which are(a) the legal sector,responsibility and nature;(b) the urban expansion process;and (c) the elements of urban infrastructure.The results show a clear shift in land use in the study area.The major changes in LULC classes occurred in urban areas and herbaceous vegetation,while the greatest loss was in arboreal/shrub vegetation.Thus,an increase in the pressure to occupy zones intended for environmental preservation could be estimated.Our results showed similar accuracies with other studies and more spatial details.The characteristics of the patterns,traces,and hotspots of urban expansion and forest cover loss were explored.We highlighted the potential use of this proposed framework to be applied and validated in other parts of the world to help better understand and quantify various aspects of ur-ban-related problems such as urban forest loss mapping using instruments for integrated urban planning and low-cost approaches.     

Regularized sparse‐grid geometric sampling for uncertainty analysis in non‐linear inverse problems

This paper introduces an efficiency improvement to the sparse‐grid geometric sampling methodology for assessing uncertainty in non‐linear geophysical inverse problems. Traditional sparse‐grid geometric sampling works by sampling in a reduced‐dimension parameter space bounded by a feasible polytope, e.g., a generalization of a polygon to dimension above two. The feasible polytope is approximated by a hypercube. When the polytope is very irregular, the hypercube can be a poor approximation leading to computational inefficiency in sampling. We show how the polytope can be regularized using a rotation and scaling based on principal component analysis. This simple regularization helps to increase the efficiency of the sampling and by extension the computational complexity of the uncertainty solution. We demonstrate this on two synthetic 1D examples related to controlled‐source electromagnetic and amplitude versus offset inversion. The results show an improvement of about 50% in the performance of the proposed methodology when compared with the traditional one. However, as the amplitude versus offset example shows, the differences in the efficiency of the proposed methodology are very likely to be dependent on the shape and complexity of the original polytope. However, it is necessary to pursue further investigations on the regularization of the original polytope in order to fully understand when a simple regularization step based on rotation and scaling is enough.