bluuebunny/crossref_metadata_2025 huggingface dataset:
Freshwater components and transports in the Fram Strait – recent observations and changes since the late 1990s
Authors: Rabe B. , Schauer U. , Mackensen A. , Karcher M. , Hansen E. , Beszczynska-Möller A.
Published: April 2010

Abstract. We present the late summer distribution and transports of freshwater components in the upper western part of the Fram Strait during 1998, 2004 and 2005. Hydrographic data and and water δ18O values are analyzed to distinguish Atlantic Water, ice-melt (IMW) and freshwater removal from ice formation (IFW), and Meteoric Water (precipitation and riverine sources; MW). Concentrations of these water masses are combined with volume transport estimates from an inverse model. The average liquid freshwater transport relative to a reference salinity of 34.92, was 2500 km3/yr or 80 mSv southward, which is at the upper end of values reported in the literature. Our results indicate that not only the region of the continental slope but also parts of the East Greenland Shelf are important for freshwater transports. The average transports of MW and IFW were 160 mSv (5000 km3/yr) and 90 mSv (2800 km3/yr) southward, respectively. The southward transport of MW was higher in 2005 than in 1998, but was compensated by a higher IFW transport. These differences in transports were associated with stronger southward velocities and the absence of northward velocities over the continental slope and the eastern East Greenland Shelf in 2005. A simulation using the North Atlantic-Arctic Ocean Sea Ice Model (NAOSIM) shows that the high transport of MW in the Fram Strait in 2005 is in agreement with the temporary storage of river water on the Siberian shelf in the mid-1990s, which reached the north of Greenland in 2003. Our results indicate that IFW follows the same pathways as MW before reaching the Fram Strait.
On the permeability of barrier layers
Authors: Mignot J. , de Boyer Montégut C. , Tomczak M.
Published: April 2010

Abstract. We present a revisited 2-degree resolution global ocean climatology of monthly mean Barrier Layer Thickness (BLT) first proposed by de Boyer Montégut et al. (2007a). In addition to using an extended data set, we present a modified computation method in order to take into account the observed permeability of Barrier Layers (BL). We name permeability the fact that barrier layers can, in some areas, be very patchy regarding the space and time scales that are considered. This patchiness may have important consequences for the climatic impact of BLs. Differences between the two computation methods are weak for robust BLs that are formed by large-scale processes. The former approach can however largely underestimate the thickness of short and/or localized barrier layers. This is especially the case for the ones formed by mesoscale mechanisms (under the ITCZ for example and along western boundary currents) and patchy barrier layers detected equatorward of the sea surface salinity subtropical maxima. Complete characterisation of regional BL dynamics therefore requires the description of BL robustness through the determination of both BLT and BL permeability.
Observing using sound and light – a short review of underwater acoustic and video-based methods
Authors: Jonsson P. , Sillitoe I. , Dushaw B. , Nystuen J. , Heltne J.
Published: April 2010

Abstract. This paper is a review which briefly describes a selection of acoustic observation techniques and certain aspects of underwater video technology suitable for observations in an underwater environment. The review is divided into two sections, one for each subject, where each section concludes with a discussion of the current challenges within the respective fields. The acoustic section of the review covers bathymetric and geometrical measurements, imaging sonars, subsurface penetrating profilers, positioning methods, acoustic underwater communication and sensor networks, and water speed measurements. The section ends by considering temperature measurements by ocean acoustic tomography and passive acoustic monitoring. The underwater video section initially deals with questions of acquisition including underwater visibility, the type of platform, and video formats, image sensors and specialized cameras. This is followed by notes on processing techniques including mosaicking, stereo video, structured light, recording and transmission, image enhancement techniques and ends with a short discussion of underwater holographic cameras.
Operational oceanography in support to indicator reporting
Authors: Johannessen J. A. , Coppini G. , Soulat F. , Larnicol G.
Published: April 2010

Abstract. Operational Oceanography (OO) has now emerged to a stage that allows the design, development and execution of marine core services tailored to user requirements. As such it is also feasible to provide routine production of environmental and climate indicators. Indicators are synthetic indices of environmental changes at various temporal and spatial scales. In this paper we outline the possible contribution and strengthening of existing indicator reporting based on OO products followed by a discussion of the relevance of such improved reporting for marine environmental policy implementation and regulation. In particular, it capitalizes on the main achievements of the Marine Environment and Security of the European Area (MERSEA) project, the outcome of a European Marine Monitoring and Assessment (EMMA) workshop on the connection between operational oceanography and the European Marine Strategy (EMS) Directive and the regular European Environmental Agency (EEA) assessment reports.
The role of continental shelves in nitrogen and carbon cycling
Authors: Fennel K.
Published: April 2010

