Analysis of soil drying incorporating a constitutive model for curling

This paper focuses on the shrinkage behavior of soil specimens involving sand, kaolinite, and kaolinite/sand mixtures subjected to desiccation under controlled conditions. Both, free and restrained shrinkage conditions are studied. The experiments show that pure soils do not curl upon unrestrained shrinkage; however, (under the same conditions) kaolinite/sand mixtures exhibited a marked curling. Furthermore, the mixture with the higher sand content broke through the middle of the sample after displaying a significant curling. Soils subjected to restricted shrinkage developed cracks with slight curling. To simulate the observed behavior, a mechanical model able to reproduce the detachment of the soil sample from the mold is proposed in this work and implemented in a fully coupled hydro-mechanical finite-element code. It is concluded that suction and differential shrinkage are key factors influencing the curling behavior of soils. The proposed framework was able to satisfactorily explain and reproduce the different stages and features of soil behavior observed in the experiments.

     

Climatic signals in lacustrine deposits of the Upper Yacoraite Formation,Western Argentina: Evidence from clay minerals,analcime, dolomite and fibrous calcite

Authigenic clays are an important control on reservoir quality in lacustrine carbonates but remain challenging to predict. Lacustrine depositional systems respond to climatic variations in rainfall, surface runoff and groundwater input, and evaporation, and result in rapid and frequent changes in lake volume; this is expressed through changing water depth and shoreline position. In the upper portion of the Early Palaeocene Yacoraite Formation of the Salta Basin in Argentina, extensive lacustrine deposits were deposited during the sag phase of rifting. Prior high-resolution stratigraphic studies have suggested that climatic factors control microbial carbonate sequences within a ‘balanced fill’ lake, with variation in the lake level having a major influence on facies association changes. This study characterizes the evolution of facies and mineralogy within the Yacoraite Formation, focusing on the distribution of clay minerals, making a link between the high, medium and low-frequency sequence stratigraphic cycles. The low-frequency transgressive hemicycle of the upper portion of the Yacoraite Formation is comprised of abundant siliciclastic facies, suggesting a wetter period. Microbialites occurring in this interval are coarse-grained and agglutinated. Detrital clay minerals such as illite and chlorite and associated siliciclastic sediments were input to the lake during high-frequency transgressive periods. During high-frequency regressive hemicycles, sedimentation was dominated by carbonate facies with Ca-rich dolomite and the authigenic clays are comprised of chlorite/smectite mixed-layers. By contrast, the low frequency regressive hemicycle records fine-grained agglutinated microbialite with horizons of fibrous calcite, more stoichiometric dolomite, barite and authigenic magnesian smectite. This indicates elevated ion concentrations in the lake under intense evaporation during an arid period. Understanding the conditions that are favourable for formation and preservation of authigenic clays within the lacustrine environment can improve understanding of reservoir quality in comparable economically important deposits.     

Separation of PP- and PS-wave reflected seismic data using two-dimensional finite offset common-reflection-surface traveltime approximation

Recently, the interest in PS-converted waves has increased for several applications, such as sub-basalt layer imaging, impedance estimates and amplitude-versus-offset analysis. In this study, we consider the problem of separation of PP- and PS-waves from pre-stacked multicomponent seismic data in two-dimensional isotropic medium. We aim to demonstrate that the finite-offset common-reflection-surface traveltime approximation is a good alternative for separating PP- and PS-converted waves in common-offset and common shot configurations by considering a two-dimensional isotropic medium. The five parameters of the finite-offset common-reflection-surface are firstly estimated through the inversion methodology called very fast simulated annealing, which estimates all parameters simultaneously. Next, the emergence angle, one of the inverted parameters, is used to build an analytical separation function of PP and PS reflection separation based on the wave polarization equations. Once the PP- and PS-converted waves were separated, the sections are stacked to increase the signal-to-noise ratio using the special curves derived from finite-offset common-reflection-surface approximation. We applied this methodology to a synthetic dataset from simple-layered to complex-structured media. The numerical results showed that the inverted parameters of the finite offset common-reflection-surface and the separation function yield good results for separating PP- and PS-converted waves in noisy common-offset and common shot gathers.     

Estimating free parameters in 2D inversion: example of gravity inversion in a rifted basement

