Glaciers in Svalbard: mass balance, runoff and freshwater flux

Gain or loss of the freshwater stored in Svalbard glaciers has both global implications for sea level and, on a more local scale, impacts upon the hydrology of rivers and the freshwater flux to fjords. This paper gives an overview of the potential runoff from the Svalbard glaciers. The freshwater flux from basins of different scales is quantified. In small basins (A < 10 km²), the extra runoff due to the negative mass balance of the glaciers is related to the proportion of glacier cover and can at present yield more than 20% higher runoff than if the glaciers were in equilibrium with the present climate. This does not apply generally to the ice masses of Svalbard, which are mostly much closer to being in balance. The total surface runoff from Svalbard glaciers due to melting of snow and ice is roughly 25 ± 5 km³ a⁻¹, which corresponds to a specific runoff of 680 ± 140 mm a⁻¹, only slightly more than the annual snow accumulation. Calving of icebergs from Svalbard glaciers currently contributes significantly to the freshwater flux and is estimated to be 4 ± 1 km³ a⁻¹ or about 110 mm a⁻¹.     

Temporal and spatial variations of post-seismic deformation following the 1999 Chi-Chi, Taiwan earthquake

We use GPS displacements collected in the 15 months after the 1999 Chi-Chi, Taiwan earthquake  ( M _w 7.6)  to evaluate whether post-seismic deformation is better explained by afterslip or viscoelastic relaxation of the lower crust and upper mantle. We find that all viscoelastic models tested fail to fit the general features in the post-seismic GPS displacements, in contrast to the satisfactory fit obtained with afterslip models. We conclude that afterslip is the dominant mechanism in the 15-month period, and invert for the space–time distribution of afterslip, using the Extended Network Inversion Filter. Our results show high slip rates surrounding the region of greatest coseismic slip. The slip-rate distribution remains roughly stationary over the 15-month period. In contrast to the limited coseismic slip on the décollement, afterslip is prominent there. Maximum afterslip of 0.57 m occurs downdip and to the east of the hypocentral region. Afterslip at hypocentral depths is limited to the southern part of the main shock rupture, with little or no slip on the northern section where coseismic slip was greatest. Whether this results from along strike variations in frictional properties or dynamic conditions that locally favour stable sliding is not clear. In general, afterslip surrounds the area of greatest coseismic slip, consistent with post-seismic slip driven by the main shock stress change. The total accumulated geodetic afterslip moment is  3.8 × 10¹⁹ N m  , significantly more than the seismic moment released by aftershocks,  6.6 × 10¹⁸ N m  . Afterslip and aftershocks appear to have different temporal evolutions and some spatial correlations, suggesting that aftershock rates may not be completely controlled by the rate of afterslip.     

Strain accumulation along an oblique plate boundary: the Reykjanes Peninsula, southwest Iceland

We use annual GPS observations on the Reykjanes Peninsula (RP) from 2000 to 2006 to generate maps of surface velocities and strain rates across the active plate boundary. We find that the surface deformation on the RP is consistent with oblique plate boundary motion on a regional scale, although considerable temporal and spatial strain rate variations are observed within the plate boundary zone. A small, but consistent increase in eastward velocity is observed at several stations on the southern part of the peninsula, compared to the 1993–1998 time period. The 2000–2006 velocities can be modelled by approximating the plate boundary as a series of vertical dislocations with left-lateral motion and opening. For the RP plate boundary we estimate left-lateral motion  18⁺⁴₋₃ mm yr⁻¹  and opening of  7⁺³₋₂ mm yr⁻¹  below a locking depth of  7⁺¹₋₂ km  . The resulting deep motion of  20⁺⁴₋₃ mm yr⁻¹  in the direction of  N(100⁺⁸₋₆)°E  agrees well with the predicted relative North America–Eurasia rate. We calculate the areal and shear strain rates using velocities from two periods: 1993–1998 and 2000–2006. The deep motion along the plate boundary results in left-lateral shear strain rates, which are perturbed by shallow deformation due to the 1994–1998 inflation and elevated seismicity in the Hengill–Hrómundartindur volcanic system, geothermal fluid extraction at the Svartsengi power plant, and possibly earthquake activity on the central part of the peninsula.     

