validation of the ice 3g model of wurm wisconsin deglaciation using a global data base of relative sea level histories a global data base consisting of radiocarbon controlled relative sea level rsl histories from 392 different geographic locations has been compiled for use in analyses of the glacial isostatic adjustment process these data may be employed to constrain both the spatiotemporal characteristics of the last deglaciation event of the current ice age and the rheological structure of the planetary interior in a previous paper in this series tushingham and peltier 1991 we employed the rsl data in this set from sites that were once ice covered to refine the apriori model of deglaciation that had been delivered by past research in the area of pleistocene geomorphology peltier and andrews 1976 wu and peltier 1983 to produce a new global model of this event that we have called ice 3g in the present paper following a description of the full data base and of the methods used to compile it we employ data from the sites that were not ice covered 200 time series to verify the plausibility of this refined model of deglaciation and to test the extent to which the totality of the observations constrain the radial viscoelastic structure of the planet the analysis of solutions to the forward problem presented here reinforce previous inferences of a preference for a thick lithosphere by rsl data from sites along the u s east coast and a preference by data from most near field sites for an increase in viscosity across the 670 km discontinuity in the mantle by a factor of approximately 2 when the viscosity of the upper mantle is fixed to the value of 1021 pa s originally deduced by haskel 1935 on the basis of his analysis of the data from fennoscandia

localization of the gravity field and the signature of glacial rebound the negative free air gravity anomaly centred on hudson bay canada shows a remarkable correlation with the location of the laurentide ice sheet suggesting that this gravity anomaly is the result of incomplete post glacial rebound this region however is also underlain by higher than average mantle seismic velocities suggesting that the gravity low might result instead from dynamic topography associated with convective downwellings here we analyse the global gravity field as a simultaneous function of geographic location and spectral content we find that the hudson bay gravity low is unique with anomalously high amplitude in the spectral band where the power from the laurentide ice load is greatest and the relaxation times predicted for viable models of viscous relaxation are longest we estimate that about half of the hudson bay gravity anomaly is the result of incomplete post glacial rebound and derive a mantle viscosity model that explains both this gravity signature and the characteristic uplift rates for the central laurentide and fennoscandian regions this model has a jump in viscosity at 670 km depth comparable to that in dynamic models of the geoid highs over subducted slabs but lacks a low viscosity asthenosphere consistent with a higher viscosity in the upper mantle beneath shields than in oceanic regions

post glacial sea level rise from a coral record at huon peninsula papua new guinea dating of coral reef terraces can provide a record of changes in sea level which should be pronounced during the transition between glacial and interglacial periods cores drilled from coral reefs at barbados in the equatorial west atlantic1have revealed the sea level changes that occurred during the younger dry as event at the end of the last glacial maximum 11 000 yr bp it has not been known however whether pacific coral reefs can grow at a rate sufficient to keep up with the rise in sea level during such a transition here we report results obtained from a 52 m drill core from the post glacial reef at huon peninsula papua new guinea spanning the interval from 7 000 to 11 000 14c yr bp which show that coral growth kept pace while the relative sea level rose by 50m although the tectonic environment is very different from that at barbados the two records compare well when corrections are made for local tectonic uplift showing that sea level rise was similar at both locations the rate of rise was greatest between 9 000 and 10 00014c yr bp corresponding to the time of the younger dryas

haskell 1935 revisited the haskell 1935 value of 1021 pa s for mantle viscosity is a classic and enduring constraint on the rheology of the earth s interior we revisit this inference using spherically symmetric self gravitating maxwell viscoelastic earth models our inference is based on both forward and inverse analyses of decay times associated with uplift at two sites considered by haskell angerman river sweden and oslo norway rather than the raw relative sea level rsl data at these sites we demonstrate that predictions of the decay time associated with the angerman river data are insensitive to variations in both the late pleistocene ice load history and the lithospheric thickness of the earth model the predictions at oslo are sensitive to both these inputs and hence decay times at this site provide a remarkably robust constraint on mantle viscosity we derive a constraint on the average viscosity of the mantle of 0 65 1 10 1021 pa s where the average resolved by the data encompasses a region which extends from the base of the lithosphere to a depth of near 1400 km this indicates that many previous analyses which have invoked the haskell value of 1021 pa s as a constraint on the average upper mantle i e above 670 km depth viscosity alone have misinterpreted the resolving power of the inference furthermore our analysis indicates that a number of apparently contradictory inferences of viscosity based on fennoscandian data satisfy the new rigorous interpretation of the haskell constraint finally we demonstrate how the ambiguity in the upper mantle lower mantle viscosity contrast associated with the haskell average may be reduced by invoking decay time constraints estimated from rsl curves in hudson bay a preliminary inversion of decay times at angerman river and richmond gulf in hudson bay suggests a contrast of approximately an order of magnitude between the average viscosities of the upper mantle and the top 1000 km of the lower mantle however a conclusive analysis in this regard must await the determination of consistent decay time estimates for the hudson bay region

ice age paleotopography a gravitationally self consistent theory of postglacial relative sea level change is used to infer the variation of surface ice and water cover since the last glacial maximum lgm the results show that lgm ice volume was approximately 35 percent lower than suggested by the climap reconstruction and the maximum heights of the main laurentian and fennoscandian ice complexes are inferred to have been commensurately lower with respect to sea level use of these ice age boundary conditions in atmospheric general circulation models will yield climates that differ significantly from those previously inferred on the basis of the climap data set

