Dynamics of the cold surface layer of polythermal Storglaciären, Sweden

Polythermal glaciers, i.e. glaciers with a combination of ice at and below the freezing point, are widespread in arctic and subarctic environments. The polythermal structure has major implications for glacier hydrology, ice flow and glacial erosion. However, the interplay of factors governing its spatial and temporal variations such as net mass balance, ice advection and water content in the ice is poorly investigated and as yet not fully understood. This study deals with a thorough investigation of the polythermal regime on Storglaciären, northern Sweden, a small valley glacier with a cold surface layer in the ablation area. Extensive field work was performed including mapping of the cold surface layer using ground-penetrating radar, ice temperature measurements, mass balance and ice velocity measurements. Analyses of these data combined with numerical modelling were used specifically to investigate the spatial and temporal variability of the cold surface layer, the spatial distribution of the water content just below the cold surface layer transition, the effect of radar frequency on the detection of the surface layer, and the sensitivity of the cold surface layer to changes in forcing. A comparison between direct temperature measurements in boreholes and ground-penetrating surveys shows that the radar-inferred cold-temperate transition depth is within ±1 m from the melting point of ice at frequencies above ~300 MHz. At frequencies below ~155 MHz, the accuracy degrades because of reduced scattering efficiency that occurs when the scatterers become much smaller compared to the wavelength. The mapped spatial pattern of the englacial cold-temperate transition boundary is complex. This pattern reflects the observed spatial variation in net loss of ice at the surface by ablation and vertical advection of ice, which is suggested to provide the predominant forcing of the cold surface layer thickness pattern. This is further supported by thermomechanical modeling of the cold surface layer, which indicates high sensitivity of the cold surface layer thickness to changes in vertical advection rates. The water content is the least investigated quantity that is relevant for the thermal regime of glaciers, but also the most difficult to assess. Spatial variability of absolute water content in the temperate ice immediately below the cold surface layer on Storglaciären was determined by combining relative estimates of water content from ground-penetrating radar data with absolute determination from temperature measurements and the thermal boundary condition at the freezing front. These measurements indicate large-scale spatial variability in the water content, which seems to arise from variations in entrapment of water at the firn-ice transition. However, this variability cannot alone explain the spatial pattern in the thermal regime on Storglaciären. Repeated surveys of the cold surface layer show a 22% average thinning of the cold surface layer on Storglaciären between 1989 and 2001. Transient thermomechanical modeling results suggest that the cold surface layer adapts to new equilibrium conditions in only a few decades after a perturbation in the forcing is introduced. An increased winter air temperature since mid-1980s seems to be the cause of the observed thinning of the cold surface layer. Over the last decades, mass balance measurements indicate that the glacier has been close to a steady state. The quasi-steady state situation is also reflected in the vertical advection, which shows no significant changes during the last decades. Increased winter temperatures at the ice surface would result in a slow-down of the formation of cold ice at the base of the cold surface layer and lead to a larger imbalance between net loss of ice at the surface and freezing of temperate ice at the cold-temperate transition.

Stability, viscous dissipation and local thermal non-equilibrium in fluid saturated porous media

In this thesis, the field of study related to the stability analysis of fluid saturated porous media is investigated. In particular the contribution of the viscous heating to the onset of convective instability in the flow through ducts is analysed. In order to evaluate the contribution of the viscous dissipation, different geometries, different models describing the balance equations and different boundary conditions are used. Moreover, the local thermal non-equilibrium model is used to study the evolution of the temperature differences between the fluid and the solid matrix in a thermal boundary layer problem. On studying the onset of instability, different techniques for eigenvalue problems has been used. Analytical solutions, asymptotic analyses and numerical solutions by means of original and commercial codes are carried out.

Granulite facies metamorphism and partial melting processes in the Ivrea Zone, Northern Italy