On the consistency of earthquake moment release and space geodetic strain rates: Europe

In this paper, approximately 100 VLBI/SLR/GPS velocities map European strain rates from <0.09 × 10⁻⁸ to >9.0 × 10⁻⁸ yr⁻¹ with regional uncertainties of 20 to 40 per cent. Kostrov's formula translates these strain-rate values into regional geodetic moment rates M¯˙ _geodetic . Two other moment rates, M¯˙ _seismic , extracted from a 100-year historical catalogue and M¯˙ _plate , taken from plate-tectonic models, contrast the geodetic rates. In Mediterranean Europe, the ratios of M¯˙ _seismic to M¯˙ _geodetic are between 0.50 and 0.71. In Turkey the ratio falls to 0.22. Although aseismic deformation may contribute to the earthquake deficit ( M¯˙ _seismic values less than M¯˙ _geodetic ), the evidence is not compelling because the magnitudes of the observed shortfalls coincide with the random variations expected in a 100-year catalogue. If the lack of aseismic deformation inferred from the 100-year catalogue holds true for longer periods, then much of Europe's strain budget would have to be accommodated by more frequent or larger earthquakes than have been experienced this century to raise the ratios of M¯˙ _seismic to M¯˙ _geodetic to unity. Improved geological fault data bases, longer historical earthquake catalogues, and densification of the continent's space geodetic network will clarify the roles of aseismic deformation versus statistical quiescence.     

Accumulation Rate Measurements at Taylor Dome, East Antarctica: Techniques and Strategies for Mass Balance Measurements in Polar Environments

Accumulation rate measurements on the East Antarctic plateau are challenging due to both spatial and temporal variability. Annual stratigraphy is often not reliably or consistently preserved in the firn, and so accumulation cannot be determined from snow pit stratigraphy. We present a suite of accumulation rate measurements collected over several seasons at Taylor Dome, East Antarctica. We compare net accumulation results from direct burial rate measurements and β -activity firn cores along a 35 km traverse. The two methods are consistent and show that the net accumulation varies from greater than 10 cm a⁻¹ to about 1 cm a⁻¹ (ice equivalent) southwest to northeast across the dome. We map the depth of shallow radar layering to interpolate and extrapolate these point-location measurements and show that considerable variations occur over kilometer scales resulting from subtle surface topography. We also present accumulation rates estimated from concentration of the cosmogenic isotope ¹⁰Be and from activity profiles of ²¹⁰Pb. Finally, satellite passive microwave data are used to estimate spatially averaged accumulation rates on the regional to continental scale to provide a context for these local observations. We show that robust mass balance measurements in this environment must rely on spatial and/or temporal averaging.     

Depositional and tectonic evolution of a supradetachment basin: 40Ar/39Ar geochronology of the Nova Formation, Panamint Range, California

The Nova Basin contains an upper Miocene to Pliocene supradetachment sedimentary succession that records the unroofing of the Panamint metamorphic core complex, west of Death Valley, California. Basin stratigraphy reflects the evolution of sedimentation processes from landslide emplacement during basin initiation to the development of alluvial fans composed of reworked, uplifted sections of the basin fill. ⁴⁰Ar/³⁹Ar geochronology of volcanic units in middle and lower parts of the sequence provide age control on the tectonic and depositional evolution of the basin and, more generally, insights regarding the rate of change of depositional environments in supradetachment basins. Our work, along with earlier research, indicate basin deposition from 11.38 Ma to 3.35 Ma. The data imply sedimentation rates, uncorrected for compaction, of ~100 m Myr⁻¹ in the lower, high-energy part to ~1000 m Myr⁻¹ in the middle part characterized by debris-flow fan deposition. The observed variation in sediment flux rate during basin evolution suggests that supradetachment basins have complex depositional histories involving rapid transitions in both the style and rate of sedimentation.     

Convergent flow of ice within the Astrolabe subglacial basin, Trre Adélie, East Antarctica: an hypothesis derived from nummerical modelling experiments