Abstract. Continental shelves play a key role in the cycling of nitrogen and carbon. Here the physical transport and biogeochemical transformation processes affecting the fluxes into and out of continental shelf systems are reviewed, and their role in the global cycling of both elements is discussed. Uncertainties in observation-based estimates of nitrogen and carbon fluxes mostly result from uncertainties in the shelf-open ocean exchange of organic and inorganic matter, which is hard to quantify based on observations alone, but can be inferred from biogeochemical models. Model-based nitrogen and carbon budgets are presented for the Northwestern North Atlantic continental shelf. Results indicate that shelves are an important sink for fixed nitrogen and a source of alkalinity, but are not much more efficient in exporting organic carbon to the deep ocean than the adjacent open ocean for the shelf region considered.
Automated gas bubble imaging at sea floor – a new method of in situ gas flux quantification
Authors: Thomanek K. , Zielinski O. , Sahling H. , Bohrmann G.
Published: April 2010

Abstract. Photo-optical systems are common in marine sciences and have been extensively used in coastal and deep-sea research. However, due to technical limitations in the past photo images had to be processed manually or semi-automatically. Recent advances in technology have rapidly improved image recording, storage and processing capabilities which are used in a new concept of automated in situ gas quantification by photo-optical detection. The design for an in situ high-speed image acquisition and automated data processing system is reported ("Bubblemeter"). New strategies have been followed with regards to back-light illumination, bubble extraction, automated image processing and data management. This paper presents the design of the novel method, its validation procedures and calibration experiments. The system will be positioned and recovered from the sea floor using a remotely operated vehicle (ROV). It is able to measure bubble flux rates up to 10 L/min with a maximum error of 33% for worst case conditions. The Bubblemeter has been successfully deployed at a water depth of 1023 m at the Makran accretionary prism offshore Pakistan during a research expedition with R/V Meteor in November 2007.
Trends in coastal upwelling intensity during the late 20th century
Authors: Narayan N. , Paul A. , Mulitza S. , Schulz M.
Published: April 2010

Abstract. This study presents linear trends of coastal upwelling intensity in the later part of the 20th century (1960–2001) employing various indices of upwelling, derived from meridional wind stress and sea surface temperature. The analysis was conducted in the four major coastal upwelling regions in the world, which are off North-West Africa, Lüderitz, California and Peru respectively. The trends in meridional wind stress showed a steady increase of intensity from 1960–2001, which was also reflected in the SST index calculated for the same time period. The steady cooling observed in the instrumental records of SST off California substantiated this observation further. Correlation analysis showed that basin-scale oscillations like the Atlantic Multidecadal Oscillation (AMO) and the Pacific Decadal Oscillation (PDO) could not be directly linked to the observed increase of upwelling intensity off NW Africa and California respectively. The relationship of the North Atlantic Oscillation (NAO) with coastal upwelling off NW Africa turned out to be ambiguous due to a negative correlation between the NAO index and the meridional wind stress and a lack of correlation with the SST index. Our results give additional support to the hypothesis that the coastal upwelling intensity increases globally because of raising greenhouse gas concentrations in the atmosphere and an associated increase of the land-sea pressure gradient and meridional wind stress.
Seismic imaging of a thermohaline staircase in the western tropical North Atlantic
Authors: Fer I. , Nandi P. , Holbrook W. S. , Schmitt R. W. , Páramo P.
Published: April 2010

Abstract. Multichannel seismic data acquired in the Lesser Antilles in the western tropical North Atlantic indicate that the seismic reflection method has imaged an oceanic thermohaline staircase. Synthetic modeling of observed density and sound speed profiles corroborates inferences from the seismic imagery. Laterally coherent, uniform layers are present at depths ranging from 550–700 m and have a separation of ~20 m, with thicknesses increasing with depth. Reflection coefficient, a measure of the acoustic impedance contrasts, associated with the interfaces is one order of magnitude greater than the background levels. Hydrography sampled in previous surveys puts a constraint on the longevity of these layers in this area to within a maximum of three years. Spectral analysis of layer horizons in the thermohaline staircase indicates that internal wave activity is anomalously low, suggesting weak internal wave-induced turbulence and mixing. Results from two independent measurements, the application of a finescale parameterization to observed high-resolution velocity profiles and direct measurements of turbulent dissipation rate, confirm the low levels of turbulence and mixing. The lack of internal wave-induced mixing allows for the maintenance of the staircase. Our observations show the potential that seismic oceanography can contribute to an improved understanding of temporal occurrence rates, and the geographical distribution of thermohaline staircases and can improve current estimates of vertical mixing rates ascribable to salt fingering in the global ocean.
Upper layer current variability in the Central Ligurian Sea
Authors: Picco P. , Cappelletti A. , Sparnocchia S. , Schiano M. E. , Pensieri S. , Bozzano R.
Published: April 2010

