Search

Understanding glaciers evolution is of major concern to evaluate their contribution to sea level rise in the context of global warming. Among the various processes involved in glacier dynamics, fractures like the calving of ice at the front of marine terminating glacier and crevasse formation affect the stress state, potentially modifying the glacier’s velocity. Crevasses also impact the melting rate. The fractures alter the roughness of the ground, increasing the amount of absorbed radiation, and open new networks in which the meltwater is likely to penetrate deeper toward the glacier bed.In this work we propose to model fractures in glacier based on finite strain elastoplasticity, using the Material Point Method (Wolper et al. 2021): we solve the classical governing equations for ice deformation in an Eulerian-Lagrangian framework and we use a strain softening Drucker-Prager constitutive model to simulate plasticity. Thanks to the Lagrangian part of the model, the fractures appear explicitly where high levels of total plastic deformation are reached. The behaviour of a glacier flowing over a step in the bedrock is investigated. The simulations show that crevasse patterns appear with regular spacing between the fractures. We perform a parametric study to determine which parameters affect the length of these patterns and potential dimensionless numbers are outlined.

Outlet glaciers and ice streams transport ice from the ice sheets to the ocean, where the glaciers lose mass by iceberg calving. Sermeq Kujalleq (SKK, Jakobshavn Isbræ) is one of the largest and most dynamic ice streams of the Greenland Ice Sheet with ice flow velocities up to 40 m/day. With extensive fieldwork and detailed repeated UAV surveys we aim at understanding the complex processes occuring at the ice stream margins and at the calving front of SKK. Such processes are often neclected in numerical models inducing uncertainties in projections of the ice sheet evolution.Within the framework of the COEBELI project we conducted drone photogrammetry surveys in July 2022 at SKK along other field measurements including in-situ GPS, GPRI, seismometers, and time-lapse imagery. Despite challenging weather conditions and constraints due to flying restrictions, we acquired 17 repeated flight surveys over the calving front of SKK during two weeks. As a result, we produced a large imagery data set, which was processed to infer high-resolution ortho-images and digital elevation models (DEM). Comparing the different products enables us to estimate changes in surface topography and ice dynamics. During the observation period several large calving events occurred allowing us to investigate the interaction between frontal processes and ice flow dynamics. With the very detailed data we can study crevasse opening, acceleration at the front, weaknesses in the ice and their origin, and the reaction of the glacier to large calving events.

Outlet glaciers and ice streams are the main channels through which ice sheets transport their mass towards the ocean. One of Greenland’s largest outlet glaciers Sermeq Kujalleq in Kangia (Jakobshavn Isbrae) has been broadly researched after experiencing a rapid retreat of the terminus and accompanying speedup to up to 40 m/day in the early 2000’s. However, such short-term ice dynamic variations remain poorly understood making numerical models difficult to constrain and predictions on future sea-level rise uncertain.The short-term ice dynamics of Sermeq Kujalleq consists in transient states and can only be captured by in-situ measurements of high spatial and temporal resolution. Glacier seismology has proven to be a valuable tool to study these dynamics, it provides data with a high temporal resolution and can provide information on processes happening below the ice surface. Within the COEBELI project we combine passive glacier seismology with global navigation satellite system (GNSS) receivers, long-range drones, time-lapse cameras and terrestrial radar interferometry to capture processes such as calving and basal sliding at their respective timescales.Here, we present results from a multi-array seismic deployment at Sermeq Kujalleq in Summer 2022. From May until September two arrays were deployed in the upstream part of the fast-flowing ice stream (>22 km from calving front) and one array on slower moving ice North of the main trunk. For a 3-week period in July, four more arrays were deployed on the fast-flowing ice stream closer to the calving front (50 Hz) icequakes. By studying these different signals, we are able to better constrain the processes and forces that control fluctuating ice-flow velocity and calving events.

Sermeq Kujalleq in Kangia (Jakobshavn Isbræ), Greenland is one of the most studied glaciers in the world mainly due to its recent retreat associated with extremely fast ice stream flow and high solid ice discharge. However, large limitations remain in the understanding of its short-term ice dynamics as the study of sub-daily variations, generally undetectable in spaceborne observations, requires high-rate field measurements that are challenging to acquire. Here, we present glacier surface velocities determined in Post-Processed Kinematic (PPK) mode from eight autonomous Global Navigation Satellite System (GNSS) stations deployed in July 2022 along the ice stream at a distance of 4 to 30 kilometers from the calving front. During this field campaign, we identified an 8-hour-long glacier speedup which was recorded at all GNSS stations and reached up to 11% of the pre-event velocity, followed by a 12-hour-long slowdown of similar magnitude. We further found the peak velocity was first measured at a GNSS station 16 kilometers away from the calving front, then recorded consecutively at each of the three other downstream GNSS stations with a positive time lag corresponding to a ~3 km/h wave propagation speed. At the station closest to the calving front, the timing of peak velocity corresponded to the occurrence of large-scale calving events. We further present line-of-sight glacier surface velocities measured along three shear margin transects with a terrestrial radar interferometer deployed simultaneously with the GNSS array. Across all profiles, we observed a widespread and simultaneous response of fast- and slow-moving ice suggesting a strong coupling between the main trunk and the shear margins of the ice stream.

Future mass loss projections of the Greenland ice sheet require understanding of the processes at a glacier terminus, especially of iceberg calving. We present detailed and high-rate terrestrial radar interferometer observations of Eqip Sermia and Bowdoin Glacier, two outlet glaciers in Greenland with comparable dimensions and investigate iceberg calving, surface elevation, velocity, strain rates and their links to air temperature, tides and topography. The results reveal that the two glaciers exhibit very different flow and calving behaviour on different timescales. Ice flow driven by a steep surface slope with several topographic steps leads to high velocities, areas of extension and intense crevassing, which triggers frequent but small calving events independent of local velocity gradients. In contrast, ice flow under smooth surface slopes leaves the ice relatively intact, such that sporadic large-scale calving events dominate, which initiate in areas with high shearing. Flow acceleration caused by enhanced meltwater input and tidal velocity variations were observed for terminus sections close to floatation. Firmly grounded terminus sections showed no tidal signal and a weak short-term reaction to air temperature. These results demonstrate reaction timescales to external forcings from hours to months, which are, however, strongly dependent on local terminus geometry., Journal of Glaciology, 69 (275), ISSN:0022-1430, ISSN:1727-5652

