Clay minerals in fluvial shaley sandstone of upper Triassic reservoir (TAGS)—Toual field SE Algeria: identification from wireline logs and core data

The creep property of rock under cyclic loading is very important in civil engineering. In order to establish a novel constitutive equation for rock under cyclic loading, a fractional-order viscoplastic body under cyclic loading was constructed based on fractional-order viscous element. A fractional-order visco-elastoplastic model (FVEPM) for rock was established by connecting constructed fractional-order viscoplastic body with Burgers model. The model was a Burgers model when the maximum value of cyclic loading was less than the critical strength of rock; otherwise, it was a FVEPM which can be used to reflect the transient, steady-state, and tertiary creep phases of rock. The cyclic loading was decomposed into a static load and a cyclic loading with a zero average stress. According to rheological mechanics theory, the rheology constitutive equation of rock under the static load can be derived. According to viscoelastic mechanics theory, the constitutive equation under cyclic loading with a zero average stress was established by introducing the variation parameters of energy storage and energy dissipation compliance caused by rock damage and fracture. Finally, a new dynamic constitutive equation of rock cyclic loading can be obtained by superimposing the constitutive equation under static load and cyclic loading with a zero average stress. Compared with existing test results of rock under cyclic loading, the proposed constitutive model can be used to describe the creep characteristics of rock under cyclic loading and reflect the presented fluctuation of strain curve of rock under cyclic loading.     

A geological record of multiple Pleistocene tsunami inundations in an oceanic island: The case of Maio,Cape Verde

In the Central Atlantic archipelagos – the Canaries, Cape Verde, Madeira and the Azores – tsunami hazard is often regarded as low, when compared with other extreme wave events such as hurricanes and storms. The geological record of many of these islands, however, suggests that tsunami hazard may be underestimated, notwithstanding being lower than in areas adjacent to subduction zones, such as the margins of the Pacific and Indian oceans. Moreover, tsunamis in oceanic islands are generally triggered by local large-scale volcanic flank collapses, for which little is known about their frequency, making it difficult to estimate the probability of a new occurrence. Part of the problem lies in the fact that tsunami deposits are usually difficult to date, and few islands in the world exhibit evidence for repeated tsunami inundation on a protracted timescale. This study reports on the presence of abundant tsunami deposits (conglomerates and sandstones) on Maio Island (Cape Verde) and discusses their stratigraphy, sedimentological characteristics, probable age and tsunamigenic source. Observations indicate that four distinct inundation events of variable magnitude took place during the Pleistocene. One of the tsunami deposits yielded a high-confidence U/Th age of 78·8 ± 0·9 ka, which overlaps within error with the 73 ± 7 ka age proposed for Fogo volcano's flank collapse, an event known to have had a significant tsunami impact on nearby Santiago Island. This shows that the Fogo tsunami also impacted Maio, resulting in runups in excess of 60 m above coeval sea-level at ca 120 km from the source. Two older deposits, possibly linked to recurrent flank collapses of the Tope de Coroa volcano in Santo Antão Island, yielded lower-confidence ages of 479 to 390 ka and 360 to 304 ka. A younger deposit (<78 ka) remains undated. In summary, the geological record of Maio exhibits well-preserved evidence of repeated tsunami inundation, reinforcing the notion that tsunami hazard is not so low at volcanic archipelagos featuring prominent and highly-active volcanoes such as in Cape Verde.     

Depositional history of the Miocene Lake Sinj (Dinaride Lake System,Croatia): a long-lived hard-water lake in a pull-apart tectonic setting

Early Miocene transpressional wrenching yielded a series of NW–SE-elongated pull-apart basins in the Dinarides of Croatia and Bosnia and Herzegovina. They accommodated a huge lake system that gave rise to spectacular endemic mollusk radiation. Lake Sinj, moderately sized at 342 km², flooded the south-westernmost basin of this system. Due to the karstic environment, the hard-water, alkaline, long-lived lake developed a sediment infill with an average thickness of 370 m, dominated by authigenic limestone. The studied section represents the upper third of the basinal infill and provides detailed insights into the critical period of the lake and of the basinal evolution during the final stages of its filling. It comprises two large-scale, shallowing-upward cycles, both starting with fossil-poor limestones, gradually passing into coal-bearing carbonate rocks and coal seams. The fossil-poor intervals are interpreted as phases of repetitive acidification events due to changing lake level, which induced periodic drying and flooding of the uppermost littoral zone inhabited by starfruit (Damasonium) meadows. The flooding of the aerated, limy mud plain introduced H+ ions from organic-matter decay reactions into the shallow lake. This decreased its pH level, with catastrophic consequences for its biota. The ecosystem then stabilized during the orbitally-forced, dry climate phases. Based on the mollusk record, streams still influenced the marginal lake environment and rich organic-matter production created swamps and mires. The onset of mollusk radiation in the section correlates with stabilized lake alkalinity, as indicated by the disappearance of starfruits, ongoing authigenic carbonate production and by coal seams representing textbook examples for coal formation in alkaline environments. The inferred basinal setting fits well with the pull-apart basin model, pointing to the presence of an extended shallow ramp in front of a steep, fault-induced hillside of the hinterland.     