The Ivrea Zone in northern Italy has been the focus of numerous petrological, geochemical and structural studies. It is commonly inferred to represent an almost complete section through the mid to lower continental crust, in which metamorphism and partial melting of the abundant metapelites was the result of magmatic underplating by a large volume of mantle-derived magma. This study concerns amphibolite and granulite facies metamorphism in the Ivrea Zone with focus on metapelites and metapsammites/metagreywackes from Val Strona di Omegna and metapelites from Val Sesia and Val Strona di Postua, with the aim to better constrain their metamorphic evolution as well as their pressure and temperature conditions via phase equilibria modelling.rnrnIn Val Strona di Omegna, the metapelites show a structural and mineralogical change from mica-schists with the common assemblage bi-mu-sill-pl-q-ilm ± liq at the lowest grades, through metatexitic migmatites (g-sill-bi-ksp-pl-q-ilm-liq) at intermediate grades, to complex diatexitic migmatites (g-sill-ru-bi-ksp-pl-q-ilm-liq) at the highest grades. Within this section several mappable isograds occur, including the first appearance of K-feldspar in the metapelites, the first appearance of orthopyroxene in the metabasites and the disappearance of prograde biotite from the metapelites. The inferred onset of partial melting in the metapelites occurs around Massiola. The prograde suprasolidus evolution of the metapelites is consistent with melting via the breakdown of first muscovite then biotite. Maximum modelled melt fractions of 30–40 % are predicted at the highest grade. The regional metamorphic field gradient in Val Strona di Omegna is constrained to range from conditions of 3.5–6.5 kbar at T = 650–730 °C to P > 9 kbar at T > 900 °C. The peak P–T estimates, particularly for granulite facies conditions, are significantly higher (around 100 °C) than those of most previous studies. In Val Sesia and Val Strona di Postua to the south the exposure is more restricted. P–T estimates for the metapelites are 750–850 °C and 5–6.5 kbar in Val Sesia and approximately 800–900 °C and 5.5–7 kbar in Val Strona di Postua. These results show similar temperatures but lower pressure than metapelites in Val Strona di Omegna. Metapelites in Val Sesia in contact with the Mafic Complex exhibit a metatexitic structure, while in Val Strona di Postua diatexitic structures occur. Further, metapelites at the contact with the Mafic Complex contain cordierite (± spinel) that overprint the regional metamorphic assemblages and are interpreted to have formed during contact metamorphism related to intrusion of the Mafic Complex. The lower pressures in the high-grade rocks in Val Sesia and Val Strona di Postua are consistent with some decompression from the regional metamorphic peak prior to the intrusion of the Mafic Complex, suggesting the rocks followed a clockwise P–T path. In contrast, the metapelites in Val Strona di Omegna, especially in the granulite facies, do not contain any cordierite or any evidence for a contact metamorphic overprint. The extrapolated granulite facies mineral isograds are cut by the rocks of the Mafic Complex to the south. Therefore, the Mafic Complex cannot have caused the regional metamorphism and it is unlikely that the Mafic Complex occurs in Val Strona di Omegna.

Biomedical and industrial applications of atmospheric pressure non-equilibrium plasmas

This dissertation will be focused on the characterization of an atmospheric pressure plasma jet source with an application oriented diagnostic approach and the description of processes supported by this plasma source. The plasma source investigated is a single electrode plasma jet. Schlieren images, optical emission spectra, temperature and heat flux profiles are analyzed to deeply investigate the fluid dynamic, the chemical composition and the thermal output of the plasma generated with a nanosecond-pulsed high voltage generator. The maximum temperature measured is about 45 °C and values close to the room temperature are reached 10 mm down the source outlet, ensuring the possibility to use the plasma jet for the treatment of thermosensitive materials, such as, for example, biological substrate or polymers. Electrospinning of polymeric solution allows the production of nanofibrous non-woven mats and the plasma pre-treatment of the solutions leads to the realization of defect free nanofibers. The use of the plasma jet allows the electrospinnability of a non-spinnable poly(L-lactic acid) (PLLA) solution, suitable for the production of biological scaffold for the wound dressing.

Modelling, diagnostics and experimental analysis of plasma assisted processes for material treatment

This work presents results from experimental investigations of several different atmospheric pressure plasmas applications, such as Metal Inert Gas (MIG) welding and Plasma Arc Cutting (PAC) and Welding (PAW) sources, as well as Inductively Coupled Plasma (ICP) torches. The main diagnostic tool that has been used is High Speed Imaging (HSI), often assisted by Schlieren imaging to analyse non-visible phenomena. Furthermore, starting from thermo-fluid-dynamic models developed by the University of Bologna group, such plasma processes have been studied also with new advanced models, focusing for instance on the interaction between a melting metal wire and a plasma, or considering non-equilibrium phenomena for diagnostics of plasma arcs. Additionally, the experimental diagnostic tools that have been developed for industrial thermal plasmas have been used also for the characterization of innovative low temperature atmospheric pressure non equilibrium plasmas, such as dielectric barrier discharges (DBD) and Plasma Jets. These sources are controlled by few kV voltage pulses with pulse rise time of few nanoseconds to avoid the formation of a plasma arc, with interesting applications in surface functionalization of thermosensitive materials. In order to investigate also bio-medical applications of thermal plasma, a self-developed quenching device has been connected to an ICP torch. Such device has allowed inactivation of several kinds of bacteria spread on petri dishes, by keeping the substrate temperature lower than 40 degrees, which is a strict requirement in order to allow the treatment of living tissues.