We advance a principle directed to assigning numerical values to free parameters usually present in inversion methods. It may be formulated as: ‘Optimum estimates of free parameters in an inversion procedure must lead, in tests using synthetic data, to solutions whose geometrical expression reflects the main qualitative or semiquantitative geological characteristic of the study area.’ To this end, the interpreter should (i) specify a typical anomalous source geometry which incorporates the most relevant geological information for the study area, (ii) compute the corresponding gravity anomaly and (iii) invert the anomaly for the source geometry finding the numerical values of the free parameters that lead to a solution closest to the typical source. Application of the above methodology to synthetic and real data from the basement relief of a rift basin has asserted its efficacy.     

Post‐Great Oxidation Event Orosirian–Statherian iron formations on the São Francisco craton: Geotectonic implications

The protocratonic core of the São Francisco craton assembled during the 2.1–2.0 Ga Transamazonian orogeny. Orosirian Fe‐rich sequences that extend from the northwestern border of the São Francisco protocraton (Colomi Group) to the southeast under the Espinhaço Belt (the < 1.99 Ga Serra da Serpentina Group) record the opening of an intracratonic basin with the episodically developed ferruginous waters prior to the initiation of the Espinhaço rift at 1.8 Ga. Ferruginous conditions developed again during deposition of the Canjica Iron Formation of the < 1.7 Ga Serra de São José Group in the Espinhaço rift (contemporaneously with felsic magmatism; Conceição do Mato Dentro Rhyolite and Borrachudos Granitic Suite) and extensive sandstones of the < (1666 ±32) Ma Itapanhoacanga and < (1683 ±11) Ma São João da Chapada Formations. In the upper São João da Chapada Formation, banded hematitic phyllite also records input of Fe‐rich fluids. The young age of these iron formations with respect to the conventionally accepted 1.88 Ga age for the youngest shallow‐marine Paleoproterozoic iron formations, the apparent absence of granular facies (granular iron formations), and yet shallow‐water (above fair‐weather base) depositional environment indicate that an unusual setting developed in a large basin after the Great Oxidation Event, in the aftermath of the Transamazonian orogeny. We propose that mantle plumes led to the opening of a previously unrecognized rift system, that could have caused the magmatism, supplied hydrothermal Fe and led to the opening of the Espinhaço, Pirapora, and Paramirim rifts, later obliterated by the Araçuaí orogenic belt during the Neoproterozoic to Early Paleozoic Brasiliano orogeny. The rift system did not develop into an open continental margin but probably evolved into a broad sag basin, stretching across the São Francisco and Congo cratons.     

Efficacy of Coagulants and Ballast Compounds in Removal of Cyanobacteria (<Emphasis Type="Italic">Microcystis</Emphasis>) from Water of the Tropical Lagoon Jacarepaguá (Rio de Janeiro,Brazil)

Eutrophication is considered the most important water quality problem in freshwaters and coastal waters worldwide promoting frequent occurrence of blooms of potentially toxic cyanobacteria. Removal of cyanobacteria from the water column using a combination of coagulant and ballast is a promising technique for mitigation and an alternative to the use of algaecides. In laboratory, we tested experimentally the efficiency of two coagulants, polyaluminium chloride (PAC) and chitosan (made of shrimp shells), alone and combined with two ballasts: red soil (RS) and the own lagoon sediment, to remove natural populations of cyanobacteria, from an urban brackish coastal lagoon. PAC was a very effective coagulant when applied at low doses (≤8 mg Al L^?1) and settled the cyanobacteria, while at high doses (≥16 mg Al L^?1) large flocks aggregated in the top of test tubes. In contrast, chitosan was not able to form flocks, even in high doses (>16 mg L^?1) and did not efficiently settle down cyanobacteria when combined with ballast. The RS itself removed 33–47 % of the cyanobacteria. This removal was strongly enhanced when combined with PAC in a dose-dependent matter; 8 mg Al L^?1 was considered the best dose to be applied. The lagoon sediment alone did not promote any settling of cyanobacteria but removal was high when combined with PAC. Combined coagulant and ballast seems a very efficient, cheap, fast and safe curative measure to lessen the harmful cyanobacteria bloom nuisance in periods when particularly needed, such as around the 2016 Olympics in Jacarepaguá Lagoon.     