The existence of a large subglacial lake beneath the antarctic Ice Sheet at Terre Adélie indicates the presence of basal ice at its pressure-melting temperature. A numerical model of the ice-sheet thermal regime is employed using the balance velocity of the ice sheet as an initial model input in order to calculate ice-sheet basal temperatures. However, the results from this model show the Terre Adélie area to be characterised by basal freezing. Heat in addition to that accounted for in the model is thus required at the ice-sheet base in order for pressure melting temperatures to be attained. The sources for such heat are (1) an enhanced geothermal heat flux and (2) an increase in frictional heating caused by the flow of ice. In this paper the latter possibility is expanded by hypothesising that subglacial topography induces convergent ice flow around Terre Adélie, causing enhanced basal ice velocities. Model experiments indicate that an increase in ice velocity (from 7 to at least 42 m yr⁻¹) is required to raise the temperature of the basal ice to the pressure melting value. Increased ice velocity, and consequent frictional heat production due to convergent ice flow, may therefore be important in explaining the location of the subglacial lake in this region. These results allow the process of convergent ice flow within a contemporary ice sheet to be quantified. A verification (or otherwise) of the model results may be possible if ice surface velocity measurements from modem GPS methods are made.     

Surface vertical displacements, potential perturbations andgravity changes of a viscoelastic earth model induced by internal point dislocations

Using the viscoelastic correspondence principle, we utilize the surface coseismic spheroidal deformation fields (i.e. vertical displacements, potential perturbations and gravity changes) of SNREI earth models caused by four typical types of point dislocation, derived by Sun & Okubo (1993 ), to deduce the fundamental formulas for spheroidal fields relevant to viscoelastic earth models. In computations, we employ a strike-slip dislocation on a vertical plane buried at the bottom of the lithosphere to estimate the maximal viscous relaxation responses to this kind of source that possibly exist on the surface of the earth. We take the seismic moment as 10²²  N  m, which is characteristic of an average large earthquake. The numerical results demonstrate that, if we take the viscosity as 10¹⁹  Pa  s in the asthenosphere, and 10²¹  Pa  s in the other mantle layers, the rates of surface vertical displacements and gravity changes within about 2.5° for the 10 postseismic years are respectively 1.5–8.1  cm  yr⁻¹ and 4.0–14.9  μgal  yr⁻¹ : the viscous relaxation for this mantle viscosity profile proceeds much faster than for a constant mantle viscosity of 10²¹  Pa  s.     

Subsurface mass redistributions at Mount Etna (Italy) during the 1995–1996 explosive activity detected by microgravity studies

Gravity changes are presented from a series of field microgravity surveys conducted at Mt Etna between August 1994 and November 1996, a period including the 1995–1996 explosive summit activity. Data were collected along a microgravity network of 69 stations at a monthly to annual sampling rate, depending on each subarray of the network.
  Results show that seasonal changes in water level within the volcano may induce gravity changes of up to 20  μgal on Etna's southern slope, and indicate that significant magma movement occurred within and below Etna's edifice between 1994 and 1996. In particular, between September 1994 and October 1995, a mass increase of 2 × 10¹⁰  kg occurred 2000  m beneath the summit craters. Between October 1995 and July 1996 this mass was lost, while another 2 × 10¹⁰  kg was injected at about 1000  m  a.s.l. into the 1989 fracture system. From the gravity data alone, it is not possible to distinguish whether the first shallow intrusion (1994–1995) was then injected laterally into the 1989 fracture, or summit activity was fed by the first shallow intrusion, while new magma entered the 1989 fracture system.
  While magma was being redistributed within the volcanic edifice, measurements along an E–W-trending profile on the southern slope of the volcano detected some 1.5 × 10¹¹  kg of magma accumulating 2–3  km below sea level between October 1995 and November 1996.     

Palaeointensity results from Ethiopian basalts: implications for the Oligocene geomagnetic field strength