Detecting and resolving position-dependent temperature effects in real-time quantitative polymerase chain reaction

Real-time quantitative polymerase chain reaction (qPCR) depends on precise temperature control of the sample during cycling. In the current study, we investigated how temperature variation in plate-based qPCR instruments influences qPCR results. Temperature variation was measured by amplicon melting analysis as a convenient means to assess well-to-well differences. Multiple technical replicates of several SYBR Green I-based qPCR assays allowed correlation of relative well temperature to quantification cycle. We found that inadequate template denaturation results in an inverse correlation and requires increasing the denaturation temperature, adding a DNA destabilizing agent, or pretreating with a restriction enzyme. In contrast, inadequate primer annealing results in a direct correlation and requires lowering the annealing temperature. Significant correlations were found in 18 of 25 assays. The critical nature of temperature-dependent effects was shown in a blinded study of 29 patients for the diagnosis of Prader-Willy and Angelman syndromes, where eight diagnoses were incorrect unless temperature-dependent effects were controlled. A method to detect temperature-dependent effects by pairwise comparisons of replicates in routine experiments is presented and applied. Systematic temperature errors in qPCR instruments can be recognized and their effects eliminated when high precision is required in quantitative genetic diagnostics and critical complementary DNA analyses.

Melting of Northern Greenland during the last interglaciation

Using simulated climate data from the comprehensive coupled climate model IPSL CM4, we simulate the Greenland ice sheet (GrIS) during the Eemian interglaciation with the three-dimensional ice sheet model SICOPOLIS. The Eemian is a period 126 000 yr before present (126 ka) with Arctic temperatures comparable to projections for the end of this century. In our simulation, the northeastern part of the GrIS is unstable and retreats significantly, despite moderate melt rates. This result is found to be robust to perturbations within a wide parameter space of key parameters of the ice sheet model, the choice of initial ice temperature, and has been reproduced with climate forcing from a second coupled climate model, the CCSM3. It is shown that the northeast GrIS is the most vulnerable. Even a small increase in melt removes many years of ice accumulation, giving a large mass imbalance and triggering the strong ice-elevation feedback. Unlike the south and west, melting in the northeast is not compensated by high accumulation. The analogy with modern warming suggests that in coming decades, positive feedbacks could increase the rate of mass loss of the northeastern GrIS, exceeding the recent observed thinning rates in the south.

Three dimensional finite element simulation of polymer melting and flow in a single-screw extruder : optimization of screw channel geometry

Single-screw extrusion is one of the widely used processing methods in plastics industry, which was the third largest manufacturing industry in the United States in 2007 [5]. In order to optimize the single-screw extrusion process, tremendous efforts have been devoted for development of accurate models in the last fifty years, especially for polymer melting in screw extruders. This has led to a good qualitative understanding of the melting process; however, quantitative predictions of melting from various models often have a large error in comparison to the experimental data. Thus, even nowadays, process parameters and the geometry of the extruder channel for the single-screw extrusion are determined by trial and error. Since new polymers are developed frequently, finding the optimum parameters to extrude these polymers by trial and error is costly and time consuming. In order to reduce the time and experimental work required for optimizing the process parameters and the geometry of the extruder channel for a given polymer, the main goal of this research was to perform a coordinated experimental and numerical investigation of melting in screw extrusion. In this work, a full three-dimensional finite element simulation of the two-phase flow in the melting and metering zones of a single-screw extruder was performed by solving the conservation equations for mass, momentum, and energy. The only attempt for such a three-dimensional simulation of melting in screw extruder was more than twenty years back. However, that work had only a limited success because of the capability of computers and mathematical algorithms available at that time. The dramatic improvement of computational power and mathematical knowledge now make it possible to run full 3-D simulations of two-phase flow in single-screw extruders on a desktop PC. In order to verify the numerical predictions from the full 3-D simulations of two-phase flow in single-screw extruders, a detailed experimental study was performed. This experimental study included Maddock screw-freezing experiments, Screw Simulator experiments and material characterization experiments. Maddock screw-freezing experiments were performed in order to visualize the melting profile along the single-screw extruder channel with different screw geometry configurations. These melting profiles were compared with the simulation results. Screw Simulator experiments were performed to collect the shear stress and melting flux data for various polymers. Cone and plate viscometer experiments were performed to obtain the shear viscosity data which is needed in the simulations. An optimization code was developed to optimize two screw geometry parameters, namely, screw lead (pitch) and depth in the metering section of a single-screw extruder, such that the output rate of the extruder was maximized without exceeding the maximum temperature value specified at the exit of the extruder. This optimization code used a mesh partitioning technique in order to obtain the flow domain. The simulations in this flow domain was performed using the code developed to simulate the two-phase flow in single-screw extruders.

