a first order global model of late cenozoic climatic change ii further analysis based on a simplification of co2 dynamics a model developed recently for the long term variations of global ice mass carbon dioxide and mean ocean temperature through the late cenozoic is simplified by hypothesizing a new equation for the co2 variations containing one less adjustable parameter but retaining the essential physical content of the previous equation including nonlinearity and the potential for instability by assuming plausible time constants for the glacial ice mass and global mean ocean temperature and setting the values of six adjustable parameters rate constants a solution for the last 5 my is obtained displaying many of the features observed over this period including the transition to the near 100 ky major ice age oscillations of the late pleistocene in obtaining this solution it is also assumed that variations in tectonic forcing lead to a reduction of the equilibrium co2 concentration perhaps due to increased weathering of rapidly uplifted mountain ranges over this period as a consequence of this co2 reduction the model dynamical system can bifurcate to a free oscillatory ice age regime that is under the pacemaker influence of earthorbital milankovitch forcing expanded discussions are given of the surface temperature variations accompanying the evolution of ice co2 and ocean temperature and of the bifurcation properties of the model from both mathematical and physical viewpoints

oxygen isotope calibration of the onset of ice rafting and history of glaciation in the north atlantic region we report here that dsdp site 552a cored with the hydraulic piston corer on the west flank of rockall bank recovered an undisturbed sequence of alternating white deep sea carbonate oozes and dark coloured layers that are rich in glacial debris oxygen isotope analysis of the sequence together with detailed nannofossil and palaeomagnetic stratigraphy shows that the first major horizon of ice rafting occurred at about 2 4 myr and was preceded by a minor pulse of ice rafting at about 2 5 myr the carbon isotope record shows that the site has been bathed by a water mass of similar characteristics to present day north atlantic deep water at least since 3 5 myr

isotopic modeling of climatic oscillations implications for a comparative study of marine and ice core records a simple climate model of quaternary glaciations is combined with a simple isotopic model in order to study phase relationships between 18o records in deep sea cores and in ice cores power spectra of the model s climatic and isotopic variables show a dominant peak near 100 kyr 1 kyr 1000 years the continuous background of all paleoclimatic spectra smaller peaks at the orbital periodicities of 41 23 and 19 kyr and additional peaks at 14 7 13 11 5 10 4 and 9 5 kyr time lags between the model s isotopic profiles are highly variable ice core events identified as sharp maxima or minima can lead apparently homologous events in the marine record or lag them by a few thousand years this time domain instability of lags is due to the complexity of the power spectrum with its continuous component and large number of lines since the complex spectrum is present in both model and paleoclimatic data we expect lag relations in the data to be quite variable as well independent dating of ice cores is preferable therefore to dating by comparison with marine records absolute dating can be based on a combination of ice flow models and of accelerator mass spectroscopy time averaged cross spectra of independently dated isotopic profiles should be stable and should help us to understand the mechanisms which lead to internal climatic variability

the influence of continental ice atmospheric co2 and land albedo on the climate of the last glacial maximum the contributions of expanded continental ice reduced atm9spheric co2 and changes in land albedo to the maintenance of the climate of the last glacial maximum lgm are examined a series of experiments is perform ed using an atmosphere mixed layer ocean model in which these changes in boundary conditions are incorporated either singly or in combination the model used has been shown to produce a reasonably reajistic simulation of the reduced temperature of the lgm manabe and broccoli 1985 by comparing the results from pairs of ex periments the effects of each of these envirornnental changes can be determined expanded continental ice and reduced atmospheric co2 are found to have a substantial impact on global mean temperature the ice sheet effect is confined almost ex clusively to the northern hemisphere while lowered co2 cools both hemispheres changes in land albedo over ice free areas have only a minor thermal effect on a global basis the reduction of co2 content in the atmosphere is the primary contributor to the cooling of the southern hemisphere the model sensitivity to both the ice sheet and co2 effects is characterized by a high latitude amplification and a late autumn and early winter maximum