deglacial sea level record from tahiti corals and the timing of global meltwater discharge the timing of the last deglaciation is important to our understanding of the dynamics of large ice sheets1 and their effects on the earth s surface2 3 moreover the disappearance of the glacial ice sheets was responsible for dramatic increases in freshwater fluxes to the oceans which probably disturbed the ocean s thermohaline circulation and hence global climate4 7 sea level increases bear witness to the melting of continental ice sheets but only two such records from barbados8 9 and new guinea10 11 corals have been accurately dated but these corals overlie active subduction zones where tectonic movements are large and often discontinuous especially in new guinea so the apparent sea level records may be contaminated by a complex tectonic component here we date fossil corals from tahiti which is far from plate boundaries and thus is likely to be tectonically relatively stable and remote from the locations of large former ice sheets the resulting record indicates a large sea level jump shortly before 13 800 calendar years bp which corresponds to meltwater pulse 1a in the barbados coral records8 9 the timing of this event is more accurately constrained in the tahiti record revealing that the meltwater pulse coincides with a short and intense climate cooling event12 15 that followed the initiation of the bolling allerod warm period12 16 but preceded the younger dryas cold event by about 1 000 years

global glacial isostatic adjustment palaeogeodetic and space geodetic tests of the ice 4g vm2 model analyses of the global process of glacial isostatic adjustment and post glacial relative sea level change continue to deliver important insights into earth system form and process one successful model of the related phenomenology is based upon a spherically symmetric internal viscoelastic structure for the solid earth which has been denoted vm2 and a model of the most recent deglaciation event of the current ice age denoted ice 4g the primary purpose of this paper is to describe several new a posteriori tests that have recently been performed to further investigate the quality of this global solution to the inverse problem for both mantle viscosity and deglaciation history that is posed by the observables associated with this large scale geodynamic phenomenon i focus especially upon the misfits of observations to the theoretical predictions of this model which i am currently using to further refine its properties and upon predictions made using it of geophysical signals that should soon become visible in the context of the gravity recovery and climate experiment grace satellite mission among the required refinements to ice 4g one that is necessary to eliminate a recently revealed misfit to space geodetic constraints on the present day rate of radial motion at the yellowknife location well to the west of hudson bay and a similar misfit to absolute gravity measurements to the southwest of the bay is the insertion of a keewatin dome of thick ice centred over yellowknife with a ridge of ice extending to the south east in the geomorphological literature the existence of such a keewatin dome previously has been hypothesised but chronological control was lacking on the surface features that suggested its former existence an important additional constraint that requires the late glacial existence of this important feature consists of new inferences of the last glacial maximum lowstand of the sea from sites in the far field of the main concentrations of land ice copyright 2002 john wiley sons ltd

holocene tectonic uplift patterns in northeastern sicily evidence from marine notches in coastal outcrops abstract marine notches are groove like features formed in rocky coastal outcrops in relation to past sea levels and are widely used as indicators of rates and patterns of tectonic activity examination of uplifted notches in the complex plate boundary setting of northeast sicily reveals the need for considerable care in their interpretation present notch forming processes are dominated by carbonate dissolution aided by microbial activity as waves surge and recede at the base of the cliffs and produce notches approximately 0 5 m above present biological mean sea level i e above an organic constructional rim of coralline algae used to define sea level at sea level dissolution is restricted through saturation and the bedrock protected by the outbuilding algal rim this assemblage varies in response to fetch exposure and uplift rate typically producing a complex band of uplifted notch features exposed above sea level correlations based on height of the shallow and discontinuous notch indentations within this band are unreliable however the upper margin of the band marked by the roof of the uppermost notch is a prominent and laterally persistent feature easily recognised as the boundary between subaerial and marine processes we propose that this is an approximately coeval datum formed when post glacial sea level rise decelerated in the mid holocene to match the rate of tectonic uplift a notion strengthened by the available dates of about 5 ka for this feature subsequent uplift outpaced later sea level rise taking the notches above sea level this proposed datum varies in height from about 2 m at capo milazzo on the north coast to about 5 m on the east coast bordered by the messina fault system the height differential is consistent with the respective elevations of nearby marine terraces dating from the last interglacial at about 125 ka but also exhibits shorter term variations suggesting holocene increases in activity rates on more localised structures affecting the localities involved regionally northeast sicily is characterised by long term differential uplift rather than behaving as a coherent block the highest uplift rate probably exceeding 2 mm a 1 calculated using the sea level curve from fairbanks r g 1989 a 17000 year glacio eustatic sea level record influence of glacial melting rates on the younger dryas event and deep ocean circulation nature 342 637 642 is suggested for a notch locality on the upthrown side of the malta escarpment rates of about 2 mm a 1 are indicated by the 5 m notch at localities bordered by the messina fault system and a rate close to 1 mm a 1 is suggested for the undated 2 m notch at capo milazzo

paleotopography of glacial age ice sheets this is technical comment and response to the subject of paleotophography of glacial age ice sheets the model presented by peltier reconstructing the paleotopography of glacial age ice sheets has implications for atmospheric general circulation models of ice age climate in addition the model suggests that the glacial age antarctic ice sheet was significantly larger than today s the commentor edwards suggests there is a discrepancy between data from papua new guinea and the model results

