past rates of accumulation in central west antarctica 1 the spatial variation in ice accumulation is critical to the form and stability of ice sheets in west antarctica knowledge of present ice accumulation rates is restricted to sparse measurements and interpolations across data free regions while information on former ice accumulation is restricted to a small number of ice cores and numerical modeling results here isochronous internal ice sheet layers recorded by airborne radio echo sounding and dated at their stratigraphic position at the byrd ice core site are used to reconstruct the spatial pattern of accumulation across the main west antarctic ice divide for the past 16 000 years mean accumulation calculated for the last 3100 years matches well to present day accumulation rates which is compatible with previous observations on siple dome accumulation between 6400 and 16 000 years ago is calculated to be around half that at present in accordance with recent ice sheet modeling results index terms 1223 geodesy and gravity ocean earth atmosphere interactions 3339 1827 hydrology glaciology 1863 9310 information related to geographic region antarctica 9604 information related to geologic time cenozoic 3349 meteorology and atmospheric dynamics polar meteorology citation siegert m j and a j payne 2004 past rates of accumulation in central west antarctica geophys res lett 31 l12403 doi 10 1029 2004gl020290

cosmic ray labeling of erosion surfaces in situ nuclide production rates and erosion models a number of in situ cosmogenic radionuclides and stable nuclides have been measured in natural exposed rock surfaces with a view to study their in situ production and rock erosion rates 1 the in situ radionuclides can be used for a high resolution tomography of the erosional history of an exposed surface two stable nuclides 3he 21ne and five radionuclides 10be 26al 36cl 14c 39ar having half lives in the range of 300 1 5 106 yr half life are measurable in many rock types a prerequisite for the application of the in situ nuclides for the study of erosional histories of surfaces is a knowledge of their production rates under different irradiation conditions altitude latitude irradiation geometry and shielding relative nuclide production rates can be determined fairly accurately using the extensive available data on cosmic ray neutrons 2 absolute nuclide production rates cannot generally be predicted with any accuracy because of lack of data on excitation functions of nuclides unless some normalization is possible as was done in the case of several cosmic ray produced isotopes in the atmosphere 3 based on a recent natural calibration experiment in which erosion free surfaces exposed to cosmic radiation for 11 000 yrs were sampled the absolute production rates of 10be and 26al in quartz have been accurately estimated for mountain altitudes in sierra nevada 4 the absolute production rates of 10be and 26al in quartz can therefore be estimated fairly accurately for any given latitude and altitude some measurements of 14c in rocks of low erosion rate 5 similarly allow an estimate of its production rate attempts made to measure the in situ production rates of 3he in rocks have not yet led to a convergent production rate in view of the importance of knowing the production rates of isotopes of he ne and ar i present here theoretical estimates of their production rates based on available cross section data i discuss the information that can be extracted from the study of the in situ nuclides in rocks useful parameters characterizing the exposure history of a rock surface are 1 the effective surface exposure age and 2 the time averaged erosion rate the implications of these parameters for single and multiple nuclide studies are discussed in terms of the erosion models considered

the rhone glacier was smaller than today for most of the holocene we present the development and application of the novel in situ cosmogenic 14 c 10 be chronometer to recently exposed proglacial bedrock of the rhone glacier switzerland results show that during the holocene the glacier was smaller than today for 6500 2000 yr and larger than today for 4500 2000 yr this pattern is consistent with limited data from other techniques for glaciers in the alps and scandinavia but in contrast to glaciers from the southern alps of new zealand emphasizing the high sensitivity of large mountain glaciers to small climate changes the 14 c 10 be chronometer also shows that abrasion rates beneath the rhone glacier increased with ice speed

