Climate variability of southern Chile since the Last Glacial Maximum: a continuous sedimentological record from Lago Puyehue (40°S)

This paper presents a multi-proxy climate record of an 11 m long core collected in Lago Puyehue (southern Chile, 40°S) and extending back to 18,000 cal yr BP. The multi-proxy analyses include sedimentology, mineralogy, grain size, geochemistry, loss-on-ignition, magnetic susceptibility and radiocarbon dating. Results demonstrate that sediment grain size is positively correlated with the biogenic sediment content and can be used as a proxy for lake paleoproductivity. On the other hand, the magnetic susceptibility signal is correlated with the aluminium and titanium concentrations and can be used as a proxy for the terrigenous supply. Temporal variations of sediment composition evidence that, since the Last Glacial Maximum, the Chilean Lake District was characterized by three abrupt climate changes superimposed on a long-term climate evolution. These rapid climate changes are: (1) an abrupt warming at the end of the Last Glacial Maximum at 17,300 cal yr BP; (2) a 13,100–12,300 cal yr BP cold event, ending rapidly and interpreted as the local counterpart of the Younger Dryas cold period, and (3) a 3,400–2,900 cal yr BP climatic instability synchronous with a period of low solar activity. The timing of the 13,100–12,300 cold event is compared with similar records in both hemispheres and demonstrates that this southern hemisphere climate change precedes the northern hemisphere Younger Dryas cold period by 500 to 1,000 years. This is the third in a series of eight papers published in this special issue dedicated to the 17,900 year multi-proxy lacustrine record of Lago Puyehue, Chilean Lake District. The papers in this special issue were collected by M. De Batist, N. Fagel, M.-F. Loutre and E. Chapron.     

Seasonal and Interannual Variability in the Distribution of Surface Nutrients and Dissolved Inorganic Carbon in the Northern North Pacific: Influence of El Niño

The seasonal and interannual changes in surface nutrients, dissolved inorganic carbon (DIC) and total alkalinity (TA) were recorded in the North Pacific (30–54°N) from 1995 to 2001. This study focuses on the region north of the subarctic boundary (∼40°N) where there was extensive monthly coverage of surface properties. The nutrient cycles showed large interannual variations in the eastern and western subarctic gyres. In the Alaska Gyre the seasonal depletion of nitrate (ΔNO₃) increased from 8–14 μmol kg⁻¹ in 1995–1999 to 21.5 μmol kg⁻¹ in 2000. In the western subarctic the shifts were similar in amplitude but more frequent. The large ΔNO₃ levels were associated with high silicate depletions, indicating enhanced diatom production. The seasonal DIC:NO₃ drawdown ratios were elevated in the eastern and central subarctic due to calcification. In the western subarctic and the central Bering Sea calcification was significant only during 1997 and/or 1998, two El Ni?o years. Regional C/N stoichiometric molar ratios of 5.7 to 7.0 (>40°N) were determined based on the years with negligible or no calcification. The annual new production (NP_a) based on ΔNO₃ and these C/N ratios showed large interannual variations. NP_a was usually higher in the western than in the eastern subarctic. However, values of 84 gC m⁻²yr⁻¹ were found in the Alaska Gyre in 2000 which is similar to that in the most productive provinces of the northern North Pacific. There were also large increases in NP_a around the Alaska Peninsula in 1997 and 1998. Finally, the net removal of carbon by the biological pump was estimated as 0.72 Gt C yr⁻¹ in the North Pacific (>30°N). This revised version was published online in August 2006 with corrections to the Cover Date.     

Nutrient dynamics at the sediment-water interface in a Mediterranean lagoon (Thau, France): influence of biodeposition by shellfish farming activities