Here, we present a compilation of 85 ice temperature profiles from 79 boreholes from the Greenland Ice Sheet and peripheral ice caps, as well as local ice caps in the Canadian Arctic. Only 25 profiles (32 %) were previously available in open-access data repositories. The remaining 54 profiles (68 %) are being made digitally available here for the first time. These newly available profiles, which are associated with pre-2010 boreholes, have been submitted by community members or digitized from published graphics and/or data tables. All 85 profiles are now made available in both absolute (meters) and normalized (0 to 1 ice thickness) depth scales, and are accompanied by extensive metadata. This metadata includes a transparent description of data provenance. The ice temperature profiles span 70 years, with the earliest profile being from 1950 at Camp VI, West Greenland. To highlight the value of this database in evaluating ice flow simulations, we compare the ice temperature profiles from the Greenland Ice Sheet with an ice flow simulation by the Parallel Ice Sheet Model (PISM). We find a cold bias in modeled near-surface ice temperatures within the ablation area, a warm bias in modeled basal ice temperatures at inland cold-bedded sites, and an apparent underestimation of deformational heating in high-strain settings. These biases provide process-level insight on simulated ice temperatures., The Cryosphere Discussions, ISSN:1994-0432, ISSN:1994-0440

Sermeq Kujalleq in Kangia (Jakobshavn Isbræ), Greenland has been extensively investigated over the past decades due to its recent retreat associated with extremely fast ice stream flow and high solid ice discharge. However, its short-term dynamics still remain poorly understood as they consist in transient states that can only be captured by high spatial and temporal in situ measurements. In the new COEBELI project, we aim at combining high resolution field data sets from seismic arrays, global navigation satellite system (GNSS) receivers, long-range uncrewed aerial vehicles and terrestrial radar interferometers (TRI) to achieve a comprehensive and detailed study of the short-term ice stream dynamics. Here, we present TRI and GNSS retrievals of surface velocity and elevation acquired during the first, exploratory field campaign of the COEBELI project in summer 2021. Seven kilometers away from the calving front, we specifically identified a slowdown of 1.12 m d-1 within a single day in the main trunk of the ice stream. While the absolute slowdown is larger in the main trunk than in the outer area of the shear margin (1.12 m d-1 versus 0.75 m d-1), it corresponds to a larger fraction of the pre-slowdown velocity in the latter zone (-4.48% versus -7.94%). We further discuss the challenges associated with the acquisition, processing and analysis of high-resolution data sets for the study of such complex and dynamic environments.

Glacier calving plays a key role in the recently observed dynamic massloss of the Greenland ice sheet. Calving waves, generated by thesudden detachment of ice from the glacier terminus, can reach tens ofmeters of height and have devastating effects upon impact onsurrounding shores. In this study, we describe a new method for thedetection of source location and timing of calving waves, and theanalysis of their magnitude and spreading properties using aterrestrial radar interferometer (TRI). This method was applied to11,500 minute-interval TRI acquisitions from Eqip Sermia, Greenland.More than 2,000 calving waves were detected within sevendays. Quantitative assessment with a Wave Power Index (WPI) showedspatially distinctive patterns: the sector of the calving front endingin deep water shows a higher wave activity (+49%) with highercumulative WPI (+34%) than the shallow sector. In combination witha detection of meltwater plume locations, we highlighted a 2.3 timeshigher occurrence of visible meltwater plumes in the deep sector than theshallow one. We found both the cumulated WPI and the number of wavesto increase by more than 80% in the presence of a meltwater plumein the deep sector while only by 30% in the shallow sector. Wetherefore explain the higher calving activity in the deep sector to be strongly related to a combination of higher occurrence of meltwater plumes and more efficient calving enhancement linked to better connections to deep warm waters.

The recent rapid retreat of many Arctic outlet glaciers has been attributed to increased oceanic melt, but the relationship between oceanic melt and iceberg calving remains poorly understood. Here, we employ a transient finite‐element model that simulates oceanic melt and ice break‐off at the terminus. The response of an idealized tidewater glacier to various submarine melt rates and seasonal variations is investigated. Our modeling shows that for zero to low oceanic melt, the rate of volume loss at the front is similar or higher than for intermediate oceanic melt rates. Only very high melt rates lead to increasing volume losses. These results highlight the complex interplay between oceanic melt and calving and question the general assumption that increased submarine melt leads to higher calving fluxes and enhanced retreat. Models for tidewater glacier evolution should therefore consider calving and oceanic melt as tightly coupled processes rather than as simple, additive parametrizations.

We observed two outlet glaciers in West- and Nordwest-Greenland with a terrestrial radar interferometer (TRI), pressure sensors and time-lapse cameras over six and two years, respectively. The resulting detailed dataset provides us with insights on the calving process and the changes in front geometry over the last years. Since the two glaciers are characterised by different geometries and velocity fields, the influence of those parameters on the calving process can be investigated. The combination of the three different observation methods enable us to overcome their individual disadvantages. With the time-lapse camera taking pictures of the glacier front every 10 seconds, we detect all calving events of different sizes and styles but cannot quantify the volume. We used the TRI to quantify the volumes of aerial calving events by DEM differentiation. Further, calving waves measured with pressure sensors are used to distinguish between different calving types. We develop a relationship between calving volumes and wave heights and use this as an additional indirect method to estimate calving volumes. We find that the calving style and size as well as the front geometry is mainly controlled by the bed topography and the presence of a subglacial discharge plume. The location of the plume is observed to migrate from year to year, which leads also to changes in the calving pattern. Calving style and pattern as well as glacier velocity fields and geometry changes are additionally compared with environmental conditions such as the air temperature and the presence of ice-mélange in the proglacial fjord. In years with an early spring we find different front characteristics and calving patterns than for years with colder conditions.

The impact of summer surface melt on the dynamics of the Greenland Ice Sheet is modulated by the state of the subglacial hydrologic system. Studies of ice motion indicate that efficiency of the subglacial system increases over the melt season, decreasing the sensitivity of ice motion to surface melt. However, these inferences are based on limited indirect observations of the subglacial hydrologic system that leave many factors poorly constrained, particularly the presence and stability of subglacial channels. Here we use observations from 11 GPS stations, from which we derive ice velocity, longitudinal strain rates, and basal uplift, alongside observations of surface ablation and supraglacial lake drainage events, to explore the coevolution of ice motion and the subglacial hydrologic system in the Pakitsoq region of western Greenland during the 2011 melt season. We observe ice acceleration after the onset of local surface melting, followed by gradual ice deceleration, consistent with the pattern expected from increased subglacial drainage efficiency. Supraglacial lake drainages appear to precipitate ice deceleration and increased basal traction, suggesting that lake drainages effectively reorganize the local subglacial hydrologic system into a more efficient state that persists through the remainder of the melt season. At high elevations, ice velocity and inferred basal uplift suggest that continued cavity growth or sediment behavior, not subglacial channelization, drive the apparent increase in subglacial efficiency. Our results further indicate that these transient perturbations are critical in the seasonal evolution of ice motion.