late holocene ice flow reconfiguration in the weddell sea sector of west antarctica abstract here we report late holocene ice sheet and grounding line changes to the weddell sea sector of west antarctica internal radio echo layering within the bungenstock ice rise which comprises very slow flowing ice separating the fast flowing institute and moller ice streams reveals ice deformed by former enhanced flow overlain by un deformed ice the ice rise surface is traversed by surface lineations explicable as diffuse ice flow generated stripes which thus capture the direction of flow immediately prior to the creation of the ice rise the arrangement of internal layers can be explained by adjustment to the flow path of the institute ice stream during either a phase of ice sheet retreat not longer than 4000 years ago or by wholesale expansion of the grounding line from an already retreated situation not sooner than 400 years ago some combination of these events involving uplift of the ice rise bed during ice stream retreat and reorganisation is also possible whichever the case the implication is that the ice sheet upstream of the bungenstock ice rise which currently grounds over a 1 5 km deep basin has been and therefore may be susceptible to significant change

data driven numerical simulations of equatorial spread f in the peruvian sector 2 autumnal equinox an ongoing effort to simulate plasma instability in the equatorial ionosphere leading to equatorial spread f esf in the american sector is described ionospheric state parameters including plasma number density and vector drift velocity profiles were measured at the jicamarca radio observatory in the period between 20 september and 3 october 2013 coherent radar backscatter from plasma irregularities was recorded simultaneously and images of the irregularities were calculated using aperture synthesis methods neutral winds were measured by the red line fabry perot interferometers at jicamarca and arequipa peru a fully 3 d numerical simulation of ionospheric irregularities initialized and forced using parameterizations derived from measurements and empirical models was used to reproduce the esf activity observed simulations were able to recover many of the features of the irregularities although some important anomalies can be noted esf events in which the first appearance of radar plumes occurred either very early or very late were not reproduced in simulation and may be indicative of nonlocal influence

deformation in the vicinity of ice divides calculs numeriques aux elements finis de la variation de la vitesse horizontale avec la profondeur au voisinage d une difference de la glace symetrique isotherme et ou glissement sur un lit plat

steep reverse bed slope at the grounding line of the weddell sea sector in west antarctica the bed of the west antarctic ice sheet is in places more than 1 5 km below sea level1 2 it has been suggested that a positive ice loss feedback may occur when an ice sheet s grounding line retreats across a deepening bed1 2 3 applied to the west antarctic ice sheet this process could potentially raise global sea level4 by more than 3 m hitherto attention has focussed on changes at the siple coast5 6 7 and amundsen sea embayment8 9 10 sectors of west antarctica here we present radio echo sounding information from the ice sheet s third sector the weddell sea embayment that reveals a large subglacial basin immediately upstream of the grounding line the reverse bed slope is steep with about 400 m of decline over 40 km the basin floor is smooth and flat with little small scale topography that would delay retreat indicating that it has been covered with marine sediment5 11 and was previously deglaciated upstream of the basin well defined glacially carved fjords with bars at their mouths testify to the position of a former ice margin about 200 km inland from the present margin evidence so far suggests that the weddell sea sector of the west antarctic ice sheet has been stable but in the light of our data we propose that the region could be near a physical threshold of substantial change

major ice sheet change in the weddell sea sector of west antarctica over the last 5 000 years until recently little was known about the weddell sea sector of the west antarctic ice sheet in the last 10 years a variety of expeditions and numerical modelling experiments have improved knowledge of its glaciology glacial geology and tectonic setting two of the sector s largest ice streams rest on a steep reverse sloping bed yet despite being vulnerable to change satellite observations show contemporary stability there is clear evidence for major ice sheet reconfiguration in the last few thousand years however knowing precisely how the ice sheet has changed in the past and when would allow us to better understand whether it is now at risk two competing hypotheses have been established for this glacial history in one the ice sheet retreated and thinned progressively from its last glacial maximum position retreat stopped at or very near the present position in the late holocene alternatively in the late holocene the ice sheet retreated significantly upstream of the present grounding line and then advanced to the present location due to glacial isostatic adjustment and ice shelf and ice rise buttressing both hypotheses point to data and theory in their support yet neither has been unequivocally tested or falsified here we review geophysical evidence to determine how each hypothesis has been formed where there are inconsistencies in the respective glacial histories how they may be tested or reconciled and what new evidence is required to reach a common model for the late holocene ice sheet history of the weddell sea sector of west antarctica

