Simulation of discrete cracks driven by nearly incompressible fluid via 2D combined finite-discrete element method

In this work, we propose a novel hydraulic solver in order to simulate key mechanisms that control fluid-driven cracks in the framework of the combined finite-discrete element method (FDEM). The main innovative aspect of the present work is the independence of the fluid's critical time step size with the fracture opening. This advantage is extremely important because it means that very fine meshes can be used around areas of interest, such as boreholes, without penalizing the computational cost as fractures propagate (ie, open) and the fluid flows through them. This is a great advantage over other recently introduced approaches that exhibit a dependency of the time step in the form of Δt_crit ∝ (l/a)² where l is the element size and a is the fracture aperture. This paper presents a series of benchmark cases for the proposed solver. The rationale adopted by the authors was to benchmark and validate the implementation of the hydraulic solver in an incremental fashion, starting from the simplest cases and building in complexity. The results shown in this work clearly demonstrate that the proposed approach is able to reproduce analytical results for fluid flow through a single crack. The results presented in this paper also demonstrate that the new approach is robust enough to deal with complex fracture patterns and complex geometries; the obtained fluid-driven fracture patterns in the vicinity of a borehole certainly stand to the scrutiny of human visual perception.     

Climatic variations in the coastal zone — Comparison between the Po River delta (Adriatic Sea,Italy) and the Huanghe River delta (Bohai Sea,China)

The aim of this work is to research the influence of natural climatic changes on the evolution of the coastal zone in modern times and the possible implication of human activities on the configuration of the present coastline. Comparison of data of two very far and different areas, the Po River delta, Adriatic Sea and the Huanghe River delta, Bohai Sea, reveals the planetary diffusion of climatic fluctuations and their effects on coastal evolution. Contribution No. 820 from the Marine Geology Institute, C. N. R., Italy. Contribution No.2207 from the Institute of Oceanology, Academia Sinica. Project supported by the National Natural Science Foundation of China.     

A formal framework for the representation of stack-based terrains

This paper presents a formal framework for the representation of three-dimensional geospatial data and the definition of common geographic information system (GIS) spatial operations. We use the compact stack-based representation of terrains (SBRT) in order to model geological volumetric data, both at the surface and subsurface levels, thus preventing the large storage requirements of regular voxel models. The main contribution of this paper is fitting the SBRT into the geo-atom theory in a seamless way, providing it with a sound formal geographic foundation. In addition we have defined a set of common spatial operations on this representation using the tools provided by map algebra. More complex geoprocessing operations or geophysical simulations using the SBRT as representation can be implemented as a composition of these fundamental operations. Finally a data model and an implementation extending the coverage concept provided by the Geography Markup Language standard are suggested. Geoscientists and GIS professionals can take advantage of this model to exchange and reuse geoinformation within a well-specified framework.     

A genetic classification of sinkholes illustrated from evaporite paleokarst exposures in Spain

This contribution analyses the processes involved in the generation of sinkholes from the study of paleokarst features exposed in four Spanish Tertiary basins. Bedrock strata are subhorizontal evaporites, and in three of the basins they include halite and glauberite in the subsurface. Our studies suggest that formation of dolines in these areas results from a wider range of subsidence processes than those included in the most recently published sinkhole classifications; a new genetic classification of sinkholes applicable to both carbonate and evaporite karst areas is thus proposed. With the exception of solution dolines, it defines the main sinkhole types by use of two terms that refer to the material affected by downward gravitational movements (cover, bedrock or caprock) and the main type of process involved (collapse, suffosion or sagging). Sinkholes that result from the combination of several subsidence processes and affect more than one type of material are described by combinations of the different terms with the dominant material or process followed by the secondary one (e.g. bedrock sagging and collapse sinkhole). The mechanism of collapse includes any brittle gravitational deformation of cover and bedrock material, such as upward stoping of cavities by roof failure, development of well-defined failure planes and rock brecciation. Suffosion is the downward migration of cover deposits through dissolutional conduits accompanied with ductile settling. Sagging is the ductile flexure of sediments caused by differential corrosional lowering of the rockhead or interstratal karstification of the soluble bedrock. The paleokarsts we analysed suggest that the sagging mechanism (not included in previous genetic classifications) plays an important role in the generation of sinkholes in evaporites. Moreover, collapse processes are more significant in extent and rate in areas underlain by evaporites than in carbonate karst, primarily due to the greater solubility of the evaporites and the lower mechanical strength and ductile rheology of gypsum and salt rocks.     

Geology,fluid inclusion and sulphur isotope characteristics of the El Cobre VHMS deposit,Southern Cuba

