Effects of crab–halophytic plant interactions on creek growth in a S.W. Atlantic salt marsh: A Cellular Automata model

The Bahía Blanca Estuary (38° 50′ S, and 62° 30′ W) presents salt marshes where interactions between the local flora (Sarcocornia perennis) and fauna (Chasmagnathus granulatus) generate some kind of salt pans that alter the normal water circulation and condition its flow and course towards tidal creeks. The crab–vegetation dynamics in the salt marsh presents variations that cannot be quantified in a reasonable period of time. The interaction between S. perennis plant and C. granulatus crab is based on simple laws, but its result is a complex biological mechanism that causes an erosive process on the salt marsh and favors the formation of tidal creeks. To study it, a Cellular Automata model is proposed, based on the laws deduced from the observation of these phenomena in the field, and then verified with measurable data within macroscale time units. Therefore, the objective of this article is to model how the interaction between C. granulatus and S. perennis modifies the landscape of the salt marsh and influences the path of tidal creeks. The model copies the basic laws that rule the problem based on purely biological factors.The Cellular Automata model proved capable of reproducing the effects of the interaction between plants and crabs in the salt marsh. A study of the water drainage of the basins showed that this interaction does indeed modify the development of tidal creeks. Model dynamics would likewise follow different laws, which would provide a different formula for the probability of patch dilation. The patch shape can be obtained changing the pattern that dilates.     

Effect of Multipoint Heterogeneity on Nonlinear Transformations for Geological Modeling： Porosity-Permeability Relations Revisited

An analysis of statistical expected values for transformations is performed in this study to quantify the effect of heterogeneity on spatial geological modeling and evaluations. Algebraic transformations are frequently applied to data from logging to allow for the modeling of geological properties. Transformations may be powers, products, and exponential operations which are commonly used in well-known relations (e.g., porosity-permeability transforms). The results of this study show that correct computations must account for residual transformation terms which arise due to lack of independence among heterogeneous geological properties. In the case of an exponential porosity-permeability transform, the values may be positive. This proves that a simple exponential model back-transformed from linear regression underestimates permeability. In the case of transformations involving two or more properties, residual terms may represent the contribution of heterogeneous components which occur when properties vary together, regardless of a pair-wise linear independence. A consequence of power- and product-transform models is that regression equationswithin those transformations need corrections via residual cumulants. A generalization of this result isthat transformations of multivariate spatial attributes require multiple-point random variable relations. This analysis provides practical solutions leading to a methodology for nonlinear modeling using correct back transformations in geology.     

Comparison of dissolved and particulate arsenic distributions in shallow aquifers of Chakdaha, India, and Araihazar, Bangladesh

Background  

The origin of the spatial variability of dissolved As concentrations in shallow aquifers of the Bengal Basin remains poorly understood. To address this, we compare here transects of simultaneously-collected groundwater and aquifer solids perpendicular to the banks of the Hooghly River in Chakdaha, India, and the Old Brahmaputra River in Araihazar, Bangladesh.     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.     

Contributions of boron isotopes to understanding the hydrogeochemistry of the coastal detritic aquifer of Castellón Plain,Spain

The Castellón Plain alluvial aquifer, Spain, is intensively exploited to meet the demand for agricultural irrigation and industrial water supply. The geochemistry of its groundwater shows complex salinization in the northern and southern parts of the aquifer, with significant pollution from human origin in the central portion. Boron content and B isotope geochemistry are useful for distinguishing between various sources of pollution and their relative importance in different parts of this aquifer. Boron concentrations in the groundwater vary between 0.01 and 0.85 mg/L. In the more saline groundwaters, found at the northern and southern ends of the study area, the presence of B is linked to inputs from seawater and water with a calcium-magnesium sulphate facies, which feed the aquifer and clearly influence the chemistry of its waters. Evidence of B adsorption processes in some samples is shown by the low B/Cl ratios and the high values of δ¹¹B. In the central portion of the aquifer, the high B/Cl ratios and the strongly negative δ¹¹B are related to pollution of human origin.     

New age constraints on emplacement of the Cévenol granitoids, South French Massif Central