A New Appraisal of Sri Lankan BB Zircon as a Reference Material for LA‐ICP‐MS U‐Pb Geochronology and Lu‐Hf Isotope Tracing

A potential zircon reference material (BB zircon) for laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) U‐Pb geochronology and Hf isotope geochemistry is described. A batch of twenty zircon megacrysts (0.5–1.5 cm³) from Sri Lanka was studied. Within‐grain rare earth element (REE) compositions are largely homogeneous, albeit with some variation seen between fractured and homogeneous domains. Excluding fractured cathodoluminescence bright domains, the variation in U content for all analysed crystals ranged from 227 to 368 μg g^?1 and the average Th/U ratios were between 0.20 and 0.47. The Hf isotope composition (0.56–0.84 g/100 g Hf) is homogeneous within and between the grains – mean ¹⁷⁶Hf/¹⁷⁷Hf of 0.281674 ± 0.000018 (2s). The calculated alpha dose of 0.59 × 10¹⁸ g^?1 for a number of BB grains falls within the trend of previously studied, untreated zircon samples from Sri Lanka. Aliquots of the same crystal (analysed by ID‐TIMS in four different laboratories) gave consistent U‐Pb ages with excellent measurement reproducibility (0.1–0.4% RSD). Interlaboratory assessment (by LA‐ICP‐MS) from individual crystals returned results that are within uncertainty equivalent to the TIMS ages. Finally, we report on within‐ and between‐grain homogeneity of the oxygen isotope systematic of four BB crystals (13.16‰ VSMOW).     

Exploiting historical rainfall and landslide data in a spatial database for the derivation of critical rainfall thresholds

Critical rainfall thresholds for landslides are powerful tools for preventing landslide hazard. The thresholds are commonly estimated empirically starting from rainfall events that triggered landslides in the past. The creation of the appropriate rainfall–landslide database is one of the main efforts in this approach. In fact, an accurate agreement between the landslide and rainfall information, in terms of location and timing, is essential in order to correctly estimate the rainfall–landslide relationships. A further issue is taking into account the average moisture conditions prior the triggering event, which reasonably may be crucial in determining the sufficient amount of precipitation. In this context, the aim of this paper is exploiting historical landslide and rainfall data in a spatial database for the derivation of critical rainfall thresholds for landslide occurrence in Sicily, southern Italy. The hourly rainfall events that caused landslides occurred in the twentieth century were specifically identified and reconstructed. A procedure was proposed to automatically convert rain guages charts recorded on paper tape into digital format and then to provide the cumulative rainfall hyetograph in digital format. This procedure is based on a segmentation followed by signal recognition techniques which allow to digitalize and to recognize the hyetograph automatically. The role of rainfall prior to the landslide events was taken into account by including in the analysis the rainfall occurred 5, 15 and 30 days before each landslide. Finally, cumulated rainfall duration thresholds for different exceedance probability levels were determined. The obtained thresholds resulted in agreement with the regional curves proposed by other authors for the same area; antecedent rainfall turned out to be particularly important in triggering landslides.     

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.     

3D Gravity Inversion of Northern Sinai Peninsula: A Case Study

The Sinai Peninsula has attracted the attention of many geological and geophysical studies as it is influenced and bounded by major tectonic events. Those are (1) the Mesozoic to Early Cenozoic tectonically active opening of the Tethys, (2) the Late Cretaceous to Early Tertiary (Laramide) Syrian arc system, due to closing of the Tethys (3) the Oligo-Miocene Gulf of Suez rifted basin, and (4) the Late Miocene to Recent transform Dead Sea–Gulf of Aqaba rift. Additionally, the shear zones inside Sinai such as the Ragabet El-Naam and Minsherah-Abu Kandu Shear Zones. Each of these major tectonic events has affected dramatically the structure evolution of the northern Sinai area. The present paper estimates the 3D density contrast model using the gravity data of northern Sinai. The estimated 3D density contrast model elucidated the peculiarities of the main structural elements in the region. The estimated 3D density contrast model showed the high and low gravity anomalies that form the main mountains and main valleys in northern Sinai. The estimated low density zones are in agreement with the inferred faults resulting from the first horizontal derivative. Comparing the 3D model with the tectonic history of the region and the results of the first horizontal derivative and least square separation increased the reliability of the model.