Search

The influence of surface melt on the flow of Greenland's largest outlet glaciers remains poorly known and in situ observations are few. We use field observations to link surface meltwater forcing to glacier-wide diurnal velocity variations on East Greenland's Helheim Glacier over two summer melt seasons. We observe diurnal variations in glacier speed that peak ~6.5 h after daily maximum insolation and extend from the terminus region to the equilibrium line. Both the amplitude of the diurnal speed variation and its sensitivity to daily melt are largest at the glacier terminus and decrease up-glacier, suggesting that the magnitude of the response is controlled not only by melt input volume and temporal variability, but also by background effective pressure, which approaches zero at the terminus. Our results provide evidence that basal lubrication by meltwater drives diurnal velocity variations at Greenland's marine-terminating glaciers in a similar manner to alpine glaciers and Greenland's land-terminating outlet glaciers.

12 pages, 13 figures, The Helheim glacier, located in southeast Greenland, has more than ten campaign-type Global Positioning System (GPS) sites; data processing led to the observation of a very rapid change in the ionospheric delay. To identify the cause of these sporadic disturbances, we analyzed the slant total electron content (STEC), single-differenced STEC (SD-STEC) and scintillation proxy index called the delta phase rate (DPR). From this analysis, the abrupt change of those ionospheric indicators was attributed to the line-of-sight direction to the satellite and the temporal sequence of the event was found to be highly correlated with the geometry of the GPS sites. In addition, the disturbance based on the result of SD-STEC occurred mostly during the night, from 17 UTC through 7 UTC, and across a band spanning the east-west direction. Based on the DPR indices obtained from GPS stations distributed across all of Greenland, Iceland, and northeastern Canada, the rapid ionospheric variation was found to be correlated with the time of the day and the geomagnetic latitude of the station. The disturbance was larger at the relatively low geomagnetic latitudes at night but was more significant at higher latitudes in the daytime. These rapid ionospheric variations tended to appear in band shapes parallel to the geomagnetic field. These results allow us to attribute such disturbance observed at the Helheim glacier to aurora-related phenomena, This work was partly supported by the basic research fund from Korea Astronomy and Space Science Institute (2020-1-8-5005). This work was also partly supported by INHA UNIVERSITY Research Grant (INHA-61597), With the funding support of the ‘Severo Ochoa Centre of Excellence’ accreditation (CEX2019-000928-S), of the Spanish Research Agency (AEI)

Glacial earthquakes are caused by large iceberg calving events, which are an important mechanism for mass loss from the Greenland ice sheet. The number of glacial earthquakes in Greenland has increased sixfold over the past two decades. We use teleseismic surface waves to analyze the 145 glacial earthquakes that occurred in Greenland from 2011 through 2013, and successfully determine source parameters for 139 events at 13 marine-terminating glaciers. Our analysis increases the number of events in the glacial-earthquake catalog by nearly 50% and extends it to 21 years. The period 2011–13 was the most prolific 3-year period of glacial earthquakes on record, with most of the increase over earlier years occurring at glaciers on Greenland's west coast. We investigate changes in earthquake productivity and geometry at several individual glaciers and link patterns in glacial-earthquake production and cessation to the absence or presence of a floating ice tongue. We attribute changes in earthquake force orientations to changes in calving-front geometry, some of which occur on timescales of days to months. Our results illustrate the utility of glacial earthquakes as a remote-sensing tool to identify the type of calving event, the grounded state of a glacier, and the orientation of an active calving front.

Glacial earthquakes are slow earthquakes of magnitude M~5 associated with major calving events at near-grounded marine-terminating glaciers. These globally detectable earthquakes provide information on the grounding state of outlet glaciers and the timing of large calving events. Seismic source modeling of glacial earthquakes provides information on the size and orientation of forces associated with calving events. We compare force orientations estimated using a centroid-single-force technique with the calving-front orientations of the source glaciers at or near the time of earthquake occurrence. We consider earthquakes recorded at four glaciers in Greenland – Kangerdlugssuaq Glacier, Helheim Glacier, Kong Oscar Glacier, and Jakobshavn Isbræ – between 1999 and 2010. We find that the estimated earthquake force orientations accurately represent the orientation of the calving front at the time of the earthquake, and that seismogenic calving events are produced by a preferred section of the calving front, which may change with time. We also find that estimated earthquake locations vary in a manner consistent with changes in calving-front position, though with large scatter. We conclude that changes in glacial-earthquake source parameters reflect true changes in the geometry of the source glaciers, providing a means for identifying changes in glacier geometry and dynamics that complements traditional remote-sensing techniques.

