Constraining magma sources using primitive olivine-hosted melt inclusions from Puñalica and Sangay volcanoes (Ecuador)

Constraining arc magma sources at continental arc settings is a delicate task, because chemical signatures from crustal processes obscure the slab and mantle signatures. Here, we present major, trace, and volatile element compositions of olivine-hosted melt inclusions (Fo_82–89) selected from the most primitive lavas (Mg# >?60) from two Ecuadorian volcanoes (Puñalica and Sangay) situated at the southern termination of the Andean Northern Volcanic Zone. Melt inclusions (MI) from Puñalica are nepheline normative and have basaltic-to-basaltic-andesite compositions (45–56 wt% SiO₂) similar to peridotite-derived melts. Sangay MI is also nepheline normative, with high CaO (up to 16 wt% and CaO/Al₂O₃?<?1) and low silica contents (41.9–44.5 wt%) pointing out an amphibole-bearing clinopyroxenite source. Both volcanoes display volatile-rich compositions (up to 6100 ppm Cl, 2200 ppm F, and 6700 ppm S). These MI cannot be related to their host lavas by fractional crystallization, implying that they represent true primitive liquids. The source of Puñalica MI was metasomatized by slab-derived melts that imprints its low Ba/Th, Sr/Th, and high Th/La (average values of 66, 129, and 0.22, respectively). On the contrary, the slab component added to the source of Sangay MI has a higher Ba/Th, Sr/Th, and low Th/La (average values of 261, 517, and 0.11, respectively) which could suggest a relative contribution of aqueous fluids. This dichotomy is related to the presence of the Grijalva Fracture Zone that separates a younger and hotter oceanic crust to the north (below Puñalica) from a colder and older oceanic crust to the south (below Sangay).     

An elasto-visco-plastic model using the finite element method for crustal and lithospheric deformation

A novel numerical model based on solid deformation is presented in this paper. This thermo-mechanical model can simulate the tectonic evolution of crust and (lithospheric and asthenospheric) mantle under different conditions. Our implementation uses the finite element method (FEM) in order to solve the equations. As a Lagrangian approach is employed, remeshing techniques are implemented to avoid distortion problems when a certain deformation threshold is reached. The translation of the state between the old and new mesh is achieved by means of the information stored on Lagrangian particles, which minimizes the diffusion. The model is able to represent elastic, viscous and plastic behaviour inside the studied domain. Three types of creep mechanism (diffusion, dislocation and Peierls) are included. Two different quadrilateral isoparametric elements were implemented and can be employed to perform the calculations. The first one is an element with 4 nodes, selective reduced integration and a stabilization operator to diminish hourglass modes, which reduces the computational time needed. The second one has 8 nodes located in standard positions, uses full integration scheme and has no hourglass modes as it satisfies the Inf-Sup condition. Several test cases with known solutions were run to validate the different aspects of the implementation.     

The Chimborazo sector collapse and debris avalanche: Deposit characteristics as evidence of emplacement mechanisms

Chimborazo is a Late Pleistocene to Holocene stratovolcano located at the southwest end of the main Ecuadorian volcanic arc. It experienced a large sector collapse and debris avalanche (DA) of the initial edifice (CH-I). This left a 4 km wide scar, removing 8.0 ± 0.5 km³ of the edifice. The debris avalanche deposit (DAD) is abundantly exposed throughout the Riobamba Basin to the Río Chambo, more than 35 km southeast of the volcano. The DAD averages a thickness of 40 m, covers about 280 km², and has a volume of > 11 km³. Two main DAD facies are recognized: block and mixed facies. The block facies is derived predominantly from edifice lava and forms > 80 vol.% of the DAD, with a probable volume increase of 15–25 vol.%. The mixed facies was essentially created by mixing brecciated edifice rock with substratum and is found mainly in distal and marginal areas. The DAD has clear surface ridges and hummocks, and internal structures such as jigsaw cracks, injections, and shear-zone features are widespread. Structures such as stretched blocks along the base contact indicate high basal shear. Substratum incorporation is directly observed at the base and is inferred from the presence of substratum-derived material in the DAD body. Based on the facies and structural interpretation, we propose an emplacement model of a lava-rich avalanche strongly cataclased before and/or during failure initiation. The flow mobilises and incorporates significant substrata (10–14 vol.%) while developing a fine lubricating basal layer. The substrata-dominated mixed facies is transported to the DAD interior and top in dykes invading previously-formed fractures.     

Time Evolution of Deformation Using Time Series of Differential Interferograms: Application to La Palma Island (Canary Islands)

Differential interferometry is a very powerful tool for detecting changes in the Earth’s crust where coherence conditions are good, but is difficult to employ in some volcanic areas due to dense vegetation. We apply two differential InSAR methods using the time series associated with the interferograms to perform a phase analysis on a data set for La Palma island (Canary Islands) from the ERS-1 and ERS-2 European Space Agency (ESA) satellites for the time period 1992 to 2000. Both methods involve choosing a master image from the database and creating a series of interferograms with respect to this image. The “Coherent Pixel Time Series” (CPTS) technique chooses pixels with good average coherence, aligns the unwrapped interferograms with a stable area and then performs an inversion to calculate the linear velocity to quantify the deformation. The Coherent Target Modeling (CTM) method calculates the temporal coherence of each pixel to identify stable targets and then determines the best velocity for each pixel by using a linear fit that maximizes the temporal coherence. Using these two methods we have been able to detect deformation on La Palma Island that has been previously undetectable by conventional InSAR methods. There is a roughly circular region on the Southern part of the island that is actively deforming at ~ −4 to −8 mm/yr. This region is located near the Teneguia valcano, the host of the last known eruption on La Palma in 1971. A thorough investigation of the possible sources for this deformation revealed that it was most likely created by a subsurface thermal source.     

Morphodynamics of alternate bars in the presence of riparian vegetation

Alpine gravel-bed rivers are dynamic systems that have been subjected to many anthropic alterations in the past centuries. Riparian vegetation development on previously bare sediment bedforms has been a common adjustment, raising important management issues in terms of flood risks and biodiversity. Many of these rivers are also channelized, and as a result present a pattern of alternate bars. Considering recent advances in numerical biomorphodynamic modeling, this study aims at exploring numerically the morphodynamics of alternate bars in the presence of riparian vegetation. To this end, a dynamic vegetation module has been implemented on top of an existing morphodynamic model, accounting for ecological processes of seed dispersal, seedling recruitment, growth, and mortality. Numerical simulations have been performed on a simplified reach of a gravel-bed river with free migrating alternate bars at initial state. In this work 96 scenarios have been simulated, each representing 50 years of channel evolution, with different flood regimes characterized by various peak discharges and flood durations. Yearly peak discharge variability is explicitly modeled in 48 scenarios. Model outcomes present two possible equilibrium biomorphodynamic behaviors: stationary vegetated bars, or free migrating bars in the case of frequent vegetation removal during floods. This binary behavior holds true when the stochasticity of annual peak discharges is represented, and for a wide range of parameter values included in vegetation dynamic modeling. Transient mobility of vegetated bars is observed in specific configurations where large sediment deposits deflect the flow field, eroding bar heads. Modeled bar wavelengths are in the range of values predicted for free bars by linear bar theory, and remain far from the theoretical values of hybrid, steady bars. The shift from unvegetated migrating bars to steady vegetated bars seems to show that in these simulations vegetation constitutes a hydraulic forcing, leading to a shift from free bars to forced bars, with a final configuration largely inherited from the initial state. © 2019 John Wiley & Sons, Ltd.     

Reducing systematic dome errors in digital elevation models through better UAV flight design

It is well established that digital elevation models (DEMs) derived from unmanned aerial vehicle (UAV) images and processed by structure from motion may contain important systematic vertical errors arising from limitations in camera geometry modelling. Even when significant, such ‘dome’-shaped errors can often remain unnoticed unless specific checks are conducted. Previous methods used to reduce these errors have involved: the addition of convergent images to supplement traditional vertical datasets, the usage of a higher number of ground control points, precise direct georeferencing techniques (RTK/PPK) or more refined camera pre-calibration. This study confirms that specific UAV flight designs can significantly reduce dome errors, particularly those that have a higher number of tie points connecting distant images, and hence contribute to a strengthened photogrammetric network. A total of 22 flight designs were tested, including vertical, convergent, point of interest (POI), multiscale and mixed imagery. Flights were carried out over a 300 × 70 m² flat test field area, where 143 ground points were accurately established. Three different UAVs and two commercial software packages were trialled, totalling 396 different tests. POI flight designs generated the smallest systematic errors. In contrast, vertical flight designs suffered from larger dome errors; unfortunately, a configuration that is ubiquitous and most often used. By using the POI flight design, the accuracy of DEMs will improve without the need to use more ground control or expensive RTK/PPK systems. Over flat terrain, the improvement is especially important in self-calibration projects without (or with just a few) ground control points. Some improvement will also be observed on those projects using camera pre-calibration or with stronger ground control. © 2020 John Wiley & Sons, Ltd.     

Mid Campanian‐Lower Maastrichtian magnetostratigraphy of the James Ross Basin,Antarctica: Chronostratigraphical implications

The James Ross Basin, in the northern Antarctic Peninsula, exposes which is probably the world thickest and most complete Late Cretaceous sedimentary succession of southern high latitudes. Despite its very good exposures and varied and abundant fossil fauna, precise chronological determination of its infill is still lacking. We report results from a magnetostratigraphic study on shelfal sedimentary rocks of the Marambio Group, southeastern James Ross Basin, Antarctica. The succession studied covers a ~1,200 m‐thick stratigraphic interval within the Hamilton Point, Sanctuary Cliffs and Karlsen Cliffs Members of the Snow Hill Island Formation, the Haslum Crag Formation, and the lower López de Bertodano Formation. The basic chronological reference framework is given by ammonite assemblages, which indicate a Late Campanian – Early Maastrichtian age for the studied units. Magnetostratigraphic samples were obtained from five partial sections located on James Ross and Snow Hill islands, the results from which agree partially with this previous biostratigraphical framework. Seven geomagnetic polarity reversals are identified in this work, allowing to identify the Chron C32/C33 boundary in Ammonite Assemblage 8‐1, confirming the Late Campanian age of the Hamilton Point Member. However, the identification of the Chron C32/C31 boundary in Ammonite Assemblage 8‐2 assigns the base of the Sanctuary Cliffs Member to the early Maastrichtian, which differs from the Late Campanian age previously assigned by ammonite biostratigraphy. This magnetostratigraphy spans ~14 Ma of sedimentary succession and together with previous partial magnetostratigraphies on Early‐Mid Campanian and Middle Maastrichtian to Danian columns permits a complete and continuous record of the Late Cretaceous distal deposits of the James Ross Basin. This provides the required chronological resolution to solve the intra‐basin and global correlation problems of the Late Cretaceous in the Southern Hemisphere in general and in the Weddellian province in particular, given by endemism and diachronic extinctions on invertebrate fossils, including ammonites. The new chronostratigraphic scheme allowed us to calculate sediment accumulation rates for almost the entire Late Cretaceous infill of the distal James Ross Basin (the Marambio Group), showing a monotonous accumulation for more than 8 Myr during the upper Campanian and a dramatic increase during the early Maastrichtian, controlled by tectonic and/or eustatic causes.     

Structural evolution and the partitioning of deformation during basin growth and inversion: A case study from the Mizen Basin Celtic Sea,offshore Ireland

The Celtic Sea basins lie on the continental shelf between Ireland and northwest France and consist of a series of ENE–WSW trending elongate basins that extend from St George’s Channel Basin in the east to the Fastnet Basin in the west. The basins, which contain Triassic to Neogene stratigraphic sequences, evolved through a complex geological history that includes multiple Mesozoic rift stages and later Cenozoic inversion. The Mizen Basin represents the NW termination of the Celtic Sea basins and consists of two NE–SW-trending half-grabens developed as a result of the reactivation of Palaeozoic (Caledonian, Lower Carboniferous and Variscan) faults. The faults bounding the Mizen Basin were active as normal faults from Early Triassic to Late Cretaceous times. Most of the fault displacement took place during Berriasian to Hauterivian (Early Cretaceous) times, with a NW–SE direction of extension. A later phase of Aptian to Cenomanian (Early to Late Cretaceous) N–S-oriented extension gave rise to E–W-striking minor normal faults and reactivation of the pre-existing basin bounding faults that propagated upwards as left-stepping arrays of segmented normal faults. In common with most of the Celtic Sea basins, the Mizen Basin experienced a period of major erosion, attributed to tectonic uplift, during the Paleocene. Approximately N–S Alpine regional compression-causing basin inversion is dated as Middle Eocene to Miocene by a well-preserved syn-inversion stratigraphy. Reverse reactivation of the basin bounding faults was broadly synchronous with the formation of a set of near-orthogonal NW–SE dextral strike-slip faults so that compression was partitioned onto two fault sets, the geometrical configuration of which is partly inherited from Palaeozoic basement structure. The segmented character of the fault forming the southern boundary of the Mizen Basin was preserved during Alpine inversion so that Cenozoic reverse displacement distribution on syn-inversion horizons mirrors the earlier extensional displacements. Segmentation of normal faults therefore controls the geometry and location of inversion structures, including inversion anticlines and the back rotation of earlier relay ramps.     

Diatom assemblages from surface sediments of the Río de la Plata estuary,Uruguay

The Río de la Plata estuary (RdlP) exhibits environmental gradients associated with the freshwater input and oceanic water intrusion. The aim of this study was to assess diatom species distribution in surface sediment samples related to such environmental gradients. The internal section of RdlP was dominated by Aulacoseira spp., Eunotia spp., Staurosirella martyi, Actinocyclus normanii and Thalassiosira baltica, indicatives of low salinity levels and high trophic conditions, associated with the riverine and estuarine regimes. The external section was dominated by Coscinodiscus radiatus, Thalassiosira spp., Paralia sulcata, Cyclotella striata, among other marine taxa, indicatives of high salinity and low trophic conditions, associated with the influence of the Southwestern Atlantic Ocean. Furthermore, the intermediate section presents a mixture of both diatom groups representing mixing conditions. The observed diatom species groups capture fairly well the RdlP environmental variability and can be reliably used for paleoenvironmental studies in this and other similar estuarine systems.