Abstract. Long-time series of surface currents and meteorological parameters were analysed to estimate the variability of the upper layer circulation as a preliminary study of the Ligurian Air-Sea Interaction Experiment (LASIE07). Current meter data were collected by an upward-looking RDI Sentinel 300 kHz ADCP deployed in the Central Ligurian Sea (43°47.77' N; 9°02.85' E) near the meteo-oceanographic buoy ODAS ITALIA1 for over eight months. The ADCP sampled the upper 50 m of water column at 8 m vertical resolution and 1 h time interval; surface marine and atmospheric hourly data were provided by the buoy. Currents were mainly barotropic and directed NW, according to the general circulation of the area, had a mean velocity of about 18 cm s−1 and hourly mean peaks up to 80 m s−1. Most of the observed variability in the upper thermocline was determined by inertial currents and mesoscale activity due to the presence of the Ligurian Front. Local wind had a minor role in the near-surface circulation but induced internal waves propagating downward in the water column.
Temporal energy partitions of Florida extreme sea level events as a function of Atlantic multidecadal oscillation
Authors: Park J. , Obeysekera J. , Barnes J.
Published: April 2010

Abstract. An energy-conservative metric based on the discrete wavelet transform is applied to assess the relative energy distribution of non-stationary extreme sea level events across different temporal scales. The metric is applied to coastal events at Key West and Pensacola Florida as a function of two Atlantic Multidecadal Oscillation (AMO) regimes. Under AMO warm conditions there is a small but significant redistribution of event energy from nearly static into more dynamic timescales at Key West, while at Pensacola the AMO-dependent changes in temporal event behaviour are less pronounced. Extreme events with increased temporal dynamics are consistent with an increase in total energy of event forcings which may be a reflection of more energetic storm events during AMO warm phases. As dynamical models mature to the point of providing regional climate index predictability, coastal planners may be able to consider such temporal change metrics in planning scenarios.
On the variability of Florida Straits and wind driven transports at 26° N in the Atlantic Ocean
Authors: Atkinson C. P. , Bryden H. L. , Hirschi J. J.-M. , Kanzow T.
Published: April 2010

Abstract. Since April 2004 the RAPID array has made continuous measurements of the Atlantic Meridional Overturning Circulation (AMOC) at 26° N. Two key components of this system are Ekman transport zonally integrated across 26° N and western boundary current transport in the Florida Straits. Whilst measurements of the AMOC as a whole are somewhat in their infancy, this study investigates what useful information can be extracted on the variability of the Ekman and Florida Straits transports using the decadal timeseries already available. Analysis is also presented for Sverdrup transports zonally integrated across 26° N. The seasonal cycles of Florida Straits, Ekman and Sverdrup transports are quantified at 26° N using harmonic analysis of annual and semi-annual constituents. Whilst Sverdrup transport shows clear semi-annual periodicity, calculations of seasonal Florida Straits and Ekman transports show substantial interannual variability due to variability at non-seasonal frequencies; the mean seasonal cycle for these transports only emerges from decadal length observations. The Florida Straits and Ekman mean seasonal cycles project on the AMOC with a combined peak-to-peak seasonal range of 3.5 Sv. The combined seasonal range for heat transport is 0.40 PW. The Florida Straits seasonal cycle possesses a smooth annual periodicity in contrast with previous studies suggesting a more asymmetric structure. No clear evidence is found to support significant changes in the Florida Straits seasonal cycle at sub-decadal periods. Whilst evidence of wind driven Florida Straits transport variability is seen at sub-seasonal and annual periods, model runs from the 1/4° eddy-permitting ocean model NEMO are used to identify an important contribution from internal oceanic variability at sub-annual and interannual periods. The Ekman transport seasonal cycle possesses less symmetric structure, due in part to different seasonal transport regimes east and west of 50 to 60° W. Around 60% of non-seasonal Ekman transport variability occurs in phase section-wide at 26° N and is related to the NAO, whilst Sverdrup transport variability is more difficult to decompose.
Earth's surface heat flux
Authors: Davies J. H. , Davies D. R.
Published: April 2010

Abstract. We present a revised estimate of Earth's surface heat flux that is based upon a heat flow data-set with 38 347 measurements, which is 55% more than used in previous estimates. Our methodology, like others, accounts for hydrothermal circulation in young oceanic crust by utilising a half-space cooling approximation. For the rest of Earth's surface, we estimate the average heat flow for different geologic domains as defined by global digital geology maps; and then produce the global estimate by multiplying it by the total global area of that geologic domain. The averaging is done on a polygon set which results from an intersection of a 1 degree equal area grid with the original geology polygons; this minimises the adverse influence of clustering. These operations and estimates are derived accurately using methodologies from Geographical Information Science. We consider the virtually un-sampled Antarctica separately and also make a small correction for hot-spots in young oceanic lithosphere. A range of analyses is presented. These, combined with statistical estimates of the error, provide a measure of robustness. Our final preferred estimate is 47±2 TW, which is greater than previous estimates.
Particle size distributions by laser diffraction – Part 1: Sensitivity of granular matter strength to analytical operating procedures
Authors: Storti F. , Balsamo F.
Published: April 2010