We analyze the first ever GPS observations of static surface deformation from a deep earthquake: the 24 May 2013 Mw 8.3 Sea of Okhotsk, 611 km-deep, event. Previous studies of deep earthquake sources relied on seismology and might have missed evidence for slow slip in the rupture. We observed coseismic static offsets on a GPS network of 20 stations over the Sea of Okhotsk region. The offsets were inverted for the best fitting double-couple source model assuming a layered spherical Earth. The seismic moment calculated from static offsets is only 7% larger than the seismological estimate from Global Centroid Moment Tensor (GCMT). Thus, GPS observations confirm shear faulting as the source model, with no significant slow-slip component. The relative locations of the U.S. Geological Survey hypocenter, GCMT centroid, and the fault from GPS indicate slip extending for tens of kilometers across most of the slab thickness.

Using a new biaxial friction apparatus, we conducted experiments of ice-on-rock friction in order to better understand basal sliding of glaciers and ice streams. A series of velocity-stepping and slide–hold–slide tests were conducted to measure friction and healing at temperatures between −20°C and melting. Experimental conditions in this study are comparable to subglacial temperatures, sliding rates and effective pressures of Antarctic ice streams and other glaciers, with load-point velocities ranging from 0.5 to 100 µm s −1 and normal stress σ n = 100 kPa. In this range of conditions, temperature dependences of both steady-state friction and frictional healing are considerable. The friction increases linearly with decreasing temperature (temperature weakening) from μ = 0.52 at −20°C to μ = 0.02 at melting. Frictional healing increases and velocity dependence shifts from velocity-strengthening to velocity-weakening behaviour with decreasing temperature. Our results indicate that the strength and stability of glaciers and ice streams may change considerably over the range of temperatures typically found at the ice–bed interface. This article is part of the themed issue ‘Microdynamics of ice’.

SUMMARY On 2010 November 14, an intense swarm of earthquakes began in the western Gulf of Aden. Within a 48-hr period, 82 earthquakes with magnitudes between 4.5 and 5.5 were reported along an ∼80-km-long segment of the east–west trending Aden Ridge, making this swarm one of the largest ever observed in an extensional oceanic setting. In this study, we calculate centroid-moment-tensor solutions for 110 earthquakes that occurred between 2010 November and 2011 April. Over 80 per cent of the cumulative seismic moment results from earthquakes that occurred within 1 week of the onset of the swarm. We find that this sequence has a b-value of ∼1.6 and is dominated by normal-faulting earthquakes that, early in the swarm, migrate westwards with time. These earthquakes are located in rhombic basins along a section of the ridge that was previously characterized by low levels of seismicity and a lack of recent volcanism on the seafloor. Body-wave modelling demonstrates that the events occur in the top 2–3 km of the crust. Nodal planes of the normal-faulting earthquakes are consistent with previously mapped faults in the axial valley. A small number of strike-slip earthquakes observed between two basins near 44°E, where the axial valley changes orientation, depth and width, likely indicate the presence of an incipient transform fault and the early stages of ridge-transform segmentation. The direction of extension accommodated by the earthquakes is intermediate between the rift orthogonal and the direction of relative motion between the Arabian and Somalian plates, consistent with the oblique style of rifting occurring along the slow-spreading Aden Ridge. The 2010 swarm shares many characteristics with dyke-induced rifting episodes from both oceanic and continental settings. We conclude that the 2010 swarm represents the seismic component of an undersea magmatic rifting episode along the nascent Aden Ridge, and attribute the large size of the earthquakes to the combined effects of the slow spreading rate, relatively thick crust and recent quiescence. We estimate that the rifting episode was caused by dyke intrusions that propagated laterally for 12–18 hr, accommodating ∼1–14 m of opening or ∼85–800 yr of spreading along this section of the ridge. Our findings demonstrate the westward propagation of active seafloor spreading into this section of the western Gulf of Aden and illustrate that deformation at the onset of seafloor spreading may be accommodated by discrete episodes of faulting and magmatism. A comparison with similar sequences on land suggests that the 2010 episode may be only the first of several dyke-induced rifting episodes to occur in the western Gulf of Aden.

