A chronology of the Little Ice Age in the tropical Andes of Bolivia (16°S) and its implications for climate reconstruction

Dating moraines by lichenometry enabled us to reconstruct glacier recession in the Bolivian Andes since the Little Ice Age maximum. On the 15 proglacial margins studied, we identified a system of ten principal moraines that marks the successive positions of glaciers over the last four centuries. Moraines were dated by performing statistical analysis of lichen measurements based on the extreme values theory. Like glaciers in many mid-latitude mountain areas, Bolivian glaciers reached their maximal extent during the second half of the 17th century. This glacier maximum coincides with the Maunder minimum of solar irradiance. By reconstructing the equilibrium-line altitude and changes in mass-balance, we think the glacier maximum may be due to a 20 to 30% increase in precipitation and a 1.1 to 1.2 °C decrease in temperature compared with present conditions. In the early 18th century, glaciers started to retreat at varying rates until the late 19th to early 20th century; this trend was generally associated with decreasing accumulation rates. By contrast, glacier recession in the 20th century was mainly the consequence of an increase in temperature and humidity. These results are consistent with observations made in the study region based on other proxies.     

Correction to: Using Causal Inference in Field Development Optimization: Application to Unconventional Plays

Mathematical Geosciences - Unfortunately, in the original version of the article the first and second name of the fourth author were wrong.     

Imprint of North-Atlantic abrupt climate changes on western European loess deposits as viewed in a dust emission model

Western European loess sequences of the last glaciation (100,000–15,000 years BP) exhibit strong, cyclic variations of the sedimentation rate, which are coeval to the Greenland stadial/interstadial cycles and the Heinrich events. These North-Atlantic rapid climate changes appear, thus, as a potential cause for the sedimentation variations, via changes in dust intensity cycle. Here we make a first step in testing this hypothesis, by modelling the impact of the North-Atlantic abrupt climate variations on dust emission. Our dust emission calculations use meteorological fields generated by the LMDZ atmospheric general circulation model at a resolution down to 60 km over Western Europe. Three numerical experiments are run, representing a Greenland stadial, an interstadial and a Heinrich event. Orbital parameters and ice-sheet configuration correspond to conditions from Marine Isotope Stage 3 (60,000–25,000 years BP), a period characterized by strong millennial-scale climate variability. The only differences we impose in the boundary conditions regard the North-Atlantic surface temperature and sea-ice cover in the latitudinal band 30°–63°N. The changes in wind, precipitation, soil moisture and snow cover from one simulated state to another result in small differences in dust emission intensity. In contrast, when the inhibition of the aeolian erosion by vegetation is taken into account, the dust fluxes for the cold climate states (Greenland stadial and Heinrich event) become generally more than twice higher than those for the relatively warmer Greenland interstadial, in agreement with the loess data. These results support the hypothesis that the North-Atlantic millennial-scale variability is imprinted in Western European loess profiles, and point to vegetation changes as the main factor responsible for millennial-scale sedimentation variations. An analysis for the English Channel and southern North Sea areas, major potential dust sources, shows that the seasonality of dust emission is not controlled by the wind speed, as in modern large deserts, but by the surface conditions. Consequently, the dusty season lasts from late winter to early summer, with maximum activity in April–May, and is shifted towards summer when the climate is colder.     

Mineralogy, geochemistry and petrogenesis of the recent magmatic formations from Mbengwi, a continental sector of the Cameroon Volcanic Line (CVL), Central Africa

The Mbengwi recent magmatic formations consist of volcanics and syenites belonging to the same magmatic episode. Lavas form a bimodal basanite-rhyolite alkaline series with a gap between 50 and 62?wt.% SiO₂. Mafic lavas (basanite-hawaiite) are sodic while felsic rocks (trachyte-rhyolite-syenites) are sodi-potassic, slightly metaluminous to peralkaline. The geochemical and isotopic characteristics (0.7031?<?(⁸⁷Sr/⁸⁶Sr)_initial?<?0.7043; 1.03?<?ε_Ndi?<?5.17) of these rocks are similar to those of other rocks from the CVL. The main differentiation process is fractional crystallization with two trends of fractionation. Their Rb/Sr isochron age of 28.2?Ma, almost similar to 27.40?±?0.6?Ma?K/Ar age obtained in a trachyte from neighboring Bamenda Mountains system, precludes any local age migration of an hypothetic hotspot. Mafic lavas have OIB features displaying an isotopic signature similar to that of HIMU mantle source different from FOZO known as source of most parental magmas along the CVL.     

Internally consistent geothermometers for garnet peridotites and pyroxenites