Factors determining subsidence in urbanized floodplains: evidence from MT‐InSAR in Seville (southern Spain)

Major rivers have traditionally been linked with important human settlements throughout history. The growth of cities over recent river deposits makes necessary the use of multidisciplinary approaches to characterize the evolution of drainage networks in urbanized areas. Since under‐consolidated fluvial sediments are especially sensitive to compaction, their spatial distribution, thickness, and mechanical behavior must be studied. Here, we report on subsidence in the city of Seville (Southern Spain) between 2003 and 2010, through the analysis of the results obtained with the Multi‐Temporal InSAR (MT‐InSAR) technique. In addition, the temporal evolution of the subsidence is correlated with the rainfall, the river water column and the piezometric level. Finally, we characterize the geotechnical parameters of the fluvial sediments and calculate the theoretical settlement in the most representative sectors. Deformation maps clearly indicate that the spatial extent of subsidence is controlled by the distribution of under‐consolidated fine‐grained fluvial sediments at heights comprised in the range of river level variation. This is clearly evident at the western margin of the river and the surroundings of its tributaries, and differs from rainfall results as consequence of the anthropic regulation of the river. On the other hand, this influence is not detected at the eastern margin due to the shallow presence of coarse‐grain consolidated sediments of different terrace levels. The derived results prove valuable for implementing urban planning strategies, and the InSAR technique can therefore be considered as a complementary tool to help unravel the subsidence tendency of cities located over under‐consolidated fluvial deposits. Copyright © 2017 John Wiley & Sons, Ltd.     

Speleothems in a north Cuban cave register sea‐level changes and Pleistocene uplift rates

A flight of marine terraces along the Cuban coast records Quaternary sea‐level highstands and a general slowly uplifting trend during the Pleistocene. U/Th dating of these limestone terraces is difficult because fossil reef corals have been affected by open system conditions. Terrace ages are thus often based on geological and geomorphological observations. In contrast, the minimum age of the terraces can be constrained by dating speleothems from coastal mixing (flank margin) caves formed during past sea‐level highstands and carving the marine limestones. Speleothems in Santa Catalina Cave have ages >360 ka and show various cycles of subaerial–subaqueous corrosion and speleothem growth. This suggests that the cave was carved during the MIS 11 sea‐level highstand or earlier. Some stalagmites grew during MIS 11 through MIS 8 and were submerged twice, once at the end of MIS 11 and then during MIS 9. Phreatic overgrowths (POS) covering the speleothems suggest anchialine conditions in the cave during MIS 5e. Their altitude at 16 m above present sea level indicates a late Pleistocene uplift rate of <0.1 mm/ka, but modelling also shows uplift to have been insignificant over a long timespan during the middle Pleistocene since the cave was carved. Our study shows that some flank margin caves in the region of Matanzas are older than commonly believed (i.e. MIS 11 rather than MIS 5). These caves not only can be preserved but are good markers of interglacial sea‐level highstands, more reliable than marine abrasion surfaces. Copyright © 2018 John Wiley & Sons, Ltd.     

Better models are more effectively connected models

Water‐ and sediment‐transfer models are commonly used to explain or predict patterns in the landscape at scales different from those at which observations are available. These patterns are often the result of emergent properties that occur because processes of water and sediment transfer are connected in different ways. Recent advances in geomorphology suggest that it is important to consider, at a specific spatio‐temporal scale, the structural connectivity of system properties that control processes, and the functional connectivity resulting from the way those processes operate and evolve through time. We argue that a more careful consideration of how structural and functional connectivity are represented in models should lead to more robust models that are appropriate for the scale of application and provide results that can be upscaled. This approach is necessary because, notwithstanding the significant advances in computer power in recent years, many geomorphic models are still unable to represent the landscape in sufficient detail to allow all connectivity to emerge. It is important to go beyond the simple representation of structural connectivity elements and allow the dynamics of processes to be represented, for example by using a connectivity function. This commentary aims to show how a better representation of connectivity in models can be achieved, by considering the sorts of landscape features present, and whether these features can be represented explicitly in the model spatial structure, or must be represented implicitly at the subgrid scale. Copyright © 2017 John Wiley & Sons, Ltd.     

Geostatistical facies simulation with geometric patterns of a petroleum reservoir

