Glacial lake depth and volume estimation based on a large bathymetric dataset from the Cordillera Blanca,Peru

Glacial lakes are most often located in remote places making it difficult to carry out detailed bathymetric surveys. Consequently, lake depths and volumes for unmeasured lakes are often estimated using empirical relationships developed mainly from small bathymetric datasets. In this study, we use the bathymetry dataset of the Cordillera Blanca, Peru comprising 121 detailed lake bathymetries, the most extensive dataset in the world. We assess the performance of the most commonly applied empirical relationships for lake mean depth and volume estimation, but also investigate relationships between different geometric lake variables. We find that lake volume estimation performs better when derived from lake mean depth, which in turn is estimated from lake width. The findings also reveal the extreme variability of lake geometry, which depends on glacio-geomorphological processes that empirical–statistical relationships cannot adequately represent. Such relationships involve characteristic uncertainty ranges of roughly ±50%. We also estimate potential peak discharges of outburst floods from these lakes by applying empirical relationships from the literature, which results in discharges varying by up to one-order of magnitude. Finally, the results are applied to the 860 lakes without bathymetric measurements from the inventory dataset of the Cordillera Blanca to estimate lake mean depth, volume and possible peak discharge for all unmeasured lakes. Estimations show that ca. 70% (610) of the lakes have a mean depth lower than 10 m and very few longer than 40 m. Lake volume of unmeasured lakes represent ca. 32% (5.18 × 10⁸ m³) of the total lake volume (1.15 × 10⁹ m³) in the Cordillera Blanca. Approximately, 50% of the lakes have potential peak discharges > 1000 m³/s in case of lake outburst floods, implying a need for additional studies for risk assessment. © 2020 John Wiley & Sons, Ltd.     

Flash Code Simulations of Rayleigh-Taylor and Richtmyer-Meshkov Instabilities in Laser-Driven Experiments

We present two- and three-dimensional simulations involving Richtmyer–Meshkov and Rayleigh-Taylor instabilities run with the adaptive mesh refinement code, flash. Variations in the rate of mixing layer growth due to dimensionality, perturbation modes, and simulation resolution are explored. These simulations are designed for detailed comparisons with experiments run on the Omega laser to gain understanding of the mixing processes and to prepare for validation of the Flash code.     

Last glacial maximum equilibrium-line altitude and paleo-temperature reconstructions for the Cordillera de Mérida, Venezuelan Andes

The pattern and magnitude of glacier equilibrium-line altitude (ELA) lowerings in the tropics during the last glacial maximum (LGM) are topics of current debate. In the northern tropics, paleo-ELA data are particularly limited, inhibiting the ability to make regional and large-scale paleoclimatic inferences. To improve these records, nine paleo-glaciers in the Venezuelan Andes were reconstructed based on field observations, aerial photographs, satellite imagery and high-resolution digital topographic data. Paleo-glacier equilibrium-line altitudes (ELAs) were estimated using the accumulation-area ratio (AAR) and the area-altitude balance ratio (AABR) methods. During the local LGM in Venezuela (∼ 22,750 to 19,960 cal yr BP), ELAs were ∼ 850 to 1420 m lower than present. Local LGM temperatures were are at least 8.8 ± 2°C cooler than today based on a combined energy and mass-balance equation to account for an ELA lowering. This is greater than estimates using an atmospheric lapse rate calculation, which yields a value of 6.4 ± 1°C cooler. The paleo-glacial data from the Venezuelan Andes support other published records that indicate the northern tropics experienced a greater ELA lowering and possibly a greater cooling than the Southern Hemisphere tropics during the LGM.     

Basalt basement from the Kerguelen Plateau and the trail of a Dupal plume

The first samples of volcanic basement recovered from the Kerguelen Plateau are Lower Cretaceous transitional tholeiites. Isotope and incompatible element abundance ratios for these rocks are similar to ocean island basalts from the southern hemisphere Dupal anomaly region, and geochemical, geological and geophysical data are consistent with volcanic activity associated with a mantle plume. A reconstruction of plate motions suggests that the Kerguelen Plateau formed above a mantle plume in the interval 118-95 Ma, during the opening of the Indian Ocean between India and Australia-Antarctica. This plume was the source of other plateaus and ridges of the eastern Indian Ocean and possibly the Bunbury Basalt of southwestern Australia, and is now beneath Heard Island.     