Mutual relationships among temperatures estimated with the most widely used geothermometers for garnet peridotites and pyroxenites demonstrate that the methods are not internally consistent and may diverge by over 200°C even in well-equilibrated mantle xenoliths. The Taylor (N Jb Min Abh 172:381–408, 1998) two-pyroxene (TA98) and the Nimis and Taylor (Contrib Mineral Petrol 139:541–554, 2000) single-clinopyroxene thermometers are shown to provide the most reliable estimates, as they reproduce the temperatures of experiments in a variety of simple and natural peridotitic systems. Discrepancies between these two thermometers are negligible in applications to a wide variety of natural samples (≤30°C). The Brey and Köhler (J Petrol 31:1353–1378, 1990) Ca-in-Opx thermometer shows good agreement with TA98 in the range 1,000–1,400°C and a positive bias at lower T (up to +90°C, on average, at T _TA98 = 700°C). The popular Brey and Köhler (J Petrol 31:1353–1378, 1990) two-pyroxene thermometer performs well on clinopyroxene with Na contents of ~0.05 atoms per 6-oxygen formula, but shows a systematic positive bias with increasing Na_Cpx (+150°C at Na_Cpx = 0.25). Among Fe–Mg exchange thermometers, the Harley (Contrib Mineral Petrol 86:359–373, 1984) orthopyroxene–garnet and the recent Wu and Zhao (J Metamorphic Geol 25:497–505, 2007) olivine–garnet formulations show the highest precision, but systematically diverge (up to ca. 150°C, on average) from TA98 estimates at T far from 1,100°C and at T < 1,200°C, respectively; these systematic errors are also evident by comparison with experimental data for natural peridotite systems. The older O’Neill and Wood (Contrib Mineral Petrol 70:59–70, 1979) version of the olivine–garnet Fe–Mg thermometer and all popular versions of the clinopyroxene–garnet Fe–Mg thermometer show unacceptably low precision, with discrepancies exceeding 200°C when compared to TA98 results for well-equilibrated xenoliths. Empirical correction to the Brey and Köhler (J Petrol 31:1353–1378, 1990) Ca-in-Opx thermometer and recalibration of the orthopyroxene–garnet thermometer, using well-equilibrated mantle xenoliths and TA98 temperatures as calibrants, are provided in this study to ensure consistency with TA98 estimates in the range 700–1,400°C. Observed discrepancies between the new orthopyroxene–garnet thermometer and TA98 for some localities can be interpreted in the light of orthopyroxene–garnet Fe³⁺ partitioning systematics and suggest localized and lateral variations in mantle redox conditions, in broad agreement with existing oxybarometric data. Kinetic decoupling of Ca–Mg and Fe–Mg exchange equilibria caused by transient heating appears to be common, but not ubiquitous, near the base of the lithosphere.     

Three-dimensional computer-assisted basin modelling to generate exploration target areas: an example from the late Archaean-early Proterozoic Transvaal Supergroup, South Africa

Second- and third-order fault-bounded Precambrian basins frequently host deposits of the sedimentary massive sulphide group. Three-dimensional geometric modelling of the thickness of preserved basin-fill successions of the Transvaal Supergroup, using DATAMINE software, and residual gravity modelling of the contemporary basement floor, help delineate areas of exploration potential in this unit. Two main depositional axes are tentatively identified for the basal volcano-sedimentary protobasinal Transvaal successions. A sheet-like geometry was indicated for the succeeding Black Reef sandstones and Chuniespoort Group chemical sedimentary rocks. The uppermost Pretoria Group thickness model delineates eastern and western second-order basins separated by a central submerged palaeohigh. A similar isopach pattern is noted for the thick shales of the Silverton Formation in this group, with, in addition, a well-defined third-order basin in the northwest of the western second-order basin. The residual gravity model indicates two linear palaeovalleys adjacent to this western basin, one coincident with one of the axes inferred for the protobasinal rocks. The fault-bounded second- and third-order basins and depositional axes postulated here are consistent with known geological data and suggested sedimentation models. Cumulative distortions implicit in the DATAMINE computer modelling technique are reduced when the method is applied on the basin-wide scale, enabling identification of regional exploration target areas rather than immediate prospecting targets. Received: 14 August 1996 / Accepted: 13 March 1997     

Climatic variations in MIS 11 recorded by stable isotopes and trace elements in a French tufa (La Celle,Seine Valley)

Marine Isotope Stage (MIS) 11 palaeoclimate has so far been documented in marine and ice sheet isotopic records. However, excepting some lacustrine pollen records, very little is known about palaeoclimatic conditions in continental areas. This study uses geochemical records in calcareous tufa deposits from rivers as a basis for reconstructing temperate palaeoclimatic conditions. Tufa deposits are now proven to record high‐quality palaeoclimatic information in recent to Holocene deposits. Work on older interglacial tufas is just starting and in this paper we present the first comprehensive results from a MIS 11 tufa. The tufa comes from the Seine Valley (La Celle, northern France). Geochemical data in the tufa calcite are interpreted to record primarily air temperature (δ¹⁸O) and humidity (δ¹³C and Mg/Ca and Sr/Ca). The combined data identify a warm and wet climatic optimum followed by two temperature decreases associated with oscillations in humidity. These marked climatic variations recorded through the La Celle profile are strongly coherent with the palaeoenvironmental reconstructions from malacological data. The abrupt climatic and environmental events recorded could be related to short‐term degradation of vegetation cover in Europe, which is itself controlled by global palaeoclimatic events. Copyright © 2012 John Wiley & Sons, Ltd.     