glacial isostatic adjustment i the forward problem summary the isostatic adjustment of a radially stratified visco elastic spheroid is treated using space time green functions for the associated surface mass load boundary value problem these impulse response functions are convolved with a heaviside function to give the time dependent deformation of the planet which would be produced by a unit point mass brought up from infinity at t 0 and allowed to remain on the surface the resulting heaviside green functions can be employed to simulate all of the important signatures of glacial isostatic adjustment given a space time dependent surface mass load consisting of ice sheet ablation histories and a model of the simultaneous filling of the ocean basins these source terms are simply convolved with the heaviside green function appropriate to a specific response signature a realistic model of the spatial distribution of the main late pleistocene ice loads and of their temporal disintegration is constructed the response of two visco elastic earth models is computed and compared to a global set of relaxation data relative sea level curves on the basis of this initial comparison of theory and observation the possibility that the lower mantle has a viscosity which is significantly in excess of the viscosity of the upper mantle is excluded in addition clear evidence of the presence of the lithosphere has been found in relaxation data from sites which were near the edge of the laurentide ice sheet such data should therefore prove useful as a basis for analysis of lateral variations in lithosphere thickness further extensions of the calculation are suggested

late quaternary variations in relative sea level due to glacial cycle polar wander growth and decay of continental ice sheets can excite significant motion of the earth s rotation pole and cause a complex spatio temporal pattern of changes in relative sea level these two effects have generally been considered separately but may interact in important ways in particular a simple model of the melting of the laurentide ice sheet causes a uniform eustatic sea level rise of 55 m and also induces a motion of the rotation pole by 0 1 to 1 degree depending on viscosity structure in the mantle this motion produces a secular pole tide which is a spherical harmonic degree 2 order 1 component of the relative sea level pattern with peak to peak amplitude of 20 to 40 m the maximum effect is along the great circle passing through the path of the pole and at latitudes of 45 this secular pole tide has been ignored in most previous attempts to estimate ice sheet loading history and mantle viscosity from global patterns of relative sea level change it has a large influence along the east coast of north america and the west coast of south america and significantly contributes to present day rates of relative sea level change

radial profile of mantle viscosity results from the joint inversion of convection and postglacial rebound observables we present new inferences of the radial profile of mantle viscosity that simultaneously fit long wavelength free air gravity harmonics associated with mantle convection and a large set of decay times estimated from the postglacial uplift of sites within previously glaciated regions hudson bay arctic canada and fennoscandia the relative sea level variation at these latter sites is constrained by age height pairs obtained by geological survey rather than the subjective trends which are commonly used in glacial isostatic adjustment gia studies our viscosity inferences are generated using two approaches first we adopt a relative viscosity profile which is known to provide a good fit to the free air gravity harmonics and determine an absolute scaling which yields a best fit to the gia decay time constraints second we perform an iterative nonlinear joint inversion of the two data sets in both cases our inferred profiles are characterized by a significant increase of viscosity 2 orders of magnitude with depth to values of 1022 pa s in the bottom half of the lower mantle the new viscosity profiles are shown to satisfy constraints based on the postglacial uplift of both fennoscandia the classic haskell 1935 number and hudson bay which have commonly been invoked to argue for an isoviscous mantle furthermore the models are used to predict a set of long wavelength signatures of the gia process these include predictions of gia induced variations in 1 the length of day over the late holocene period 2 the earth s precession constant and obliquity over the last 2 6 myr and 3 the present day zonal harmonics of the geopotential jl l 7 the predictions 1 and 3 bound the late holocene and ongoing mass flux between the large polar ice sheets greenland and antarctic and the global oceans to small values 0 4 mm yr equivalent eustatic sea level rise

rheology of the upper mantle a synthesis rheological properties of the upper mantle of the earth play an important role in the dynamics of the lithosphere and asthenosphere however such fundamental issues as the dominant mechanisms of flow have not been well resolved a synthesis of laboratory studies and geophysical and geological observations shows that transitions between diffusion and dislocation creep likely occur in the earth s upper mantle the hot and shallow upper mantle flows by dislocation creep whereas cold and shallow or deep upper mantle may flow by diffusion creep when the stress increases grain size is reduced and the upper mantle near the transition between these two regimes is weakened consequently deformation is localized and the upper mantle is decoupled mechanically near these depths

correction to glacial isostatic adjustment and earth rotation refined constraints on the viscosity of the deepest mantle by w r peltier and xianhua jiang

on postglacial geoid subsidence over the equatorial oceans we develop two new spectral formalisms for the gravitationally self consistent solution of the sea level equation which governs the redistribution of glacial meltwater on a visco elastic earth the first is a purely spectral technique based on an extension of theory outlined by dahlen 1976 for determining the equilibrium oceanic tide on an elastic planet using this technique with a feasible spherical harmonic truncation level l 30 we can obtain convergence in the computed sea level variations within the degree range l 10 we also develop a second pseudospectral technique however which permits the construction of gravitationally self consistent solutions to much higher degree and order l 100 the pseudospectral formalism is employed here to develop a comprehensive physical explanation for the global pattern of present day sea level variations due to ongoing glacial isostatic adjustment in particular we focus upon a mechanism which we term equatorial ocean syphoning that acts to draw water toward the oceanic portion of the collapsing peripheral bulge that encircles previously glaciated regions the collapse of the forebulge induces a flow of water which is required to maintain hydrostatic equilibrium the syphoning mechanism dominates the relative sea level rsl variation in oceans in the far field of the ice sheets that is beyond the peripheral bulge during periods such as the interglacial of the past 4000 years during which the volume of the ice sheets has not apparently changed appreciably while at the same time isostatic adjustment persists the rsl change in the near field is in contrast dominated by the local vertical displacement of the solid surface finally comparison of gravitationally self consistent predictions of rsl change in the far field with the observational data at small pacific island sites has provided an upper bound of 1 to 2 m for the eustatic sea level rise produced by any recent last 3000 years melting of the antarctic ice sheet