ice age storm trajectories inferred from radar stratigraphy at taylor dome antarctica the taylor dome ice core east antarctica provides a climate record that extends through the previous interglacial we use ice dynamics modeling and geochemical records to infer the local accumulation rate history a pronounced accumulation reduction occurred during the last glacial maximum lgm the spatial gradient of accumulation across the surface of taylor dome indicates the modern local prevailing storm direction radar internal layering shows that this gradient was reversed during the lgm suggesting a reorganization of synoptic weather systems occurred possibly linked to the advance and retreat of the west antarctic ice sheet

antarctic ice rise formation evolution and stability antarctic ice rises originate from the contact between ice shelves and one of the numerous topographic highs emerging from the edge of the continental shelf while investigations of the raymond effect indicate their millennial scale stability little is known about their formation and their role in ice shelf stability here we present for the first time the simulation of an ice rise using the bisicles model the numerical results successfully reproduce several field observable features such as the substantial thinning downstream of the ice rise and the successive formation of a promontory and ice rise with stable radial ice flow center showing that ice rises are formed during the ice sheet deglaciation we quantify the ice rise buttressing effect found to be mostly transient delaying grounding line retreat significantly but resulting in comparable steady state positions we demonstrate that ice rises are key in controlling simulations of antarctic deglaciation

sea level as a stabilizing factor for marine ice sheet grounding lines climate change could potentially destabilize marine ice sheets which would affect projections of future sea level rise1 2 3 4 specifically an instability mechanism5 6 7 8 has been predicted for marine ice sheets such as the west antarctic ice sheet that rest on reversed bed slopes whereby ice sheet thinning or rising sea level leads to irreversible retreat of the grounding line however existing analyses of this instability mechanism have not accounted for deformational and gravitational effects that lead to a sea level fall at the margin of a rapidly shrinking ice sheet9 10 11 here we present a suite of predictions of gravitationally self consistent sea level change following grounding line migration our predictions vary the initial ice sheet size and also consider the contribution to sea level change from various subregions of the simulated ice sheet using these results we revisit a canonical analysis of marine ice sheet stability5 and demonstrate that gravity and deformation induced sea level changes local to the grounding line contribute a stabilizing influence on ice sheets grounded on reversed bed slopes we conclude that accurate treatments of sea level change should be incorporated into analyses of past and future marine ice sheet dynamics

computations of the viscoelastic response of a 3 d compressible earth to surface loading an application to glacial isostatic adjustment in antarctica and canada s u m m a r y we develop a 3 d finite element model to study the viscoelastic response of a compressible earth to surface loads the effects of centre of mass motion polar wander feedback and selfconsistent ocean loading are implemented to assess the model s accuracy we benchmark the numerical results against a semi analytic solution for spherically symmetric structure we force our model with the ice 5g global ice loading history to study the effects of laterally varying viscosity structure on several glacial isostatic adjustment gia observables including relative sea level rsl measurements incanada and present day time variable gravity and uplift rates in antarctica canadian rsl observations have been used to determine the earth s globally averaged viscosity profile antarcticgps uplift rates have been used to constrainantarcticgia models and gia time variable gravity and uplift signals are error sources for grace and altimeter estimates of present day antarctic ice mass loss and must be modelled and removed from those estimates computinggia results for a 3 d viscosity profile derived from a realistic seismic tomography model and comparing with results computed for 1 d averages of that 3 d profile we conclude that 1 a gia viscosity model based on canadian relative sea level data is more likely to represent a canadian average than a true global average 2 the effects of 3 d viscosity structure on grace estimates of present day antarctic mass loss are probably smaller than the difference between gia models based on different antarctic deglaciation histories and 3 the effects of 3 d viscosity structure on antarctic gps observations of present day uplift rate can be significant and can complicate efforts to use gps observations to constrain 1 d gia models

