bluuebunny/crossref_metadata_2025 huggingface dataset:
Exchange across the shelf break at high southern latitudes
Authors: Klinck J. M. , Dinniman M. S.
Published: April 2010

Abstract. Exchange of water across the Antarctic shelf break has considerable scientific and societal importance due to its effects on circulation and biology of the region, conversion of water masses as part of the global overturning circulation and basal melt of glacial ice and the consequent effect on sea level rise. The focus in this paper is the onshore transport of warm, oceanic Circumpolar Deep Water (CDW); export of dense water from these shelves is equally important, but has been the focus of other recent papers and will not be considered here. A variety of physical mechanisms are described which could play a role in this onshore flux. The relative importance of some processes are evaluated by simple calculations. A numerical model for the Ross Sea continental shelf is used as an example of a more comprehensive evaluation of the details of cross-shelf break exchange. In order for an ocean circulation model simulate these processes at high southern latitudes, it needs to have high spatial resolution, realistic geometry and bathymetry. Grid spacing smaller than the first baroclinic radius deformation (a few km) is required to adequately represent the circulation. Because of flow-topography interactions, bathymetry needs to be represented at these same small scales. Atmospheric conditions used to force these circulation models also need to be known at a similar small spatial resolution (a few km) in order to represent orographically controlled winds (coastal jets) and katabatic winds. Significantly, time variability of surface winds strongly influences the structure of the mixed layer. Daily, if not more frequent, surface fluxes must be imposed for a realistic surface mixed layer. Sea ice and ice shelves are important components of the coastal circulation. Ice isolates the ocean from exchange with the atmosphere, especially in the winter. Melting and freezing of both sea ice and glacial ice influence salinity and thereby the character of shelf water. These water mass conversions are known to have an important effect on export of dense water from many Antarctic coastal areas. An artificial dye, as well as temperature, is used to diagnose the flux of CDW onto the shelf. Model results for the Ross Sea show a vigorous onshore flux of oceanic water across the shelf break both at depth and at the surface as well as creation of dense water (High Salinity Shelf Water) created by coastal polynyas in the western Ross Sea.
A numerical scheme to calculate temperature and salinity dependent air-water transfer velocities for any gas
Authors: Johnson M. T.
Published: April 2010