The El Cobre deposit is located in eastern Cuba within the volcanosedimentary sequence of the Sierra Maestra Paleogene arc. The deposit is hosted by tholeiitic basalts, andesites and tuffs and comprises thick stratiform barite and anhydrite bodies, three stratabound disseminated up to massive sulphide bodies produced by silicification and sulphidation of limestones or sulphates, an anhydrite stockwork and a siliceous stockwork, grading downwards to quartz veins. Sulphides are mainly pyrite, chalcopyrite and sphalerite; gold occurs in the stratabound ores. Fluid inclusions measured in sphalerite, quartz, anhydrite and calcite show salinities between 2.3 and 5.7 wt% NaCl eq. and homogenisation temperatures between 177 and 300°C. Sulphides from the stratabound mineralisation display δ ³⁴S values of 0‰ to +6.0‰, whilst those from the feeder zone lie between −1.4‰ and +7.3‰. Sulphides show an intra-grain sulphur isotope zonation of about 2‰; usually, δ ³⁴S values increase towards the rims. Sulphate sulphur has δ ³⁴S in the range of +17‰ to +21‰, except two samples with values of +5.9‰ and +7.7‰. Sulphur isotope data indicate that the thermochemical reduction of sulphate from a hydrothermal fluid of seawater origin was the main source of sulphide sulphur and that most of the sulphates precipitated by heating of seawater. The structure of the deposit, mineralogy, fluid inclusion and isotope data suggest that the deposit formed from seawater-derived fluids with probably minor supply of magmatic fluids.     

Climatic zoning for the calculation of the thermal demand of buildings in Extremadura (Spain)

The present work reports on a methodology to assess the climatic severity of a particular geographic region as compared to specific information available in the current regulations. The viability for each of the 387 municipalities in the Autonomous Community of Extremadura (Spain) is analysed, making a distinction between those with reliable climate reports and those for which no such information is available. In the case study, although the weather conditions in Extremadura are quite homogeneous according to the Spanish Technical Building Code (STBC 2015) classification and most areas are associated to zone C4 (soft winters and hot summers), the southern area in the region is associated to zone D1, similar to the north of Spain, where winters and summers are cool, which does not coincide with the actual climate in the south of Extremadura. The general climatic homogeneity in Extremadura was also highlighted with the new procedure, predominating zone C4, but unexpected or unreal climatic zoning was not generated, giving place to a consistent spatial distribution of zones throughout the region. Consequently, the proposed method allows a more accurate climatic zoning of any region in agreement with the Spanish legislation on energy efficiency in buildings, which would enhance the setting of thermal demand rates according to the actual climatic characterisation of the area in which a particular municipality is located.     

Linking sedimentation rates and large‐scale architecture for facies prediction in nonmarine basins (Paleogene,Almazán Basin,Spain)

This article focuses on the relationships between the large‐scale stratigraphic architecture of the Almazán basin infill and the sedimentation rates (SR) calculated for precise time intervals. Our aim was to improve the understanding of the timing and causes of the architectural changes, their significance in terms of accommodation space and sediment supply and their relationship with climate and tectonics. The study area includes the Gómara fluvial fan, the main sediment transfer system of the Almazán basin during Paleogene times. Its large‐scale architecture shifted through time between a stacking pattern of low density ribbon‐like and high density sheet‐like channel fills. Laterally to the fluvial system, mudstone and evaporitic mudstone units represented evaporitic mudflats which passed laterally into palustrine/lacustrine limestone units interpreted as lakes and ponds. Stacked calcretes occurred in distal alluvial and distal floodplain settings. A magnetostratigraphy encompassing 2600 m guided by available fossil mammal biochronology has provided a temporal framework that spans the complete Paleogene infill of the basin, from Late Lutetian to Late Oligocene, filling a gap in the Cenozoic chronostratigraphy of Spanish basins. This permits to constrain the kinematics of the structures both in the basin and in its margins, and to provide the timing for the depositional sequences. These data, combined with a magnetostratigraphic map, where magnetic reversals were traced through the Gómara monocline, allow a detailed analysis of the SR variability across the fluvial system and its adjacent depositional environments. The results show that high sedimentation rates (around 30–40 cm kyr^?1) are related to fluvial environments with low density ribbon‐shaped channels, while low SR (around or below 10 cm kyr^?1) are related to high density sheet‐like channels. Laterally, mud dominated environments with high SR (15–20 cm kyr^?1) grade into palustrine/lacustrine carbonated environments with low SR (around 9 cm ky^?1). The lowest SR (about 3 cm kyr^?1) are related to the development of stacked calcrete profiles in distal floodplain and in the connection of distal alluvial and palustrine/lacustrine units.