a model for subglacial flooding through a preexisting hydrological network during the rapid drainage of supraglacial lakes increasingly large numbers of supraglacial lakes form and drain every summer on the greenland ice sheet some of these lakes drain rapidly within the timescale of a few hours and the vertical discharge of water during these events may find a preexisting film of water potentially within a distributed drainage system of linked cavities here we present a model for subglacial flooding applied specifically to such circumstances one of many interesting results we find is that water flows in from the far field prior to the arrival of flooding we systematically evaluate the effect of initial ice bed opening on the degree of perturbation to the subglacial system of particular importance we find that floods propagate much faster and vertical displacements are much greater for larger openings for 10 and 1 cm of initial opening for example floods travel about 68 and 50 farther than in the fully coupled ice bed scenario after 2 h of drainage respectively for the same choices of initial opening the elastostatic displacement at the injection point is about 1 39 and 1 26 times that of the fully coupled scenario using the framework with a preexisting water film results in avoiding the pressure singularity that is inherent to classical hydrofracture models thus opening an avenue for integrating the likes of our model within continuum subglacial hydrological models furthermore we foresee that the theory presented can be used to potentially infer subglacial hydrological conditions from surface observables

contributions and unrealized potential contributions of cosmogenic nuclide exposure dating to glacier chronology 1990 2010 abstract this paper reviews the application of cosmogenic nuclide exposure dating to glacier chronology exposure dating of glacial landforms has made an outsize impact on this field because the technique filled an obvious need that had already been recognized by glacial geologists by now hundreds of studies have used cosmogenic nuclide exposure dating to date glacial deposits and in fact it is rare to find a study of glacial geology or glacier chronology or any paleoclimate synthesis that makes use of such studies that does not involve exposure dating these developments have resulted in major contributions to glacier chronology and paleoclimate in particular i reconstructing antarctic ice sheet change ii establishing the chronology of late pleistocene and holocene glacier change in mountain regions where it was previously unknown iii establishing the broad chronological outlines of mountain glaciations prior to the last glacial maximum and iv gaining insight into subglacial erosional processes through the observation that many glaciated surfaces preserve cosmogenic nuclide inventories from long past ice free periods as well as the present one an important potential future contribution will be the application of the large data set of exposure dated glacier chronologies to better understand global and regional climate dynamics during lateglacial and holocene millennial scale climate changes however this contribution cannot be realized without significant progress in two areas i understanding and accounting for geologic processes that cause apparent exposure ages on glacial landforms to differ from the true age of the landform and ii minimizing systematic uncertainties in exposure ages that stem from cosmogenic nuclide production rate estimates and scaling schemes at present there exists an enormous data set of exposure ages on glacial deposits but these data cannot be used to their full potential in paleoclimate syntheses due to an inadequate understanding of geologic scatter and production rate uncertainties the intent of this paper is to highlight this situation and suggest some strategies for realizing this potential

holocene deglaciation of marie byrd land west antarctica surface exposure ages of glacial deposits in the ford ranges of western marie byrd land indicate continuous thinning of the west antarctic ice sheet by more than 700 meters near the coast throughout the past 10 000 years deglaciation lagged the disappearance of ice sheets in the northern hemisphere by thousands of years and may still be under way these results provide further evidence that parts of the west antarctic ice sheet are on a long term trajectory of decline west antarctic melting contributed water to the oceans in the late holocene and may continue to do so in the future

a deglacial model for antarctica geological constraints and glaciological modelling as a basis for a new model of antarctic glacial isostatic adjustment abstract we present a new reconstruction of the antarctic ice sheets between 20 ka bp and the present day our reconstruction is derived using a numerical model to generate a physically consistent ice surface across the whole of the continent we define the extent of the ice sheet at five time slices 20 15 10 5 and 0 ka bp assuming an equilibrium state for the 20 ka bp time slice and a transient state for the deglacial time slices the evolution of the ice sheet within the numerical model is driven by variations in temperature accumulation rate and relative sea level in order to reconstruct the concave profile of the ice sheet in marine grounded regions such as the ross and weddell seas we force our model to develop channels of faster flow by defining greater basal sliding along the trajectory of former ice streams we find a strong dependence upon the basal sliding parameters and also the position of the grounding line we use an extensive data base of geological and glaciological data to constrain our ice sheet reconstruction grounding line extent is prescribed from marine geological data and we test ice sheet thickness against onshore geological data at 62 sites of the five time slices considered our 20 ka bp reconstruction is the best constrained by data and has an rms misfit of 147 6 m when compared to observations of ice thickness change between 20 ka bp and the present day across all time slices there are large regions of the ice sheet which are poorly constrained especially after 20 ka bp we estimate the spatial distribution of uncertainty in our ice sheet reconstruction and note that the solutions are least reliable in regions of complex topography we predict that the antarctic ice sheets contributed 9 1 5 m of eustatic sea level to the global ocean between 20 ka bp and the present and our reconstruction with minimum misfit contributes 8 m eustatic sea level during this period these values which we argue are an upper bound are lower than many previous estimates the reconstructed pattern of ice unloading can serve as a new input for glacial isostatic models

