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...     

Impact of irrigated agriculture on groundwater-recharge salinity: a major sustainability concern in semi-arid regions

Intensive irrigated agriculture substantially modifies the hydrological cycle and often has major environmental impacts. The article focuses upon a specific concern—the tendency for progressive long-term increases in the salinity of groundwater recharge derived from irrigated permeable soils and replenishment of unconfined aquifers in more arid regions. This process has received only scant attention in the water-resource literature and has not been considered by agricultural science. This work makes an original contribution by analysing, from scientific principles, how the salinisation of groundwater recharge arises and identifies the factors affecting its severity. If not proactively managed, the process eventually will impact irrigation waterwell salinity, the productivity of agriculture itself, and can even lead to land abandonment. The types of management measure required for mitigation are discussed through three detailed case histories of areas with high-value groundwater-irrigated agriculture (in Spain, Argentina and Pakistan), which provide a long-term perspective on the evolution of the problem over various decades.     

Rainfall thresholds for the triggering of landslides considering previous soil moisture conditions (Asturias,NW Spain)

Given its geological and climatic conditions and its rugged orography, Asturias is one of the most landslide prone areas in the North of Spain. Most of the landslides occur during intense rainfall episodes. Thus, precipitation is considered the main triggering factor in the study area, reaching average annual values of 960 mm. Two main precipitation patterns are frequent: (i) long-lasting periods of moderate rainfall during autumn and winter and (ii) heavy short rainfall episodes during spring and early summer. In the present work, soil moisture conditions in the locations of 84 landslides are analysed during two rainfall episodes, which represent the most common precipitation patterns: October–November 2008 and June 2010. Empirical data allowed the definition of available water capacity percentages of 99–100% as critical soil moisture conditions for the landslide triggering. Intensity-duration rainfall thresholds were calculated for each episode, considering the periods with sustained high soil moisture levels before the occurrence of each analysed landslide event. For this purpose, data from daily water balance models and weather stations were used. An inverse relationship between the duration of the precipitation and its intensity, consistent with published intensity-duration thresholds, was observed, showing relevant seasonal differences.     

An adjoint-based optimization method for gas production in shale reservoirs

In this work, we consider a new model for flow in a multiporosity shale gas reservoir constructed within the framework of an upscaling procedure where hydraulic fractures are treated as (\(n-1\)) interfaces (\(n=2,3\)). Within this framework, the hydrodynamics is governed by a new pressure equation in the shale matrix which is treated as a homogenized porous medium composed of organic matter (kerogen aggregates with nanopores) and inorganic impermeable solid (clay, calcite, quartz) separated from each other by a network of interparticle pores of micrometer size. The solution of the pressure equation is strongly influenced by the constitutive response of the retardation parameter and effective hydraulic conductivity where the former incorporates gas adsorption/desorption in the nanopores of the kerogen. By focusing our analyses on this nonlinear diffusion equation in the domain occupied by the shale matrix, an optimization strategy seated on the adjoint sensitivity method is developed to minimize a cost functional related to gas production and net present value in a single hydraulic fracture. The gradient of the objective functional computed with the adjoint formulation is explored to update the controlled pressure drop aiming to optimize production in a given window of time. The combination of the direct approach and gradient-based optimization using the adjoint formulation leads to the construction of optimal production scenarios under controlled pressure decline in the well. Numerical simulations illustrate the potential of the methodology proposed herein in optimizing gas production.     

Multiple timescale constraints for high-flux magma chamber assembly prior to the Late Bronze Age eruption of Santorini (Greece)