During exploration and pre-feasibility studies of a typical petroleum project many analyses are required to support decision making. Among them is reservoir lithofacies modeling, preferably using uncertainty assessment, which can be carried out with geostatistical simulation. The resulting multiple equally probable facies models can be used, for instance, in flow simulations. This allows assessing uncertainties in reservoir flow behavior during its production lifetime, which is useful for injector and producer well planning. Flow, among other factors, is controlled by elements that act as flow corridors and barriers. Clean sand channels and shale layers are examples of such reservoir elements that have specific geometries. Besides simulating the necessary facies, it is also important to simulate their shapes. Object-based and process-based simulations excel in geometry reproduction, while variogram-based simulations perform very well at data conditioning. Multiple-point geostatistics (MPS) combines both characteristics, consequently it was employed in this study to produce models of a real-world reservoir that are both data adherent and geologically realistic. This work aims at illustrating how subsurface information typically available in petroleum projects can be used with MPS to generate realistic reservoir models. A workflow using the SNESIM algorithm is demonstrated incorporating various sources of information. Results show that complex structures (e.g. channel networks) emerged from a simple model (e.g. single branch) and the reservoir facies models produced with MPS were judged suitable for geometry-sensitive applications such as flow simulations.     

Sangay volcano, Ecuador: structural development, present activity and petrology

Sangay (5230 m), the southernmost active volcano of the Andean Northern Volcanic Zone (NVZ), sits 130 km above a >32-Ma-old slab, close to a major tear that separates two distinct subducting oceanic crusts. Southwards, Quaternary volcanism is absent along a 1600-km-long segment of the Andes. Three successive edifices of decreasing volume have formed the Sangay volcanic complex during the last 500 ka. Two former cones (Sangay I and II) have been largely destroyed by sector collapses that resulted in large debris avalanches that flowed out upon the Amazon plain. Sangay III, being constructed within the last avalanche amphitheater, has been active at least since 14 ka BP. Only the largest eruptions with unusually high Plinian columns are likely to represent a major hazard for the inhabited areas located 30 to 100 km west of the volcano. However, given the volcano's relief and unbuttressed eastern side, a future collapse must be considered, that would seriously affect an area of present-day colonization in the Amazon plain, 30 km east of the summit. Andesites greatly predominate at Sangay, there being few dacites and basalts. In order to explain the unusual characteristics of the Sangay suite—highest content of incompatible elements (except Y and HREE) of any NVZ suite, low Y and HREE values in the andesites and dacites, and high Nb/La of the only basalt found—a preliminary five-step model is proposed: (1) an enriched mantle (in comparison with an MORB source), or maybe a variably enriched mantle, at the site of the Sangay, prior to Quaternary volcanism; (2) metasomatism of this mantle by important volumes of slab-derived fluids enriched in soluble incompatible elements, due to the subduction of major oceanic fracture zones; (3) partial melting of this metasomatized mantle and generation of primitive basaltic melts with Nb/La values typical of the NVZ, which are parental to the entire Sangay suite but apparently never reach the surface and subordinate production of high Nb/La basaltic melts, maybe by lower degrees of melting at the periphery of the main site of magma formation, that only infrequently reach the surface; (4) AFC processes at the base of a 50-km-thick crust, where parental melts pond and fractionate while assimilating remelts of similar basaltic material previously underplated, producing andesites with low Y and HREE contents, due to garnet stability at this depth; (5) low-pressure fractionation and mixing processes higher in the crust. Both an enriched mantle under Sangay prior to volcanism and an important slab-derived input of fluids enriched in soluble incompatible elements, two parameters certainly related to the unique setting of the volcano at the southern termination of the NVZ, apparently account for the exceptionally high contents of incompatible elements of the Sangay suite. In addition, the low Cr/Ni values of the entire suite—another unique characteristic of the NVZ—also requires unusual fractionation processes involving Cr-spinel and/or clinopyroxene, either in the upper mantle or at the base of the crust.     

An evaluation of wave propagation simulations over a barred beach with a Boussinesq-type model

The FUNWAVE model is used for simulating simulation of monochromatic and irregular wave propagation in a channel with a bar-trough profile. FUNWAVE is based upon the extended Boussinesq equations. The study aims to analyze the model's performance when simulating shoaling, wave breaking and nonlinear interactions that are present in nearshore wave propagation. For that, high-order time domain statistics (root mean-square wave height, skewness, asymmetry and the kurtosis) of the model simulations and of the observations were compared along the whole channel. Also, a frequency domain analysis including standard spectral analysis and the bispectrum was carried out in selected points of the flume. The evaluation included the role of the wave breaking internal model parameters. The main conclusion is that, in general, the one-dimensional version of FUNWAVE simulates quite well the nonlinear transformation of a wave over a bottom with a bar-tough profile, for both regular and irregular wave conditions. The model reproduces the transformation of the wave shape, specially the increasing sharper wave crests and flatter troughs and also the lack of vertical symmetry with crests pitching forward, as it propagates along the domain. However, some differences persist after wave breaking, mainly due to the nature of the wave-breaking module. In this module, the energy dissipation is induced by the increase of viscosity, a rather simple mechanism, without the modification of the wave shape. Also, the energy dissipation develops in a smooth way which is appropriated for spilling breaking waves, but not for plunging breaking waves where the dissipation starts more abruptly.