the mid quaternary climatic transition as the free response of a three variable dynamical model abstract a simplified version of a previously described dynamical model governing global ice mass atmospheric carbon dioxide and mean ocean temperature that may also be a proxy for some other co2 controlling oceanic variable e g nutrient supply is shown to possess solutions in a realistic parameter range that can replicate the main features of the climatic variations implied by the full two million years quarternary 18o record these variations include a major transition between a low ice 18o low variance mode before roughly 900 kyr bp to a high ice near 100 kyr period variance mode after that time to the present the model contains only three free parameters in this simplified form no external earth orbital forcing is prescribed i e the model represents only internal dynamics from the previous studies it seems clear that additional variance representing such features as the rapid deglaciation and the phase of the major quaternary oscillation can be largely explained with no more

bore like surf beat on reef coasts reef coasts because the wave height of incoming stormy high waves are extremely dampened by coral reef have been believed as calm sea areas however many coastal structures have actually been damaged by stormy waves dampened by the reef surf beat with large wave height was discovered through field observations the surf beat has been named as bore like surf beat the wave height and the wave velocity of the bore like surf beat is lager than that of individual waves numerical simulation taking into account of non linear effect of surf beat phenomenon on the reef shows good agreement with laboratory data the main reason of the disasters of coastal structures on the reef coasts is existence of the bore like surf beat resonantly excited incoming wave groups

influence of late cenozoic mountain building on ocean geochemical cycles in a steady state ocean input fluxes of dissolved salts to the sea must be balanced in mass and isotopic value by output fluxes for the elements strontium calcium and carbon rivers provide the primary input whereas marine biogenic sedimentation dominates removal dissolved fluxes in rivers are related to rates of continental weathering which in turn are strongly dependent on rates of uplift the largest dissolved fluxes today arise in the himalayan and andean mountain ranges and the tibetan plateau during the past 5 m y uplift rates in these areas have increased significantly this suggests that weathering rates and river fluxes may have increased also the oceanic records of carbonate sedimentation level of the calcite compensation depth and delta sup 13 c and delta sup 87 sr in biogenic sediments are consistent with a global increase in river fluxes since the late miocene the cooling of global climate over the past few million years may be linked to a decrease in atmospheric co sub 2 driven by enhanced continental weathering in these tectonically active regions

evolution of atlantic pacific δ13c gradients over the last 2 5 m y the evolution of interocean carbon isotopic gradients over the last 2 5 m y is examined using high resolution 13c records from deep sea cores in the atlantic and pacific oceans over much of the northern hemisphere ice ages relative reductions in north atlantic deep water production occur during ice maxima from 2 5 to 1 5 ma glacial reductions in nadw are less than those observed in the late pleistocene glacial suppression of nadw intensified after 1 5 ma earlier than the transition to larger ice sheets around 0 7 ma at a number of times during the pleistocene 13c values at dsdp site 607 in the north atlantic were indistinguishable from eastern equatorial pacific 13c values from approximately the same depth odp site 677 indicating significant incursions of low 13c water into the deep north atlantic atlantic pacific 13c values converge during glaciations between 1 13 1 05 m y 0 83 0 70 m y and 0 46 0 43 m y this represents a pseudo periodicity of approximately 300 kyr which cannot easily be ascribed to global ice volume or orbital forcing this partial decoupling at low frequencies of the 18o and 13c signals at site 607 indicates that variations in north atlantic deep water circulation cannot be viewed simply as a linear response to ice sheet forcing