From a large collection of Ethiopian flood basalts (~30  Myr old) sampled for magnetostratigraphy, ⁴⁰Ar/³⁹Ar geochronology and geochemical investigations, 47 samples were selected in order to test their suitability for Thellier palaeointensity experiments. Only 17 samples from eight individual flows yielded reliable palaeointensity estimates, with flow-mean virtual dipole moments ranging from 3.0 to 10.5 × 10²²  A  m² .
  A critical review of the Oligocene palaeointensity data set, including these new Ethiopian data, indicates an Oligocene mean virtual dipole moment of 5.1 ± 2.5 × 10²²  A  m² for the complete data set. After applying mild selection criteria, the reduced data set yields a mean value of only 4.6 ± 1.9 × 10²²  A  m² . This value is significantly lower than the present-day field strength but is higher than the Mesozoic dipole low mean field. This low Oligocene field might be in agreement with the high palaeosecular variation and rather high non-dipole field invoked around 30  Ma. However, the Oligocene data set is largely dependent on the palaeointensity determinations from Armenia, obtained mainly from baked contacts, which show amazingly low dispersion at both flow and between-flow levels. More data are needed to reduce the weight of these determinations on the mean value and avoid a possible bias.     

40Ar-39Ar incremental heating studies on the Tudor Gabbro, Grenville Province, Ontario: its bearing on the North American apparent polar wander path in late Proterozoic times

Summary. ⁴⁰Ar-³⁹Ar incremental heating studies have been carried out on samples taken from the Tudor Gabbro, Grenville Province, Ontario. In an earlier K-Ar study, these rocks have yielded an isochron age of ∼700 Ma together with very high initial argon ratios. Age spectrum plots on whole-rock samples, in general, display a saddle-shaped character, with two of them exhibiting minima close to 700 Ma. No clear plateaus are observed for these rocks. A hornblende separate records the time ∼1110 Ma at which the stock finally cooled through the ∼590°C isotherm. The Tudor Gabbro was probably intruded into an area undergoing middle-amphibolite facies meta-morphism about 1180 Ma ago. The age spectra of two whole-rock samples together with that of their plagioclase separates, suggest that the stock cooled to ∼200–250°C at about 720 Ma. Slow cooling, averaging about 1°C Ma⁻¹ is indicated for this section of the Grenville Province for the period 1100–700 Ma. If the age of the Tudor Gabbro's palaeomagnetic pole position is taken to be 720 Ma, the Hadrynian Track Hypothesis leads to very high polar wander rates of > 20cm a⁻¹ for the period 820–720 Ma. If this hypothesis is rejected, the average drift rate for this period would be ∼4 cm a⁻¹, in much better agreement with published values of ∼5 cm a⁻¹ for the period 1400–820 Ma.     

Continuous GPS measurements across the Western Alps, 1996–1998

Global Positioning System (GPS) data are analysed from three permanent sites, two in the Western Alps (Grasse, France, and Zimmerwald, Switzerland) and the third in the Po Basin (Torino, Italy), for the 2.5-year period from 1 January 1996 to 1 July 1998. An analysis of the stochastic properties of the position time series reveals a significant amount of spatially and temporally correlated noise, which best fits a model combining white noise and flicker noise. The coloured noise is drastically reduced by spatially filtering the time series, suggesting that it is not due to site-specific effects such as monument motion, but rather to noise sources common to the three sites, such as reference frame, satellite orbit or Earth orientation errors. We find velocity uncertainties (95 per cent confidence interval) of 2.7  mm  yr⁻¹ at GRAS, 2.3  mm  yr⁻¹ at ZIMM, and 2.9  mm  yr⁻¹ at TORI. The residual velocity in the Eurasian reference frame is statistically greater than zero at ZIMM only, with 2.4 ± 2.3  mm  yr⁻¹ of motion in a NW direction. These uncertainties place an upper bound on the expected deformation in the Western Alps. Velocities relative to stable Eurasia do not exceed 3  mm  yr⁻¹, corresponding to a maximum strain rate of 0.03  μstrain  yr⁻¹ over the GRAS–TORI–ZIMM triangle.     

The crustal structure of the eastern Marmara region, Turkey by teleseismic receiver functions