Computational simulations of latent heat thermal energy storage systems - with innovative and first-principles based simulation for the underlying unsteady melting (and solidification) processes

This thesis develops an effective modeling and simulation procedure for a specific thermal energy storage system commonly used and recommended for various applications (such as an auxiliary energy storage system for solar heating based Rankine cycle power plant). This thermal energy storage system transfers heat from a hot fluid (termed as heat transfer fluid - HTF) flowing in a tube to the surrounding phase change material (PCM). Through unsteady melting or freezing process, the PCM absorbs or releases thermal energy in the form of latent heat. Both scientific and engineering information is obtained by the proposed first-principle based modeling and simulation procedure. On the scientific side, the approach accurately tracks the moving melt-front (modeled as a sharp liquid-solid interface) and provides all necessary information about the time-varying heat-flow rates, temperature profiles, stored thermal energy, etc. On the engineering side, the proposed approach is unique in its ability to accurately solve – both individually and collectively – all the conjugate unsteady heat transfer problems for each of the components of the thermal storage system. This yields critical system level information on the various time-varying effectiveness and efficiency parameters for the thermal storage system.

High-resolution summer temperature reconstruction from Lake Silvaplana based on in-situ reflectance spectroscopy

Annually laminated (varved) sediments of proglacial Lake Silvaplana (46 ̊27’N, 9 ̊48’E, 1791 m a.s.l., Engadine, eastern Swiss Alps) provide an excellent archive for quantitative high-resolution (seasonal – annual) reconstruction of high- and lowfrequency climate signals back to AD 1580. The chronology of the core is based on varve counting, Cs-137, Pb-210 and event stratigraphy. In this study we present a reconstruction based on in-situ reflectance spectroscopy. In situ reflectance spectroscopy is known as a cost- and time-effective non destructtive method for semi-quantitative analysis of pigments (e.g., chlorines and carotenoids) and of lithoclastic sediment fractions. Reflectance-dependent absorption (RDA) was measured with a Gretac Macbeth spectrolino at 2 mm resolution. The spectral coverage ranges from 380 nm to 730 nm at 10 nm band resolution. In proglacial Lake Silvaplana, 99% of the sediment is lithoclastic prior to AD 1950. Therefore, we concentrate on absorption features that are characteristic for lithoclastic sediment fractions. In Lake Silvaplana, two significant correlations that are stable in time were found between RDA typical for lithoclastics and meteorological data: (1) the time series R 570 /R 630 (ratio between RDA at 570 nm and 630 nm) of varves in Lake Silvaplana and May to October temperatures at nearby station of Sils correlate highly significantly (calibration period AD 1864 – 1951, r = 0.74, p < 0.01 for 5ptsmoothed series; RMSE is 0.28 ̊C, RE = 0.41 and CE = 0.38), and (2) the minimum reflectance within the 690nm band (min690) data correlate with May to October (calibration period AD 1864 – 1951, r = 0.68, p < 0.01 for 5pt-smoothed series; RMSE = 0.22 ̊C, RE = 0.5, CE = 0.31). Both proxy series (min690nm and R 570 /R 630 values) are internally highly consistent (r = 0.8, p < 0.001). In proglacial Lake Silvaplana the largest amount of sediment is transported by glacial meltwater. The melting season spans approximately from May to October, which gives us a good understanding of the geophysical processes explaining the correlations between lithoclastic proxies and the meteorological data. The reconstructions were extended back to AD 1580 and show a broad corresponddence with fully independent reconstructions from tree rings and documentary data.

Metamorphic and geochronogical study of the Triassic El Oro metamorphic complex, Ecuador: Implications for high-temperature metamorphism in a forearc zone