The Thau Lagoon, a French Mediterranean shallow lagoon, is a site where extensive shellfish farming occurs. The aim of the present work is to evaluate the role of this activity on nutrient exchange at the sediment-water interface in relation to organic matter (OM) sedimentation and degradation. Two stations inside (C5) and outside (C4) of the shellfish farming areas were sampled at three seasons. Porewater chemistry surveys and calculated diffusive fluxes were used to evaluate the trophic status of the Thau lagoon. Quantitative (Particulate Organic Carbon) as well as qualitative OM (Hydrogen Index, Carbohydrates) analyses were performed on sediments to assess OM characteristics. Results emphasized that surficial sediments at C5 are always more enriched in OM. Porewater nutrient concentrations are 10-20 times higher at C5 than at C4. In June 2003, the porewater profiles exhibit a sharp gradient at the bottom waters, indicating a hypereutrophic status, leading to an anoxic crisis.     

Multivariate Analysis of Incomplete Mapped Data

Classical multivariate analyses are based on matrix algebra and enable the analysis of a table containing measurements of a set of variables for a set of sites. Incomplete mapped data consist of measurements of a set of variables recorded for the same geographical region but for different zonal systems and with only a partial sampling of this zone. This kind of data cannot be analysed with usual multivariate methods because there is no common system of sites for all variables. We propose a new approach using GIS technology and NIPALS, an iterative multivariate method, to analyse the spatial patterns of this kind of data. Moreover, an extension of our method is that it can be used for areal interpolation purposes. We illustrate the method in analysing data concerning the distribution of roe deer weights over several years in a reserve.     

Performance of the LandSoil expert-based model to map erosion and sedimentation: application to a cultivated catchment in central Belgium

Intensive agricultural practices on sensitive soils induce high erosion rates in central Belgium. Expert-rules models quantify runoff and erosion at catchment scale, avoiding over-parameterization, and can include some direct or indirect connectivity features. The aim of this article is to test the ability of an expert-based model, LandSoil, to quantify runoff and to locate erosion and sedimentation areas in a small cultivated loamy catchment in Belgium during the years 2014, 2015 and 2016. Spatialized data are important for assessing model outputs and the erosive response. Measurements of runoff and observation of spatial erosion/deposition patterns, especially around major connectivity points, permitted an assessment of the reliability of the model results. Runoff modelling gave contrasting results (good linear adjustment at the outlet of the 83 ha sub-catchment (point 1): r² of 0.96, Nash–Sutcliffe criterion of 0.95; less good at the outlet of the 3.9 ha sub-catchment (point 2): r² of 0.28, Nash–Sutcliffe criterion of –0.47). For point 2 the poor results are explained by the very few runoff events observed, a scaling effect and the small area with a single land use. Graduated rulers demonstrate that the model is able to provide a coherent pattern of erosion/deposition. The study highlights great sensitivity to the effect of land use, land allocation, landscape design and slope gradients. Grass strips induce deposition of eroded particles when slopes are gentle (< 2%). Woodland strips decrease connectivity by being in the stream but deposit thinner sediment layers. Field boundaries have a role in the transport, but not really the quantity, of sediments. This model validation in the Belgian loess context allows us to use LandSoil in other similar environments in order to estimate the effects of landscape management scenarios. © 2020 John Wiley & Sons, Ltd.     

Discrimination between different water bodies from a multilayered aquifer (Kaluvelly watershed,India): Trace element signature