We study the crustal structure of eastern Marmara region by applying the receiver function method to the data obtained from the 11 broad-band stations that have been in operation since the 1999 İzmit earthquake. The stacked single-event receiver functions were modelled by an inversion algorithm based on a five-layered crustal velocity model to reveal the first-order shear-velocity discontinuities with a minimum degree of trade-off. We observe crustal thickening from west (29–32 km) to east (34–35 km) along the North Anatolian Fault Zone (NAFZ), but we observe no obvious crustal thickness variation from north to south while crossing the NAFZ. The crust is thinnest beneath station TER (29 km), located near the Black Sea coast in the west and thickest beneath station TAR (35 km), located inland in the southeast. The average crustal thickness and S -wave velocity for the whole regions are  31 ± 2  km and  3.64 ± 0.15 km s⁻¹  , respectively. The eastern Marmara region with its average crustal thickness, high heat flow value (101 ± 11 mW m⁻²) and with its remarkable extensional features seems to have a Basin and Range type characteristics, but the higher average shear velocities (∼3.64 km s⁻¹) and crustal thickening from 29 to 35 km towards the easternmost stations indicate that the crustal structure shows a transitional tectonic regime. Therefore, we conclude that the eastern Marmara region seems to be a transition zone between the Marmara Sea extensional domain and the continental Anatolian inland region.     

A surge of Perseibreen, Svalbard, examined using aerial photography and ASTER high resolution satellite imagery

The identification of surge activity is important in assessing the duration of the active and quiescent phases of the surge cycle of Svalbard glaciers. Satellite and aerial photographic images are used to identify and describe the form and flow of Perseibreen, a valley glacier of 59 km² on the east coast of Spitsbergen. Heavy surface crevassing and a steep ice front, indicative of surge activity, were first observed on Perseibreen in April 2002. Examination of high resolution (15 m) Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite imagery confirmed this surge activity. Perseibreen retreated by almost 750 m between 1961 and 1990. Between 1990 and the summer of 2000, Perseibreen switched from retreat and its front began to advance. Rapid advance was underway during the period June 2000 to May 2001, with terminus advance at over 400 m yr⁻¹. Between May and August 2001 the rate increased to over 750 m yr⁻¹. The observed crevasse orientation indicates that ice was in longitudinal tension, suggesting the down-glacier transfer of mass. Ice surface velocities, derived from image correlation between ASTER images, were 2-2.5 m d⁻¹ between May and August 2001. The glacier was flowing at a relatively uniform speed with sharp velocity gradients located close to its lateral margins, a velocity structure typical of ice masses in the active phase of the surge cycle. The stress regime is extensional throughout and the surge appears to be initiated low on the glacier. This is similar to the active-phase dynamics of other Svalbard tidewater glaciers. Perseibreen has probably been inactive since at least 1870, a period of about 130 years to the present surge which defines a minimum length for the quiescent phase.     

A comparison of annual and seasonal carbon dioxide effluxes between sub-Arctic Sweden and High-Arctic Svalbard

Recent climate change predictions suggest altered patterns of winter precipitation across the Arctic. It has been suggested that the presence, timing and quantity of snow all affect microbial activity, thus influencing CO₂ production in soil. In this study annual and seasonal emissions of CO₂ were estimated in High-Arctic Adventdalen, Svalbard, and sub-Arctic Latnjajaure, Sweden, using a new trace gas-based method to track real-time diffusion rates through the snow. Summer measurements from snow-free soils were made using a chamber-based method. Measurements were obtained from different snow regimes in order to evaluate the effect of snow depth on winter CO₂ effluxes. Total annual emissions of CO₂ from the sub-Arctic site (0.662–1.487 kg CO₂ m^–2 yr^–1) were found to be more than double the emissions from the High-Arctic site (0.369–0.591 kg CO₂ m^–2 yr^–1). There were no significant differences in winter effluxes between snow regimes or vegetation types, indicating that spatial variability in winter soil CO₂ effluxes are not directly linked to snow cover thickness or soil temperatures. Total winter emissions (0.004–0.248 kg CO₂ m^–2) were found to be in the lower range of those previously described in the literature. Winter emissions varied in their contribution to total annual production between 1 and 18%. Artificial snow drifts shortened the snow-free period by 2 weeks and decreased the annual CO₂ emission by up to 20%. This study suggests that future shifts in vegetation zones may increase soil respiration from Arctic tundra regions.     

Modelling of short-interval silent slip events in deeper subduction interfaces considering the frictional properties at the unstable–stable transition regime