In the forearc of the Andean active margin in southwest Ecuador, the El Oro metamorphic complex exhibits a well exposed tilted forearc section partially migmatized. We used Raman spectroscopy on carbonaceous matter (RSCM) thermometry and pseudosections coupled with mineralogical and textural studies to constrain the pressure–temperature (P–T) evolution of the El Oro metamorphic complex during Triassic times. Our results show that anatexis of the continental crust occurred by white-mica and biotite dehydration melting along a 10 km thick crustal domain (from 4.5 to 8 kbar) with increasing temperature from 650 to 700 °C. In the biotite dehydration melting zone, temperature was buffered at 750–820 °C in a 5 km thick layer. The estimated average thermal gradient during peak metamorphism is of 30 °C/km within the migmatitic domain can be partitioned into two apparent gradients parts. The upper part from surface to 7 km depth records a 40–45 °C/km gradient. The lower part records a quasi-adiabatic geotherm with a 10 °C/km gradient consistent with an isothermal melting zone. Migmatites U–Th–Pb geochronology yielded zircon and monazite ages of 229.3 ± 2.1 Ma and 224.5 ± 2.3 Ma, respectively. This thermal event generated S-type magmatism (the Marcabeli granitoid) and was immediately followed by underplating of the high-pressure low-temperature (HP-LT) Arenillas–Panupalí unit at 225.8 ± 1.8 Ma. The association of high-temperature low-pressure (HT-LP) migmatites with HP-LT unit constitutes a new example of a paired metamorphic belt along the South American margin. We propose that in addition to crustal thinning, underplating of the Piedras gabbroic unit before 230 Ma provided the heat source necessary to foster crustal anatexis. Furthermore, its MORB signature shows that the asthenosphere was involved as the source of the heat anomaly. S-type felsic magmatism is widespread during this time and suggests that a large-scale thermal anomaly affected a large part of the South American margin during the late Triassic. We propose that crustal anatexis is related to an anomaly that arose during subduction of the Panthalassa ocean under the South American margin. Slab verticalization or slab break-off can be invoked as the origin of the upwelling of the asthenosphere.

Experimental Determination of Salinity, Temperature, Growth, and Metabolic Effects on Shell Isotope Chemistry of Mytilus edulis Collected from Maine and Greenland

To study the effects of temperature, salinity, and life processes (growth rates, size, metabolic effects, and physiological/ genetic effects) on newly precipitated bivalve carbonate, we quantified shell isotopic chemistry of adult and juvenile animals of the intertidal bivalve Mytilus edulis (Blue mussel) collected alive from western Greenland and the central Gulf of Maine and cultured them under controlled conditions. Data for juvenile and adult M. edulis bivalves cultured in this study, and previously by Wanamaker et al. (2006), yielded statistically identical paleotemperature relationships. On the basis of these experiments we have developed a species-specific paleotemperature equation for the bivalve M. edulis [T degrees C = 16.28 (+/- 0.10) -4.57 (+/- 0.15) {delta(18)O(c) VPBD - delta(18)O(w) VSMOW} + 0.06 (+/- 0.06) {delta(18)O(c) VPBD - delta(18)O(w) VSMOW}(2); r(2) = 0.99; N = 323; p < 0.0001]. Compared to the Kim and O'Neil (1997) inorganic calcite equation, M. edulis deposits its shell in isotope equilibrium (delta(18)O(calcite)) with ambient water. Carbon isotopes (delta(13)C(calcite)) from sampled shells were substantially more negative than predicted values, indicating an uptake of metabolic carbon into shell carbonate, and delta(13)C(calcite) disequilibrium increased with increasing salinity. Sampled shells of M. edulis showed no significant trends in delta(18)O(calcite) based on size, cultured growth rates, or geographic collection location, suggesting that vital effects do not affect delta(18)O(calcite) in M. edulis. The broad modern and paleogeographic distribution of this bivalve, its abundance during the Holocene, and the lack of an intraspecies physiologic isotope effect demonstrated here make it an ideal nearshore paleoceanographic proxy throughout much of the North Atlantic Ocean.

Experimental Determination of Salinity, Temperature, Growth, and Metabolic Effects on Shell Isotope Chemistry of Mytilus Edulis Collected from Maine and Greenland