Moulins permit access of surface meltwater to the glacier bed, causing basal lubrication and ice speedup in the ablation zone of western Greenland during summer. Despite the substantial impact of moulins on ice dynamics, the conditions under which they form are poorly understood. We assimilate a time-series of ice surface velocity from a network of eleven Global Positioning System receivers into an ice sheet model to estimate ice sheet stresses during winter, spring, and summer in an approx. 30 x 10 km region. Surface-parallel von Mises stress increases slightly during spring speedup and early summer, sufficient to allow formation of 16% of moulins mapped in the study area. In contrast, 63% of moulins experience stresses over the tensile strength of ice during a short (hours) supraglacial lake drainage event. Lake drainages appear to control moulin density, which is itself a control on subglacial drainage efficiency and summer ice velocities.

In previous studies, the incorporation of polychlorinated biphenyls (PCBs) has been quantified in the accumulation areas of Alpine glaciers. Here, we introduce a model framework that quantifies mass fluxes of PCBs in glaciers and apply it to the Silvretta glacier (Switzerland). The models include PCB incorporation into the entire surface of the glacier, downhill transport with the flow of the glacier ice, and chemical fate in the glacial lake. The models are run for the years 1900-2100 and validated by comparing modeled and measured PCB concentrations in an ice core, a lake sediment core, and the glacial streamwater. The incorporation and release fluxes, as well as the storage of PCBs in the glacier increase until the 1980s and decrease thereafter. After a temporary increase in the 2000s, the future PCB release and the PCB concentrations in the glacial stream are estimated to be small but persistent throughout the 21st century. This study quantifies all relevant PCB fluxes in and from a temperate Alpine glacier over two centuries, and concludes that Alpine glaciers are a small secondary source of PCBs, but that the aftermath of environmental pollution by persistent and toxic chemicals can endure for decades.

Our ability to measure, quantify and assimilate hydrological properties and processes of snow in operational models is disproportionally poor compared to the significance of seasonal snowmelt as a global water resource and major risk factor in flood and avalanche forecasting. We show here that strong electrical self-potential fields are generated in melting in situ snowpacks at Rhone Glacier and Jungfraujoch Glacier, Switzerland. In agreement with theory, the diurnal evolution of self-potential magnitudes ( ∼ 60–250mV) relates to those of bulk meltwater fluxes (0–1.2 × 10−6m3s−1) principally through the permeability and the content, electrical conductivity and pH of liquid water. Previous work revealed that when fresh snow melts, ions are eluted in sequence and electrical conductivity, pH and self-potential data change diagnostically. Our snowpacks had experienced earlier stages of melt, and complementary snow pit measurements revealed that electrical conductivity ( ∼ 1–5 × 10−6Sm−1) and pH ( ∼ 6.5–6.7) as well as permeabilities (respectively ∼ 9.7 × 10−5 and ∼ 4.3 × 10−5m2 at Rhone Glacier and Jungfraujoch Glacier) were invariant. This implies, first, that preferential elution of ions was complete and, second, that our self-potential measurements reflect daily changes in liquid water contents. These were calculated to increase within the pendular regime from ∼ 1 to 5 and ∼ 3 to 5.5% respectively at Rhone Glacier and Jungfraujoch Glacier, as confirmed by ground truth measurements. We conclude that the electrical self-potential method is a promising snow and firn hydrology sensor owing to its suitability for (1) sensing lateral and vertical liquid water flows directly and minimally invasively, (2) complementing established observational programs through multidimensional spatial mapping of meltwater fluxes or liquid water content and (3) monitoring autonomously at a low cost. Future work should focus on the development of self-potential sensor arrays compatible with existing weather and snow monitoring technology and observational programs, and the integration of self-potential data into analytical frameworks., The Cryosphere, 10 (1), ISSN:1994-0416, ISSN:1994-0424

Penetration of surface meltwater to the bed of the Greenland Ice Sheet each summer causes an initial increase in ice speed due to elevated basal water pressure, followed by slowdown in late summer that continues into fall and winter. While this seasonal pattern is commonly explained by an evolution of the subglacial drainage system from an inefficient distributed to efficient channelized configuration, mounting evidence indicates that subglacial channels are unable to explain important aspects of hydrodynamic coupling in late summer and fall. Here we use numerical models of subglacial drainage and ice flow to show that limited, gradual leakage of water and lowering of water pressure in weakly connected regions of the bed can explain the dominant features in late and post melt season ice dynamics. These results suggest that a third weakly connected drainage component should be included in the conceptual model of subglacial hydrology., Surface meltwater draining to the bed of the Greenland Ice Sheet each summer causes ice flow changes inconsistent with the prevailing theory of channelizing subglacial drainage. Here, the authors show this is caused by limited, gradual leakage of water from previously ignored weakly connected regions of the bed.

Current understanding of ice dynamics predicts that increasing availability and variability of meltwater will have an impact on basal motion and therefore on the evolution and future behavior of the Greenland ice sheet. We present measurements of ice deformation, subglacial water pressure, and surface velocity that show periodic and episodic variations on several time scales (seasonal, multiday, and diurnal). These variations, observed with GPS and sensors at different depths throughout the ice column, are not synchronous but show delayed responses of ice deformation with increasing depth and basal water pressure in antiphase with surface velocity. With the help of a Full-Stokes ice flow model, these observations are explained as ice motion in a caterpillar-like fashion. Caused by patches of different basal slipperiness, horizontal stress transfer through the stiff central part of the ice body leads to spatially varying surface velocities and ice deformation patterns. Variation of this basal slipperiness induces characteristic patterns of ice deformation variability that explain the observed behavior. Ice flow in the ablation zone of the Greenland ice sheet is therefore controlled by activation of basal patches by varying slipperiness in the course of a melt season, leading to caterpillar-like ice motion superposed on the classical shear deformation.