reconstruction of changes in the weddell sea sector of the antarctic ice sheet since the last glacial maximum the weddell sea sector is one of the main formation sites for antarctic bottom water and an outlet for about one fifth of antarctica s continental ice volume over the last few decades studies on glacial geological records in this sector have provided conflicting reconstructions of changes in ice sheet extent and ice sheet thickness since the last glacial maximum lgm at ca 23 19 calibrated kiloyears before present cal ka bp terrestrial geomorphological records and exposure ages obtained from rocks in the hinterland of the weddell sea ice sheet thickness constraints from ice cores and some radiocarbon dates on offshore sediments were interpreted to indicate no significant ice thickening and locally restricted grounding line advance at the lgm other marine geological and geophysical studies concluded that subglacial bedforms mapped on the weddell sea continental shelf subglacial deposits and sediments over compacted by overriding ice recovered in cores and the few available radiocarbon ages from marine sediments are consistent with major ice sheet advance at the lgm reflecting the geological interpretations different ice sheet models have reconstructed conflicting lgm ice sheet configurations for the weddell sea sector consequently the estimated contributions of ice sheet build up in the weddell sea sector to the lgm sea level low stand of 130 m vary considerably in this paper we summarise and review the geological records of past ice sheet margins and past ice sheet elevations in the weddell sea sector we compile marine and terrestrial chronological data constraining former ice sheet size thereby highlighting different levels of certainty and present two alternative scenarios of the lgm ice sheet configuration including time slice reconstructions for post lgm grounding line retreat moreover we discuss consistencies and possible reasons for inconsistencies between the various reconstructions and propose objectives for future research the aim of our study is to provide two alternative interpretations of glacial geological datasets on antarctic ice sheet history for the weddell sea sector which can be utilised to test and improve numerical ice sheet models

sequential stagnation of kamb ice stream west antarctica 1 ice penetrating radar data confirm the presence of two relict northern margins of kamb ice stream kis west antarctica in an area called the duckfoot between the two relict margins deep layers are preserved and buried crevasses remain at a constant depth until roughly 1 5 km from the innermost youngest margin suggesting that this outer portion of the ice stream stagnated suddenly 10 years stagnation was likely accomplished through changes in basal conditions i e reduced lubrication beneath this outer portion of the ice stream and resulted in a narrowing of the ice stream trunk by 27 and a reduction in flow speed 200 years prior to complete stagnation flowstripe deformation as a result of duckfoot stagnation provides an estimate of kis paleo velocity of 210 m a after duckfoot stagnation

airborne radar evidence for tributary flow switching in institute ice stream west antarctica implications for ice sheet configuration and dynamics despite the importance of ice streaming to the evaluation of west antarctic ice sheet wais stability we know little about mid to long term dynamic changes within the institute ice stream iis catchment here we use airborne radio echo sounding to investigate the subglacial topography internal stratigraphy and holocene flow regime of the upper iis catchment near the ellsworth mountains internal layer buckling within three discrete topographically confined tributaries through ellsworth independence and horseshoe valley troughs provides evidence for former enhanced ice sheet flow we suggest that enhanced ice flow through independence and ellsworth troughs during the mid to late holocene was the source of ice streaming over the region now occupied by the slow flowing bungenstock ice rise although buckled layers also exist within the slow flowing ice of horseshoe valley trough a thicker sequence of surface conformable layers in the upper ice column suggests slowdown more than 4000 years ago so we do not attribute enhanced flow switch off here to the late holocene ice flow reorganization intensely buckled englacial layers within horseshoe valley and independence troughs cannot be accounted for under present day flow speeds the dynamic nature of ice flow in iis and its tributaries suggests that recent ice stream switching and mass changes in the siple coast and amundsen sea sectors are not unique to these sectors and that they may have been regular during the holocene and may characterize the decline of the wais