Abstract. We tested laser diffraction particle size analysis in poorly coherent carbonate platform cataclastic breccias and unfaulted quartz-rich eolian sands, representing low- and high-strength granular materials, respectively. We used two different instruments with different sample dispersion and pumping systems and several wet and dry analytical procedures that included different pump speeds, measure precision tests with and without sample ultrasonication, and different dispersant liquids. Results of our work indicate that high strength material is not strongly affected by analytical operating procedures, whereas low strength material is very sensitive to the pump speed, ultrasonication intensity, and measurement run time. To reduce such a data variability, we propose a workflow for analytical tests preliminary to the set up of the most appropriate SOP.
Candidates for multiple impact craters: popigai and chicxulub as seen by EGM08, a global 5'×5' gravitational model
Authors: Klokočník J. , Kostelecký J. , Pešek I. , Novák P. , Wagner C. A. , Sebera J.
Published: April 2010

Abstract. In 2008 the new Earth Gravitational Model (EGM08) was released. It contains a complete set of spherical harmonic coefficients of the Earth's gravitational potential (Stokes parameters) to degree 2190 and order 2159 that can be used for evaluation of various potential quantities with both the unprecedented accuracy and high spatial resolution. Two such quantities, the gravity anomaly and second-order radial derivative of the disturbing potential, were computed over selected areas with known impact craters. The displays of these derivatives for two such sites clearly show not only the strong circular-like features known to be associated with them but also other symmetrical structures which appear to make them multiple impact sites. At Popigai, Siberia, the secondary circular features fall in a line from the primary in the SE direction. At Chicxulub, Yucatán, there appears to be one secondary crater close to the primary in the NE direction, as well as possibly others in the vicinity of the main crater. Gravity information alone is not proof of the impact craters but it is useful in identifying candidate sites for further study, for future examination by geologists and geophysicists.
Detrital zircon geochronologic tests of the SE Siberia-SW Laurentia paleocontinental connection
Authors: MacLean J. S. , Sears J. W. , Chamberlain K. R. , Khudoley A. K. , Prokopiev A. V. , Kropachev A. P. , Serkina G. G.
Published: April 2010

Abstract. Strikingly similar Late Mesoproterozoic stratigraphic sequences and correlative U-Pb detrital-zircon ages may indicate that the Sette Daban region of southeastern Siberia and the Death Valley region of southwestern North America were formerly contiguous parts of a Grenville foreland basin. The Siberian section contains large numbers of detrital zircons that correlate with Grenville, Granite-Rhyolite, and Yavapai basement provinces of North America. The sections in both Siberia and Death Valley exhibit west-directed thrust faults that may represent remnants of a Grenville foreland thrust belt. North American detrital-zircon components do not occur in Siberian samples above a ~600 Ma breakup unconformity, suggesting that rifting and continental separation blocked transfer of clastic sediment between the cratons by 600 Ma. Faunal similarities suggest, however, that the two cratons remained within the breeding ranges of Early Cambrian trilobites and archeocyathans.
Glacier annual balance measurement, 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, 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, correlation coefficient 0.91, using 1 April snowpack water equivalent and ablation season temperature. Utilizing ba from 10 glaciers 1984–2006 and net balance (bn) from South Cascade 1960–2005, a set of forecast rules for glacier mass balance were derived utilizing October–April Pacific Decadal Oscillation and Multivariate El Nino Southern Oscillation index values. The forecast rules provide a correct assessment in 41 of the 46 years for South Cascade Glacier and 20 of 23 years for NCGCP glaciers. 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 ba.
Landsat-derived glacier inventory for Jotunheimen, Norway, and deduced glacier changes since the 1930s
Authors: Andreassen L. M. , Paul F. , Kääb A. , Hausberg J. E.
Published: April 2010

Abstract. A Landsat Thematic Mapper (TM) scene from 2003 covering the Jotunheimen and Breheimen region has been used to map the recent glacier extents using thresholded ratio images (TM3/TM5). Orthoprojected aerial photographs and glacier outlines from digital maps have been used to validate the method and control the results. We further calculated glacier changes by comparing the Landsat-derived 2003 glacier outlines with previous maps and inventories from the 1930s, 1960s and 1980s. Our results confirm that the applied automatic mapping method is robust and agrees with the reference data used. Some manual editing was necessary to correct the outline at ice-lake contacts and at debris covered glaciers. However, for most of the glaciers no corrections were required. The most laborious task has been to assign ID numbers and couple the new Landsat inventory to previous inventories to assess area changes. The glaciers investigated shrank since the 1930s with an overall area reduction of about 23% for 38 glaciers. Since the 1960s the area reduction was 12% for 164 glaciers. Although the general trend is glacier recession, some glaciers have increased their size or remained nearly unchanged over these decades.
The equilibrium flow and mass balance of the Taku Glacier, Alaska 1950–2006
Authors: Pelto M. S. , Miller M. M. , Adema G. W. , Beedle M. J. , McGee S. R. , Sprenke K. F. , Lang M.
Published: April 2010