The rhyodacitic magma discharged during the 30–80 km³ DRE (dense rock equivalent) Late Bronze Age (LBA; also called ‘Minoan’) eruption of Santorini caldera is known from previous studies to have had a complex history of polybaric ascent and storage prior to eruption. We refine the timescales of these processes by modelling Mg–Fe diffusion profiles in orthopyroxene and clinopyroxene crystals. The data are integrated with previously published information on the LBA eruption (phase equilibria studies, melt inclusion volatile barometry, Mg-in-plagioclase diffusion chronometry), as well as new plagioclase crystal size distributions and the established pre-LBA history of the volcano, to reconstruct the events that led up to the assembly and discharge of the LBA magma chamber. Orthopyroxene, clinopyroxene and plagioclase crystals in the rhyodacite have compositionally distinct rims, overgrowing relict, probably source-derived, more magnesian (or calcic) cores, and record one or more crystallization (plag???opx?>?cpx) events during the few centuries to years prior to eruption. The crystallization event(s) can be explained by the rapid transfer of rhyodacitic melt from a dioritic/gabbroic region of the subcaldera pluton (mostly in the 8–12 km depth range), followed by injection, cooling and mixing in a large melt lens at 4–6 km depth (the pre-eruptive magma chamber). Since crystals from all eruptive phases yield similar timescales, the melt transfer event(s), the last of which took place less than 2 years before the eruption, must have involved most of the magma that subsequently erupted. The data are consistent with a model in which prolonged generation, storage and segregation of silicic melts were followed by gravitational instability in the subcaldera pluton, causing the rapid interconnection and amalgamation of melt-rich domains. The melts then drained to the top of the pluton, at fluxes of up to 0.1–1 km³ year^??1, where steep vertical gradients of density and rheology probably caused them to inject laterally, forming a short-lived holding chamber prior to eruption. This interpretation is consistent with growing evidence that some large silicic magma chambers are transient features on geological timescales. A similar process preceded at least one earlier caldera-forming eruption on Santorini, suggesting that it may be a general feature of this rift-hosted magmatic system.     

An ecosystem-based composite spatial model for floodplain vulnerability assessment: a case study of Artigas,Uruguay

Floods are natural processes that constitute a hazard to society when associated to improper land use. Anthropic activities in floodplains are a factor of vulnerability that converts a natural hazard into a threat factor, eventually leading to disaster. Nowadays, natural and social complex processes demand integrated assessments in order to improve their understanding, helping decision making over sustainable use of territory, as well as integrating society’s activity in ecosystems and potentials, restrictions and benefits that society obtain from them. In this context, the objective of this work was to build a composite vulnerability model for a floodplain under urban influence, using an integrated assessment approach. This model was based on three dimensions; threat, fragility and an ecosystem services provision. These dimensions were calculated using both primary and secondary information, and weights by specialists. Main results show that the area presents high vulnerability with an increasing gradient towards high and urbanized areas, associated with an important number and relevant ecosystem services. Also, a spatial heterogeneity of the three dimensions emerged, making evident this area’s complexity and the need of integrated assessments to approach it. The composite vulnerability model proposed presents an elevated potential for natural and social processes analysis in floodplains, which is crucial for these territory management. Moreover, these integrated dimensions could contribute to decision making in different levels, as well as generating important supplies for environmental management and land planning.

     

The P-y Response of Laterally Loaded Flexible Piles in Residual Soil

This paper describes the main features related to lateral displacements with depth after successive lateral loading–unloading cycles applied to the top of reinforced-concrete flexible bored piles embedded in naturally bonded residual soil. The bored piles under study have a cylindrical shape, with 0.40-m in diameter and 8.0-m in length. Both bored piles types (P1 and P2) include an embedded steel pipe section in their center as longitudinal steel reinforcements: pile type P1 has another 16 steel rods as steel reinforcement to concrete while pile type P2 has no further steel reinforcement. Pile type P1 has three times as much stiffness (EI) and four and a half times the plastic moment (M_y) than pile type P2. A similar load–displacement performance was observed at initial loads as for small displacements of both piles. At this initial loading stage, the response of the reinforced concrete piles is a function of the soil characteristics and of a linear elastic pile deformation. During this stage, piles can even be understood as probes for evaluating soil reactions. For larger horizontal displacements, after the concrete section starts undergoing large deformations, approaching the ultimate bending moment, pile behavior and consequently the load–displacement relation starts to diverge for both piles. For pile P1 the values of relevant lateral displacements are extended to about 2.5-m in depth, while for pile P2 lateral displacements are mostly constrained to about 2.0-m in depth. Measurements of horizontal displacements of pile P1 against depth recorded with a slope indicator show that, after unloading, lateral loads at distinct stages (small and near failure loads), exhibits a much higher elastic phase of the system response. An analytical fitting model of soil reaction is proposed based on the measured displacements from slope indicator. The integration of a continuous model proposed for the soil reaction agrees fairly well with the measured displacements up to moments close to plastic limit. Results of load–displacement show that the stiffer pile (P1) was able to mobilize twice as much lateral load compared to pile P2 for a service limit displacement of about 20 mm. The paper shows results that enable the isolation of the structural variable through real scale pile load tests, thus granting understanding of its importance and enabling its quantitative visualization in examples of piles embedded in residual soil sites.

     