Geochemical variations in Andean basaltic and silicic lavas from the Villarrica-Lanin volcanic chain (39.5° S): an evaluation of source heterogeneity,fractional crystallization and crustal assimilation

At 39.5° S in the southern volcanic zone of the Andes three Pleistocene-recent stratovolcanoes, Villarrica, Quetrupillan and Lanin, form a trend perpendicular to the strike of the Andes, 275 to 325 km from the Peru-Chile trench. Basalts from Villarrica and Lanin are geochemically distinct; the latter have higher incompatible element abundances and La/Sm but lower Ba/La and alkali metal/La ratios. These differences are consistent with our previously proposed models involving: a) a west to east decrease in an alkali metal-rich, high Ba/La slab-derived component which causes an across strike decrease in degree of melting; or b) a west to east increase in the contamination of subduction-related magma by enriched subcontinental lithospheric mantle. Silicic and mafic lavas from the stratovolcanoes have overlapping Sr, Nd and O isotopic ratios. Silicic lavas also have geochemical differences that parallel those of their associated basalts, e.g., rhyolite from Villarrica has lower La/Sm and incompatible element contents than high-SiO₂ andesite from Lanin. At each volcano the most silicic lavas can be modelled by closed system fractional crystallization while andesites are best explained by magma mixing. Apparently crustal contamination was not an important process in deriving the evolved lavas. Basaltic flows from small scoria cones, 20–35 km from Villarrica volcano have high incompatible element contents and low Ba/La, like Lanin basalts, but trend to higher K/Rb (356–855) and lower ⁸⁷Sr/ ⁸⁶Sr (0.70361–0.70400) than basalts from either stratovolcano. However all basalts have similar Nd, Pb and O isotope ratios. The best explanation for the unique features of the cones is that the sources of SVZ magmas, e.g., slab-derived fluids or melts of the subcontinental lithospheric mantle, have varying alkali metal and radiogenic Sr contents. These heterogeneities are not manifested in stratovolcano basalts because of extensive subcrustal pooling and mixing. This model is preferable to one involving crustal contamination because it can account for variable Sr isotope ratios and uniform Nd and Pb isotope ratios among the basalts, and the divergence of the cones from across-strike geochemical trends defined by the stratovolcanoes.     

The evolution of illite to muscovite: mineralogical and isotopic data from the Glarus Alps,Switzerland

Thirty-five illite and muscovite concentrates were extracted from Triassic and Permian claystones, shales, slates and phyllites along a cross-section from the diagenetic Alpine foreland (Tabular Jura and borehole samples beneath the Molasse Basin) to the anchi- and epimetamorphic Helvetic Zone of the Central Alps. Concentrates and thin sections were investigated by microscopic, X-ray, infrared, Mössbauer, thermal (DTA and TG), wet chemical, electron microprobe, K-Ar, Rb-Sr, ⁴⁰Ar/³⁹Ar and stable isotope methods.With increasing metamorphic grade based on illite crystallinity data (XRD and IR) the following continuous changes are observed: (i) the 1Md2M₁ polymorph transformation is completed in the higher grade anchizone; (ii) K₂O increases from 6–8 wt. % (diagenetic zone) to 8.5–10% (anchizone) to 10–11.5% (epizone), reflecting an increase in the total negative layer charge from 1.2 to 2.0; (iii) a decrease of the chemical variation of the mica population with detrital muscovite surviving up to the anchizone/ epizone boundary; iv) a shift of an endothermic peak in differential thermal curves from 500 to 750° C; (v) K-Ar and Rb-Sr apparent ages of the fraction <2 m decrease from the diagenetic zone to the epizone, K-Ar ages being generally lower than Rb-Sr ages. The critical temperature for total Ar resetting is estimated to be 260±30° C. K-Ar and Rb-Sr ages become concordant when the anchizone/ epizone boundary is approached. The stable isotope data, on the other hand, show no change with metamorphic grade but are dependent on stratigraphic age.These results suggest that the prograde evolution from 1 Md illite to 2M₁ muscovite involves a continuous lattice restructuration without rupture of the tetrahedral and octahedral bonds and change of the hydroxyl radicals, however this is not a recrystallization process. This restructuration is completed approximately at the anchizone/epizone boundary. The isotopic data indicate significant diffusive loss of ⁴⁰Ar and ⁸⁷Sr prior to any observable lattice reorganization. The restructuration progressively introduces a consistent repartition of Ar and K in the mineral lattices and is outlined by the ⁴⁰Ar/³⁹Ar age spectra.Concordant K-Ar and Rb-Sr ages of around 35-30 Ma. with concomitant concordant ⁴⁰Ar/³⁹Ar release spectra are representative for the main phase of Alpine metamorphism (Calanda phase) in the Glarus Alps. A second age group between 25 and 20 Ma. can probably be attributed to movements along the Glarus thrust (Ruchi phase), while values down to 9 Ma., in regions with higher metamorphic conditions, suggest thermal conditions persisting at least until the middle Tortonian.     