Abstract. The Taku Glacier, Alaska has advanced 7.5 km since the late nineteenth century, while all other primary outlet glaciers of the Juneau Icefield are in retreat. The Juneau Icefield Research Program has completed field work on the Taku Glacier annually since 1946. The collected observations of surface mass balance, glacier velocity and glacier thickness at Profile IV 29 km above the terminus and 4 km above the equilibrium line provide a means to assess the equilibrium nature of the Taku Glacier. Annual velocity measured and summer velocity measurements completed at a Profile IV from 1950–2006 indicate insignificant variations in velocity seasonally or from year to year. The consistency of velocity over the 56-year period indicates that in the vicinity of the equilibrium line, the flow of the Taku Glacier has been in an equilibrium state. Surface mass balance was positive from 1946–1988 averaging +0.42 m a−1. This led to glacier thickening. From 1988–2006 an important change has occurred and annual balance has been −0.14 m a−1, and the glacier thickness has ceased increasing along Profile IV. Field measurements of ice depth and surface velocity allow calculation of the volume flux at Profile IV. Volume flux is then compared with the surface balance flux from the region of the glacier above Profile IV, determined annually in the field. Above Profile IV the observed mean surface flux from 1950–2006 is 5.50×108 m3 a−1 (±5%), while the calculated volume flux range for the same period flowing through profile IV is 5.00–5.47×108 m3 a−1. The mean surface flux has been greater than the volume flux, which has led to slow thickening of the Taku Glacier up to 1988. The thickening has not led to a change in the flow of Taku Glacier at Profile IV.
On the limit to resolution and information on basal properties obtainable from surface data on ice streams
Authors: Gudmundsson G. H. , Raymond M.
Published: April 2010

Abstract. An optimal estimation method for simultaneously determining both basal slipperiness and basal topography from variations in surface flow velocity and topography along a flow line on ice streams and ice sheets is presented. We use Bayesian inference to update prior statistical estimates for basal topography and slipperiness using surface measurements along a flow line. Our main focus here is on how errors and spacing of surface data affect estimates of basal quantities and on possibly aliasing/mixing between basal slipperiness and basal topography. We find that the effects of spatial variations in basal topography and basal slipperiness on surface data can be accurately separated from each other, and mixing in retrieval does not pose a serious problem. For realistic surface data errors and density, small-amplitude perturbations in basal slipperiness can only be resolved for wavelengths larger than about 50 times the mean ice thickness. Bedrock topography is well resolved down to horizontal scale equal to about one ice thickness. Estimates of basal slipperiness are not significantly improved by accurate prior estimates of basal topography. However, retrieval of basal slipperiness is found to be highly sensitive to unmodelled errors in basal topography.
Partitioning of melt energy and meltwater fluxes in the ablation zone of the west Greenland ice sheet
Authors: van den Broeke M. , Smeets P. , Ettema J. , van der Veen C. , van de Wal R. , Oerlemans J.
Published: April 2010

Abstract. We present four years (August 2003–August 2007) of surface mass balance data from the ablation zone of the west Greenland ice sheet along the 67° N latitude circle. Sonic height rangers and automatic weather stations continuously measured accumulation/ablation and near-surface climate at distances of 6, 38 and 88 km from the ice sheet margin at elevations of 490, 1020 and 1520 m a.s.l. Using a melt model and reasonable assumptions about snow density and percolation characteristics, these data are used to quantify the partitioning of energy and mass fluxes during melt episodes. The lowest site receives very little winter accumulation, and ice melting is nearly continuous in June, July and August. Due to the lack of snow accumulation, little refreezing occurs and virtually all melt energy is invested in runoff. Higher up the ice sheet, the ice sheet surface freezes up during the night, making summer melting intermittent. At the intermediate site, refreezing in snow consumes about 10% of the melt energy, increasing to 40% at the highest site. The sum of these effects is that total melt and runoff increase exponentially towards the ice sheet margin, each time doubling between the stations. At the two lower sites, we estimate that radiation penetration causes 20–30% of the ice melt to occur below the surface.
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. To illustrate the importance of accurate glacier outlines and the effectiveness of GLIMS standards we conducted a case study on Bering Glacier System (BGS), Alaska. BGS is a complex glacier system aggregated from multiple drainage basins, numerous tributaries, 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 (3630 km2) and analysis of the area-altitude distribution (hypsometry) of BGS using Landsat images from 2000 and 2001 and a US Geological Survey 15-min digital elevation model. We compared this new hypsometry with three different hypsometries to illustrate the errors that result from the widely varying estimates of BGS extent. The use of different BGS hypsometries results in highly variable measures of volume change and net balance (bn). Applying a simple hypsometry-dependent mass-balance model to different hypsometries results in a bn rate range of −1.0 to −3.1 m a−1 water equivalent (W.E.), a volume change range of −3.8 to −6.7 km3 a−1 W.E., and a near doubling in contributions to sea level equivalent, 0.011 mm a−1 to 0.019 mm a−1. Current inaccuracies in glacier outlines hinder our ability to correctly quantify glacier change. Understanding of glacier extents can become comprehensive and accurate. Such accuracy is possible with the increasing volume of satellite imagery of glacierized regions, recent advances in tools and standards, and dedication to this important task.
Antarctic summer sea ice concentration and extent: comparison of ODEN 2006 ship observations, satellite passive microwave and NIC sea ice charts
Authors: Ozsoy-Cicek B. , Xie H. , Ackley S. F. , Ye K.
Published: April 2010