the lageos constraint on deep mantle viscosity results from a new normal mode method for the inversion of viscoelastic relaxation spectra the recently observed acceleration in the node of the orbit of the lageos satellite has previously been shown to provide a good constraint on the viscosity of the earth s mantle beneath the 670 km seismic discontinuity this paper explores the extent to which the previously inferred lower mantle viscosity may be traded off against variations of other properties of the radial viscoelastic structure we show explicitly that the nontidal acceleration of rotation inferred from the lageos data and predicted as a consequence of the earth s response to the last deglaciation event of the current ice age is relatively insensitive both to variations of lithospheric thickness and to the presence or absence of plausible chemical heterogeneity across the 670 km interface it therefore provides a high quality constraint on the viscosity of the deep mantle with the upper mantle viscosity fixed at the nominal value of 1021 pa s the bounds on lower mantle viscosity lm are such that 2 1021 pa s lm 6 1021 pa s even when allowance is made for the presence of a viscous boundary layer at the base of the upper mantle such as would exist if the mantle convective circulation were layered

the viscoelastic relaxation of a realistically stratified earth and a further analysis of postglacial rebound summary we discuss the isostatic adjustment of the earth in response to pleistocene deglaciation and derive constraints on the earth s viscosity profile we model the earth as viscoelastic self gravitating realistically stratified and spherically symmetric it is shown that the earth s relaxation eigenspectrum besides possessing a small set of well known discrete modes also includes infinitely dense sets of modes a continuous spectrum which separate into classes that depend on the principal restoring force of the mode methods for estimating the continuous spectrum are discussed and the effects on the predicted uplift are assessed we find that the effects of the continuous spectrum are far less important than those of the discrete modes though the presence of the continuous spectrum can make the discrete modes more difficult to identify the continuous spectrum does become important when computing the rebound due to more recent loading such as might be associated with present day changes in polar ice we use our model to refine predictions of the earth s response to pleistocene glacial loading with particular emphasis on implications for mantle viscosity we have used time histories that are consistent with the recently recalibrated 14c time scale a special effort is made to match the geologic records of sea level at sites in the vicinity of hudson bay we find consistent with the results of other authors that the curvatures of the records are sensitive indicators of lower mantle viscosity but that there are inconsistencies between the records at different sites however we find we are also able to change the predicted curvatures by modifying the temporal and spatial distribution of the ice loads if we assume a factor of 50 increase in viscosity between the upper and lower mantles and if we assume that the ice on the eastern side of hudson bay melted a few thousand years later than the ice on the western side a result consistent with recent geologic evidence we are able to match the observed sea level data about as well as we are able to if we use instead a uniform mantle viscosity profile a model with a factor of 50 jump in viscosity predicts a large free air gravity anomaly over northern canada that is reasonably consistent with what is observed whereas a uniform viscosity model does not though some authors have argued that the observed north american gravity anomaly may be mostly a consequence of mantle convection and largely unrelated to postglacial rebound our objective here is not to present a preferred ice model or viscosity profile rather we conclude only that the interpretation of the data is ambiguous and that the mantle viscosity profile can not yet be uniquely inferred from them

glacial isostatic adjustment and earth rotation refined constraints on the viscosity of the deepest mantle we explore the extent to which two properties of earth s present day rotational state may be applied to constrain the viscosity of the lower regions of the lower mantle the analysis builds upon recent advances in understanding the radial resolving power of the relative sea level data of postglacial rebound which have demonstrated anew that no significant increase of viscosity across the 660 km seismic discontinuity appears to be allowed by these data such observations do not constrain the viscosity of the mantle below about 1400 km depth however and here we estimate the average viscosity in this deepest region by invoking observations of the present day rate and direction of polar wander and the present day magnitude of the nontidal acceleration of the rate of planetary rotation since the rotational response to the glaciation cycle depends only upon the degree 2 spherical harmonic components of the induced deformation this response is expected to provide the most specific sampling possible of deep mantle properties our analysis is based upon the use of glacial excitation functions calculated from complete gravitationally self consistent solutions of the sea level equation for the most recent refinement of the history of ice sheet thickness variations across the last glacial interglacial transition namely the ice 4g model this analysis shows that when earth rotation observations are combined with relative sea level constraints the viscosity of the deepest mantle is required to be 1 order of magnitude higher than the viscosity of the upper part of the lower mantle such significant elevation of viscosity in the lowermost region of the lower mantle is also in accord with recent inferences based upon analysis of the nonhydrostatic geoid the analyses presented herein therefore suffice to effect a reconciliation between the requirements of these distinct geophysical observations