During the development of the Variscan orogeny, large amounts of granitic melt were produced, giving rise to the intrusion of granitoids at different structural levels. Despite numerous studies, ages available from previous work on the Cévennes granites remain largely imprecise. In order to better constrain the age and emplacement mode of these granites, we have combined U–Pb dating on monazites and zircons and ⁴⁰Ar/³⁹Ar dating on biotites with petrological observations, major element chemical analysis and SEM zircon imaging on five samples from the Aigoual–St Guiral–Liron and Mont Lozère granitic massifs. The results revealed that granitic intrusions and cooling in Southern Cévennes occurred in a short time span at ∼306 Ma after the main episode of regional metamorphism. Petrological and chemical data suggest that they result from a mixing between mantle-derived basic magmas (lamprophyres) and lower crust acid magmas. At a regional scale the production of these melts occurred at the end of crustal thickening induced by nappe stacking, at the same time as the late anatectic events recorded further north in the Velay dome and the granulite facies metamorphism recorded in metasedimentary granulite enclaves brought up by Tertiary volcanoes of the Velay area (Bournac).     

Giant ancient landslide in the Alma water gap (Crimean Mountains,Ukraine): notes to the predisposition,structure, and chronology

Large-scale ancient landslides of the area of more than 5 km² and volume exceeding 200 × 10⁶ m³ are characteristic features of the valleys incised in the northern periphery of the Crimean Mountains (Ukraine). The largely affected area is located in the outermost cuesta range of the Crimean Mountains which consists of rigid Sarmatian limestones overlying weak Middle Miocene and Upper Palaeogene deposits. A giant landslide arose in the Alma water gap as a reflection of several coincident preparatory factors such as suitable bedrock stratification, smectite-rich bedrock exposed to swelling activity, presence of faults parallel to the valley trend, and river capture event which preceded the landslide event. The occurrence of such ancient megaslides is particularly interesting in the area which is characterized by low precipitation (<500 mm/year) and weak contemporary seismicity. It probably reflects a more dynamic environment in humid phases of the Holocene; however, seismic triggering along the Mesozoic suture zone cannot be rejected. Compressional features such as gravitational folds in the central and distal parts of the landslide, which probably correlate with the whole landslide genesis or its significant reactivation, arose, according to the radiocarbon dating, during the Holocene climatic optimum in the Atlantic period. The slope deformation has been relatively quiescent since that time, except minor historic reactivization which took place in the frontal part of the landslide. We suppose that the studied landslide could be classified as a transitional type of slope deformation with some signs of spreading and translational block slides.     

Provenance of the Bosnian Flysch

Sandwiched between the Adriatic Carbonate Platform and the Dinaride Ophiolite Zone, the Bosnian Flysch forms a c. 3000 m thick, intensely folded stack of Upper Jurassic to Cretaceous mixed carbonate and siliciclastic sediments in the Dinarides. New petrographic, heavy mineral, zircon U/Pb and fission-track data as well as biostratigraphic evidence allow us to reconstruct the palaeogeology of the source areas of the Bosnian Flysch basin in late Mesozoic times. Middle Jurassic intraoceanic subduction of the Neotethys was shortly followed by exhumation of the overriding oceanic plate. Trench sedimentation was controlled by a dual sediment supply from the sub-ophiolitic high-grade metamorphic soles and from the distal continental margin of the Adriatic plate. Following obduction onto Adria, from the Jurassic–Cretaceous transition onwards a vast clastic wedge (Vranduk Formation) was developed in front of the leading edge, fed by continental basement units of Adria that experienced Early Cretaceous synsedimentary cooling, by the overlying ophiolitic thrust sheets and by redeposited elements of coeval Urgonian facies reefs grown on the thrust wedge complex. Following mid-Cretaceous deformation and thermal overprint of the Vranduk Formation, the depozone migrated further towards SW and received increasing amounts of redeposited carbonate detritus released from the Adriatic Carbonate Platform margin (Ugar Formation). Subordinate siliciclastic source components indicate changing source rocks on the upper plate, with ophiolites becoming subordinate. The zone of the continental basement previously affected by the Late Jurassic–Early Cretaceous thermal imprint has been removed; instead, the basement mostly supplied detritus with a wide range of pre-Jurassic cooling ages. However, a c. 80 Ma, largely synsedimentary cooling event is also recorded by the Ugar Formation, that contrasts the predominantly Early Cretaceous cooling of the Adriatic basement and suggests, at least locally, a fast exhumation.