12 pages, 10 figures., The flow speed of Greenland outlet glaciers is governed by several factors, the relative importance of which is poorly understood. The delivery of surface-generated meltwater to the bed of alpine glaciers has been shown to influence glacier flow speed when the volume of water is sufficient to increase basal fluid pressure and hence basal lubrication. While this effect has also been demonstrated on the Greenland ice-sheet margin, little is known about the influence of surface melting on the large, marine-terminating outlet glaciers that drain the ice sheet. We use a validated model of meltwater input and GPS-derived surface velocities to quantify the sensitivity of glacier flow speed to changes in surface melt at Helheim Glacier during two summer seasons (2007-08). Our observations span ∼55 days near the middle of each melt season. We find that relative changes in glacier speed due to meltwater input are small, with variations of ∼45% in melt producing changes in velocity of ∼2-4%. These velocity variations are, however, of similar absolute magnitude to those observed at smaller glaciers and on the ice-sheet margin. We find that the glacier's sensitivity to variations in meltwater input decreases approximately exponentially with distance from the calving front. Sensitivity to melt varies with time, but generally increases as the melt season progresses. We interpret the time-varying sensitivity of glacier flow to meltwater input as resulting from changes in subglacial hydraulic routing caused by the changing volume of meltwater input., The Helheim 2007 and 2008 project was supported by the Gary Comer Science and Education Foundation, the US National Science Foundation, the Danish Commission for Scientific Research in Greenland (KVUG), the Spanish Ministry of Science and Innovation (MICINN), the Geological Survey of Denmark and Greenland (GEUS), Geocenter Copenhagen, the Danish National Space Center, NASA, the Lamont–Doherty Climate Center, and the Dan and Betty Churchill Exploration Fund. GPS equipment and technical support were provided by UNAVCO, Inc.

SUMMARY Previous analyses of geological and geodetic data suggest that the obliquely compressive relative motion across the Caribbean–North America plate boundary in Hispaniola is accommodated through strain partitioning between near-vertical transcurrent faults on land and low-angle thrust faults offshore. In the Dominican Republic, earthquake focal-mechanism geometries generally support this interpretation. Little information has been available about patterns of seismic strain release in Haiti, however, due to the small numbers of moderate-to-large earthquakes occurring in western Hispaniola during the modern instrumental era. Here, we analyse the damaging MW= 7.0 earthquake that occurred near Port au Prince on 2010 January 12 and aftershocks occurring in the four months following this event, to obtain centroid–moment-tensor (CMT) solutions for 50 earthquakes with magnitudes as small as MW= 4.0. While the 2010 January main shock exhibited primarily strike-slip motion on a steeply dipping nodal plane (strike=250°, dip=71° and rake=22°), we find that nearly all of the aftershocks show reverse-faulting motion, typically on high-angle (30°–45°) nodal planes. Two small aftershocks (MW 4.5 and 4.6), located very close to the main shock epicentre, show strike-slip faulting with geometries similar to the main shock. One aftershock located off the south coast of Haiti shows low-angle thrust faulting. We also examine earthquakes occurring in this region from 1977–2009; successful analysis of four such events provides evidence for both strike-slip and reverse faulting. The pattern of seismic strain release in southern Haiti thus indicates that partitioning of plate motion between transcurrent and reverse structures extends far west within Hispaniola. While we see limited evidence for low-angle underthrusting offshore, most reverse motion appears to occur on high-angle fault structures adjacent to the Enriquillo fault. Our results highlight the need to incorporate seismogenic slip on compressional structures into hazard assessments for southern Haiti.

Glacial earthquakes are a new class of seismic events, first discovered as signals in long-period seismograms recorded on the Global Seismographic Network. Most of these events occur along the coasts of Greenland, where they are spatially related to large outlet glaciers. Glacial earthquakes show a strong seasonality, with most events occurring during the late summer. The rate of glacial-earthquake occurrence increased between 2000 and 2005, with a stabilization of earthquake frequency at 2003–2004 levels in 2006–2008. Recent observations establish a strong temporal correlation between the distinct seismic signals of glacial earthquakes and large ice-loss events in which icebergs of cubic-kilometer scale collapse against the calving face, linking the seismogenic process to the force exerted by these icebergs on the glacier and the underlying solid earth. A sudden change in glacier speed results from these glacial-earthquake calving events. Seasonal and interannual variations in glacier-terminus position account for general characteristics of the temporal variation in earthquake occurrence. Glacial earthquakes in Antarctica are less well studied, but they exhibit several characteristics similar to glacial earthquakes in Greenland.

Movers and shakers When the edge of an ice sheet breaks off and falls into the sea (calves), the remaining section of the ice sheet moves backward and down and can suffer a glacial earthquake. Murray et al. studied calving from Greenland's Helheim Glacier. The forces that cause the change in the motion of the ice sheet at its terminus also trigger the accompanying earthquakes. Because these seismic signals can be detected by instruments located all over the globe, it should be possible to use these glacial earthquakes as proxies for glacier calving. Science , this issue p. 305