holocene sea level determination relative to the australian continent u th tims and 14c ams dating of coral cores from the abrolhos islands u th tims and 14 c ams measurements are presented from two coral cores from the easter group of the houtman abrolhos islands between 28 s and 29 s on the western continental margin of australia the u th measurements on the morley core from morley island cover a depth interval from 0 2 m above present sea level to 24 4 m below present sea level and comprise eleven samples the ages vary between 6320 50 a at 0 2 m above sea level and 9809 95 a at 24 4 m below sea level all errors are 2 the mean growth rate is 7 1 0 9 m ka the 14 c dates of selected morley core corals show that the 14 c ages are 1000 a younger than their corresponding u th ages which agrees with previous results the main purpose of our 14 c measurements is to be able to compare them precisely with other coral cores where no u th measurements are available the u th measurements of the suomi core from suomi island cover a depth interval from 0 05 m to 14 2 m below present sea level and consist of four samples the ages vary between 4671 40 a at 0 05 m below sea level and 7102 82 a at 14 2 m below sea level with a mean growth rate of 5 8 0 2 m ka the growth history of both cores is explained by a simple model in which the growth rates of the morley core can be interpreted as reflecting local rates of sea level rise whereas the suomi core is interpreted as reflecting lateral growth during the past 6000 a our results indicate that sea level relative to the western margins of the australian continent was about 24 m lower than present at about 9800 a b p 14 c gives a date of 8500 a b p sea level then rose and reached a highstand slightly higher than the present position at about 6300 a b p 14c date 5500 a this highstand declined but was still higher than present at 4600 a b p this is in agreement with previous observations along the australian coastal margins and with observations from the huon peninsula papua new guinea our results are very similar to theoretical numerical models which take into consideration water loading and isostatic compensation and viscous mantle flow in contrast coral cores from barbados show that corals with a 14 c age of 5500 a b p are some 10 m b p s l we interpret the difference between the barbados core and the morley core as resulting from additional flooding of barbados by water redistribution due to changes in the earth s geoid but not reflecting global sea level rise or major addition of melt waters over the past 6000 a the difference in the geoid at barbados between 6000 a b p and the present will require a refinement in the geophysical models precise 230 th tims measurements on continental coasts will be required to provide an adequate data base for modelling deformation flow of mantle material and sea level height

earth s gravitational field seismic tomography resolves the enigma of the laurentian anomaly we demonstrate that the negative free air gravity anomaly over hudson bay is not due to glacial isostatic disequilibrium but rather to the influence of large scale mantle convection

glacial isostatic adjustment and the free air gravity anomaly as a constraint on deep mantle viscosity summary a gravitationally self consistent theory for relative sea level variations forced by pleistocene deglaciation events is employed to explore the extent to which rsl and free air gravity observations together constrain the viscosity of the mantle beneath the seismic discontinuity at 670 km depth a trade off is revealed between errors in the assumed deglaciation history and errors in the inferred value of the viscosity of the lower mantle taking full account of such uncertainty plausible values of the viscosity beneath the transition region are bounded above by 1023 poise cgs units the preferred value is at least a factor of 2 lower than this and is strongly constrained by the observed free air gravity anomaly over hudson bay the calculations described in detail here show for the first time that the relatively large gravity anomalies observed over sites of wurm wisconsin deglaciation do not require any pronounced increase of mantle viscosity with depth in order to explain them

terminating the 100 kyr ice age cycle we report a simulation of the most recent 100 000 year glaciationdeglaciation cycle of the late pleistocene ice age a simulation that delivers an ice sheet chronology that is in close accord with that inferred from the geological record our analyses are performed with a reduced model of the climate system that incorporates significant improvements to the representation of both climate forcing and mass balance response in a previously described theory based upon a coupled one level energy balance model ebm and vertically integrated ice sheet model ism the theory fully incorporates the influences of orbital insolation forcing glacial isostatic adjustment and variations in the atmospheric concentrations of greenhouse gases it correctly predicts the main geographical regions of the northern hemisphere that were glaciated at last glacial maximum 21 000 years ago as well as the abrupt termination of the glacial epoch that occurred subsequently the latter feature of the ice age cycle is obtained without the need to incorporate unconstrained and therefore controversial physical processes into the model a limitation of all previous attempts to understand this global scale climate cycle our analyses suggest that the radiative impact on surface glaciation due to the changing atmospheric concentration of co2 is critical to the ability of the model to deliver a synthetic history of glaciation and deglaciation that is in accord with inferences based upon surface geological and geomorphological evidence with the incorporation of this influence model predicted ice thickness distributions at last glacial maximum lgm are very similar to those of the recently described ice 4g reconstruction that was based upon the inversion of postglacial relative sea level histories

calibration of the c 14 timescale over the past 30 000 years using mass spectrometric u th ages from barbados corals uranium thorium ages obtained by mass spectrometry from corals raised off the island of barbados confirm the high precision of this technique over at least the past 30 000 years comparison of the u th ages with 14c ages obtained on the holocene samples shows that the u th ages are accurate because they accord with the dendrochronological calibration before 9 000 yr bp the 14c ages are systematically younger than the u th ages with a maximum difference of 3 500 yr at 20 000 yrbp the u th technique thus provides a way of calibrating the radiocarbon timescale beyond the range of dendrochronological calibration

mantle viscosity and ice age ice sheet topography ice age paleotopography and mantle viscosity can both be inferred from observations of earth s response to the most recent deglaciation event of the current ice age this procedure requires iterative application of a theoretical model of the global process of glacial isostatic adjustment results demonstrate that the iterative inversion procedure converges to a paleotopography that is extremely close to that from the ice 4g model the accompanying mantle viscosity profile is furthermore shown to reconcile the requirements of aspherical geoid anomalies related to the mantle convection process thus resolving a fundamental issue concerning mantle rheology the combined model also explains postglacial sea level histories for the east coast of the united states

dynamic surface topography a new interpretation based upon mantle flow models derived from seismic tomography the very long wavelength i e in the degree range l 1 8 dynamic topography at the earth s surface estimated by isostatic reduction of the observed topography is shown to be in good agreement with the topography predicted by a model of viscous flow in the mantle this flow model was previously derived by independent consideration of the observed nonhydrostatic geoid and is based upon a recent model of large scale seismic shear velocity heterogeneity in the mantle our estimate and prediction of the long wavelength dynamic surface topography reveals significant depressions 2 3 km of the continental shields relative to their hydrostatic reference positions