9000 years of changes in peat organic matter composition in Store Mosse (Sweden) traced using FTIR-ATR

Store Mosse (the ‘Great Bog’ in Swedish) is one of the most extensive bog complexes in southern Sweden (~77 km²), where pioneering palaeoenvironmental research has been carried out since the early 20th century. This includes, for example, vegetation changes, carbon and nitrogen dynamics, peat decomposition, atmospheric metal pollution, mineral dust deposition, dendrochronology, and tephrochronology. Even though organic matter (OM) represents the bulk of the peat mass and its compositional change has the potential to provide crucial ecological information on bog responses to environmental factors, peat OM molecular composition has not been addressed in detail. Here, a 568-cm-deep peat sequence was studied at high resolution, by attenuated reflectance Fourier-transform infrared spectroscopy (FTIR-ATR) in the mid-infrared region (4000–400 cm^–1). Principal components analysis was performed on selected absorbances and change-point modelling was applied to the records to determine the timing of changes. Four components accounted for peat composition: (i) depletion/accumulation of labile (i.e. carbohydrates) and recalcitrant (i.e. lignin and other aromatics, aliphatics, organic acids and some N compounds) compounds, due to peat decomposition; (ii) variations in N compounds and carbohydrates; (iii) residual variation of lignin and organic acids; and (iv) residual variation of aliphatic structures. Peat decomposition showed two main patterns: a long-term trend highly correlated to peat age (r = 0.87), and a short-term trend, which showed five main phases of increased decomposition (at ~8.4–8.1, ~7.0–5.6, ~3.5–3.1, ~2.7–2.1 and ~1.6–1.3 ka) – mostly corresponding to drier climate and its effect on bog hydrology. The high peat accumulation event (~5.6–3.9 ka), described in earlier studies, is characterized by the lowest degree of peat decomposition of the whole record. Given that FTIR-ATR is a quick, non-destructive, cost-effective technique, our results indicate that it can be applied in a systematic way (including multicore studies) to peat research and provide relevant information on the evolution of peatlands.     

The Colider and Roosevelt Volcanic rocks (sw amazonian craton): geochemistry and sm-nd isotope characteristics of a silicic large igneous province

The volcanic rocks of the Colíder and Roosevelt formations are extensively exposed in the south-central portion of the Amazonian Craton where effusive and pyroclastic rocks have been mapped. Both units, topped by chemical sediments and oceanic facies as rhyolite and andesite lavas, rhyodacite, and porphyritic dacite, with frequent intercalations of pyroclastic and epiclastic deposits. Whole-rock geochemistry for 55 samples of rhyolitic to andesitic composition suggests the involvement of fertile mantle-derived components with E-MORB to OIB compositions. The analyzed rocks display calc-alkaline to shoshonitic affinity consistent with generation related to an active continental margin. The whole-rock Sm-Nd isotope data from selected felsic volcanic rocks of the Colíder and Roosevelt formations yield negative initial ε_Nd values between –3 and –9, indicating the predominantly crustal nature of the parental magmas with early Archean to late Paleoproterozoic (ca. 2.5–2.0 Ga) depleted mantle model ages.