Observations of glacier melt and runoff are of fundamental interest in the study of glaciers and their interactions with their environment. Considerable recent interest has developed around distributed acoustic sensing (DAS), a sensing technique which utilizes Rayleigh backscatter in fiber optic cables to measure the seismo-acoustic wavefield in high spatial and temporal resolution. Here, we present data from a month-long, 9 km DAS deployment extending through the ablation and accumulation zones on Rhonegletscher, Switzerland, during the 2020 melt season. While testing several types of machine learning (ML) models, we establish a regression problem, using the DAS data as the dependent variable, to infer the glacier discharge observed at a proglacial stream gauge. We also compare two predictive models that only depend on meteorological station data. We find that the seismo-acoustic wavefield recorded by DAS can be utilized to infer proglacial discharge. Models using DAS data outperform the two models trained on meteorological data with mean absolute errors of 0.64, 2.25 and 2.72 m3 s−1, respectively. This study demonstrates the ability of in situ glacier DAS to be used for quantifying proglacial discharge and points the way to a new approach to measuring glacier runoff.

Glacial hydraulic tremor (GHT) can be monitored to observe changes in location and distribution of water flow beneath glacial ice, allowing the spatiotemporal evolution of subglacial hydrology to be studied continuously and remotely. We use frequency-dependent polarization analysis (FDPA) to classify types of GHT and assess its spatio-temporal extent beneath Rhonegletscher, Switzerland, in a continuous seismic record through 2018 and 2019 at three ice-proximal bedrock-based seismometers. We determine the frequency bands composing the GHT and calculate back azimuth angles pointing to a previously known subglacial channel. We also inspect the relationship between GHT seismic power and water discharge from the glacier to compare daily and seasonal shifts with the observed GHT. We observed the seasonal shift from a distributed system to a channelized system, and our dataset allowed comparison of channel locations within and across seasons, with implications for sediment evacuation and bed erosion. The successful use of this method to assess GHT previously on Taku glacier (the methods of which this project follows) and now Rhonegletscher shows that existing ice-proximal passive seismic installations can be used to easily and continuously monitor subglacial hydraulic activity.

Retreat of calving glaciers worldwide has contributed substantially to sea-level rise in recent decades. Mass loss by calving contributes significantly to the uncertainty of sea-level rise projections. At Bowdoin Glacier, Northwest Greenland, most calving occurs by a few large events resulting from kilometre-scale fractures forming parallel to the calving front. High-resolution terrestrial radar interferometry data of such an event reveal that crevasse opening is fastest at low tide and accelerates during the final 36 h before calving. Using the ice flow model Elmer/Ice, we identify the crevasse water level as a key driver of modelled opening rates. Sea water-level variations in the range of local tidal amplitude (1 m) can reproduce observed opening rate fluctuations, provided crevasse water level is at least 4 m above the low-tide sea level. The accelerated opening rates within the final 36 h before calving can be modelled by additional meltwater input into the crevasse, enhanced ice cliff undercutting by submarine melt, ice damage increase due to tidal cyclic fatigue, crevasse deepening or a combination of these processes. Our results highlight the influence of surface meltwater and tides on crevasse opening leading to major calving events at grounded tidewater glaciers such as Bowdoin.

Luminous quasars are powered by accretion onto supermassive black holes. Such luminous quasars have been discovered up to the highest redshifts, z > 7. Here we discuss recent observations of the host galaxies of luminous quasars at z ≳ 6. We do not find a correlation between ongoing black hole growth and star-formation rate in the high redshift quasars, possibly indicating that black holes and their hosts do not co-evolve. We further show that even with high spatial resolution observations of the gas kinematics, dynamical mass estimates remain highly uncertain and should be used with caution.

Hydraulic processes within and beneath glacial bodies exert a far-reaching control on ice flow through their influence on basal sliding. Within the subglacial system, rapid changes in these processes may excite resonances whose interpretation requires an understanding of the underlying wave mechanics. Here, we explore these mechanics using observations from a kHz-sampled pressure sensor installed in a borehole directly above the hard granite bedrock of a temperate mountain glacier in Switzerland. We apply a previously established theory of wave propagation along thin, water-filled structures such as water-filled voids, basal water layers, or hydraulic fractures. Within such structures, short-wavelength waves experience restoring forces due to compressibility and are composed of sound waves. Long-wavelength resonances, in contrast, experience restoring forces due to elasticity and are composed of anomalously dispersed crack waves or Krauklis waves. Our borehole observations confirm the occurrence of both sound and crack waves within the basal water layer. Using both the resonance frequencies and attenuation of recorded crack waves we estimate thickness, aperture and length of the resonating basal water layer patch into which we drilled. We demonstrate that high-frequency observations of subglacial hydraulic processes provide new insights into this evolving dynamic system.