the evolution of pleistocene climate a time series approach abstract in the later part of the pleistocene variations in global ice volume have been dominated by an approximate 100 000 year cycle analysis of 2 myr long oxygen isotope record from an equatorial pacific core indicates that this is true only for the last 900 000 years prior to this time the amplitude of the 100 000 year cycle is much reduced as is the variance of all oscillations with periods greater than 60 000 years based on results of time series analysis of this 2 myr long record the pleistocene glacial cycles can be divided into three sections 1 the late pleistocene 0 900 kyr b p where the variations in the isotope record are dominated by the 100 000 year cycle 2 the middle pleistocene 900 1450 kyr b p in which low frequency components are not as important as in the later period of the pleistocene and 3 the early pleistocene late pliocene 1450 2000 kyr b p where general reductions of importance at all frequencies is seen as compared to the later intervals recent modeling efforts which describe variations in global ice volume show that the dominant low frequency component observed in the late pleistocene can result from different time constants for the rate of glacial growth and decay in response to variations in the earth s orbital parameters it is hypothesized that during the early pleistocene the rate of growth and decay of glaciers were more similar and that continental erosion by successive glacial advances lowered the land surface in areas of ice cap formation to below sea level when the ice caps became marine based more rapid decay of the ice became possible

sensitivity experiments with a model of the ice age cycle the response to milankovitch forcing abstract we address the question of the extent to which a simple one dimensional time dependent climate model is able to explain the fluctuations of continental ice volume which have occurred during the last quarter of the pleistocene epoch the model includes explicit ice physics based on the glen flow law and explicit glacial isostatic adjustment based on a realistic viscoelastic field theory when the model is forced by the milankovitch variations of insulation it delivers ice volume fluctuations which are very similar to those implied by the specmap 18o record in particular the model appears capable of explaining the observed correlation between eccentricity and ice volume variations the dominance of the 105 year period in the power spectrum of the ice volume time series and the thickness and latitudinal extents of the continental ice sheets which existed at wurm wisconsin maximum 18 000 years ago

a one dimensional model of continental ice volume fluctuations through the pleistocene implications for the origin of the mid pleistocene climate transition abstract a new one dimensional zonally averaged high resolution spherical model of the accumulation and flow of continental ice sheets in response to variations in insolation forcing has been constructed the model includes a detailed representation of the influence of the isostatic adjustment of the earth under the weight of the ice at a rate determined by the viscosity of the planetary mantle with this model we address the issue as to whether the combination of physical processes that it includes is sufficient to allow us to reconcile the variations of continental ice volume that are inferred to have dominated the pleistocene epoch on the basis of the analysis of 18o data from deep sea sedimentary cores we use high resolution data from the newly analyzed odp site 677 and from dsdp site 607 as principal observational constraints these datasets reveal a marked transition in climate to have occurred at mid pleistocene time about 900 000 years bp in order to understand this transition we find that we

sediment thicknesses and holocene glacial marine sedimentation rates in three east greenland fjords ca 68 n we compared measured and estimated sediment budgets in heavily glaciated fjords in east greenland mass balance calculations and regional glacio climatic conditions suggest that the sediment flux to the seafloor in kangerdlugssuaq and nansen fjords should be dominated by iceberg rafting and not by the rain out of suspended particulates in meltwater as the glacier calving flux is estimated at 15 and 2 km 3 yr compared to meltwater volumes of 4 4 and 1 7 km 3 yr respectively gravity cores in the three fjords indicate that the uppermost 1 2 5 m of sediment consists of diamictons or fine grained laminated muds ams radiocarbon dates on calcareous foramininfera or shells 16 total indicate sedimentation rates of 110 to 340 cm ka within the fjords over the last 1 ka and 10 20 cm ka during the holecene on the inner and middle shelf annual sediment discharge is around 0 67 times 10 6 tonnes yr within the kangerdlugssuaq fjord and trough system which translates into an average basin wide rate