Abstract. The transfer velocity determines the rate of exchange of a gas across the air-water interface for a given deviation from Henry's law equilibrium between the two phases. In the thin film model of gas exchange, which is commonly used for calculating gas exchange rates from measured concentrations of trace gases in the atmosphere and ocean/freshwaters, the overall transfer is controlled by diffusion-mediated films on either side of the air-water interface. Calculating the total transfer velocity (i.e. including the influence from both molecular layers) requires the Henry's law constant and the Schmidt number of the gas in question, the latter being the ratio of the viscosity of the medium and the molecular diffusivity of the gas in the medium. All of these properties are both temperature and (on the water side) salinity dependent and extensive calculation is required to estimate these properties where not otherwise available. The aim of this work is to standardize the application of the thin film approach to flux calculation from measured and modelled data, to improve comparability, and to provide a numerical framework into which future parameter improvements can be integrated. A detailed numerical scheme is presented for the calculation of the gas and liquid phase transfer velocities (ka and kw respectively) and the total transfer velocity, K. The scheme requires only basic physical chemistry data for any gas of interest and calculates K over the full range of temperatures, salinities and wind-speeds observed in and over the ocean. Improved relationships for the wind-speed dependence of ka and for the salinity-dependence of the gas solubility (Henry's law) are derived. Comparison with alternative schemes and methods for calculating air-sea flux parameters shows good agreement in general but significant improvements under certain conditions. The scheme is provided as a downloadable program in the supplementary material, along with input files containing molecular weight, solubility and structural data for 80 gases of general interest, enabling calculation of the total transfer velocity over ranges of temperature and salinity for each gas.
On the numerical resolution of the bottom layer in simulations of oceanic gravity currents
Authors: Laanaia N. , Wirth A. , Molines J. M. , Barnier B. , Verron J.
Published: April 2010

Abstract. The role of an increased numerical vertical resolution, leading to an explicit resolution of the bottom Ekman layer dynamics, is investigated. Using the hydrostatic ocean model NEMO-OPA9, we demonstrate that the dynamics of an idealised gravity current (on an inclined plane), is well captured when a few (around five) sigma-coordinate levels are added near the ocean floor. Such resolution allows to considerably improve the representation of the descent and transport of the gravity current and the Ekman dynamics near the ocean floor, including the important effect of Ekman veering, which is usually neglected in today's simulations of the ocean dynamics. Results from high resolution simulations (with σ and z-coordinates) are compared to simulations with a vertical resolution commonly employed in today's ocean models. The latter show a downslope transport that is reduced by almost an order of magnitude and the decrease in the along slope transport is reduced six-fold. We strongly advocate for an increase of the numerical resolution at the ocean floor, similar to the way it is done at the ocean surface and at the lower boundary in atmospheric models.
Importance of the variability of hydrographic preconditioning for deep convection in the Gulf of Lion, NW Mediterranean
Authors: Grignon L. , Smeed D. A. , Bryden H. L. , Schroeder K.
Published: April 2010

Abstract. We study the variability of hydrographic preconditioning defined as the heat and salt contents in the Ligurian Sea before convection. The stratification is found to reach a maximum in the intermediate layer in December, whose causes and consequences for the interannual variability of convection are investigated. Further study of the interannual variability and correlation tests between the properties of the deep water formed and the winter surface fluxes support the description of convection as a process that transfers the heat and salt contents from the top and intermediate layers to the deep layer. A proxy for the rate of transfer is given by the final convective mixed layer depth, that is shown to depend equally on the surface fluxes and on the preconditioning. In particular, it was found that deep convection in winter 2004–2005 would have happened even with normal winter conditions, due to low pre-winter stratification.
Technical Note: A low cost Unmanned Aerial Vehicle for ship based science missions
Authors: Waugh E. , Mowlem M.
Published: April 2010

Abstract. A low-cost Unmanned Aerial Vehicle is compared with those already available and the motivation for its development is established. It is targeted at ship-based science missions and potential applications are described including a specific science case to measure white capping in the deep ocean. The current vehicle includes a range of more than 1000 km, carrying a payload of 2 kg and it can be launched and recovered from a coastal research vessel. The vehicle has flown successfully in Force 4 gusting Force 6–7 wind conditions, an important requirement for operation at sea. Data analysis is performed on images captured by the vehicle to provide a measurement of wave period and white capping fraction. The next stage of the project is to develop a suitable payload and perform a demonstration science mission.
Geochronology and thermochronology of Cretaceous plutons and metamorphic country rocks, Anyui-Chukotka fold belt, North East Arctic Russia
Authors: Miller E. L. , Katkov S. M. , Strickland A. , Toro J. , Akinin V. V. , Dumitru T. A.
Published: April 2010

Abstract. U-Pb isotopic dating of seven granitoid plutons and associated intrusions from the Bilibino region (Arctic Chukotka, Russia) was carried out using the SHRIMP-RG. The crystallization ages of these granitoids, which range from approximately 116.9±2.5 to 108.5±2.7 Ma, bracket two regionally significant deformational events. The plutons cut folds, steep foliations and thrust-related structures related to sub-horizontal shortening at lower greenschist facies conditions (D1), believed to be the result of the collision of the Arctic Alaska-Chukotka microplate with Eurasia along the South Anyui Zone (SAZ). Deformation began in the Late Jurassic, based on fossil ages of syn-orogenic clastic strata, and involves strata as young as early Cretaceous (Valanginian) north of Bilibino and as young as Hauterivian-Barremian, in the SAZ. The second phase of deformation (D2) is developed across a broad region around and to the east of the Lupveem batholith of the Alarmaut massif and is interpreted to be coeval with magmatism. D2 formed gently-dipping, high-strain foliations (S2). Growth of biotite, muscovite and actinolite define S2 adjacent to the batholith, while chlorite and white mica define S2 away from the batholith. Sillimanite (± andalusite) at the southeastern edge the Lupveem batholith represent the highest grade metamorphic minerals associated with D2. D2 is interpreted to have developed during regional extension and crustal thinning. Extension directions as measured by stretching lineations, quartz veins, boudinaged quartz veins is NE-SW to NW-SE. Mapped dikes associated with the plutons trend mostly NW-SE and indicate NE-SW directed extension. 40Ar/39Ar ages from S2 micas range from 109.3±1.2 to 103.0±1.8 Ma and are interpreted as post-crystallization cooling ages following a protracted period of magmatism and high heat flow. Regional uplift and erosion of many kilometers of cover produced a subdued erosional surface prior to the eruption of volcanic rocks of the Okhotsk-Chukotka volcanic belt (OCVB) whose basal units (~87 Ma) overlie this profound regional unconformity. A single fission track age on apatite from granite in the Alarmaut massif yielded an age of 90±11 Ma, in good agreement with this inference.
Structural studies near Pevek, Russia: implications for formation of the East Siberian Shelf and Makarov Basin of the Arctic Ocean
Authors: Miller E. L. , Verzhbitsky V. E.
Published: April 2010

Abstract. The Pevek region of Arctic Russia provides excellent beach cliff exposure of sedimentary and igneous rocks that yield detailed information on the nature, progression and timing of structural events in this region. Regional folding and thrust faulting, with the development of a south-dipping axial plane cleavage/foliation developed during N-S to NE-SW directed shortening and formation of the Chukotka-Anyui fold belt. This deformation involves strata as young as Valanginian (136–140 Ma, Gradstein et al., 2004). Fold-related structures are cut by intermediate to silicic batholiths, plutons and dikes of Cretaceous age. Reported K-Ar whole rock and mineral ages on the granitoids range from 144 to 85 Ma, but to the south, more reliable U-Pb zircon ages on compositionally similar plutons yield a much narrower age range of ~120–105 Ma (Miller et al., this volume) and a pluton in Pevek yields a U-Pb age on zircon of 108.1±1.1 Ma with evidence for inheritance of slightly older 115 Ma zircons. Magmas were intruded during an episode of E-W to ENE-WSW directed regional extension based on the consistent N-S to NNW-SSE orientation of over 800 mapped dikes and quartz veins. Analysis of small-offset faults and slickensides yield results compatible with those inferred from the dikes. Younger tectonic activity across this region is minor and the locus of magmatic activity moved southward towards the Pacific margin as represented by the <90 Ma Okhotsk-Chukotsk volcanic belt (OCVB). A lengthy period of uplift and erosion occurred after emplacement of Cretaceous plutons and produced the peneplain beneath the younger OCVB. Based on our studies, we speculate that ~120–105 Ma magmatism, which heralds a change in tectonic regime from compression to extension, could represent one of the consequences of the inception of rifting in the Amerasian Basin of the Arctic, forming the Makarov Basin north of the Siberian shelf at this longitude. A synthesis of available seismic reflection, gravity and magnetic data for the offshore Siberian Shelf reveals a widespread, seismically mappable basement-sedimentary cover contact that deepens northward towards the edge of the shelf with few other significant basins. Various ages have been assigned to the oldest strata above the unconformity, ranging from Cretaceous (Albian – 112–100 Ma) to Tertiary (Paleocene–Eocene – ~60–50 Ma). The period of uplift and erosion documented along the Arctic coast of Russia at this longitude could represent the landward equivalent of the (yet undrilled) offshore basement-sedimentary cover contact, thus overlying sedimentary sequences could be as old as early Late Cretaceous. Although quite speculative, these conclusions suggest that land-based geologic, structural, petrologic and geochronologic studies could provide useful constraints to help resolve the plate tectonic history of the Arctic Ocean.
Age and paleomagnetism of the Okhotsk-Chukotka Volcanic Belt (OCVB) near Lake El'gygytgyn, Chukotka, Russia
Authors: Stone D. B. , Layer P. W. , Raikevich M. I.
Published: April 2010

Abstract. Paleomagnetic results from the upper two thirds of the whole section of the Okhotsk-Chukotka Volcanic Belt (OCVB) volcanics exposed in the area around Lake El'gygytgyn, Chukotka yield stable, consistent magnetic vectors and well-preserved reversed directions. The magnetostratigraphy and 40Ar/39Ar geochronologic data reported here indicate that the sampled OCVB volcanics were erupted between about 90 and 67 Ma, and show no significant change in the apparent pole position over that time. The OCVB extends from northeast China, across Northeast Russia to the Bering Straight. This belt is made up of both extrusive and intrusive rocks, with the extrusive rocks and their associated sediments being dominant. The whole belt important in interpreting the paleogeography of the region because it overlies many of the accreted terranes of Northeast Russia. Most importantly, it overlies parts of the Chukotka-Alaska block which is thought to have moved out of the Arctic Ocean region, as well as terranes accreted from the south. These latter terranes have been rafted northwards on the paleo-plates of the Pacific, implying that the present relative paleogeography of all of the terranes overlain by the OCVB were essentially in place by 67 Ma, and possibly as early as 90 Ma. However, comparing our paleomagnetic pole position for the OCVB with those for North America and Eurasia (a proxy for Siberia) shows a statistically significant displacement of the OCVB pole to the south west. This implies that not only the OCVB, but the underlying terranes of northeast Russia, experienced southerly displacement with respect to the Siberian and North American platforms since the Late Cretaceous.
Geotectonic setting of the Tertiary Uyandina and Indigirka-Zyryanka basins, Republic Sakha (Yakutia), Northeast Russia, using coal rank data
Authors: Paech H.-J.
Published: April 2010

Abstract. Outcrops along the Inach River in the Uyandina basin and those along the Myatis' River in the Indigirka-Zyryanka basin were studied in detail and sampled for coal rank determinations. The Uyandina basin is an intramontane pull-apart basin characterized by extensional structures within the Moma rift system. The coal rank is below 0.3% vitrinite reflectance (Rr), which indicates shallow, immature conditions of basin formation and very low subsidence. The Myatis' River coal-bearing outcrops in the Indigirka-Zyryanka basin reveal compression induced by continent collision. The compressive deformation includes also lowermost Pliocene strata. Due to the position in the Verkhoyansk-Chersky fold belt adjacent to the Kolyma-Omolon microcontinent the Indigirka-Zyryanka basin has much in common with a foredeep, i.e. the asymmetry in thickness and tectonic structure. The vitrinite reflectance data (Rr) which range from 0.25% to more than 5% reinforce the accepted models that describe basin subsidence and geothermal history and the tectonic deformation.
Direct evidence for continuous radar reflector originating from changes in crystal-orientation fabric
Authors: Eisen O. , Hamann I. , Kipfstuhl S. , Steinhage D. , Wilhelms F.
Published: April 2010

Abstract. The origin of a strong continuous radar reflector observed with airborne radio-echo sounding (RES) at the EPICA deep-drilling site in Dronning Maud Land, Antarctica, is identified as a transition in crystal fabric orientation from a vertical girdle to an increased single-pole orientation seen along the ice core. The reflector is observed with a 60 ns and 600 ns long pulse at a frequency of 150 MHz, spans one pulse length, is continuous over 5 km, and occurs at a depth of about 2025–2045 m at the drill site. Changes in conductivity as reflector origin are excluded by investigating the ice-core profile, synthetic RES data, and a RES profile with different electromagnetic polarisation azimuths. The reflector's magnitude shows maximum values for polarisation parallel to the nearby ice divide and disappears for polarisation perpendicular to it, identifying the orientation of the girdle to lie in the vertical plane parallel to the ice divide. Observations allow us to extrapolate the crystal orientation feature along the reflector in space, with implications for ice-sheet dynamics and modeling.
Reconstruction of the 1979–2006 Greenland ice sheet surface mass balance using the regional climate model MAR
Authors: Fettweis X.
Published: April 2010

Abstract. Results from a 28-year simulation (1979–2006) over the Greenland ice sheet (GrIS) reveal an increase of solid precipitation (+0.4±2.5 km3 yr−2) and run-off (+7.9±3.3 km3 yr−2) of surface meltwater. The net effect of these competing factors is a significant Surface Mass Balance (SMB) loss of −7.2±5.1 km3 yr−2. The contribution of changes in the net water vapour flux (+0.02±0.09 km3 yr−2) and rainfall (+0.2±0.2 km3 yr−2) to the SMB variability is negligible. The meltwater supply has increased because the GrIS surface has been warming up +2.4°C since 1979. Sensible heat flux, latent heat flux and net solar radiation have not varied significantly over the last three decades. However, the simulated downward infrared flux has increased by 9.3 W m−2 since 1979. The natural climate variability (e.g. the North Atlantic Oscillation) does not explain these changes. The recent global warming, due to the greenhouse gas concentration increase induced by human activities, could be a cause of these changes. The doubling of surface meltwater flux into the ocean over the period 1979–2006 suggests that the overall ice sheet mass balance has been increasingly negative, given the likely meltwater-induced acceleration of outlet glaciers. This study suggests that increased melting overshadows over an increased accumulation in a warming scenario and that the GrIS is likely to keep losing mass in the future. An enduring GrIS melting will probably affect in the future an certain effect on the stability of the thermohaline circulation and the global sea level rise.
Spatial structures in the heat budget of the Antarctic atmospheric boundary layer
Authors: van de Berg W. J. , van den Broeke M. R. , van Meijgaard E.
Published: April 2010

Abstract. Output from the regional climate model RACMO2/ANT is used to calculate the heat budget of the Antarctic atmospheric boundary layer (ABL). The main feature of the wintertime Antarctic ABL is a persistent temperature deficit compared to the free atmosphere. The magnitude of this deficit is controlled by the heat budget. During winter, transport of heat towards the surface by turbulence and net longwave emission are the primary ABL cooling terms. These processes show horizontal spatial variability only on continental scales. Vertical and horizontal, i.e. along-slope, advection of heat are the main warming terms. Over regions with convex ice sheet topography, i.e. domes and ridges, warming by downward vertical advection is enhanced due to divergence of the ABL wind field. Horizontal advection balances excess warming caused by vertical advection, hence the temperature deficit in the ABL weakens over domes and ridges along the prevailing katabatic wind. Conversely, vertical advection is reduced in regions with concave topography, i.e. valleys, where the ABL temperature deficit enlarges along the katabatic wind. Along the coast, horizontal and vertical advection is governed by the inability of the large-scale circulation to adapt to small scale topographic features. Meso-scale topographic structures have thus a strong impact on the ABL winter temperature, besides latitude and surface elevation. During summer, this mechanism is much weaker, and the horizontal variability of ABL temperatures is smaller.
Diagnosing the extreme surface melt event over southwestern Greenland in 2007
Authors: Tedesco M. , Serreze M. , Fettweis X.
Published: April 2010

Abstract. Analysis of passive microwave brightness temperatures from the space-borne Special Sensor Microwave Imager (SSM/I) documents a record surface snowmelt over high elevations (above 2000 m) of the Greenland ice sheet during summer of 2007. To interpret this record, results from the SSM/I are examined in conjunction with fields from the National Centers for Environmental Prediction/National Center for Atmospheric Research reanalysis and output from a regional climate model. The record surface melt reflects unusually warm conditions, seen in positive summertime anomalies of surface air temperatures, downwelling longwave radiation, 1000–500 hPa atmospheric thickness, and the net surface energy flux, linked in turn to southerly airflow over the ice sheet. Low snow accumulation may have contributed to the record through promoting anomalously low surface albedo.
Exploring uncertainty in glacier mass balance modelling with Monte Carlo simulation
Authors: Machguth H. , Purves R. S. , Oerlemans J. , Hoelzle M. , Paul F.
Published: April 2010

Abstract. By means of Monte Carlo simulations we calculated uncertainty in modelled cumulative mass balance over 400 days at one particular point on the tongue of Morteratsch Glacier, Switzerland, using a glacier energy balance model of intermediate complexity. Before uncertainty assessment, the model was tuned to observed mass balance for the investigated time period and its robustness was tested by comparing observed and modelled mass balance over 11 years, yielding very small deviations. Both systematic and random uncertainties are assigned to twelve input parameters and their respective values estimated from the literature or from available meteorological data sets. The calculated overall uncertainty in the model output is dominated by systematic errors and amounts to 0.7 m w.e. or approximately 10% of total melt over the investigated time span. In order to provide a first order estimate on variability in uncertainty depending on the quality of input data, we conducted a further experiment, calculating overall uncertainty for different levels of uncertainty in measured global radiation and air temperature. Our results show that the output of a well calibrated model is subject to considerable uncertainties, in particular when applied for extrapolation in time and space where systematic errors are likely to be an important issue.
The ISMIP-HOM benchmark experiments performed using the Finite-Element code Elmer
Authors: Gagliardini O. , Zwinger T.
Published: April 2010

Abstract. The aim of this paper is to describe in detail how the benchmark tests ISMIP-HOM (Ice Sheet Model Intercomparison Project-Higher-Order ice-sheet Model) have been performed using the open source finite element (FE) code Elmer (http://www.csc.fi/elmer). The ISMIP-HOM setup consists of five diagnostic and one prognostic experiments, for both 2-D and 3-D geometries. For all the tests, the full-Stokes equations are solved. Some technical points concerning FE, such as mesh characteristics, stabilisation methods, numerical methods used to solve the linear system and parallel performance are discussed. For all these setups, the CPU time consumption in relation to the accuracy of the solution is analysed. Based on these findings, some general rules on optimising the computing time versus the accuracy of the results are deduced.
Analytical solutions for the surface response to small amplitude perturbations in boundary data in the shallow-ice-stream approximation
Authors: Gudmundsson G. H.
Published: April 2010

Abstract. New analytical solutions describing the effects of small-amplitude perturbations in boundary data on flow in the shallow-ice-stream approximation are presented. These solutions are valid for a non-linear Weertman-type sliding law and for Newtonian ice rheology. Comparison is made with corresponding solutions of the shallow-ice-sheet approximation, and with solutions of the full Stokes equations. The shallow-ice-stream approximation is commonly used to describe large-scale ice stream flow over a weak bed, while the shallow-ice-sheet approximation forms the basis of most current large-scale ice sheet models. It is found that the shallow-ice-stream approximation overestimates the effects of bed topography perturbations on surface profile for wavelengths less than about 5 to 10 ice thicknesses, the exact number depending on values of surface slope and slip ratio. For high slip ratios, the shallow-ice-stream approximation gives a very simple description of the relationship between bed and surface topography, with the corresponding transfer amplitudes being close to unity for any given wavelength. The shallow-ice-stream estimates for the timescales that govern the transient response of ice streams to external perturbations are considerably more accurate than those based on the shallow-ice-sheet approximation. In particular, in contrast to the shallow-ice-sheet approximation, the shallow-ice-stream approximation correctly reproduces the short-wavelength limit of the kinematic phase speed given by solving a linearised version of the full Stokes system. In accordance with the full Stokes solutions, the shallow-ice-sheet approximation predicts surface fields to react weakly to spatial variations in basal slipperiness with wavelengths less than about 10 to 20 ice thicknesses.
Diagnostic and prognostic simulations with a full Stokes model accounting for superimposed ice of Midtre Lovénbreen, Svalbard
Authors: Zwinger T. , Moore J. C.
Published: April 2010

Abstract. We present steady state (diagnostic) and transient (prognostic) simulations of Midtre Lovénbreen, Svalbard performed with the thermo-mechanically coupled full-Stokes code Elmer. This glacier has an extensive data set of geophysical measurements available spanning several decades, that allow for constraints on model descriptions. Consistent with this data set, we included a simple model accounting for the formation of superimposed ice. Diagnostic results indicated that a dynamic adaptation of the free surface is necessary, to prevent non-physically high velocities in a region of under determined bedrock depths. Observations from ground penetrating radar of the basal thermal state agree very well with model predictions, while the dip angles of isochrones in radar data also match reasonably well with modelled isochrones, despite the numerical deficiencies of estimating ages with a steady state model. Prognostic runs for 53 years, using a constant accumulation/ablation pattern starting from the steady state solution obtained from the configuration of the 1977 DEM show that: 1 the unrealistic velocities in the under determined parts of the DEM quickly damp out; 2 the free surface evolution matches well measured elevation changes; 3 the retreat of the glacier under this scenario continues with the glacier tongue in a projection to 2030 being situated ≈500 m behind the position in 1977.
Changes of Wilkins Ice Shelf over the past 15 years and inferences on its stability
Authors: Braun M. , Humbert A. , Moll A.
Published: April 2010

Abstract. The Wilkins Ice Shelf is situated on the Antarctic Peninsula, a region where seven ice shelves disintegrated or retreated between 1995 and 2002. This study combines various remote sensing datasets from Wilkins Ice Shelf, with the aim of detecting its present and recent dynamics as well as recent changes. The survey includes structural mapping, ERS-1/2 SAR interferometry and analysis of ICESat GLAS ice surface elevation data. Ice front retreat rates from 1986 to 2008 showed several distinct break-up events, including one in February 2008, when 40% of a part of the ice shelf that connected two islands broke off. Surface elevations have been used to study tidal effects, crack formation and to estimate the ice thickness over the floating area. The derived interferometric velocities cover the south-eastern part of the ice shelf as well as major tributaries and reveal maximum inflow speeds of up to 330 m a−1. We show that drainage of melt ponds into crevasses were of no relevance for the break-up at Wilkins Ice Shelf. Buoyancy forces caused rift formation before the break-up in February 2008. Additionally, the evolution of failure zones of the order of tenths of kilometres in length in pre-conditioned locations at ice rises is shown. Investigation of the current (February 2009) situation shows that about 3100 km2 at the Northern Wilkins Ice Shelf are endangered, however, there is no visible signature that the remaining 8000 km2 are at risk.
Comparison of the meteorology and surface energy balance at Storbreen and Midtdalsbreen, two glaciers in southern Norway
Authors: Giesen R. H. , Andreassen L. M. , van den Broeke M. R. , Oerlemans J.
Published: April 2010

Abstract. We compare 5 years of meteorological records from automatic weather stations (AWSs) on Storbreen and Midtdalsbreen, two glaciers in southern Norway, located approximately 120 km apart. The records are obtained from identical AWSs with an altitude difference of 120 m and cover the period September 2001 to September 2006. Air temperature at the AWS locations is found to be highly correlated, even with the seasonal cycle removed. The most striking difference between the two sites is the difference in wind climate. Midtdalsbreen is much more under influence of the large-scale circulation with wind speeds on average a factor 1.75 higher. On Storbreen, weaker katabatic winds are dominant. The main melt season is from May to September at both locations. During the melt season, incoming and net solar radiation are larger on Midtdalsbreen, whereas incoming and net longwave radiation are larger on Storbreen, primarily caused by thicker clouds on the latter. The turbulent fluxes are a factor 1.7 larger on Midtdalsbreen, mainly due to the higher wind speeds. Inter-daily fluctuations in the surface energy fluxes are very similar at the AWS sites. On average, melt energy is a factor 1.3 larger on Midtdalsbreen, a result of both larger net radiation and larger turbulent fluxes. The relative contribution of net radiation to surface melt is larger on Storbreen (76%) than on Midtdalsbreen (66%). As winter snow depth at the two locations is comparable in most years, the larger amount of melt energy results in an earlier disappearance of the snowpack on Midtdalsbreen and 70% more ice melt than on Storbreen. We compare the relative and absolute values of the energy fluxes on Storbreen and Midtdalsbreen with reported values for glaciers at similar latitudes. Furthermore, a comparison is made with meteorological variables measured at two nearby weather stations, showing that on-site measurements are essential for an accurate calculation of the surface energy balance and melt rate.
Testing longwave radiation parameterizations under clear and overcast skies at Storglaciären, Sweden
Authors: Sedlar J. , Hock R.
Published: April 2010

Abstract. Energy balance based glacier melt models require accurate estimates of incoming longwave radiation but direct measurements are often not available. Multi-year near-surface meteorological data from Storglaciären, Northern Sweden, were used to evaluate commonly used longwave radiation parameterizations in a glacier environment under clear-sky and all-sky conditions. Parameterizations depending solely on air temperature performed worse than those which include water vapor pressure. All models tended to overestimate incoming longwave radiation during periods of low longwave radiation, while incoming longwave was underestimated when radiation was high. Under all-sky conditions root mean square error (RMSE) and mean bias error (MBE) were 17 to 20 W m−2 and −5 to 1 W m−2, respectively. Two attempts were made to circumvent the need of cloud cover data. First cloud fraction was parameterized as a function of the ratio, τ, of measured incoming shortwave radiation and calculated top of atmosphere radiation. Second, τ was related directly to the cloud factor (i.e. the increase in sky emissivity due to clouds). Despite large scatter between τ and both cloud fraction and the cloud factor, resulting calculations of hourly incoming longwave radiation for both approaches were only slightly more variable with RMSE roughly 3 W m−2 larger compared to using cloud observations as input. This is promising for longwave radiation modeling in areas where shortwave radiation data are available but cloud observations are not.
Role of glaciers in watershed hydrology: a preliminary study of a "Himalayan catchment"
Authors: Thayyen R. J. , Gergan J. T.
Published: April 2010

Abstract. A large number of Himalayan glacier catchments are under the influence of humid climate with snowfall in winter (November–April) and south-west monsoon in summer (June–September) dominating the regional hydrology. Such catchments are defined as "Himalayan catchment", where the glacier meltwater contributes to the river flow during the period of annual high flows produced by the monsoon. The winter snow dominated Alpine catchments of the Kashmir and Karakoram region and cold-arid regions of the Ladakh mountain range are the other major glacio-hydrological regimes identified in the region. Factors influencing the river flow variations in a "Himalayan catchment" were studied in a micro-scale glacier catchment in the Garhwal Himalaya, covering an area of 77.8 km2. Three hydrometric stations were established at different altitudes along the Din Gad stream and discharge was monitored during the summer ablation period from 1998 to 2004, with an exception in 2002. These data have been analysed along with winter/summer precipitation, temperature and mass balance data of the Dokriani glacier to study the role of glacier and precipitation in determining runoff variations along the stream continuum from the glacier snout to 2360 m a.s.l. The study shows that the inter-annual runoff variation in a "Himalayan catchment" is linked with precipitation rather than mass balance changes of the glacier. This study also indicates that the warming induced an initial increase of glacier runoff and subsequent decline as suggested by the IPCC (2007) is restricted to the glacier degradation-derived component in a precipitation dominant Himalayan catchment and cannot be translated as river flow response. The preliminary assessment suggests that the "Himalayan catchment" could experience higher river flows and positive glacier mass balance regime together in association with strong monsoon. The important role of glaciers in this precipitation dominant system is to augment stream runoff during the years of low summer discharge. This paper intends to highlight the importance of creating credible knowledge on the Himalayan cryospheric processes to develop a more representative global view on river flow response to cryospheric changes and locally sustainable water resources management strategies.
Reduced glacier sliding caused by persistent drainage from a subglacial lake
Authors: Magnússon E. , Björnsson H. , Rott H. , Pálsson F.
Published: April 2010

Abstract. We present velocity observations of a glacier outlet in Vatnajökull, Iceland, deduced from interferometric SAR (InSAR) data obtained during the ERS1/2 tandem mission in 1995–2000. More than a 50% decrease in glacier velocity was observed subsequent to a large jökulhlaup from the subglacial lake Grímsvötn in 1996. The glacier had not reached its former flow rate in 2000. The jökulhlaup damaged the lake's ice-dam causing persistent drainage from the lake. InSAR based studies of water accumulation within Grímsvötn suggest that a leakage of >3 m3 s−1 prevailed throughout our study period. We suggest that the lake leakage kept open a tunnel at low water pressure underneath the whole length of the glacier. The tunnel flow drained water from its surroundings, hence lowering the water pressure of a distributed drainage system, underneath the upper and centre parts of the glacier, which prior to the jökulhlaup sustained significant basal sliding. This is in accordance with theoretical prediction that tunnel flow in a steady state may cause slow-down in glacier velocity by reducing the subglacial water pressure. The width of the affected areas was ~5 km on the upper part of the glacier and ~8 km on the centre part of the glacier. This indicates that the water pressure reduction propagates laterally from the tunnel over a distance of a few km.
Simulation of the specific surface area of snow using a one-dimensional physical snowpack model: implementation and evaluation for subarctic snow in Alaska
Authors: Jacobi H.-W. , Domine F. , Simpson W. R. , Douglas T. A. , Sturm M.
Published: April 2010

Abstract. The specific surface area (SSA) of the snow constitutes a powerful parameter to quantify the exchange of matter and energy between the snow and the atmosphere. However, currently no snow physics model can simulate the SSA. Therefore, two different types of empirical parameterizations of the specific surface area (SSA) of snow are implemented into the existing one-dimensional snow physics model CROCUS. The parameterizations are either based on diagnostic equations relating the SSA to parameters like snow type and density or on prognostic equations that describe the change of SSA depending on snow age, snowpack temperature, and the temperature gradient within the snowpack. Simulations with the upgraded CROCUS model were performed for a subarctic snowpack, for which an extensive data set including SSA measurements is available at Fairbanks, Alaska for the winter season 2003/2004. While a reasonable agreement between simulated and observed SSA values is obtained using both parameterizations, the model tends to overestimate the SSA. This overestimation is more pronounced using the diagnostic equations compared to the results of the prognostic equations. Parts of the SSA deviations using both parameterizations can be attributed to differences between simulated and observed snow heights, densities, and temperatures. Therefore, further sensitivity studies regarding the thermal budget of the snowpack were performed. They revealed that reducing the thermal conductivity of the snow or increasing the turbulent fluxes at the snow surfaces leads to a slight improvement of the simulated thermal budget of the snowpack compared to the observations. However, their impact on further simulated parameters like snow height and SSA remains small. Including additional physical processes in the snow model may have the potential to advance the simulations of the thermal budget of the snowpack and, thus, the SSA simulations.
Reconstruction of the 1979–2006 Greenland ice sheet surface mass balance using the regional climate model MAR
Authors: Fettweis X.
Published: April 2010

Abstract. Results from a 28-year simulation (1979–2006) over the Greenland ice sheet (GIS) reveal an increase of the solid precipitation (+0.4±2.5 km3 yr−2) and the run-off (+7.9±3.3 km3 yr−2) of surface melt water. The net effect of these competing factors leads to a significant Surface Mass Balance (SMB) loss rate of −7.2±5.1 km3 yr−2. The contribution of changes in the net water vapour fluxes (+0.02±0.09 km3 yr−2) and rainfall (+0.2±0.2 km3 yr−2) to the SMB variability is negligible. The melt water supply has increased because the GIS surface has been warming up +2.4°C since 1979. Latent heat flux, sensible heat flux and net solar radiation have not varied significantly over the last three decades. However, the simulated downward infra-red flux has increased by 9.3 W m−2 since 1979. The natural climate variability (e.g. the North Atlantic Oscillation) does not explain these changes on the GIS. The recent global warming, due to the greenhouse gas concentration increase induced by the human activities, could be a cause of these changes. The doubling of the surface melt water flux into the ocean over the period 1979–2006 suggests that the overall ice sheet mass balance has been increasingly negative, given the probable meltwater-induced outlet glacier acceleration. This study suggests that an increased melting dominates over an increased accumulation in a warming scenario and that the GIS would likely continue to loose mass in the future. A GIS melting would have an effect on the stability of the thermohaline circulation (THC) and the global sea level rise.
Improving estimation of glacier volume change: a GLIMS case study of Bering Glacier System, Alaska
Authors: Beedle M. J. , Dyurgerov M. , Tangborn W. , Khalsa S. J. S. , Helm C. , Raup B. , Armstrong R. , Barry R. G.
Published: April 2010

Abstract. The Global Land Ice Measurements from Space (GLIMS) project has developed tools and methods that can be employed by analysts to create accurate glacier outlines and resultant measures of glacier extent. To illustrate the importance of accurate glacier outlines and the effectiveness of GLIMS standards we have conducted a case study on Bering Glacier System (BGS), Alaska. BGS is a complex glacier system aggregated from multiple drainage basins, numerous individual ice streams, and many accumulation areas. Published measurements of BGS surface area vary from 1740 to 6200 km2, depending on how the boundaries of this system have been defined. Utilizing GLIMS tools and standards we have completed a new outline and analysis of the area-altitude distribution (hypsometry) of BGS using Landsat images from 2000 and 2001. We compared this new outline (3632 km2) with three previous outlines to illustrate the errors that result from the widely varying estimates used in previous analysis of BGS area. The use of different BGS outlines results in highly variable measures of volume change and net balance (bn). Outline variability alone results in a net balance rate range of –1.0 to –3.2 m/yr water equivalent (W.E.), a volume change range of –4.2 to –8.2 km3/yr, and a near doubling in contributions to sea level equivalent (SLE), 0.0122 mm/yr to 0.0236 mm/yr. A study of three different models of BGS net balance leads us to favor estimates of bn of –1.2 m/yr W.E. and total volume change of –4.2 km3/yr for the period 1950–2004. These estimates result in a near doubling of contributions to sea level equivalent when compared with previous studies. While current inaccuracies in glacier outlines hinder our ability to fully understand glacier change, there is no reason why our understanding of glacier extents should not be comprehensive and accurate. Such accuracy is possible with the increasing volume of satellite imagery of glacierized regions, and recent advances in tools and standards.
Glacier annual balance measurement, prediction, forecasting and climate correlations, North Cascades, Washington 1984–2006
Authors: Pelto M. S.
Published: April 2010

Abstract. North Cascade glacier annual balance measured on 10 glaciers from 1984–2006 yielded mean annual balance (ba) of –0.54 m/a, and –12.38 m cumulatively. This is a significant loss for glaciers that average 30–60 m in thickness, representing 20–40% of their entire volume. Two observed glaciers, Lewis Glacier and Spider Glacier, no longer exist. The ba of North Cascade glaciers is reliably calculated based on 1 April snowpack water equivalent and ablation season temperature. 1 May forecasting of ba using the Pacific Decadal Oscillation and the Multivariate El Nino Southern Oscillation circulation indices correctly determined the sign of mass balance in 42 of 47 years. Glacier annual balance forecasting is an important step for summer water resource management in glacier runoff dominated stream systems. The forecast for North Cascade glaciers in 2007 is for a negative annual balance.
Modelling historical and recent mass loss of McCall Glacier, Alaska, USA
Authors: Delcourt C. , Pattyn F. , Nolan M.
Published: April 2010

Abstract. Volume loss of valley glaciers is now considered to be a significant contribution to sea level rise. Understanding and identifying the processes involved in accelerated mass loss are necessary to determine their impact on the global system. Here we present results from a series of model experiments with a higher-order thermomechanically coupled flowline model (Pattyn, 2002). Boundary conditions to the model are parameterizations of surface mass balance, geothermal heating, observed surface and 10 m ice depth temperatures. The time-dependent experiments aim at simulating the glacier retreat from its LIA expansion to present according to different scenarios and model parameters. Model output was validated against measurements of ice velocity, ice surface elevation and terminus position at different stages. Results demonstrate that a key factor in determining the glacier retreat history is the importance of internal accumulation (>50%) in the total mass balance. The persistence of a basal temperate zone characteristic for this polythermal glacier depends largely on its contribution. Accelerated glacier retreat since the early nineties seems directly related to the increase in ELA and the sudden reduction in AAR due to the fact that a large lower elevation cirque – previously an important accumulation area – became part of the ablation zone.
Reconstructing the glacier contribution to sea-level rise back to 1850
Authors: Oerlemans J. , Dyurgerov M. , van de Wal R. S. W.
Published: April 2010

Abstract. We present a method to estimate the glacier contribution to sea-level rise from glacier length records. These records form the only direct evidence of glacier changes prior to 1946, when the first systematic mass-balance observations began. A globally representative length signal is calculated from 197 length records from all continents by normalisation and averaging of 14 different regions. Next, the resulting signal is calibrated with mass-balance observations for the period 1961–2000. We find that the glacier contribution to sea level rise was 5.5±1.0 cm during the period 1850–2000 and 4.5±0.7 cm during the period 1900–2000.