In the multilayered aquifer of Kaluvelly (India), comprising various sedimentary layers overlying a charnockitic basement, concentrations of trace elements were measured in aquifer formations and in groundwaters to identify geochemical tracers for water bodies. The two main sandstones (Cuddalore and Vanur) originate from the charnockites and the Cuddalore sandstone has experienced lateritization. In the studied area, two charnockite end-members were identified: a dioritic and a granitic one. Mineralogical composition and whole-rock Ti concentrations confirmed the charnockite which displayed the granitic composition as the parent rock of the two sandstones. Titanium distribution indicates that the Cuddalore sandstone originates from a more intense weathering of the parent material than the Vanur sandstone. Despite extensive differences in trace element contents recorded in aquifer formations, only a few trace elements were suitable to distinguish the water bodies. Among soluble elements, Li (in the Vanur aquifer) and Ba (in the charnockite and carbonaceous aquifers) can be used as tracers. As the input of these elements in solution is mainly regulated by the available stock, for a given mineralogical origin there is a direct link between the aquifer formation composition and water signature. With the exception of As, concentrations of redox-sensitive elements were not preserved during pumping because of oxidation, preventing their use as tracers. Low-mobility elements such as La, Ce, Th, Zr, Nb, Hf, or Ta were too insoluble to be detected in waters and/or to record the aquifer formation signature. Their input in solution was not regulated by the available stock but by the dissolution rate of rock-forming minerals. Only Ti can be used to distinguish between two out of the three aquifers (charnockite and Vanur). The specific behavior of Ti recorded in these waters may be linked to rutile inclusions within plagioclases and to the influence of climate on Ti solubility.     

Integrated stochastic modeling of pharmaceuticals in sewage networks

Pharmaceuticals and particularly antibiotics can harm sensitive aquatic species. Their occurrence in urban wastewater systems is the consequence of five successive processes: (i) ingestion of the substance, (ii) accumulation in the urine, (iii) excretion, (iv) degradation in the sewer system and (v) transport to the wastewater treatment plant (WTP). These processes were included in an integrated model that can be used to assess the dynamics of pharmaceuticals at a WTP inlet. First, information on sales data, posology, pharmacokinetics and toilet flushing frequency were combined to create a source model of pharmaceuticals entering a sewer system. This production function was then coupled with a transport/degradation model to simulate concentrations of pharmaceuticals at a WTP inlet. In an example application, the full model was applied to simulate the concentration of the antibiotic ciprofloxacin on an hourly time scale. In this application, the model was calibrated and validated for a case study at a WTP in Lausanne, Switzerland. Validation of the integrated approach was successful despite the high variability evident in the model results. This modeling approach has potential use in pollution management and epidemiology related to wastewater.     

Diatom assemblage changes in lacustrine sediments from Isla de los Estados, southernmost South America, in response to shifts in the southwesterly wind belt during the last deglaciation

Isla de los Estados (54° 45′S, 63° 10′–64° 46′W) lies east of the main island of Tierra del Fuego and is the southeastern-most point in Argentina. Because of its geographic position near the latitudes of the Southern Hemisphere Westerlies and the strong influence of the Antarctic Circumpolar Current (ACC), the area is suitable for paleoecological and paleoclimate research. The island is not far north of the Subantarctic Front, which limits the northern boundary of the ACC. Paleoenvironmental study in this geographic location can shed light on past changes in atmospheric and marine circulation patterns. Diatom analysis of the lower part of a sediment sequence from Laguna Cascada (54° 45′ 51.3′′S, 64° 20′ 20.07′′W) enabled inference of changing lake conditions between 16 and 11.1 cal ka BP. Between 16 and 14.4 cal ka BP fragilarioid diatom species, often a pioneer group, dominated the record. Their presence shows seasonally open-water conditions from the onset of sedimentation. In zone II (14.4–12.8 cal ka BP), the dominance of planktonic/tychoplanktonic Aulacoseira spp. might represent longer ice-free periods and windier conditions, which would have kept this heavy species suspended in the water column. This period corresponds to the Antarctic Cold Reversal, when the Southern Hemisphere Westerlies were possibly centered on the latitudes of Tierra del Fuego, resulting in windy and wet conditions. Zone III (12.8–11.1 cal ka BP) is dominated by benthic diatom taxa that are mainly associated with peat and wetland vegetation. This suggests that climate conditions had become milder and less windy, favoring aquatic productivity and terrestrial vegetation development. This change in environmental conditions may have been a consequence of the southward movement of the Southern Hemisphere Westerlies at the start of the Antarctic Holocene thermal optimum.     

Sensitivity of simulated Asian and African summer monsoons to orbitally induced variations in insolation 126, 115 and 6 kBP