To study the effects of temperature, salinity, and life processes (growth rates, size, metabolic effects, and physiological/ genetic effects) on newly precipitated bivalve carbonate, we quantified shell isotopic chemistry of adult and juvenile animals of the intertidal bivalve Mytilus edulis (Blue mussel) collected alive from western Greenland and the central Gulf of Maine and cultured them under controlled conditions. Data for juvenile and adult M. edulis bivalves cultured in this study, and previously by Wanamaker et al. (2006), yielded statistically identical paleotemperature relationships. On the basis of these experiments we have developed a species-specific paleotemperature equation for the bivalve M. edulis [T degrees C = 16.28 (+/- 0.10) -4.57 (+/- 0.15) {delta(18)O(c) VPBD - delta(18)O(w) VSMOW} + 0.06 (+/- 0.06) {delta(18)O(c) VPBD - delta(18)O(w) VSMOW}(2); r(2) = 0.99; N = 323; p < 0.0001]. Compared to the Kim and O'Neil (1997) inorganic calcite equation, M. edulis deposits its shell in isotope equilibrium (delta(18)O(calcite)) with ambient water. Carbon isotopes (delta(13)C(calcite)) from sampled shells were substantially more negative than predicted values, indicating an uptake of metabolic carbon into shell carbonate, and delta(13)C(calcite) disequilibrium increased with increasing salinity. Sampled shells of M. edulis showed no significant trends in delta(18)O(calcite) based on size, cultured growth rates, or geographic collection location, suggesting that vital effects do not affect delta(18)O(calcite) in M. edulis. The broad modern and paleogeographic distribution of this bivalve, its abundance during the Holocene, and the lack of an intraspecies physiologic isotope effect demonstrated here make it an ideal nearshore paleoceanographic proxy throughout much of the North Atlantic Ocean.

Dependence of Eemian Greenland temperature reconstructions on the ice sheet topography

The influence of a reduced Greenland Ice Sheet (GrIS) on Greenland's surface climate during the Eemian interglacial is studied using a set of simulations with different GrIS realizations performed with a comprehensive climate model. We find a distinct impact of changes in the GrIS topography on Greenland's surface air temperatures (SAT) even when correcting for changes in surface elevation, which influences SAT through the lapse rate effect. The resulting lapse-rate-corrected SAT anomalies are thermodynamically driven by changes in the local surface energy balance rather than dynamically caused through anomalous advection of warm/cold air masses. The large-scale circulation is indeed very stable among all sensitivity experiments and the Northern Hemisphere (NH) flow pattern does not depend on Greenland's topography in the Eemian. In contrast, Greenland's surface energy balance is clearly influenced by changes in the GrIS topography and this impact is seasonally diverse. In winter, the variable reacting strongest to changes in the topography is the sensible heat flux (SHF). The reason is its dependence on surface winds, which themselves are controlled to a large extent by the shape of the GrIS. Hence, regions where a receding GrIS causes higher surface wind velocities also experience anomalous warming through SHF. Vice-versa, regions that become flat and ice-free are characterized by low wind speeds, low SHF, and anomalous low winter temperatures. In summer, we find surface warming induced by a decrease in surface albedo in deglaciated areas and regions which experience surface melting. The Eemian temperature records derived from Greenland proxies, thus, likely include a temperature signal arising from changes in the GrIS topography. For the Eemian ice found in the NEEM core, our model suggests that up to 3.1 °C of the annual mean Eemian warming can be attributed to these topography-related processes and hence is not necessarily linked to large-scale climate variations.

Alpine lacustrine varved record reveals summer temperature as main control of glacier fluctuations over the past 2250 years

Glacier fluctuations are a key indicator of changing climate. Their reconstruction beyond historical times unravels glacier variability and its forcing factors on long time scales, which can considerably improve our understanding of the climate–glacier relationship. Here, we present a 2250-year-long reconstruction of particle-mass accumulation rates recorded in the lacustrine sediments of Lake Trüebsee (Central Swiss Alps) that are directly related to glacier extent, thus reflecting a continuous record of fluctuations of the upstream-located Titlis Glacier. Mass accumulation rate values show strong centennial to multi-centennial fluctuations and reveal 12 well-pronounced periods of enhanced values corresponding to times of maximum extent of the neighboring Lower Grindelwald Glacier. This result supports previous studies of proglacial lake sediments that documented high mass accumulation rate values during glacier advances. The strong variability in the Lake Trüebsee mass accumulation rate record thus represents a highly sensitive paleoclimatic archive, which mirrors rapid and pronounced feedbacks of Titlis Glacier to climatic changes over the past 2250years. The comparison of our data with independent paleo-temperature reconstructions from tree rings suggests that variations in mean summer temperature were the primary driving factor of fluctuations of Titlis Glacier. Also, advances of Titlis Glacier occurred during the grand solar minima (Dalton, Maunder, Spörer, Wolf) of the last millennium. This relation of glacier extent with summer temperature reveals strong evidence that the mass balance of this Alpine glacier is primarily controlled by the intensity of glacier melting during summer.