Simultaneous observations of moulins and boreholes in western Greenland show that water delivery to the base of the ice sheet via moulins affects short-term ice velocity fluctuations, but not late-season ice velocity decelerations. Increased meltwater delivery to the base of the Greenland Ice Sheet will increase the ice sheet velocity, accelerating its inevitable rush towards the ocean and subsequent sea level rise. Or will it? Debate on this topic is at the forefront of cryospheric research, but has been hampered by the lack of simultaneous observations of hydraulic head in moulins, vertical shafts that deliver water to the base of the ice sheet, and boreholes, which monitor basal water pressure. Lauren Andrews and colleagues now provide these observations from a small region in western Greenland and reveal that water delivery by moulins into channelized basal flow is indeed linked to short-term fluctuations in ice velocity. Late season decelerations in ice velocity, however, seem to be controlled by changes in unchannelized flow, rather than any shifts in the moulin system. Seasonal acceleration of the Greenland Ice Sheet is influenced by the dynamic response of the subglacial hydrologic system to variability in meltwater delivery to the bed1,2 via crevasses and moulins (vertical conduits connecting supraglacial water to the bed of the ice sheet). As the melt season progresses, the subglacial hydrologic system drains supraglacial meltwater more efficiently1,2,3,4, decreasing basal water pressure4 and moderating the ice velocity response to surface melting1,2. However, limited direct observations of subglacial water pressure4,5,6,7 mean that the spatiotemporal evolution of the subglacial hydrologic system remains poorly understood. Here we show that ice velocity is well correlated with moulin hydraulic head but is out of phase with that of nearby (0.3–2 kilometres away) boreholes, indicating that moulins connect to an efficient, channelized component of the subglacial hydrologic system, which exerts the primary control on diurnal and multi-day changes in ice velocity. Our simultaneous measurements of moulin and borehole hydraulic head and ice velocity in the Paakitsoq region of western Greenland show that decreasing trends in ice velocity during the latter part of the melt season cannot be explained by changes in the ability of moulin-connected channels to convey supraglacial melt. Instead, these observations suggest that decreasing late-season ice velocity may be caused by changes in connectivity in unchannelized regions of the subglacial hydrologic system. Understanding this spatiotemporal variability in subglacial pressures is increasingly important because melt-season dynamics affect ice velocity beyond the conclusion of the melt season8,9,10.

The geological nature of glacier beds plays a key role in ice sheet dynamics. Whereas little is known about Greenland’s subglacial geology, the presence of basal sediments is a necessary condition for fast-flowing Antarctic ice streams. Such sediments sustain subglacial till layers, which if water saturated and under high pore-water pressure provide little resistance to ice flow. Using receiver function modeling of teleseismic P-waves, we report a thick (at least tens of meters) sediment layer beneath a site in Greenland’s ablation zone, ∼15 km away from the western ice sheet margin. Although we do not discuss the origin or detailed nature of these subglacial sediments, we suggest that they are capable of sustaining a till layer. Due to the prevalence of an inefficient, pressurized subglacial drainage system, this till layer would typically be under high pore pressures and fail at the shear stresses transmitted from the overlaying ice. Ice flow resistance is thus focused on regions where till is consolidated or absent, or where form drag over obstacles takes place. In contrast to Antarctica, Greenland’s surface melt now affects practically all latitudes, and efficient hydraulic connections between ice sheet surface and bed are pervasive throughout the ablation zone. This implies that rapid, melt-induced changes in subglacial water pressure are possible. The time required for basal till strength to react to these changes depends on the till’s properties. We estimate that formation and destruction of flow resistance can occur on time scales of less than a few years. This could lead to changes in ice flow that are currently difficult to predict.

Mountain glaciers sample a combination of climate fields – temperature, precipitation and radiation – by accumulation and melting of ice. Flow dynamics acts as a transfer function that maps volume changes to a length response of the glacier terminus. Long histories of terminus positions have been assembled for several glaciers in the Alps. Here I analyze terminus position histories from an ensemble of seven glaciers in the Alps with a macroscopic model of glacier dynamics to derive a history of glacier equilibrium line altitude (ELA) for the time span 400–2010 C.E. The resulting climatic reconstruction depends only on records of glacier variations. The reconstructed ELA history is similar to recent reconstructions of Alpine summer temperature and Atlantic Multidecadal Oscillation (AMO) index, but bears little resemblance to reconstructed precipitation variations. Most reconstructed low-ELA periods coincide with large explosive volcano eruptions, hinting at a direct effect of volcanic radiative cooling on mass balance. The glacier advances during the LIA, and the retreat after 1860, can thus be mainly attributed to temperature and volcanic radiative cooling., The Cryosphere, 8 (2), ISSN:1994-0416, ISSN:1994-0424

Dynamic changes on Greenland outlet glaciers are a primary driver for increases in ice-sheet mass loss and its contribution to sea-level rise. One dramatic example of such change has been observed at Jakobshavn Isbræ, which has thinned, retreated and doubled in speed since the early 2000s. Complementary to large changes on decadal scales, we observe the glacier response on shorter timescales, driven by tidal forcing and calving events. During a 14 day period in August 2009, we documented changes in geometry and speed near the terminus. On this timescale, ice flow responds to forcing at the front from iceberg calving and ocean tides. We observe a step-increase in velocity near the terminus during a large calving event, with transient deceleration in the days following the event. A simple calving-response model explains 94–99% of variations in detrended positions at the four sites considered. During each day, variability due to tidal forcing covers 10–90% of the variability that remains after removing effects accounted for by the calving-response model. The influence of the tidal forcing on flow decays upstream with a characteristic length scale of 2 km, comparable with about two ice thicknesses.

During summer 2011, seismic activity in the ablation zone of the western Greenland ice sheet (GrIS) was monitored using a network of three-component seismometers. The seismic record includes a large variety of icequakes and seismic tremors that demonstrate a clear correlation with subglacial water flow. We verified the existence of well-known shallow icequakes (related to surface crevasse formation), deep icequakes (located at 100–160 m depth) and narrow-banded short-term seismic tremors (tens of seconds in duration). In addition, we present previously unreported long-term tremors lasting several hours. Using attenuation of the measured tremor amplitude, we locate the epicentre of this long-term tremor to a large moulin within our study area. Between 3 and 11 Hz, our continuous seismic record is dominated by this ‘moulin tremor’ and shows strong correlation with the water level of the generating moulin. We argue that monitoring of icequake and glacial tremor sources bears high potential for investigating glacier hydraulics and dynamics, and is thus an ideal supplement to traditional glaciological measurements.