An anisotropy of magnetic susceptibility (AMS) fabric record of till kinematics within a Late Weichselian low Baltic till,southern Sweden

Herein we report on the results of an anisotropy of magnetic susceptibility (AMS) fabric case‐study of two Late Weichselian tills exposed in a bedrock quarry in Dalby, Skåne, southern Sweden. The region possesses a complex glacial history, reflecting alternating and interacting advances of the main body of the Scandinavian Ice Sheet (SIS) and its ice lobes from the Baltic basin, perhaps driven by streaming ice. AMS till fabrics are robust indicators of ice‐flow history and till kinematics, and provide a unique tool to investigate till kinematics within and amongst till units. The till section investigated here contains ~8 m of the Dalby Till – a dark grey silt‐clay rich till deposited during one or more Baltic advance – overlain by ~1.5 m of the regional surface diamicton. AMS fabrics within the lower part of the Dalby Till conform to the regional surface fluting, and reflect sustained flow from the ENE with progressive increases in basal strain. A boulder‐rich horizon approximately 3 m from the base of the till marks a restricted excursion in till fabric direction, fabric strength and style of strain. Ice flow is from the SW and W in the upper section. We interpret these fabrics to record shifting ice flow and bed conditions at the margins of the Young Baltic Advance ice lobe in southern Sweden, prior to a short‐lived re‐advance of the main body of the SIS over mainland Sweden recorded by the surface diamicton.     

Geochemical, mineralogical and isotopic signatures of the Semikan, West Kasrık “Turkish” phosphorites from the Derik–Mazıdağı–Mardin area, SE Anatolia

The geochemistry and mineralogy of a condensed section, ~6 m thick of the West Kasrik member (Coniacian–Santonian) near the uplifted northern flank of the Mardin–Derik anticline (south-eastern Turkey) was studied. The only deposit exploited in Turkey is found in this area. The sediment textures as well as the mineralogical and the geochemical results collectively suggest that these recycled phosphorites accumulated in areas of intensive very early diagenesis of the sediments in highly oxic bottom waters, and almost no detrital apport. The total concentration of redox-sensitive trace metals is very low (<600 μg/g); in addition the structural CO₂ and F/P₂O₅ are unusually high (~5 wt % and 0.14, respectively). REE distribution shows a clear “seawater” pattern with a strongly negative Ce-anomaly (0.20 ± 0.02) and heavy REE enrichment (Lu_N/La_N = 1.50 ± 0.12), however their total concentration is very low. The low REE contents, quite unexpected in recycled phosphorites, are explained by the scarcity of terrigenous components which leads to minimal incorporation of REEs from detrital clastic phases in the CFA fraction. Weathering previously suggested as responsible for the high P enrichment of the Mazidagi phosphorites is rather improbable in the samples we studied, in view of their high structural CO₂, high F/P₂O₅, high (La/Nd)_N and (La/Sm)_N ratios, and high Sr/P and Ca/P of separated CFA fractions that all negate post-depositional weathering. More likely, the high P enrichment of some rocks (P₂O₅ content reaches 34%) was produced by sedimentary and early diagenetic processes acting in oxygenated areas of starved sedimentation on and around tectonic highs.