marine ice sheet collapse potentially under way for the thwaites glacier basin west antarctica resting atop a deep marine basin the west antarctic ice sheet has long been considered prone to instability using a numerical model we investigated the sensitivity of thwaites glacier to ocean melt and whether its unstable retreat is already under way our model reproduces observed losses when forced with ocean melt comparable to estimates simulated losses are moderate 1 mm per year of sea level rise collapse in the different simulations within the range of 200 to 900 years

a 3 d coupled ice sheet sea level model applied to antarctica through the last 40 ky abstract we present results from a three dimensional ice sheet shelf model of antarctica coupled to a gravitationally self consistent global sea level model that incorporates maxwell viscoelastic deformation of the solid earth the coupled model captures complex post glacial changes in sea level associated with the gravitational deformational and rotational effects of the evolving surface mass ice plus ocean load over the global ocean including at the grounding lines of marine based ice the simulations are initiated at 40 ka and we focus on ice distributions and sea levels from the last glacial maximum to present our results extend and confirm the key conclusions of earlier work using a simplified 1 d ice sheet model by demonstrating that the sea level coupling has a significant stabilizing influence on marine ice sheet grounding lines the feedback of sea level changes into the ice sheet model acts to slow down the retreat and advance of the grounding line relative to simulations in which the full coupling is not incorporated differences in ice thickness between these simulations can reach 1 km close to the grounding line finally we perform preliminary comparisons of our results to relative sea level rsl histories and gps derived present day uplift rates at sites near the margins of antarctica we find that the coupling improves fits to uplift rates in several regions and that the rsl predictions of the coupled model yield a fit to the observations that is comparable to recent uncoupled simulations in which the underlying earth model was varied to obtain a best fit to the rsl histories

traveling slippery patches produce thickness scale folds in ice sheets large complex stratigraphic folds that rise as high as 60 of the local ice thickness have been observed in ice sheets on antarctica and greenland here we show that ice deformation caused by heterogeneous and time variable basal sliding can produce the observed structures we do this using a thermomechanical ice sheet model in which sliding occurs when the base approaches the melting point and slippery patches develop these slippery patches emerge and travel downstream because of a feedback between ice deformation vertical flow and temperature our model produces the largest overturned structures comparable to observations when the patches move at about the ice column velocity we conclude that the history of basal slip conditions is recorded in the ice sheet strata these basal conditions appear to be dynamic and heterogeneous even in the slow flowing interior regions of large ice sheets

exposure ages from mountain dipsticks in mac robertson land east antarctica indicate little change in ice sheet thickness since the last glacial maximum past changes in east antarctic ice sheet eais volume are poorly known and diffi cult to measure yet are critical for predicting the response of the ice sheet to modern climate change in particular it is important to identify the sources of sea level rise since the last glacial maximum lgm and ascertain the present day stability of the world s largest ice sheet we present altitudinal transects of 10 be and 26 al exposure ages across the framnes mountains in mac robertson land that allow the magnitude and timing of eais retreat to be quantifi ed our data show that the coastal eais thinned by at most 350 m in this region during the past 13 k y this reduction in ice sheet volume occurred over a 5 k y period and the present icesheet profi le was attained ca 7 ka in contrast to the west antarctic ice sheet which continues to retreat today combined with regional offshore and terrestrial geologic evidence our data suggest that the reduction in eais volume since the lgm was smaller than that indicated by contemporary ice sheet models and added little meltwater to the global oceans stability of the ice margin since the middle holocene provides support for the hypothesis that eais volume changes are controlled by growth and decay of northern hemisphere ice sheets and associated global sea level changes