holocene coral reef terraces and coseismic uplift of huon peninsula papua new guinea abstract as many as six levels of emerged holocene coral terraces occur along 40 km of coastline on the huon peninsula papua new guinea recording uplift history since culmination of the postglacial transgression the holocene reef crest ca 6000 yr b p is tilted down to the northwest parallel to the coast and concordant with the deformation of the last interglacial coral reef terrace and descends from 23 to 12 m in the study area the pattern and rate of deformation have been uniform in the late quaternary because average uplift rates have remained the same since the last interglaciation the holocene terraces described here are erosional features with regressive encrusting corals developed upon the holocene transgressive reef the multiple levels represent episodic probably coseismic uplift which has occurred repeatedly in the last ca 6000 yr significant longshore variation in the age of the lowest terrace from 1700 to 2500 yr b p suggests independent coseismic uplift on different sectors of the coast this is supported by age height relationships of the higher holocene terraces nonlinear uplift during the holocene with recurrence intervals increasing toward the present is clearly recorded by the regressive terraces in each subregion some of the holocene regressive terraces grade laterally into fluvial terraces capped with debris flow deposits probably reflecting seismically triggered mass movement

ice 3g a new global model of late pleistocene deglaciation based upon geophysical predictions of post glacial relative sea level change a new high resolution global model of late pleistocene deglaciation is inferred on the basis of geophysical predictions of postglacial relative sea level variations in which the ice ocean solid earth interaction is treated in a gravitationally self consistent fashion for the purpose of these analyses the radial viscoelastic structure of the planet is assumed known on the basis of previously published sensitivity tests on solutions of the forward problem only radiocarbon controlled relative sea level histories from sites that were actually ice covered with one or two additions are employed to constrain the model leaving relative sea level rsl data from sites that were not ice covered to be employed to confirm its consistency results for these confirmatory analyses are reported elsewhere here the new deglaciation model referred to as ice 3g is compared to previous models derived by several independent means and tested against a number of additional observations other than sea level histories including geologically controlled retreat isochrones oxygen isotope data from deep sea sedimentary cores and coral terrace elevations the latter two observations strongly constrain the net sea level rise that has occurred since the onset of deglaciation and therefore the mass of ice that melted during the last glacial interglacial transition

mantle viscosity from the simultaneous inversion of multiple data sets pertaining to postglacial rebound we employ the formalism of bayesian inference to invert three distinct data sets related to the phenomenon of glacial isostatic adjustment these data consist respectively of the wavenumber spectrum of relaxation times for fennoscandian rebound inferred by mcconnell 1968 the characteristic relaxation times for a large collection of distinct 14c dated relative sea level curves from both canada and fennoscandia and two data related to the rotational response of the earth to the glaciation deglaciation cycle of the late pleistocene ice age our analyses show that the spherically symmetric viscosity structure objectively inferred on the basis of these data from the once glaciated regions of the northern hemisphere allows a fully satisfactory reconciliation of both us east coast relative sea level observations as well as providing a close fit to the requirements of aspherical geoid anomalies

pleistocene deglaciation and the global gravity field the present day global geoid and free air gravity anomaly signals induced by pleistocene deglaciation have been computed along with their secular variations a revised green function has been derived for the free air gravity anomaly and we find that an earth model with a lower mantle viscosity moderately larger than the upper mantle value satisfies observations of the anomaly over both hudson s bay and fennoscandia using this preferred earth model the predicted global geoid anomaly map is characterized by peak negative values of 36 7 9 4 and 21 9 m over hudson s bay northern europe and antarctica respectively and by a smooth small amplitude 2 5 m upwarping over the major ocean basins although the field has an order of magnitude less power than the observed there is a significant spatial correlation between the two inferences of mantle viscosity based upon comparisons between the observed geoid and the computed deglaciation induced anomaly are hampered by the large contribution from the convective circulation in the earth s mantle we have subtracted a predicted convection signal based on seismic tomography from the observed field however the residual geoid thus obtained is still too inaccurate to allow a useful comparison in contrast a filtered version of the observed geoid over canada which retains harmonic coefficients in the degree range 10 22 provides a relatively uncontaminated by convection datum that is indeed satisfied by the preferred earth model secular variations in the earth s gravitational field are naturally dominated by the signal from pleistocene deglaciation in this respect we compute zonal harmonics in the expansion of the present day rate of change of the geopotential and in particular we predict amplitudes comparable to that for the previously analysed j2 for j4 j5 j7 and j9 all of which reflect in large part the response due to the melting of the massive laurentide ice sheet this analysis suggests that measurement of these higher order secular variations of the zonal harmonics based upon analysis of long time series of lageos ranging data could provide very useful additional constraints on the radial variation of mantle viscosity these constraints along with all others described in this paper are shown to be relatively insensitive to assumptions concerning the deflection during the glacial isostatic adjustment process of the 400 and 670 km density discontinuities within the mantle

pleistocene deglaciation and the earth s rotation a new analysis summary this paper is concerned with an analysis of the effect of the pleistocene glacial cycle upon the earth s rotation we demonstrate that two important geophysical observables may be explained as aspects of the rotational response to surface mass loading by ice sheets these are the astronomically observed non tidal component of the acceleration of planetary rotation and the secular drift of the rotation pole relative to the surface geography which is evident in the ils pole path the former observation is shown to provide an unambiguous constraint upon the viscosity of the planetary mantle and requires a depth dependence of this parameter which is the same as that which has previously been inferred in studies of the relative sea level variations and free air gravity anomalies associated with post glacial rebound the observed secular drift of the pole cannot provide an independent estimate of the viscosity of the mantle because it depends jointly upon mantle viscosity and lithospheric thickness with the viscosity profile fixed by the rebound data the observed speed of polar wander suggests a continental lithospheric thickness which is in excess of that appropriate for old ocean basins but nevertheless in accord with independent constraints