12 pages, 9 figures, 1 table, 1 appendix, We have used tide-gauge data from near Helheim Glacier, East Greenland, and GPS data acquired on the glacier to investigate the spectra of tidal forcing and flow response. For both the tidegauge and GPS time series, we calculated amplitudes and phases for a harmonic expansion using a limited set of harmonic constituents. We find that the semidiurnal constituents of the glacier flow are well modeled using a single admittance and lag with respect to the tide-gauge data. However, diurnal variations in the glacier flow cannot be simply described using this model. We find an additional signal in glacier position, in phase with the S1 solar tide, with some modulations at other frequencies. These non-tidal variations account for a peak-to-peak variation in glacier flow speed at a site close to the terminus of ~0.7md-1, compared with a mean flow rate at this location of ~22md-1. The speed variations reach their daily maximum value ~6 hours after local noon. We hypothesize that these additional diurnal variations are associated with peaks in lubrication of the glacier bed due to surface melting driven by diurnal solar heating, We thank members of the Helheim Project for collecting the datasets used in this analysis. This work was supported by US National Science Foundation grants ARC-0713970, ARC-07131749 and ARC-1110322, the Gary Comer Science and Education Foundation, the Spanish Ministry of Science and Innovation (MICINN; now MINECO), the Danish Commission for Scientific Research in Greenland (KVUG), the Geological Survey of Denmark and Greenland (GEUS), Geocenter Copenhagen, the Danish National Space Center and NASA. GPS equipment and technical support were provided by UNAVCO, Inc.

The March 25, 1998, Antarctic plate earthquake (Mw=8.1) occurred ∼250 km from the nearest plate boundary, in oceanic lithosphere with an age of 35–55 my. Analysis of aftershock patterns shows that the earthquake ruptured a fault, or series of strike-slip fault segments, nearly 300 km long. The strike of the fault(s) is nearly perpendicular to the north-south trending fossil fracture zones which are the most marked bathymetric features of this region. Moment release during the mainshock was concentrated in two large subevents, clearly visible in the teleseismic body wave waveforms. Modeling of these body waves using a finite fault source places the first of the two subevents near the point of rupture initiation, on the eastern end of the fault, with the second large subevent occurring 220–280 km to the west. The two pulses of moment release are found to be separated in time by ∼65 s. Comparison of the relative S wave amplitudes of the first and second pulses suggests that a rotation of the focal mechanism by ∼10° occurred between the two subevents.

The source characteristics of 10 shallow earthquakes of moderate size (M>5) beneath Bardarbunga Volcano in Iceland are investigated using long-period and broadband teleseismic data. The highly non-double-couple nature of the focal mechanisms for several of the events previously reported in the Harvard centroid moment tensor (CMT) catalog is confirmed by analysis of an expanded long-period body wave data set, as well as by an independent inversion of intermediate period surface wave data. Modeling of broadband P waveforms constrains the depths of the events to lie between 3 km and 8 km beneath the volcano. Significant complexity in the P waves observed for some events is best explained by a sequence of subevents of varying geometry, rather than by wave propagation effects caused by shallow structure near the earthquake. Inversion of P wave data for two of the earthquakes using a parameterization with subsources results in sequences of events with primarily thrust motion on planes of varying strike. The seismological results are interpreted as faulting on an outward dipping cone-shaped ring fault beneath the Bardarbunga caldera. The association in time and space between the September 29, 1996, earthquake and the September 30 eruption near Bardarbunga and Grimsvotn Volcanos suggests that inflation of a shallow magma chamber, and the associated stress loading of the deeper ring fault, may be the tectonic process which generates the earthquakes.

Some glaciers and ice streams periodically lurch forward with sufficient force to generate emissions of elastic waves that are recorded on seismometers worldwide. Such glacial earthquakes on Greenland show a strong seasonality as well as a doubling of their rate of occurrence over the past 5 years. These temporal patterns suggest a link to the hydrological cycle and are indicative of a dynamic glacial response to changing climate conditions.