Since the 2000s, Greenland ice sheet mass loss has been accelerating, followed by increasing numbers of glacial earthquakes (GEs) at near-grounded glaciers. GEs are caused by calving of km-scale icebergs which capsize against the terminus. Seismic record inversion allows a reconstruction of the history of GE sources which captures capsize dynamics through iceberg-to-terminus contact. When compared with a catalog of contact forces from an iceberg capsize model, seismic force history accurately computes calving volumes while the earthquake magnitude fails to uniquely characterize iceberg size, giving errors up to 1 km3. Calving determined from GEs recorded ateight glaciers in 1993–2013 accounts for up to 21% of the associated discharge and 6% of the Greenland mass loss. The proportion of discharge attributed to capsizing calving may be underestimated by at least 10% as numerous events could not be identified by standard seismic detections (Olsen and Nettles, 2018). While calving production tends to stabilize in East Greenland, Western glaciers have released more and larger icebergs since 2010 and have become major contributors to Greenland dynamic discharge. Production of GEs and calving behavior are controlled by glacier geometry with bigger icebergs being produced when the terminus advances in deepening water. We illustrate how GEs can help in partitioning and monitoring Greenland mass loss and characterizing capsize dynamics.

Seasonal lake Gornersee forms at the confluence of Gornergletscher and Grenzgletscher, Switzerland, and experiences outburst floods annually in midsummer. To study the interplay between lake drainage, glacier movement and crevasse activity, high-frequency seismometers and GPS receivers were deployed in networks near Gornersee during the summer ablation seasons of 2004, 2006 and 2007. We use a Rayleigh wave coherence method to locate 3289, 7939 and 4087 icequakes, respectively, primarily along well-defined surface crevasses. We calculate two-dimensional strains from triads of GPS stations and find mean differential strain rates of ~300 × 10−6 d−1 with diurnal variations up to 800 × 10−6 d−1. Crevasse icequake activity and glacial velocity are highest during early season, then decrease as meltwater channels erode and subglacial water pressure decreases. Glacial response to Gornersee drainage varied year-to-year, with icequake activity promoted at some crevasses and inhibited at others, suggesting syn-drainage icequakes may be indicative of local drainage patterns and small-scale features of the stress field. Diurnal pulses in icequake activity exhibit peak activity at different times of day in different locations, coincident with a southeast-to-northwest trending concentrated shear zone near the Gornergletscher–Grenzgletscher confluence, likely due to differences in the timing of peak strain rate in these regions.

Passive seismology allows measurement of the structure of glaciers and ice sheets. However, most techniques used so far in this context are based on horizontally homogeneous media where parameters vary only with depth (1-D approximations), which are appropriate only for a subset of glaciers. Here, we analyze seismic noise records from three different types of glaciers (plateau, valley and avalanching glacier) to characterize the influence of the glacier geometry on the seismic wavefield. Using horizontal-to-vertical spectral ratios, polarization analysis and modal analysis, we show that the plateau glacier and the valley glacier can be seen as 1-D, whereas the relatively small avalanching glacier shows 3-D effects due to its bed topography and the deep crevasses. In principle, the techniques proposed here might allow monitoring such crevasses and their depth, and thus to constrain a key parameter of avalanching and calving glacier fronts.

Crevasses and englacial fracture networks route meltwater from a glacier's surface to the subglacial drainage system and thus influence glacial hydraulics. However, rapid fracture growth may also lead to sudden and potentially hazardous structural failure of unstable glaciers and ice dams, rifting of ice shelves, or iceberg calving. Here, we use passive seismic recordings to investigate the englacial fracture network on Glacier de la Plaine Morte, Switzerland. Glacier dynamics and the drainage of an ice-marginal lake give rise to numerous icequakes, the majority of which generate dispersed, high-frequency Rayleigh waves. A wide distribution of events allows us to study azimuthal anisotropy between 10 and 30 Hz in order to extract englacial seismic velocities in regions of preferentially oriented crevasses. Beamforming applied to a 100-m-aperture array reveals azimuthal anisotropy of Rayleigh-wave phase velocities reaching a strength of 8% at high frequencies. In addition, we find that the fast direction of wave propagation coincides with the observed surface strike of the narrow crevasses. Forward modeling and inversion of dispersion curves suggest that the azimuthal anisotropy is induced by a 40-m-thick crevassed layer at the surface of the glacier with 8% anisotropy in shear-wave velocity.