[1] Cold ice within a polythermal ice body controls its flow dynamics through the temperature dependence of viscosity, and affects glacier hydrology by blocking water flow paths. Lakes on the surface, linked by persistent, deeply incised meltwater streams, are hallmark features of cold ice in the ablation zone of a glacier or ice sheet. Ice radar is a convenient method to map scattering from internal water bodies present in ice at the pressure melting temperature (PMT). Consequently, lack of internal scatters is indicative of cold ice. We use a helicopter-borne 30 MHz ice radar to delineate the extent of cold ice within Grenzgletscher (Zermatt, Swiss Alps). The inferred thermal structure is validated with temperature measurements in 15 deep boreholes, showing excellent agreement. The cold ice occupies 80–90 % of the total ice thickness in a 400 m wide flow band along the central flow line. Quantitative interpretation of ice radar scattering power indicates a decrease of ice water content between PMT and 0.5 K below PMT, as predicted by theory, and observed in the laboratory. The cold ice which emerges at the surface in the lower ablation zone is impermeable to water, and is thus devoid of moulins if not crevassed. The surface water from melt and rain is routed through deeply incised, persistent streams and lakes, and cryoconite holes are frequent, in stark contrast to the adjacent temperate ice from other tributaries. The cold ice thus has a strong control on glacier hydrology, but is likely to change due to continued warming.

We use time-lapse photography, MODIS satellite imagery, ocean wave measurements and regional broadband seismic data to demonstrate that icebergs that calve from Jakobshavn Isbræ, Greenland, can generate ocean waves that are detectable over 150 km from their source. The waves, which are recorded seismically, have distinct spectral peaks, are not dispersive and persist for several hours. On the basis of these observations, we suggest that calving events at Jakobshavn Isbræ can stimulate seiches, or basin eigenmodes, in both Ilulissat Icefjord and Disko Bay. Our observations furthermore indicate that coastal, land-based seismometers located near calving termini (e.g. as part of the new Greenland Ice Sheet Monitoring Network (GLISN)) can aid investigations into the largely unexplored, oceanographic consequences of iceberg calving.

The length change record of 91 glaciers in the Swiss Alps was analyzed with a novel macroscopic glacier model (LV-model). Based on a history of equilibrium line variations, synthetic length change data were calculated. From the LV-models matching best the measured length changes, characteristic parameters were obtained. The volume time scale thus determined ranges from 5 to 170 years for glaciers of different slope and length. The analysis shows that the observed glacier length changes cannot be reproduced with an equilibrium line variation based on temperature and precipitation alone. The equilibrium line has to be lowered by 100 to 200 meters during several phases of the Little Ice Age (in the time span 1650 to 1850) to obtain observed glacier responses. Such an effect might be attributable to either higher winter precipitation in the Alps, or to radiation forcing.

Recently, increasing concentrations of persistent organic pollutants (POPs) have been observed in the sediment of glacier-fed Lake Oberaar, Switzerland. Melting glaciers have been suggested as a secondary source of POPs released to Alpine lakes. Here we further investigate whether climate warming may accelerate the release of POPs previously deposited to Alpine glaciers ("glacier hypothesis"). To test this hypothesis, a dynamic multimedia mass balance model is developed for the catchment area of Lake Oberaar and is applied to polychlorinated biphenyls (PCBs), polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs), and dichlorodiphenyl trichloroethane (DDT). This lake model is combined with two other models. The first is a dynamic multimedia mass balance model parametrized for the Swiss lowlands that is used to calculate (on the basis of historical emission data) the atmospheric concentrations that are an advective input into the model of the lake catchment. The second is a flow model of Oberaar Glacier that determines the residence time of persistent chemicals in the glacier after their deposition to the glacier surface. According to results from these three models in combination, the release of POPs by the glacier is currently increasing and accounts for the observed increase in concentrations in the lake sediment. The models indicate that approximately half of the amount of PCBs, PCDD/Fs, and DDT initially incorporated into the glacier ice is still stored in the glacier. Under the assumption that the climate is warming, accelerated release of POPs is to be expected for the future; in a model run where no climate warming is assumed, the period of time required for release of the same amount of chemicals is longer by several decades than in the scenario with a changing climate.

The transient response of glaciers to climate variations is investigated with a novel low-order model of an idealized glacier resting on uniformly inclined bedrock. The model consists of two volumes representing the accumulation and ablation areas, which are joined by a flux gate controlling mass flow according to the shallow-ice approximation. Under the assumption of a constant vertical mass-balance gradient, a volume–length scaling relation is derived which depends explicitly on mass-balance gradient and bedrock slope. Analytic expressions for the volume and area timescales are given, which are inversely proportional to the mass-balance gradient and a geometric factor which is the ratio between the vertical extent of the ablation area and the ice thickness at the equilibrium line. From the low-order model, a dynamical system in length and volume is obtained. Results from this system are in good agreement with solutions obtained from a transient finite-element model solving the full force-balance and mass-conservation equations. Under periodic forcing there are significant deviations of the length response, which show that the usual relaxation-type parameterization of length change is not well suited for short-term reactions. Switching on and off periodic climate forcing, the model glaciers show surprisingly large initial and final transient responses that have not been investigated before. These results are of significance for the interpretation of length and thickness changes observed on glaciers.

The puzzling fact that Jakobshavn Isbræ, West Greenland, is flowing very fast but without any significant seasonal velocity changes, despite big amounts of surface-derived meltwater entering the ice stream (Echelmeyer and Harrison, 1990), has been explained by a combination of different types of measurements. It is now well established from seismic measurements and radio-echo sounding that Jakobshavn Isbræ flows through a deeply eroded subglacial trench that, even 50 km inland of the grounding line, extends as far as 1500 m below sea level (Clarke and Echelmeyer, 1996; Legarsky and Huang, 2006). Temperature measurements in boreholes down to 65% of the 2500 m thick ice stream at site B, some 50 km upstream of the calving front (Fig. 1b; Iken and others, 1993), were used to infer the presence of a substantial layer of temperate ice, the thickness of which was estimated to be at least 300 m by modeling and matching internal layering structures (Funk and others, 1994; Lüthi and others, 2002). The presence of a thick layer of temperate ice under very high driving stress allows for high ice-deformation rates, which contribute substantially to the observed fast flow velocities. Basal motion, while certainly important, seems to be barely influenced by the seasonal meltwater input (Echelmeyer and Harrison, 1990).

We present the thermal distribution in the confluence area of Gorner- and Grenzgletscher, Valais, Switzerland. The area was mapped by ice-penetrating radar at 1–5 and 40 MHz. The higher-frequency data reveal a thick surface layer of low backscatter in the center of the Grenzgletscher branch. Based on datasets of borehole-temperature measurements and flow velocity, we interpret this as a thick layer of cold ice, advected from the accumulation region of Grenzgletscher. Along seven profiles the base of the low-backscatter zone can be found at a maximum depths between approximately 100 and 200 m. Laterally, the layer extends some 400 m, ∼1/3 of the width of the Grenzgletscher branch. The lower boundary of the low-backscatter zone is systematically higher than the cold–temperate transition surface found in the boreholes. This discrepancy is attributed to the direct sensitivity of radar backscatter to liquid-water inclusions, rather than to the temperature distributions as observed in boreholes. We present the current state of the cold layer and discuss its influence on other glacier characteristics.