forcing of late cenozoic northern hemisphere climate by plateau uplift in southern asia and the american west geologic evidence indicates that net vertical uplift occurred on a large kilometer scale and at accelerating rates during the middle and late cenozoic in plateaus of southern asia and the american west based on this evidence general circulation model sensitivity tests were run to isolate the unique effects of plateau uplift on climate the experiments simulated significant climatic changes in many places some far from the uplifted regions the basic direction of most of these simulated responses to progressive uplift is borne out by changes found in the geologic record winter cooling of north america northern europe northern asia and the arctic ocean summer drying of the north american west coast the eurasian interior and the mediterranean winter drying of the north american northern plains and the interior of asia and changes over the north atlantic ocean conducive to increased formation of deep water the modeled changes result from increased orographic diversion of westerly winds from cyclonic and anticyclonic surface flow induced by summer heating and winter cooling of the uplifted plateaus and from the intensification of vertical circulation cells in the atmosphere caused by exchanges of mass between the summer heated and winter cooled plateaus and the mid latitude oceans disagreements between the geologic record and the model simulations in alaska and the southern rockies and plains may be related mainly to the lack of narrow mountain barriers in the model orography taken together the observed regional trends comprise much of the pattern of late cenozoic climatic deterioration in the northern hemisphere that culminated in the plio pleistocene ice ages the success of the uplift sensitivity experiment in simulating the correct pattern and sign of most of the observed regional climatic trends points to uplift as an important forcing function of late cenozoic climatic change in the northern hemisphere at time scales longer than orbital variations however the modest amplitude of the uplift induced cooling simulated at high latitudes indicates a probable need for additional climatic forcing

stability of the astronomical frequencies over the earth s history for paleoclimate studies the expected changes over the past 500 million years in the principal astronomical frequencies influencing the earth s climate may be strong enough to be detectable in the geological records and such effects have been inferred in several cases calculations suggest that the shortening of the earth moon distance and of the length of the day back in time induced a shortening of the fundamental periods for the obliquity and climatic precession from 54 to 35 41 to 29 23 to 19 and 19 to 16 thousand years over the last half billion years at the same time the precessional constant increased from 50 to 61 arc seconds per year the changes in the frequencies of the planetary system due to its chaotic motion are much smaller their influence on the changes of the periods of climatic precession obliquity and eccentricity of the earth s orbit around the sun can be neglected eccentricity periods used for quaternary climate studies may therefore be considered to have been more or less constant for pre quaternary times

climate evolution in the pliocene and pleistocene from marine sediment records and simulations internal variability versus orbital forcing three deep sea 18o records are compared in the temporal and spectral domain with each other a simple ice sheet and bedrock climate model simulation and insolation anomalies at 65 n the model is one dimensional it represents the accumulation and ablation processes of ice its viscoplastic flow and the isostatic adjustment of bedrock under it being forced solely by orbital insolation anomalies the objective of the comparison is to determine the importance of internal variability arising from feedback between the ice sheets and bedrock through the late pliocene and pleistocene in the simulations and in the climate record it is found that a considerable part of the spectral variability at suborbital frequencies i e below the precessional and obliquity frequencies as seen in the sediment cores is simulated by the model simulated variability near 100 kyr has its period determined by the difference between two precessional frequencies the model s internal variability mechanism has a period near 30 50 kyr but is highly damped still it enhances suborbital variability near 100 kyr and 300 kyr and appears to induce nonstationarity of the orbital and suborbital peaks on even longer time scales since the insolation forcing exhibits no noticeable peaks at suborbital or lower frequencies it appears that the nonstationarity in the model simulations as well as the suborbital peaks themselves are due to coupling between model nonlinearities and the orbital forcing since the variability at 100 kyr in the deep sea records is correlated with the simulated variability such a control by the orbital forcing may be also operating in the climate system the absence from the simulation of a trend toward larger ice sheets and colder climate seen throughout the pleistocene in the sediment record is possibly due to slow tectonics and co2 forcing changes not included in the model