[1] Glacial earthquakes are anomalous earthquakes associated with large ice-loss events occurring at marine-terminating glaciers, primarily in Greenland. They are detectable teleseismically, and a proper understanding of the source mechanism may provide a remote-sensing tool to complement glaciological observations of these large outlet glaciers. We model teleseismic surface-wave waveforms to obtain locations and centroid–single-force source parameters for 121 glacial earthquakes occurring in Greenland during the period 2006–2010. We combine these results with those obtained by previous workers to analyze spatial and temporal trends in glacial-earthquake occurrence over the 18-year period from 1993–2010. We also examine earthquake occurrence at six individual glaciers, comparing the earthquake record to independently obtained observations of glacier change. Our findings confirm the inference that glacial-earthquake seismogenesis occurs through the capsize of large, newly calved icebergs. We find a close correspondence between episodes of glacier retreat, thinning, and acceleration and the timing of glacial earthquakes, and document the northward progression of glacial earthquakes on Greenland's west coast over the 18-year observing period. Our results also show that glacial earthquakes occur when the termini of the source glaciers are very close to the glacier grounding line, i.e., when the glaciers are grounded or nearly grounded.

Some of the most dramatic effects of climate change have been observed in the Earth's polar regions. In Greenland, ice loss from the Greenland ice sheet has accelerated in recent years [Shepherd et al ., 2012]. Outlet glaciers are changing their behavior rapidly, with many thinning, retreating, and accelerating [Joughin et al ., 2004]. The loss of ice weighing on the crust and mantle below has allowed both to rebound, resulting in high rock uplift rates [Bevis et al ., 2012]. Changes in ice cover and meltwater production influence sea level and climate feedbacks; they are expected to contribute to increasing vulnerability to geohazards such as landslides, flooding, and extreme weather. ISSN:0096-3941 ISSN:2324-9250

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.

[1] Large calving events at Greenland’s largest outlet glaciers are associated with glacial earthquakes and near‐ instantaneous increases in glacier flow speed. At some glaciers and ice streams, flow is also modulated in a regular way by ocean tidal forcing at the terminus. At Helheim Glacier, analysis of geodetic data shows decimeter‐level periodic position variations in response to tidal forcing. However, we also observe transient increases of more than 100% in the glacier’s responsiveness to such tidal forcing following glacial‐earthquake calving events. The timing and amplitude of the changes correlate strongly with the step‐like increases in glacier speed and longitudinal strain rate associated with glacial earthquakes. The enhanced response to the ocean tides may be explained by a temporary disruption of the subglacial drainage system and a concomitant reduction of the friction at the ice‐bedrock interface, and suggests a new means by which geodetic data may be used to infer glacier properties. Citation: de Juan, J., et al. (2010), Sudden increase in tidal response linked to calving and acceleration at a large Greenland outlet glacier, Geophys. Res. Lett., 37, L12501

We used satellite images to examine the calving behavior of Helheim and Kangerdlugssuaq Glaciers, Greenland, from 2001 to 2006, a period in which they retreated and sped up. These data show that many large iceberg-calving episodes coincided with teleseismically detected glacial earthquakes, suggesting that calving-related processes are the source of the seismicity. For each of several events for which we have observations, the ice front calved back to a large, pre-existing rift. These rifts form where the ice has thinned to near flotation as the ice front retreats down the back side of a bathymetric high, which agrees well with earlier theoretical predictions. In addition to the recent retreat in a period of higher temperatures, analysis of several images shows that Helheim retreated in the 20th Century during a warmer period and then re-advanced during a subsequent cooler period. This apparent sensitivity to warming suggests that higher temperatures may promote an initial retreat off a bathymetric high that is then sustained by tidewater dynamics as the ice front retreats into deeper water. The cycle of frontal advance and retreat in less than a century indicates that tidewater glaciers in Greenland can advance rapidly. Greenland's larger reservoir of inland ice and conditions that favor the formation of ice shelves likely contribute to the rapid rates of advance.

[1] High precision relocation of earthquakes recorded by the Hawaiian Volcano Observatory (HVO) seismic network provides new information on the characteristics of seismic faulting at this oceanic hot spot. Using waveform cross correlation, we have measured correlation coefficients and travel time differences for a set of 14,605 deep (≥13 km) earthquakes recorded from 1988 to 1998. We find that about half of the analyzed earthquakes are in similar event clusters that delineate fault zones in the lower crust and upper mantle. We suggest that much of this deep seismicity reflects rupture in the brittle lithosphere away from the magma pathways, although at Kilauea the stresses from magma movement may additionally help trigger mantle earthquakes on preexisting faults in regions with high differential ambient stresses. Focal mechanisms of similar event clusters throughout Hawaii display characteristic patterns and appear consistent with the hypothesis that deep earthquakes on preexisting faults reflect the stresses due to volcano loading and flexure. We also present the results of applying cross correlation analyses and relocation to ∼7000 earthquakes at all depths located west of 155.55°W and recorded from 1988 to 1998. The pattern of relocated earthquakes at the Kealakekua fault zone is consistent with the presence of a low-angle detachment on the west flank of Mauna Loa.