Glaciological and related observations from 1961 to 2005 at the summit of Mt Wrangell (62.00° N, 144.02°W; 4317m a.s.l.), a massive glacier-covered shield volcano in south-central Alaska, show marked changes that appear to have been initiated by the Great Alaska Earthquake (Mw = 9.2) of 27 March 1964. The 4×6 km diameter, ice-filled Summit Caldera with several post-caldera craters on its rim, comprises the summit region where annual snow accumulation is 1–2m of water equivalent and the mean annual temperature, measured 10 m below the snow surface, is –20°C. Precision surveying, aerial photogrammetry and measurements of temperature and snow accumulation were used to measure the loss of glacier ice equivalent to about 0.03 km3 of water from the North Crater in a decade. Glacier calorimetry was used to calculate the associated heat flux, which varied within the range 20–140Wm–2; total heat flow was in the range 20–100MW. Seismicity data from the crater’s rim show two distinct responses to large earthquakes at time scales from minutes to months. Chemistry of water and gas from fumaroles indicates a shallow magma heat source and seismicity data are consistent with this interpretation.

We present results from a chemical fate model quantifying incorporation of polychlorinated biphenyls (PCBs) into the Silvretta glacier, a temperate Alpine glacier located in Switzerland. Temperate glaciers, in contrast to cold glaciers, are glaciers where melt processes are prevalent. Incorporation of PCBs into cold glaciers has been quantified in previous studies. However, the fate of PCBs in temperate glaciers has never been investigated. In the model, we include melt processes, inducing elution of water-soluble substances and, conversely, enrichment of particles and particle-bound chemicals. The model is validated by comparing modeled and measured PCB concentrations in an ice core collected in the Silvretta accumulation area. We quantify PCB incorporation between 1900 and 2010, and discuss the fate of six PCB congeners. PCB concentrations in the ice core peak in the period of high PCB emissions, as well as in years with strong melt. While for lower-chlorinated PCB congeners revolatilization is important, for higher-chlorinated congeners, the main processes are storage in glacier ice and removal by particle runoff. This study gives insight into PCB fate and dynamics and reveals the effect of snow accumulation and melt processes on the fate of semivolatile organic chemicals in a temperate Alpine glacier.

Our ability to measure, quantify and assimilate hydrological properties and processes of snow in operational models is disproportionally poor compared to the significance of seasonal snowmelt as a global water resource and major risk factor in flood and avalanche forecasting. Encouraged by recent theoretical, modelling and laboratory work, we show here that the diurnal evolution of aerially-distributed self-potential magnitudes closely track those of bulk meltwater fluxes in melting in-situ snowpacks at Rhone and Jungfraujoch glaciers, Switzerland. Numerical modelling infers temporally-evolving liquid water contents in the snowpacks on successive days in close agreement with snow-pit measurements. Muting previous concerns, the governing physical and chemical properties of snow and meltwater became temporally invariant for modelling purposes. Because measurement procedure is straightforward and readily automated for continuous monitoring over significant spatial scales, we conclude that the self-potential geophysical method is a highly-promising non-intrusive snow-hydrological sensor for measurement practice, modelling and operational snow forecasting.

Observations of surface motion and ice deformation from 2002–03 were used to infer mean stress fields in a cross-section of Black Rapids Glacier, Alaska, USA, over seasonal timescales. Basal shear stresses in a well-defined zone north of the center line (orographic left) were approximately 7% and 16% lower in spring and summer, respectively, than in winter. Correspondingly higher stresses were found near the margins. These changes in the basal shear stress distribution were sufficiently large to cause mean surface velocities to be 1.2 and 1.5 times larger in spring and summer than in winter. These results were inferred with a simple inverse finite-element flow model that can successfully reproduce bulk surface velocities and tiltmeter data. Stress redistribution between the well-defined zone and the margins may also occur over much shorter time periods as a result of rapidly changing basal conditions (ice–bed decoupling or enhanced till deformation), thereby causing large variations in surface velocity and strongly influencing the glacier’s net motion during summer.

In the frame of the Swiss Alpine Airborne SAR Experiment 2003 (SASARE), multitemporal/multibaseline polarimetric synthetic aperture radar (SAR) interferometry (Pol-InSAR) datasets were acquired with the German Experimental SAR (E-SAR) platform at L- and P-band in the region of the Aletsch glacier in Switzerland. For SASARE, several sub-test-sites had been selected below and significantly above the equilibrium line. We provide results from the analysis of the multibaseline Pol-InSAR L- and P-band measurements of the polarization-dependent complex interferometric coherence of an Alpine glacier in the cold infiltration zone. Interferometric volume decorrelation is demonstrated as a function of polarization. Conventional polarimetric decomposition techniques confirm the polarization-dependent backscattering properties of the analyzed snow and ice surfaces.

Continuum damage mechanics describes the progressive deterioration of material subjected to loading. Jointly used with a level-set method, it proves to be a promising approach to computing the interface motion of a damaged material. For polycrystalline ice, a local isotropic damage evolution law (generalized Kachanow’s law) applied to Glen’s flow law allows the description of tertiary creep and facilitates the modeling of crevasse opening using a failure criterion based on damage accumulation. The use of a level-set method permits the description, in a continuum approach, of the motion of a fractured glacier surface. Using these methods, a model is developed. The ability of this model to describe phenomena connected to crevasse opening is presented. The rupture of a large ice block from a hanging glacier is computed and analyzed. The regular acceleration of such an unstable ice block prior to its collapse is calculated and compared to the acceleration function obtained from observations. A good agreement between the two acceleration functions was found.

Ice temperature profiles from the Greenland Ice Sheet contain information on the deformation history, past climates and recent warming. We present full-depth temperature profiles from two drill sites on a flowline passing through Swiss Camp, West Greenland. Numerical modeling reveals that ice temperatures are considerably higher than would be expected from heat diffusion and dissipation alone. The possible causes for this excess heat are evaluated using a Lagrangian heat flow model. The model results reveal that the observations can be explained with a combination of different processes: enhanced dissipation (strain heating) in ice-age ice, temperate paleo-firn, and cryo-hydrologic warming in deep crevasses.