a model study of the role of high latitude topography in the climatic response to orbital insolation anomalies abstract the response of a combined global zonal averaged energy balance model and a continental ice sheet model to insulation anomalies produced by orbital perturbations is examined although the role of ice sheet dynamics is included in the model this paper focuses mostly on the response of the atmosphere hydrosphere to the perturbations the model developed from that of suarez and held 1979 explicitly incorporates surface topography surface energy balances snow and sea ice budgets and variable static stability the ice sheet model is based on a power law rheology it is found that the presence of a high latitude plateau significantly increases the sensitivity of the climate to the insolation perturbation the sensitivity is maximized when the elevation of the summer snowline is near the elevation of the plateau as appears to be the case with baffin island today there appears to be a rapid response of the seasonal cycle to an insolation anomaly through albedo temperature feedback and with the

slow physics of large continental ice sheets and underlying bedrock and its relation to the pleistocene ice ages very simple two component models of the earth s cryosphere lithosphere system are analyzed the models consist of a continental ice sheet coupled to a responding earth model because of relatively long response times of the bedrock at least over shorter wavelengths long period of the order of 50 000 years internal or free oscillations in the system occur over a wide range of parameters and parameterizations the oscillatory mechanism operates primarily in the ablation zone and in particular at the firn line these free oscillations could play a dynamic role in the response of the climate system to eccentricity forcing a simple analysis presented strongly suggests that for the simplest rheological models of the earth response to pleistocence ice sheets elastic response is confined to a layer considerably less than the thermal thickness of the lithosphere furthermore because of the striking lack of geological evidence for large depressions of bedrock at the leading edge of growing continental ice sheets predicted by viscoelastic models of the earth more sophisticated rheological models of the earth s uppermost layer appear necessary

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

a climate model with cryodynamics and geodynamics we present a highly simplified zero dimensional model of the climatic system the model attempts to incorporate mechanisms important on the time scale of glaciation cycles namely 104 105 years in particular the radiation balance of the ocean atmosphere the plastic flow of continental ice sheets and the viscous flow of the upper mantle are taken into account the stress is on the interaction between the cryosphere ice sheets and the asthenosphere upper mantle which had not been taken into account in previous energy balance or ice sheet models of climate weertman 1976 kallen et al 1979 the model exhibits free self sustained oscillations of an amplitude and period comparable to those found in the paleoclimatic record of glaciations viz o 10 deg k and o 104 years respectively such oscillations had already been found in the simpler model of kallen et al 1978 1979 the fact that nonlinear self sustained oscillations also obtain for the present slightly more complex model is interesting it offers mild support to the idea that unforced oscillations can actually exist in the real climatic system itself these oscillations are characterized by a quarter phase lag between ice extent and global temperature given that the climatic system is a nonlinear oscillator with a period of o 104 years it is expected that astronomical forcing at 19 000 23 000 and 41 000 years will lead to subharmonic oscillations with a period close to 100 000 years the dominant period of glaciation cycles such effects have already been noticed in the work of birchfield and weertman 1978 with an even simpler model including only glacial dynamics the careful study of the interplay between internal mechanisms and external forcing with comparable time scales represents an interesting challenge to the theory of ice ages the results of such theoretical studies can also provide guidance for observational work in particular periodicities of o 104 years in the climate record and phase lags between ice volume and temperature appear to be supported by the most recent paleogeochemical investigations ruddiman and mclntyre 1981

high frequency paleovariability in climate and co2 levels from vostok ice core records the high resolution of the vostok records provides a unique look at the causes of paleoclimatic variability during the last complete glacial cycle the records present strong evidence for the interaction between orbital forcing and internal physico chemical mechanisms of variability this interaction appears to account for the great wealth of spectral features found in the records

a model for atmospheric co 2 over phanerozoic time a new model has been constructed for calculating the level of atmospheric co sub 2 over phanerozoic time which is much simpler mathematically than the blag model but more complex geologically and biologically results suggest that there has been a notable pattern of varying atmospheric co sub 2 level over the past 570 my with high levels during the mesozoic and early paleozoic and low levels during the permo carboniferous and late cenozoic sensitivity analysis shows that within reasonable limits of the principal governing parameters this qualitative trend of varying co sub 2 is relatively insensitive to the actual values chosen for these parameters causes of the co sub 2 variation are multiple and no single geological or biological process can be called upon to explain all co sub 2 variation for all time the calculated trend of co sub 2 over time agrees well with independent estimations of paleoclimates thus the greenhouse theory of paleoclimate on a long geological time scale is supported by the results of the present study