Abstract. Antarctic sea ice cover has shown a slight increase (<1%/decade) in overall observed ice extent as derived from satellite mapping from 1979 to 2008, contrary to the decline observed in the Arctic regions. Spatial and temporal variations of the Antarctic sea ice however remain a significant problem to monitor and understand, primarily due to the vastness and remoteness of the region. While satellite remote sensing has provided and has great future potential to monitor the variations and changes of sea ice, uncertainties remain unresolved. In this study, the National Ice Center (NIC) ice edge and the AMSR-E (Advanced Microwave Scanning Radiometer-Earth Observing System) ice extent are examined, while the ASPeCt (Antarctic Sea Ice Process and Climate) ship observations from the Oden expedition in December 2006 are used as ground truth to verify the two products during Antarctic summer. While there is a general linear trend between ASPeCt and AMSR-E ice concentration estimates, there is poor correlation (R2=0.41) and AMSR-E tends to underestimate the low ice concentrations. We also found that the NIC sea ice edge agrees well with ship observations, while the AMSR-E shows the ice edge further south, consistent with its poorer detection of low ice concentrations. The northward extent of the ice edge at the time of observation (NIC) had mean values varying from 38 km to 102 km greater on different days for the area as compared with the AMSR-E sea ice extent. For the circumpolar area as a whole in the December period examined, AMSR-E therefore may underestimate the area inside the ice edge at this time by up to 14% or, 1.5 million km2 less area, compared to the NIC ice charts. Preliminary comparison of satellite scatterometer data however, suggests better resolution of low concentrations than passive microwave, and therefore better agreement with ship observations and NIC charts of the area inside the ice edge during Antarctic summer. A reanalysis data set for Antarctic sea ice extent that relies on the decade long scatterometer and high resolution satellite data set, instead of passive microwave, may therefore give better fidelity for the recent sea ice climatology.
The emergence of surface-based Arctic amplification
Authors: Serreze M. C. , Barrett A. P. , Stroeve J. C. , Kindig D. N. , Holland M. M.
Published: April 2010

Abstract. Rises in surface and lower troposphere air temperatures through the 21st century are projected to be especially pronounced over the Arctic Ocean during the cold season. This Arctic amplification is largely driven by loss of the sea ice cover, allowing for strong heat transfers from the ocean to the atmosphere. Consistent with observed reductions in sea ice extent, fields from both the NCEP/NCAR and JRA-25 reanalyses point to emergence of surface-based Arctic amplification in the last decade.
The role of radiation penetration in the energy budget of the snowpack at Summit, Greenland
Authors: Kuipers Munneke P. , van den Broeke M. R. , Reijmer C. H. , Helsen M. M. , Boot W. , Schneebeli M. , Steffen K.
Published: April 2010

Abstract. Measurements of the summer surface energy balance at Summit, Greenland, are presented (8 June–20 July 2007). These measurements serve as input to an energy balance model that searches for a surface temperature for which closure of all energy terms is achieved. A good agreement between observed and modelled surface temperatures was found, with an average difference of 0.45°C and an RMSE of 0.85°C. It turns out that penetration of shortwave radiation into the snowpack plays a small but important role in correctly simulating snow temperatures. After 42 days, snow temperatures in the first meter are 3.6–4.0°C higher compared to a model simulation without radiation penetration. Sensitivity experiments show that these results cannot be reproduced by tuning the heat conduction process alone, by varying snow density or snow diffusivity. We compared the two-stream radiation penetration calculations with a sophisticated radiative transfer model and discuss the differences. The average diurnal cycle shows that net shortwave radiation is the largest energy source (diurnal average of +61 W m−2), net longwave radiation the largest energy sink (−42 W m−2). On average, subsurface heat flux, sensible and latent heat fluxes are the remaining, small heat sinks (−5, −5 and −7 W m−2, respectively), although these are more important on a subdaily timescale.
Measurement of the specific surface area of snow using infrared reflectance in an integrating sphere at 1310 and 1550 nm
Authors: Gallet J.-C. , Domine F. , Zender C. S. , Picard G.
Published: April 2010