Using dense networks of three-component seismometers installed in direct contact with the ice, the seismic activity of Gornergletscher, Switzerland, was investigated during the summers of 2004 and 2006, as subglacial water pressures varied drastically. These pressure variations are due to the diurnal cycle of meltwater input as well as the subglacial drainage of Gornersee, a nearby marginal ice-dammed lake. Up to several thousand seismic signals per day were recorded. Whereas most icequakes are due to surface crevasse openings, about 200 events have been reliably located close to the glacier bed. These basal events tend to occur in clusters and have signals with impulsive first arrivals. At the same time, basal water pressures and ice-surface velocities were measured to capture the impact of the lake drainage on the subglacial hydrological system and the ice-flow dynamics. Contrary to our expectations, we did not observe an increase of basal icequake activity as the lake emptied, thereby raising the subglacial water pressures close to the flotation level for several days. In fact, the basal icequakes were usually recorded during the morning hours, when the basal water pressure was either low or decreasing. During the high-pressure period caused by the drainage of the lake, no basal icequakes were observed. Furthermore, GPS measurements showed that the glacier surface was lowering during the basal seismic activity. These observations lead us to conclude that such icequakes are connected to the diurnal variation in glacier sliding across the glacier bed.

The tongue of Triftgletscher, Switzerland, is particularly susceptible to major break-off events due to its steep slope. In order to detect precursors of such an event, we monitored the local seismic activity and detected 2426 icequakes with sources located in an area ranging between 2050 and 2350ma.s.l. Events triggered by cracks and icefalls were recorded, but no precise distinction between the two sources was possible using duration or frequency criteria. Clusters of seismic activity were located and confirmed by visual observations. We performed a surface motion analysis and found that surface motion was driven significantly by runoff changes at a timescale of 2-3 days. By means of a statistical analysis, a power-law behaviour of the released seismic energy distribution was discerned at certain times during the period investigated. Variations in power-law exponent values indicated that low- and high-energy events predominantly occurred during phases of enhanced and reduced surface motion, respectively. Substantial releases of seismic energy likely to signal the glacier recoupling were detected during phases of decreasing runoff. Clues to potential seismic precursors of break-off events are discussed.

With mass loss from the Greenland ice sheet accelerating and spreading to higher latitudes, the quantification of mass discharge in the form of icebergs has recently received much scientific attention. Here we make use of very low-frequency (0.001–0.01 Hz) seismic data from three permanent broadband stations installed in the summers of 2009–10 in northwest Greenland in order to monitor local calving activity. At these frequencies, calving seismograms are dominated by a tilt signal produced by local ground flexure in response to fjord seiching generated by major iceberg calving events. A simple triggering algorithm is proposed to detect calving events from large calving fronts with potentially no user interaction. Our calving catalogue identifies spatial and temporal differences in calving activity between Jakobshavn Isbræ and glaciers in the Uummannaq district ∼200 km to the north. The Uummannaq glaciers show clear seasonal fluctuations in seiche-based calving detections as well as seiche amplitudes. In contrast, the detections at Jakobshavn Isbræ show little seasonal variation, which may be evidence for an ongoing transition to winter calving activity. The results offer further evidence that seismometers can provide efficient and inexpensive monitoring of calving fronts.

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.

We infer the absolute time dependence of kinematic gas temperature along a proposed orbit of molecular clouds in the Central Molecular Zone (CMZ) of the Galactic Center (GC). Ammonia gas temperature maps are one of the results of the “Survey of Water and Ammonia in the Galactic Center” (SWAG, PI: J. Ott); the dynamical model of molecular clouds in the CMZ was taken from Kruijssen et al. (2015). We find that gas temperatures increase as a function of time in both regimes before and after the cloud passes pericenter on its orbit in the GC potential. This is consistent with the recent proposal that pericenter passage triggers gravitational collapse. Other investigated quantities (line width, column density, opacity) show no strong sign of time dependence but are likely dominated by cloud-to-cloud variations.