eccentricity forcing of pliocene early pleistocene climate revealed in a marine oxygen isotope record milankovitch theory that climate is controlled by variations in the earth s orbital parameters has gained wide acceptance for its ability to account for two climate cycles a 23 kyr cycle that is phase locked to the precession driven insolation cycle and a 41 kyr cycle that is phase locked to the obliquity driven insolation cycle1 6 but explaining the observed 100 kyr climate cycle in terms of milankovitch theory especially for the late pleistocene ice age cycle remains controversial in spite of a strong correlation with the 100 kyr cycle in the earth s orbital eccentricity5 one problem is that eccentricity affects insolation mainly by modulating the precession cycle its direct contribution to radiation change is too small 0 1 to effect the observed climate change directly5 7 another is the absence of a late pleistocene ice volume cycle in oxygen isotope records to match the 404 kyr component of the eccentricity cycle5 8 here we examine an oxygen isotope record spanning the interval 1 2 to 5 2 million years ago before the late pleistocene ice age regime we find 404 kyr and 100 kyr climate cycles which are coherent with eccentricity and which have amplitudes that are similar to the coexisting 23 kyr cycle analysis of these low frequency cycles suggests that they originate through an asymmetrical response mechanism that preferentially introduces variance into the climate system from the warmer portions of the eccentricity modulated precession cycle our data thus support eccentricity s role in the origin of low frequency oxygen isotope cycles before the late pleistocene ice age

exploring the structure of the mid pleistocene revolution with advanced methods of time series analysis the mid pleistocene climate transition is a complex global change leading to late pleistocene ice ages with increased mean ice volume and dominant 100 ka cycle a thorough understanding of this transition demands quantitative investigations in the time and frequency domains and in the stochastic domain three methods of time series analysis are presented which have been adapted for this purpose they are tested by means of predefined artificial time series and applied to benthic oxygen isotope d18o records which serve as ice volume indicator results are as follows a the time dependent mean shows an increase of 0 35 vs pdb from 942 to 892 ka b evolutionary spectral analysis reveals an abrupt in crease of 100 ka cycle amplitude at approximately 650 ka c probability density function exhibits a bifur cation behavior at approximately 725 ka these find ings point to a multiple transition from a more linear climate system to a strong nonlinear system the signif icant lead of the transition in mean in relation to the 100 ka amplitude change and bifurcation is left open for explanation

synchrony and causal relations between permian triassic boundary crises and siberian flood volcanism the permian triassic boundary records the most severe mass extinctions in earth9s history siberian flood volcanism the most profuse known such subaerial event produced 2 million to 3 million cubic kilometers of volcanic ejecta in approximately 1 million years or less analysis of 40 ar 39 ar data from two tuffs in southern china yielded a date of 250 0 0 2 million years ago for the permian triassic boundary which is comparable to the inception of main stage siberian flood volcanism at 250 0 0 3 million years ago volcanogenic sulfate aerosols and the dynamic effects of the siberian plume likely contributed to environmental extrema that led to the mass extinctions

sedimentary record of rapid climatic variability in the north atlantic ocean during the last glacial cycle comparisons between a greenland ice core isotopic record and marine sediment grey level records of the north atlantic ocean show that rapid temperature variations witnessed by the ice core dansgaard oeschger events during the last glacial period can also be detected in marine sediments this shows that the resolution obtained in marine cores can be sufficiently high to record rapid climatic fluctuations several spectral analyses of those grey level marine records consistently indicate that significative fluctuations of the climatic response exist with periodicities of 5 to 1 kyr in addition to orbital forcing these high frequencies are close to those predicted by various dynamic ocean models