Abstract. Even though the specific surface area (SSA) and the snow area index (SAI) of snow are crucial variables to determine the chemical and climatic impact of the snow cover, few data are available on the subject. We propose here a novel method to measure snow SSA and SAI. It is based on the measurement of the hemispherical infrared reflectance of snow samples using the DUFISSS instrument (DUal Frequency Integrating Sphere for Snow SSA measurement). DUFISSS uses the 1310 or 1550 nm radiation of laser diodes, an integrating sphere 15 cm in diameter, and InGaAs photodiodes. For SSA<60 m2 kg−1, we use the 1310 nm radiation, reflectance is between 15 and 50% and the accuracy of SSA determination is 10%. For SSA>60 m2 kg−1, snow is usually of low density (typically 30 to 100 kg m−3), resulting in insufficient optical depth and 1310 nm radiation reaches the bottom of the sample, causing artifacts. The 1550 nm radiation is therefore used for SSA>60 m2 kg−1. Reflectance is then in the range 5 to 12% and the accuracy on SSA is 12%. We propose empirical equations to determine SSA from reflectance at both wavelengths, with that for 1310 nm taking into account the snow density. DUFISSS has been used to measure the SSA of snow and the SAI of snowpacks in polar and Alpine regions.
Measured and modelled sublimation on the tropical Glaciar Artesonraju, Perú
Authors: Winkler M. , Juen I. , Mölg T. , Wagnon P. , Gómez J. , Kaser G.
Published: April 2010

Abstract. Sublimation plays a decisive role in the surface energy and mass balance of tropical glaciers. During the dry season (May–September) low specific humidity and high surface roughness favour the direct transition from ice to vapour and drastically reduce the energy available for melting. However, field measurements are scarce and little is known about the performance of sublimation parameterisations in glacier mass balance and runoff models. During 15 days in August 2005 sublimation was measured on the tongue of Glaciar Artesonraju (8°58' S, 77°38' W) in the Cordillera Blanca, Perú, using simple lysimeters. Indicating a strong dependence on surface roughness, daily totals of sublimation range from 1–3 kg m−2 for smooth to 2–5 kg m−2 for rough conditions. (The 15-day means at that time of wind speed and specific humidity were 4.3 m s−1 and 3.8 g kg−1, respectively.) Measured sublimation was related to characteristic surface roughness lengths for momentum (zm) and for the scalar quantities of temperature and water vapour (zs), using a process-based mass balance model. Input data were provided by automatic weather stations, situated on the glacier tongue at 4750 m a.s.l. and 4810 m a.s.l., respectively. Under smooth conditions the combination zm=2.0 mm and zs=1.0 mm appeared to be most appropriate, for rough conditions zm=20.0 mm and zs=10.0 mm fitted best. Extending the sublimation record from April 2004 to December 2005 with the process-based model confirms, that sublimation shows a clear seasonality. 60–90% of the energy available for ablation is consumed by sublimation in the dry season, but only 10–15% in the wet season (October–April). The findings are finally used to evaluate the parameterisation of sublimation in the lower-complexity mass balance model ITGG, which has the advantage of requiring precipitation and air temperature as only input data. It turns out that the implementation of mean wind speed is a possible improvement for the representation of sublimation in the ITGG model.
The annual surface energy budget of a high-arctic permafrost site on Svalbard, Norway
Authors: Westermann S. , Lüers J. , Langer M. , Piel K. , Boike J.
Published: April 2010

Abstract. Independent measurements of radiation, sensible and latent heat fluxes and the ground heat flux are used to describe the annual cycle of the surface energy budget at a high-arctic permafrost site on Svalbard. During summer, the net short-wave radiation is the dominant energy source, while well developed turbulent processes and the heat flux in the ground lead to a cooling of the surface. About 15% of the net radiation is consumed by the seasonal thawing of the active layer in July and August. The Bowen ratio is found to vary between 0.25 and 2, depending on water content of the uppermost soil layer. During the polar night in winter, the net long-wave radiation is the dominant energy loss channel for the surface, which is mainly compensated by the sensible heat flux and, to a lesser extent, by the ground heat flux, which originates from the refreezing of the active layer. The average annual sensible heat flux of −6.9 Wm−2 is composed of strong positive fluxes in July and August, while negative fluxes dominate during the rest of the year. With 6.8 Wm−2, the latent heat flux more or less compensates the sensible heat flux in the annual average. Strong evaporation occurs during the snow melt period and particularly during the snow-free period in summer and fall. When the ground is covered by snow, latent heat fluxes through sublimation of snow are recorded, but are insignificant for the average surface energy budget. The near-surface atmospheric stratification is found to be predominantly unstable to neutral, when the ground is snow-free, and stable to neutral for snow-covered ground. Due to long-lasting near-surface inversions in winter, an average temperature difference of approximately 3 K exists between the air temperature at 10 m height and the surface temperature of the snow. As such comprehensive data sets are sparse for the Arctic, they are of great value to improve process understanding and support modeling efforts on the present-day and future arctic climate and permafrost conditions.
Estimating basal properties of ice streams from surface measurements: a non-linear Bayesian inverse approach applied to synthetic data
Authors: Raymond M. J. , Gudmundsson G. H.
Published: April 2010