Fluid flow through the sedimentary cover in northern Switzerland recorded by calcite–celestite veins (Oftringen borehole,Olten)

Abundant veins filled by calcite, celestite and pyrite were found in the core of a 719 m deep borehole drilled in Oftringen near Olten, located in the north-western Molasse basin, close to the thrust of the Folded Jura. Host rocks are calcareous marl, argillaceous limestone and limestone of the Dogger and Malm. The δ¹⁸O values of vein calcite are lower than in host rock carbonate and, together with microthermometric data from fluid inclusions in vein calcite, indicate precipitation from a seawater-dominated fluid at average temperatures of 56–68°C. Such temperatures were reached at the time of maximum burial of the sedimentary pile in the late Miocene. The depth profile of δ¹³C and ⁸⁷Sr/⁸⁶Sr values and Sr content of both whole-rock carbonate and vein calcite show marked trends towards negative δ¹³C, high ⁸⁷Sr/⁸⁶Sr, and low Sr content in the uppermost 50–150 m of the Jurassic profile (upper Oxfordian). The ⁸⁷Sr/⁸⁶Sr of vein minerals is generally higher than that of host rock carbonate, up to very high values corresponding to Burdigalian seawater (Upper Marine Molasse, Miocene), which represents the last marine incursion in the region. No evidence for internally derived radiogenic Sr (clay minerals) has been found and so an external source is required. S and O isotope composition of vein celestite and pyrite can be explained by bacterial reduction of Miocene seawater sulphate. The available data set suggests the vein mineralization precipitated from descending Burdigalian seawater and not from a fluid originating in the underlying Triassic evaporites.     

New insights into Late Pleistocene glacial and postglacial history of northernmost Ungava (Canada) from Pingualuit Crater Lake sediments

The Pingualuit Crater was formed by a meteoritic impact ca. 1.4 million years ago in northernmost Ungava (Canada). Due to its geographical position near the center of successive North American ice sheets and its favorable morphometry, the Pingualuit Crater Lake (water depth = 246 m) promises to yield a unique continuous sedimentary sequence covering several glacial/interglacial cycles in the terrestrial Canadian Arctic. In this paper, we suggest the existence of a subglacial lake at least during the Last Glacial Maximum (LGM) by hydraulic potential modeling using LGM ice-surface elevation and bed topography derived from a digital elevation model. These results support the hypothesis that the bottom sediments of the Crater Lake escaped glacial erosion and may contain a long-term continental sedimentary sequence. We also present the stratigraphy of a 9 m-long core retrieved from the deep basin of the lake as well as a multiproxy reconstruction of its deglacial and postglacial history. The base of the core is formed by very dense diamicton reflecting basal melt-out environments marking the end of subglacial conditions at the coring site. The overlying finely laminated silt are related to the onset of proglacial conditions characterized by extremely low lacustrine productivity. Infra Red Stimulated Luminescence and AMS ¹⁴C dating, as well as biostratigraphic data indicate sediment mixing between recent (e.g. Holocene) and much older (pre- to mid-Wisconsinan) material reworked by glacier activity. This process prevents the precise dating of these sediments that we interpret as being deposited just before the final deglaciation of the lake. Two finer grained and organic-rich intervals reflect the inception of lacustrine productivity resulting from the cessation of glacial meltwater inputs and ice-free periods. The lower organic interval corresponds to the early postglacial period (6850–5750 cal BP) and marks the transition between proglacial and postglacial conditions during the Holocene Thermal Maximum, while the uppermost organic-rich core section represents late Holocene sediments (～4200–600 cal BP). The organic intervals are separated by a basin-scale erosive slide occurring around 4200 cal BP and likely related to 1) a seismic event due to the glacio-isostatic rebound following the last deglaciation or 2) slope instabilities associated with rapid discharge events of the lake.     

Micromechanical aspects of the shear strength of wet granular soils

This paper presents a micromechanical model for the analysis of wet granular soils at low saturation (below 30%). The discrete element method is employed to model the solid particles. The capillary water is assumed to be in a pendular state and thus exists in the form of liquid bridges at the particle‐to‐particle contacts. The resulting inter‐particle adhesion is accounted for using the toroidal approximation of the bridge. Hydraulic hysteresis is accounted for based on the possible mechanism of the formation and breakage of the liquid bridges during wetting and drying phases. Shear test computational simulations were conducted at different water contents under relatively low net normal stresses. The results of these simulations suggest that capillary‐induced attractive forces and hydraulic hysteresis play an important role in affecting the shear strength of the soil. These attractive forces produce a tensile stress that contributes to the apparent cohesion of the soil and increases its stiffness. During a drying phase, capillary‐induced tensile stresses, and hence shear strength, tend to be larger than those during a wetting phase. The proposed model appears to capture the macroscopic response of wet granular materials and revealed a number of salient micromechanical mechanisms and response patterns consistent with theoretical considerations. Copyright © 2008 John Wiley & Sons, Ltd.