Recent high-resolution observations of crustal movements have revealed silent slip events (SSEs) with propagation velocities of around 10–15 km d⁻¹ and with intervals of 3–14 months along the deeper parts of the Cascadia and Nankai subduction zones. This study develops 2-D and 3-D models of these short-interval SSEs considering the frictional behaviour that was confirmed experimentally by Shimamoto for the unstable–stable transition regime. To represent this frictional behaviour, a small cut-off velocity to an evolution effect is introduced in a rate- and state-dependent friction law. When the cut-off velocity to the evolution effect is significantly smaller than that to a direct effect, steady-state friction exhibits velocity weakening at low slip velocities and velocity strengthening at high slip velocities. At the deeper Cascadia and Nankai subduction interfaces, the pore pressure is inferred to be high because of the dehydration of materials in the descending plate. Under conditions where the pore-fluid pressure is nearly equal to the lithostatic pressure and the critical weakening displacement is very small, short-interval SSEs with propagation velocities and slip velocities of 4–8 km d⁻¹ and  2 − 4 × 10⁻⁷  m s⁻¹, respectively, can be reproduced. The propagation velocity of short-interval SSEs is in proportion to the slip velocity. The results also show that during the nucleation process of large earthquakes, the occurrence of short-interval SSEs becomes irregular because of the accelerated slips that occur at the bottom of the seismogenic zone. Our results suggest that monitoring of short-interval SSEs might be useful for forecasting the main earthquakes.     

Chemical denudation rates in Kärkevagge, Swedish Lapland

This study examines the spatial and temporal variability of chemical denudation rates in Kärkevagge, northern Sweden. The chemical flux rates within the valley are strongly influenced by the local geology. Chemical denudation rates determined for the study period are more than double those previously reported in the literature for this valley. Rates of greater than 46t km⁻² a⁻¹ were measured at the valley mouth over the course of the melt season. This difference is likely due to differences in measurement technique compared to that used by past researchers. This rate is also much higher than for other arctic and alpine watersheds. Chemical denudation in Kärkevagge is comparable to larger temperate rivers. The rapid chemical denudation in Kärkevagge is likely due to sulfide weathering creating acid solutions.     

Ice-free conditions in Novaya Zemlya 35 000–30 000 cal years B.P., as indicated by radiocarbon ages and amino acid racemization evidence from marine molluscs

Novaya Zemlya was covered by the eastern part of the Barents–Kara ice sheet during the glacial maximum of marine isotope stage 2 (MIS 2). We obtained ¹⁴C ages on 37 samples of mollusc shells from various sites on the islands. Most samples yielded ages in the range of 48–26 ¹⁴C Ky. Such old samples are sensitive to contamination by young ¹⁴C, and therefore their reliability was assessed using replicate analyses and amino acid geochronology. The extent of aspartic acid racemization (Asp D/L) indicates that many of the ¹⁴C ages are correct, whereas some are minimum ages only. The results indicate that a substantial part of Novaya Zemlya was ice-free about 35–27 ¹⁴C Kya, and probably even earlier. Corresponding shorelines up to >140 m a.s.l. indicate a large Barents–Kara ice sheet during early MIS 3. These results are consistent with findings from Svalbard and northern Russia: in both places a large MIS 4/3 Barents–Kara ice sheet is postulated to have retreated about 50 Kya, followed by an ice-free interstadial that lasted until up to ca. 25 Kya. The duration of the MIS 2 glaciation in Novaya Zemlya was calculated by applying the D/L values to a kinetic equation for Asp racemization. This indicates that the islands were ice covered for less than 3000 years if the basal temperature was 0^oC, and for less than 10 000 years if it was −5^oC.     

Nitrogen uptake rates in phytoplankton and ice algae in the Barents Sea

Uptake rates of NH₄⁺, NO₃⁻ and dissolved organic nitrogen (urea) were measured in phytoplankton and in ice algae in the Barents Sea using a ¹⁵N-technique. NO₃⁻ was the most important nitrogen source for the ice algae (f-ratio = 0.92). The in situ irradiances in the subsurface chlorophyll maximum and in the ice algal communities were low. The in situ NO₃⁻ uptake rate in the ice algal communities was light-limited The in situ NO₃⁻ and NH₄⁻ uptake rates in the subsurface chlorophyll maximum were at times light-limited. It is hypothesised that NH₄⁺ may accumulate in low light in the bottom of the euphotic zone and inhibit the in situ NO₃⁻ uptake rate.