Abstract. We propose a new approach to indirectly estimate basal properties of ice streams, i.e. bedrock topography and basal slipperiness, from observations of surface topography and surface velocities. We demonstrate how a maximum a posteriori estimate of basal conditions can be determined using a Bayesian inference approach in a combination with an analytical linearisation of the forward model. Using synthetic data we show that for non-linear media and non-linear sliding law only a few forward-step model evaluations are needed for convergence. The forward step is solved with a numerical finite-element model using the full Stokes equations. The Fréchet derivative of the forward function is approximated through analytical small-perturbation solutions. This approximation is a key feature of the method and the effects of this approximation on model performance are analyzed. The number of iterations needed for convergence increases with the amplitude of the basal perturbations, but generally less than ten iterations are needed.
An efficient regional energy-moisture balance model for simulation of the Greenland Ice Sheet response to climate change
Authors: Robinson A. , Calov R. , Ganopolski A.
Published: April 2010

Abstract. In order to explore the response of the Greenland ice sheet (GIS) to climate change on long (centennial to multi-millennial) time scales, a regional energy-moisture balance model has been developed. This model simulates seasonal variations of temperature and precipitation over Greenland and explicitly accounts for elevation and albedo feedbacks. From these fields, the annual mean surface temperature and surface mass balance can be determined and used to force an ice sheet model. The melt component of the surface mass balance is computed here using both a positive degree day approach and a more physically-based alternative that includes insolation and albedo explicitly. As a validation of the climate model, we first simulated temperature and precipitation over Greenland for the prescribed, present-day topography. Our simulated climatology compares well to observations and does not differ significantly from that of a simple parameterization used in many previous simulations. Furthermore, the calculated surface mass balance using both melt schemes falls within the range of recent regional climate model results. For a prescribed, ice-free state, the differences in simulated climatology between the regional energy-moisture balance model and the simple parameterization become significant, with our model showing much stronger summer warming. When coupled to a three-dimensional ice sheet model and initialized with present-day conditions, the two melt schemes both allow realistic simulations of the present-day GIS.
Geometric changes and mass balance of the Austfonna ice cap, Svalbard
Authors: Moholdt G. , Hagen J. O. , Eiken T. , Schuler T. V.
Published: April 2010

Abstract. The dynamics and mass balance regime of the Austfonna ice cap, the largest glacier on Svalbard, deviates significantly from most other glaciers in the region and is not fully understood. We have compared ICESat laser altimetry, airborne laser altimetry, GNSS surface profiles and radio echo-sounding data to estimate elevation change rates for the periods 1983–2007 and 2002–2008. The data sets indicate a pronounced interior thickening of up to 0.5 m y−1, at the same time as the margins are thinning at a rate of 1–3 m y−1. The southern basins are thickening at a higher rate than the northern basins due to a higher accumulation rate. The overall volume change in the 2002–2008 period is estimated to be −1.3±0.5 km3 w.e. y−1 (or −0.16±0.06 m w.e. y−1) where the entire net loss is due to a rapid retreat of the calving fronts. Since most of the marine ice loss occurs below sea level, Austfonna's current contribution to sea level change is close to zero. The geodetic results are compared to in-situ mass balance measurements which indicate that the 2004–2008 surface net mass balance has been slightly positive (0.05 m w.e. y−1) though with large annual variations. Similarities between local net mass balances and local elevation changes indicate that most of the ice cap is slow-moving and not in dynamic equilibrium with the current climate. More knowledge is needed about century-scale dynamic processes in order to predict the future evolution of Austfonna based on climate scenarios.
High resolution modelling of snow transport in complex terrain using downscaled MM5 wind fields
Authors: Bernhardt M. , Liston G. E. , Strasser U. , Zängl G. , Schulz K.
Published: April 2010

Abstract. Snow transport is one of the most dominant processes influencing the snow cover accumulation and ablation in high mountain environments. Hence, the spatial and temporal variability of the snow cover is significantly modified with respective consequences on the total amount of water in the snow pack, on the temporal dynamics of the runoff and on the energy balance of the surface. For the present study we used the snow transport model SnowModel in combination with MM5 (Penn State University – National Center for Atmospheric Research MM5 model) generated wind fields. In a first step the MM5 wind fields were downscaled by using a semi-empirical approach which accounts for the elevation difference of model and real topography, and vegetation. The target resolution of 30 m corresponds to the resolution of the best available DEM and land cover map of the test site Berchtesgaden National Park. For the numerical modelling, data of six automatic meteorological stations were used, comprising the winter season (September–August) of 2003/04 and 2004/05. In addition we had automatic snow depth measurements and periodic manual measurements of snow courses available for the validation of the results. It could be shown that the model performance of SnowModel could be improved by using downscaled MM5 wind fields for the test site. Furthermore, it was shown that an estimation of snow transport from surrounding areas to glaciers becomes possible by using downscaled MM5 wind fields.