glacial rebound of the british isles iii constraints on mantle viscosity summary observations of sea level change since the time of the last glacial maximum provide important constraints on the response of the earth to changes in surface loading on time scales of 103 104 years this response is conveniently described by an effective elastic lithospheric thickness and effective viscosities for one or more mantle layers considerable trade off between the parameters describing these layers can occur and different combinations can give rise to comparable predictions of sea level change in particular the trade off between lithospheric thickness and upper mantle viscosity can be important and for any reasonable value for the lithospheric thickness a corresponding mantle viscosity structure can be found that gives a plausible comparison of sealevel predictions with observations in particular thin lithosphere models will lead to low estimates for the upper mantle viscosity while thick lithosphere models lead to high viscosity values however either solution may represent only a local minimum in the model parameter space and may not correspond to the optimum solution it becomes important therefore that in the inversion of observational data a comprehensive search is conducted throughout the entire model parameter space to ensure that the solution identified does indeed correspond to the optimum solution the sea level data for the british isles lend themselves well to such an inversion because of the relatively high quality of the data the good geographic distribution of the data relative to the former ice sheet and reasonable observational constraints on the dimensions of the former ice sheet and on its retreat furthermore because of the contribution to the sea level signal from the distant ice sheets as well as from the melt water load the observational data base for the region also has some resolving power for the viscosity of the deeper mantle the parameter space explored is defined by up to five mantle layers the lithosphere of effective elastic thickness d and a series of upper mantle layers i 2 4 extending down to depths of 200 400 and 670 km respectively each of viscosity qi and a lower mantle layer of viscosity qlm extending down to the coremantle boundary the range of parameters explored is 30 d i 120 km 3 x loi9 iqi

on postglacial sea level summary an exact method is presented for calculating the changes in sea level that occur when ice and water masses are rearranged on the surface of elastic and viscoelastic non rotating earth models the method is used to calculate the instantaneous elastic and delayed vi scoelastic sea level changes following the partial melting of late quaternary ice sheets we find that there can be large errors in the usual assumption that changes in sea level are uniform over the ocean basins if a quantity of ice equivalent to a uniform 100 m rise in sea level melts from the laurentide and fennoscandian ice sheets then in the south pacific the instantaneous rise in sea level can be as large as 120m in the north atlantic the instantaneous rise is always less than 100 m there is a zone in the north atlantic with almost no sea level change and near greenland and norway the sea level falls rather than rises by over 100 m one thousand years after the melting a forebulge migrating towards the ice loads causes water to flow from the south pacific into the north pacific suggesting that raised beaches should occur in the south pacific the gravitational attraction of an ice mass upon a nearby ocean tends to hold sea level high in the vicinity of the ice this extra load near the ice may have a significant influence on postglacial isostatic adjustment

constraints on mantle viscosity based upon the inversion of post glacial uplift data from the hudson bay region summary we outline two parametrizations for post glacial relative sea level rsl histories associated with previously glaciated regions the first parametrization is based on a site dependent normalization of the rsl history while the second involves the estimate of a site dependent logarithm of the inverse decay time for the exponential like form which characterizes these histories both parametrizations are shown to yield data sets which are relatively insensitive to the details of the late pleistocene surface load history and therefore inferences of mantle viscosity based upon them will be particularly robust we apply the two parametrizations to consider the rsl record at a number of sites across the hudson bay region in this regard our inferences which are derived from both forward and inverse calculations are based upon the actual rsl age height pairs obtained by survey rather than the highly subjective set of rsl trends adopted in previous studies one of the main goals of the analysis is to assess the validity of a set of previously published and highly contradictory inferences of the radial profile of mantle viscosity based on the hudson bay rsl record forward analyses using models with isoviscous upper and lower mantle regions as adopted in the vast majority of previous analyses indicate that the parametrized versions of the rsl record in hudson bay excluding data from the cape henrietta maria site are best fitted by a lower mantle viscosity near 1021 pa s the same conclusion holds when data from only northern hudson bay are considered the rsl record in southern hudson bay is not self consistent if the error bars adopted herein are reasonable however the parametrized versions of the rsl curves from each site in this region can be reconciled by a model with a lower mantle viscosity somewhere in the rather moderate range 0 5 3 0 1021 pa s the value of 3 1021 pa s represents a lower bound on the lower mantle viscosity required to fit the rsl records at cape henrietta maria this record is characterized by a longer decay time than those associated with other sites in the data base this lower bound is in contrast with previous suggestions that the rsl record at cape henrietta maria requires a lower mantle viscosity of 30 1021 pa s inverse analyses described herein indicate that the rsl record from the entire hudson bay constrains the average viscosity in the radial region extending from the lower reaches of the upper mantle to mid lower mantle depths the inferences listed above are therefore more properly ascribed to this region rather than the entire lower mantle we have found that the data are sensitive to moderately deeper variations in the radial viscosity profile as one considers sites situated further north in hudson bay this provides one means for a spherically symmetric model to yield decay times which vary across the hudson bay however the data do not provide an unambiguous requirement for such a variation