Fractured systems developing at the base of temperate glaciers play a significant role in the transport of subglacial water. Their influence on dynamics and hydraulics of glaciers and the seasonal alterations they experience during the hydrological year, remain however, poorly understood, as they lay beyond the reach of most conventional glaciological techniques. In this work, we investigate the seasonal evolution of glacier seismicity in the ablation zone of Rhonegletscher (Switzerland) and focus on basal fracture processes. To do this, we use seismic measurements performed between June 2012 and July 2013 and complementary measurements of surface motion and changes in melt water input inferred from degree day modeling. Results of a spectrogram analysis highlight predominant seismic emissions until end of October 2012 and from mid April 2013 in line with measured surface velocities, and indicate a rapid adaptation of the basal drainage system at initiation and end of the winter period. The analysis of icequakes emitted from two active regions shows that the basal events, due to opening and closing of cracks within fracture networks, persist after winter initiation, and points out the year-long preservation of the associated sources. In this regard, we deduce that basal crack networks aseismically evolve during the winter period.

We present the first results from SWAN: “Survey of Water and Ammonia in Nearby galaxies”. Nearby galaxies are conveniently located to probe molecular gas properties on scales of 10 to 200 pc, which are appropriate for the study of Giant Molecular Clouds (GMCs). The resolution of the Very Large Array in D and C configurations corresponds to a few 10s of parsecs in these galaxies. To advance studies of galaxy evolution it is paramount to understand how processes in the molecular Interstellar Medium(ISM) and star formation are linked on these scales. We have observed the metastable transitions of ammonia and the 22GHz water maser line in four nearby galaxies: NGC 253, IC 342, NGC 2146, and NGC 6946 using the VLA. These galaxies were chosen to span an order of magnitude in star formation rate, and a range of galactic ecosystems. We use the ammonia transitions to derive kinetic temperatures, which exposes the heating and cooling balance of the ISM. We then aim to relate these conditions to energetic feedback from star formation as indicated by water masers.

Currently, our analysis is focused on NGC 253. NGC 253 is a barred spiral starburst galaxy with a nucleated star formation rate (SFR) of ≈ 3M⊙ per year. We use a distance of 3.9Mpc for analysis. We have observed ammonia transitions (1,1) to (5,5) and the 22GHz water maser line with a resolution of ≈63pc. We have identified nine regions across the nucleated starburst for study. The ammonia (3,3) line appears to be masing in the centermost 200pc. We have identified two regions of water maser emission. The first region is a minor axis extension, about the center of the galaxy, and is very close to the outflow. By means of these measurements, we gain an understanding of the molecular ISM associated with the nucleated starburst environment in NGC 253.

Our understanding of galaxy evolution has traditionally been driven by pre-selection of galaxies based on their broad-band continuum emission. This approach is potentially biased, in particular against gas-rich systems at high-redshift which may be dust-obscured. To overcome this limitation, we have recently concluded a blind CO survey at 3mm in a region of the Hubble Deep Field North using the IRAM Plateau de Bure Interferometer. Our study resulted in 1) the discovery of the redshift of the bright SMG HDF850.1 (z = 5.183); 2) the discovery of a bright line identified as CO(2-1) arising from a BzK galaxy at z = 1.785, and of other 6 CO lines associated with various galaxies in the field; 3) the detection of a few lines (presumably CO(3-2) at z ∼ 2) with no optical/NIR/MIR counterparts. These observational results allowed us to expand the parameter space of galaxy properties probed so far in high-z molecular gas studies. Most importantly, we could set first direct constraints on the cosmic evolution of the molecular gas content of the universe. The present study represents a first, fundamental step towards an unbiased census of molecular gas in ‘normal’ galaxies at high-z, a crucial goal of extragalactic astronomy in the ALMA era.

High-redshift quasars are thought to reside in the most massive halos in the early universe and should therefore be located in fields with overdensities of galaxies, which are expected to evolve into galaxy clusters seen in the local Universe. In Bañados et al. (2013), we used deep narrow-band imaging to study the environment of the z=5.72 quasar ULAS J0203 + 0012. The redshift range probed by our narrow band selection is Δz ~ 0.1. This was the first time that Lyman alpha emitters (LAEs) were searched for near a z ~ 6 quasar, to provide clues on the surroundings of quasars at the end of the epoch of reionization. The main result of this work is that no enhancement of LAEs has been found in the surroundings of ULAS J0203 + 0012. We discuss different explanations and interpretations for this non-detection of a galaxy overdensity.