Modelling the heat flow in small, cold high-altitude glaciers is important for the interpretation of paleoclimatic data from ice cores. Coupled glacier-flow and heat-flow models are presented that incorporate the densification, heat advection and possible phase transitions at the permafrost boundaries within the bedrock. Marked bends observed in the temperature profiles from two recent boreholes on Colle Gnifetti, Swiss Alps, are interpreted with the help of a transient heat-flow model, driven with a temperature history. The conclusion is that substantial warming of the mean firn temperature at shallow depths has taken place over the last few decades. This has not been observed before in cold-firn regions of the Alps. Modelled heat fluxes in the Monte Rosa massif are strongly influenced by the mountain topography. This leads to a spatial variability of the temperature gradient near the glacier base which has been observed in boreholes to the bedrock. In order to match the measured temperature profiles in the glacier, the vertical heat flux at great depth must be set to an extremely low value. It is shown with the help of the transient heat-flow model that this is a paleoclimatic effect, possibly enhanced by a degrading permafrost base.

A flow model for cold firn proves to be very successful in calculating the age-depth relation of several deep ice cores drilled on Colle Gnifetti, Monte Rosa, Swiss Alps. The compressibility of firn is taken into account by an appropriate constitutive equation, first employed in glaciology by Gagliardini and Meyssonnier (1997), which is implemented in a finite-element code. Flow models of the Colle Gnifetti saddle glaciation in two and three dimensions are based on digital elevation maps of the surface σnd the bedrock, based on radio-echo soundings of the ice thickness. Firn density and the englacial temperature fields are either prescribed or calculated in coupled models. Measured surface velocities, density profiles, the ages of chemically dated layers in ice cores and the closure of a 100 m deep borehole provide benchmarks for the models. The good agreement of modeled and measured quantities confirms that the model includes the relevant physical processes and particularly that the firn flow law is well suited for this type of glacier. The study provides new constraints on the age of the ice near the base as well as the source regions of the ice in the cores.An exceptional flow behavior of the basal ice layer was detected in measurements of borehole closure and inclination. Measurtxl6ed deformation rates exceed upper bounds derived from the flow models, and are thus attributed to altered rheological properties.

Over a two-year period, continuous measurements of temporal changes in tilt, conducted with a string of tilt meters in a borehole on Unteraargletscher, Bernese Alps, Switzerland, have been used to estimate the basal-motion component. This estimation is based on a comparison of the measurements with synthetic tilt curves, computed using a parameterization of a simplified flow field. The best agreement is found for a ratio of basal motion to forward motion due to ice deformation (slip ratio) equal to about 1.2. Measured tilt curves exhibit a number of different transient features. While an overall increase in tilt angle is observed at every tilt-meter location, two of the sensors recorded anomalous tilt behaviour. These anomalies are characterized by sudden and drastic variations in tilt. A particularly intriguing example of such short-term tilt variations was recorded with a tilt meter positioned 40 m above the bed during the 1997 summer melt season.

[1] A novel method to reconstruct glacier volume changes from a measured length record is presented and tested for 13 glaciers in the Swiss Alps. The response of a glacier to changes in climate is modeled with a two-parameter dynamical system in the variables “length” and “volume”. Driven by a history of equilibrium line altitude (ELA), the model yields variations of glacier length and volume. A dynamically equivalent simple model (DESM) is determined for each glacier by matching modeled and measured length changes. The volume changes predicted with the DESM agree well with measurements for 12 glaciers, whereas agreement is poor for one glacier with topographic breaks in the terminus area. For all glaciers, which are located in different climate regions, the length and volume changes are reproduced with the same ELA history. This agreement shows that the macroscopic glacier response to the climate history is well correlated over a whole mountain range. Modeling the future evolution of the glaciers under a constant present-day climate reveals that fast-reacting glaciers are close to equilibrium, whereas length and volume of the large valley glaciers would be reduced during the next century by an amount similar to the volume lost during the last 150 years.

[1] We used time-lapse imagery, seismic and audio recordings, iceberg and glacier velocities, ocean wave measurements, and simple theoretical considerations to investigate the interactions between Jakobshavn Isbrae and its proglacial ice melange. The melange behaves as a weak, granular ice shelf whose rheology varies seasonally. Sea ice growth in winter stiffens the melange matrix by binding iceberg clasts together, ultimately preventing the calving of full-glacier-thickness icebergs (the dominant style of calving) and enabling a several kilometer terminus advance. Each summer the melange weakens and the terminus retreats. The melange remains strong enough, however, to be largely unaffected by ocean currents (except during calving events) and to influence the timing and sequence of calving events. Furthermore, motion of the melange is highly episodic: between calving events, including the entire winter, it is pushed down fjord by the advancing terminus (at ∼40 m d−1), whereas during calving events it can move in excess of 50 × 103 m d−1 for more than 10 min. By influencing the timing of calving events, the melange contributes to the glacier's several kilometer seasonal advance and retreat; the associated geometric changes of the terminus area affect glacier flow. Furthermore, a force balance analysis shows that large-scale calving is only possible from a terminus that is near floatation, especially in the presence of a resistive ice melange. The net annual retreat of the glacier is therefore limited by its proximity to floatation, potentially providing a physical mechanism for a previously described near-floatation criterion for calving.

Objective Dental caries is one of the most common diseases in Western countries. Its pathoetiology is multifactorial, however, bacteria including Streptococcus mutans and the closely related Streptococcus sobrinus are regarded as key factors involved in this process. The fact that therapeutic approaches to eradicate these microorganisms are still limited prompted us to investigate the treatment potential of photodynamic therapy with the photoactive compounds hypericin (HYP) and meso-tetra(hydroxyphenyl)chlorin (mTHPC) in vitro. Material and methods S. mutans and S. sobrinus were cultivated under standard conditions and incubated with HYP (Invitrogen, Basel, Switzerland), the liposomal mTHPC derivative Foslipos (FOS, Biolitec, Jena, Germany), or a mixture of both at concentrations ranging between 0.625 and 10 μg/ml for various time points. Following a thorough washing step, bacteria were irradiated with a dental polymerization instrument (400–505 nm). All samples were subjected to serial dilutions and spiral plating on blood agar plates. Viable colony counts were determined after 48 h in culture. Photosensitizer fluorescence of bacteria was visualized by confocal microscopic techniques. Results One hundred percent of S. sobrinus could be killed by a 15 min incubation with as little as 2.5 μg/ml HYP, 5 μg/ml FOS or a mixture of 1.25 μg/ml of each photosensitizer followed by light activation of 120 s. In contrast to S. sobrinus, S. mutans displayed a significant dark toxicity for FOS (10–1.25 μg/ml) and no relevant PDT effects using HYP (10–0.625 μg/ml) under these conditions. HYP-mediated PDT effects (10 μg/ml) could be enhanced to more than 99.9% by prolonging photosensitizer incubation to 30 min and fractional illumination (2×120 s). Complete eradication of S. mutans was achieved by incubation for 15 min with a mixture of 0.625 μg/ml each of FOS and HYP and illumination for 120 s. Conclusion For both S. mutans and S. sobrinus, short PDT protocols with FOS and/or HYP could be established that completely eradicated these cariogenic bacteria in suspension. Our study, however, indicated that careful optimization of PDT conditions may be necessary for successful treatment of even closely related bacterial species. In multispecies microbial populations, the application of photosensitizer combinations for PDT may be useful.