pleistocene sea levels and tectonic uplift based on dating of corals from sumba island indonesia we collected and dated fossil corals from the uplifted reef tracts of sumba island mass spectrometry tims and ams was used to obtain reliable and precise radiochronological results it is possible to identify outcrops corresponding to major sea level high stands isotope stage 1 starting at 7 1 kyr bp 5a c a 86 kyr bp 5c c a 105 kyr bp 5e between 119 and 132 kyr bp and 9 c a 305 kyr bp the radiochronological data are thus broadly consistent with the timing of sea level fluctuations predicted by the astronomical theory of paleoclimates comparison with similar data from other terrace sequences and from benthic 18o records indicates that the average tectonic uplift rate at cape laundi sumba is on the order of 0 2 0 5m kyr the lowest values 0 3m kyr lead to the build up of polycyclic topographic terraces in which it is possible to find corals grown during different interstadial or even interglacial times

holocene glacial rebound and sea level change in nw europe summary observations of late pleistocene and holocene sea level change relative to the crust exhibit very considerable variations across nw europe in consequence of the response of the earth s crust to the deglaciation of fennoscandia and of the water added to the oceans from the melting of all late pleistocene ice sheets inversion of sea level observations from a site near the centre of the fennoscandian ice sheet and from three sites located beyond the margin of the ice sheet at the time of maximum glaciation yield a range of plausible models for the earth s response and for the ice models further constraints on this range of models is placed by a comparison of observed sea levels with predicted values at other sites near the former ice sheet margins the resulting mantle parameters are upper mantle viscosity 3 5 x 10 pa s lower mantle viscosity 2 7 x loz1 pa s lithospheric thickness 100150 km these values represent effective parameters that describe the response of the earth to surface loading of short to intermediate wavelengths on a time scale of 104yr the lower mantle viscosity is poorly constrained but the marked increase from upper to lower mantle is a characteristic of all plausible solutions the inversion places a constraint on the total volume of ice in the fennoscandian ice sheet such that the equivalent sea level rise from this contribution is about 13 14 m a less well determined constraint of about 10 m equivalent sea level rise is suggested for the barents kara ice sheet the inversion also indicates that a small amount of melt water from ice sheets far away from europe continued to be added into the oceans during late holocene time so as to raise the equivalent sea level by about 3 m during the past 6000 yr consistent with similar inversions of data from sites in the australian and pacific regions

variations in the isotopic composition of molybdenum in freshwater lake systems variations in molybdenum isotopic composition spanning the range of similar to 2 3 parts per thousand in the terms of mo 97 mo 95 ratio have been measured in sediment cores from three lakes in northern sweden and north western russia these variations have been produced by both isotopically variable input of mo into the lakes due to mo isotopic heterogeneity of bedrock in the drainage basins and fractionation in the lake systems due to temporal variations in limnological conditions mo isotope abundances of bedrock in the lake drainage basins have been documented by analysis of mo isotope ratios of a suite of molybdenite occurrences collected in the studied area and of detrital fractions of the lake sediment cores the median delta mo 97 value of the investigated molybdenites is 0 26 parts per thousand with standard deviation of 0 43 parts per thousand n 19 whereas the median delta mo 97 value of detrital sediment fractions from two lakes is 0 40 parts per thousand with standard deviation of 0 36 parts per thousand n 15 the isotopic compositiori of mo in the sediment cores has been found to be dependent on redox conditions of the water columns and the dominant type of scavenging phases hydrous fe oxides have been shown to be an efficient scavenger of mo from porewater under oxic conditions oxidative precipitation of fe ii in the sediments resulted in co precipitation of mo and significant authigenic enrichment at the redox boundary in spite of a pronounced increase in mo concentration associated with fe oxides at the redox boundary the isotopic composition of mo in this zone varies insignificantly suggesting little or no isotope fractionation during scavenging of mo by hydrous fe oxides in a lake with anoxic bottom water a chironomid inferred reconstruction of o 2 conditions in the bottom water through the holocene indicates that increased o 2 concentrations are generally associated with low delta mo 97 mo 95 values of the sediments whereas lowered o 2 contents of the bottom water are accompanied by relatively high delta mo 97 mo 95 values thus confirming the potential of mo isotope data to be a proxy for redox conditions of overlying waters however it is pointed out that other processes including input of isotopically heterogeneous mo and mn cycling in the redox stratified water column can be a primary cause of variations in mo isotopic compositions of lake sediments less

some examples of the middle triassic marine transgression in south western mediterranean europe in southwestern mediterranean europe ne spain sardinia and nw apennines e g monti pisani and punta bianca the middle triassic transgression on a margin of an intracratonic basin appears to be controlled by the different morphology and tectonic activity of the margin itself

postglacial rebound from vlbi geodesy on establishing vertical reference difficulty in establishing a reference frame fixed to the earth s interior complicates the measurement of the vertical radial motions of the surface i propose that a useful reference frame for vertical motions is that found by minimizing differences between vertical motions observed with vlbi ma and ryan 1995 and predictions from postglacial rebound predictions peltier 1995 the optimal translation of the geocenter is 1 7 mm yr toward 36 n 11l e when determined from the motions of 10 vlbi sites this translation gives a better fit of observations to predictions than does the vlbi reference frame used by ma and ryan 1995 but the improvement is statistically insignificant the root mean square of differences decreases 20 to 0 73 mm yr and the correlation coefficient increases from 0 76 to 0 87 postglacial rebound is evident in the uplift of points in sweden and ontario that were beneath the ancient ice sheets of fennoscandia and canada and in the subsidence of points in the northeastern u s germany and alaska that were around the periphery of the ancient ice sheets