We have conducted four numerical experiments with an atmospheric general circulation model (AGCM) to investigate the sensitivity of Asian and African monsoons to small changes (–5 to +12%), with respect to present-day, in incoming solar radiation at the top of the atmosphere. We show that, during the mid-Holocene (6 kBP where kBP means thousands of years before present-day) and the last interglacial (126 kBP), the Northern Hemisphere seasonal contrast was increased, with warmer summers and colder winters. At the time of glacial inception (115 kBP) however, summers were cooler and winters milder. As a consequence, Asia and tropical North Africa experienced stronger (weaker) summer monsoons 6 and 126 kBP (115 kBP), in agreement with previous numerical studies. This present study shows that summer warming/cooling of Eurasia and North Africa induced a shift of the main low-level convergence cell along a northwest/southeast transect. When land was warmer (during the summer months 6 and 126 kBP), the monsoon winds converged further inland bringing more moisture into northern India, western China and the southern Sahara. The southern tips of India, Indochina and southeastern China, as well as equatorial North Africa became drier. When land was cooler (during the summer 115 kBP), the main convergence zone was located over the west Pacific and the wet (dry) areas were those that were dry (wet) 6 and 126 kBP. The location and intensity of the simulated precipitation maxima were therefore very sensitive to changes in insolation. However the total amount of monsoon rain in Asia as well as in Africa remained remarkably stable through the time periods studied. These simulated migrations of convective activities were accompanied by changes in the nature of precipitation events: increased monsoon rains in these experiments were always associated with more high precipitation events (> 5 mm day ^–1), and fewer light showers (1 mm day^–). Rainy days with rates between 1 and 5 mm day^–1 were almost unchanged.     

Timing of the late glacial and Younger Dryas cold reversal in southern Chile varved sediments

The timing of the last deglaciation in southern Chile is re-evaluated from a calendar varve chronology (Lago Puyehue, 40° S). The climate shifts are analysed by continuous annual varve-thickness measurements through the ∼17,100 cal. year to 10,800 cal. year BP time window (∼3.5 m sediment core). The varve years are determined by the alternation of light (phytoplankton-rich) and dark (terrigenous and organic-rich) layers forming graded annual couplets (∼0.2 to 0.8 mm/year). The varve chronology is constructed by conventional varve-counting methods on thin sections after correction for instantaneous volcanic and/or seismic events detected in the thin sections. The calibrated varve-age model derived from the manual varve counting is constrained by high-resolution grey-scale (GS) semi-automatic counts of the annual light phytoplankton-rich layers (∼120 μm to 300 μm thick). Due to physical sediment properties the GS constitutes a proxy record for the phytoplankton/terrigenous varve-thickness variations through the sediment record. The varve couplets are thicker/thinner during humid/dry phases and darker/lighter (negative/positive annual grey-scale index) during cold/warm phases. Our results show that at 40° S the last deglaciation took place in two phases between ∼17,100 cal. year and ∼15,500 cal. year BP. We note a climate instability between ∼15,500 cal. year and 13,300 cal. year BP and a significant dry phase between ∼15,000 and 14,500 cal. year BP. We evidence a cold event in two phases between ∼13,300 and 12,200 cal. year BP interrupted by a dry event between ∼12,800 and 12,600 cal. year BP. The onset of a significant warmer period is observed after ∼11,500 cal. year BP. Our results provide new evidence of a Younger Dryas cool reversal in southern Chile, i.e., the Huelmo/Mascardi event Hajdas et al. (2003) associated with an abrupt dry pulse at ∼12,800–12,600 cal. year BP. The high-resolution grey-scale measurements performed on the biogenic varves from Lago Puyehue provide a reliable calibrated chronology of the regional environmental and climate shifts during the last deglaciation. This is eighth in a series of eight papers published in this special issue dedicated to the 17,900 year multi-proxy lacustrine record of Lago Puyehue, Chilean Lake District. The papers in this issue were collected by M. De Batist, N. Fagel, M.-F. Loutre and E. Chapron.