[1] The recent loss of Jakobshavn Isbrae's extensive floating ice tongue has been accompanied by a change in near terminus behavior. Calving currently occurs primarily in summer from a grounded terminus, involves the detachment and overturning of several icebergs within 30–60 min, and produces long-lasting and far-reaching ocean waves and seismic signals, including “glacial earthquakes”. Calving also increases near-terminus glacier velocities by ∼3% but does not cause episodic rapid glacier slip, thereby contradicting the originally proposed glacial earthquake mechanism. We propose that the earthquakes are instead caused by icebergs scraping the fjord bottom during calving.

The supply, storage, and flow of water under an ice sheet are crucial for its overall dynamics. The density difference between ice and water limits subglacial water pressure to the pressure from overlying ice. Near this limit, the ice gradually becomes decoupled from the base (see the figure), water-filled cavities form, and sediments become weakened. The velocity of ice moving over the bedrock increases due to reduced friction and faster sediment deformation. Subglacial water pressure is therefore the key parameter that controls basal motion, ice flux, and the future evolution of an ice sheet. On page 777 of this issue ( 1 ), Meierbachtol et al. present a detailed record of water pressure variability under one region of the Greenland ice sheet.

[1] During summer 2001, a 140 m long ice core was recovered from the Belukha glacier (49°48′26″N, 86°34′43″E, 4062 m a.s.l.) in the Siberian Altai. The englacial temperature of −17°C indicates that this unique glaciochemical record is well preserved and suitable for the reconstruction of air pollution levels in this previously unexplored region. The upper 86 m were dated by 210Pb and cover the period 1862–2001. A lack of strong winter minima was observed in the δ18O record and attributed to the small amount of precipitation during that season and to wind erosion. The ion concentrations are comparable to those observed in Swiss glaciers, except for ammonium and formate, where enhanced concentrations indicate biogenic emissions from Siberian forests. Sulfate, ammonium and nitrate records all show anthropogenic impacts despite the remoteness of this site.

[1] A zone of highly irregular ice motion has been detected within the ice sheet adjacent to Jakobshavn Isbrae, a fast moving ice stream in West Greenland. Two tilt sensors in different boreholes and frozen into the ice show angle variations of up to 0.1° and 0.3° within a 1200 s interval. The amplitude of these flutter oscillations increases toward the transition from Holocene to Wisconsin ice at 682 m depth, 150 m above bedrock. All other tilt sensors above and below show low background noise. A numerical model confirms that the stress configuration is exceptional in the neighborhood of the ice stream. Stress transfer from the ice stream to the margins causes an extra forcing of the basal ice at the drill site. Several possible explanations for the observed flutter oscillations are given, the most likely being dynamic rupture events or stick-slip motion along the Holocene-Wisconsin transition. If true, this would correspond to an active overthrust fault with very high loading rates.

Sliding of the Greenland ice sheet is affected by the production of surface meltwater. A new theory shows that whether the result is a long-term speed-up or slow-down of ice motion depends on the variability in melt input. See Letter p.803 Recent observations in Greenland have led to the notion that surface meltwater supply to the ice-sheet bed lubricates ice flow, suggesting that climate warming could lead to runaway glacial acceleration. Christian Schoof uses a physically based model that captures drainage channelization under the ice to challenge this view. The model shows that increased melt supply leads to channelization and a drop in water pressure. This causes ice flow to slow down rather than speed up. However, this effect can be overcome if melt supply is variable over short timescales, when temporary pressure spikes can lead to accelerated flow. The positive melt/dynamic thinning feedback is therefore still viable, but is heavily dependent on daily temperature variations and rain events that are largely ignored in current ice sheet models.

Objective Dental caries is one of the most common diseases in Western countries. Its pathoetiology is multifactorial; however, the Gram-positive, facultatively anaerobic bacterium Streptococcus mutans may be one of the key factors involved in this process. The fact that therapeutic approaches to eradicate S. mutans are still limited prompted us to investigate the potential effects of photodynamic therapy with hypericin and Foslipos in vitro . Material and methods S. mutans (OMZ 176) was cultivated in fluid universal medium (FUM) under standard conditions. In total 100 μl bacteria (OD1 at 550 nm) were incubated with the photosensitizers hypericin (Invitrogen, Basel, Switzerland), the meso-tetrahydroxyphenylchlorin derivative Foslipos (Biolitec, Jena, Germany) or a mixture of both at concentrations ranging between 1.25 and 10 μg/ml for 1–90 min in FUM, followed by a thorough washing step in buffer. Bacteria were then irradiated for 60–120 s with a dental polymerization instrument (450 nm). All samples were subjected to serial dilutions and spiral plating on blood agar plates. Viable colony counts were determined after 48 h in culture. Adhesion and uptake of photosensitizers was analyzed by confocal and TIRF microscopy. Results Using the mix of photosensitizers (10 μg/ml total), a photosensitizer incubation time of 90 min and irradiation for 120 s, 100% of S. mutans could be killed. However, further optimization studies showed that irradiation time could be shortened to 60 s without loss of efficacy. In addition, the incubation time of the photosensitizers could be scaled down to 5 min, still resulting in 99.9999994% bacterial death. Comparisons of the concentration range indicated that a total photosensitizer concentration of as low as 2.5 μg/ml is sufficient for complete mortification of S. mutans . Under the same conditions, use of single photosensitizers instead of the mix appears to be less effective. Microscopic analyses confirmed the adherence of photosensitizers to bacteria. Conclusion Using a short protocol with a novel combination of low-dose photosensitizers, a complete eradication of S. mutans could be achieved in vitro . Future studies will include investigations on the efficacy of our protocol for other oral pathogens in vitro and in biofilms. Our data may set the stage for more effective treatment modalities not only within the oral cavity, but also for other systems prone to bacterial infections, such as the lower urinary tract.