eolian deposition in the northeast pacific ocean cenozoic history of atmospheric circulation the mass accumulation rate and grain size distribution of the total eolian component of north pacific pelagic clays have been used to evaluate changes in eolian sedimentation and in the intensity of atmospheric circulation that have occurred during the past 70 m y eolian deposition was relatively low approximately 12 mg cm2 1 000 yr in the late cretaceous and paleocene reflecting an increasingly humid environment a peak in accumulation rates across the paleocene eocene boundary may be the result of episodic volcanism in mexico a sharp decrease in grain size across the paleocene eocene boundary is not readily interpreted fine eolian grain size and the low accumulation rates averaging 8 mg cm2 1 000 yr in the eocene and early oligocene are in agreement with low early tertiary thermal gradients and less vigorous atmospheric circulation large increases in grain size between 37 and 30 m y ago 20 and 10 m y ago and 5 m y ago to the present suggest that atmospheric circulation intensified by a factor of 1 5 1 8 and 1 4 respectively at these times eolian accumulation rates increased from 8 mg cm2 1 000 yr at the eocene oligocene boundary to approximately 25 mg cm2 1 000 yr at the end of the miocene with significant changes at 25 and 13 m y ago increased rates about 25 m y ago probably reflect movement of the core site into the influence of the westerlies and the increase 13 m y ago appears to be a response to steepening thermal gradients and intensifying atmospheric circulation resulting from major expansion of antarctic ice in the mid miocene along with a general increase in global aridity in the late cenozoic an order of magnitude increase in eolian accumulation rates 25 to 226 mg cm2 1 000 yr in pliocene pleistocene time coincides with the onset of northern hemisphere glaciation and consequent climatic deterioration the large increase in accumulation rates in pliocene pleistocene time coincides with a relatively smal increase in the intensity of atmospheric circulation compared to that in oligocene and miocene time suggesting that the increased deposition rates were the result of an increase in global aridity in glacial times rather than a further increase in the vigor of the global wind systems

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

back arc extension in the scotia sea the nature of back arc extension in the east scotia sea is re examined with the use of an enlarged geophysical data set well developed oceanic magnetic lineations confirm that the present spreading episode started about 8 ma ago that spreading is asymmetric and that the total rate increased from 50 to 70 m m a about 1 5 ma ago most of the currently active south sandwich volcanic island arc lies upon ocean floor only 6 8 m a old and generated at the current spreading ridge subsequent extension has not modified the curvature of the arc east west magnetic lineations of miocene age in the central scotia sea and contemporaneous low k arc tholeiites dredged from the eastern south scotia ridge discovery bank indicate a regime of coupled subduction and back arc extension preceding that occurring now a speculative model involving a series of collisions of parts of this earlier discovery trench with ridge crest sections of the south american antarctic plate boundary explains the transformation of this earlier regime into the present self contained sandwich plate regime the considerable small scale variability observed in the back arc region may be seen as an inevitable consequence of the action of the ridge trench collision mechanism the entire scotia sea could have formed by a similar kind of back arc extensional modification of the south american antarctic plate boundary

relations among subduction parameters clues to the dynamics of the subduction process are found in the many measurable parameters of modern subduction zones based on a critical appraisal of the geophysical and geological literature 26 parameters are estimated for each of 39 modern subduction zones to isolate causal relationships among these parameters multivariate analysis is applied to this data set this analysis yields empirical quantitative relations that predict strain regime and strike slip faulting in the overriding plate maximum earthquake magnitude benioff zone length slab dip arc trench gap and maximum trench depth excellent correlation is found between length of the benioff zone and the product of convergence rate and age of the downgoing slab this relationship is consistent with the conductive heating model of molnar et al 1979 if the model is modified in one respect the rate of heating of the slab is not constant it is substantially slower during passage of the slab beneath the accretionary prism and overriding plate the structural style in the overriding plate is determined by its stress state though the stress state of overriding plates cannot be quantified one can classify each individual subduction zone into one of seven semiquantitative strain classes that form a continuum from strongly extensional class 1 back arc spreading to strongly compressional class 7 active folding and thrusting this analysis indicates that strain class is probably determined by a linear combination of convergence rate slab age and shallow slab dip interplate coupling controlled by convergence rate and slab age is an important control on strain regime and the primary control on earthquake magnitude arc parallel strike slip faulting is a common feature of convergent margins forming a forearc sliver between the strike slip fault and trench optimum conditions for the development of forearc slivers are oblique convergence a compressional environment and a continental overriding plate the primary factor controlling presence of strike slip faulting is coupling strongly oblique convergence is not required the rate of strike slip faulting is affected by both convergence obliquity and convergence rate maximum trench depth is a response to flexure of the underthrusting plate the amount of flexural deflection at the trench depends on the vertical component of slab pull force which is very sensitive to slab age and shallow slab dip shallow slab dip conforms to the cross sectional shape of the overriding plate which is controlled by width of the accretionary prism and duration of subduction deep slab dip is affected by the mantle trajectory established at shallow depth but may be modified by mantle flow much of the structural complexity of convergent margins is probably attributable to terrane juxtaposition associated with temporal changes in both forearc strike slip faulting and strain regime empirical equations relating subduction parameters can provide both a focus for future theoretical studies and a conceptual and kinematic link between plate tectonics and the geology of subduction zones

trench migration and subduction zone geometry seismic tomography results show a large variety in upper mantle structure along convergent plate boundaries we numerically investigate the effect of trench migration on the evolution of a slab with temperature dependent viscosity encountering a viscosity interface we find that subduction zone geometry is sensitive to even small rates 1 cm yr of retrograde motion increase in trench migration rate decreases the dip angle of the slab and its ability to penetrate the lower mantle upon including a background mantle flow it turns out that trench migration relative to the upper mantle flow is more decisive than the absolute plate velocities subduction zone geometry appears to be strongly time dependent we conclude that individual tectonic setting and time dependent slab behaviour can account for many different types of observable subduction zone anomalies our model temperature can adequately account for magnitudes and patterns of seismic anomalies as obtained from seismic tomography

ablative subduction a two sided alternative to the conventional subduction model subduction is conventionally modelled as a convergence and overlapping of two semirigid plates such a convergence is one sided in that one lithospheric plate descends while the other remains on the surface we suggest an alternative model one based on the dynamics of fluids in this model the viscous lower lithosphere flows downward and the brittle upper lithosphere deforms in passive response this process is potentially double sided for we find that even a buoyant plate can be dragged downward by a dense descending neighbor thus an apparent overriding plate may be worn away by a process of viscous ablation with the rate of ablation a function of plate buoyancy we call this process ablative subduction ablative subduction allows us to simply interpret observations concerning slab profiles interplate seismicity back arc tectonics and complex processes such as double subduction and subduction polarity reversal in performing experiments modelling the evolution of simple fluid slabs we find that slab profile is strongly influenced by ablation in the overriding plate when ablation is weak as when a buoyant continent borders the trench deformable slabs adopt shallow andean style profiles these profiles develop over time from an initially steep shape more vigorous ablation rates as might occur in an ocean ocean convergence yield steeper marianas style profiles thus differences between marianas style and andean style slabs may result from differences in ablation rate and subduction duration the occurrence of ablation might not be easily detected at depth because material from subducting and ablating plates adhere closely as a single slab these slabs show downdip deformations that are consistent with seismic focal mechanisms plate deformations associated with ablative subduction are also consistent with observed patterns of seismicity in particular the relative aseismicity of marianas style plate boundaries is consistent with the ablative model we may thus explain how plates can converge aseismically without requiring that these plates be decoupled ablation should be weakest in andean style subduction although pulses of rapid ablation might lead to episodes of crustal shortening in oceanic regions the manifestation of ablation may depend upon tectonic regime in an extensional area like the marianas ablation may be temporally self limiting and ablative cycles may explain conflicting observations of tectonic erosion and accretion in a compressional regime vigorous ablation might lead to observations of double subduction and subduction reversal in each case the behavior of the lower lithosphere is relatively simple what varies is the rate of ablation and the response of the upper lithosphere

late pliocene onset of glaciation ice rafting and diatom stratigraphy of north pacific dsdp cores abstract the revised diatom biostratigraphy for the north pacific allows the timing of ice rafting as recorded in deep sea drilling project cores to be determined significant ice rafting began at the time of the matuyama gauss reversal boundary about 2 48 ma other bio and lithostratigraphic data from the pacific atlantic and arctic oceans and the black sea all indicate that northern hemispherical cooling and the accumulation of continental ice began at 2 4 2 5 ma we interpret all these data to indicate northern hemisphere ice cap formation at that time if so the closing of the isthmus of panama approximately 3 2 m y ago was not the immediate cause of the plio pleistocene glaciations

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

evolving pliocene pleistocene climate a north pacific perspective abstract siliceous microfauna in sediments from north pacific deep sea drilling project sites record the response of the regional surface and near surface waters to global climate change over the past 3 0 ma during the late pliocene 2 7 3 0 ma the northern subtropical and subpolar pacific experienced warmer climate conditions than those of today the first predominant polar assemblage maximum occurs in the north pacific at 2 45 ma signifying this region s response to the initiation of significant northern hemisphere glaciation north pacific faunal data show that late pliocene 1 8 2 5 ma glacial maxima were only about half to two thirds as severe as those which occurred during the late pleistocene low faunal assemblage variability during a 200 ka interval in the latest pliocene 1 6 1 82 ma in the north pacific indicates minimal variations in the surface water properties and climate conditions most like those which existed just prior to 2 5 ma a distinct change in faunal assemblage registered in sediments from site 580 between 1 2 and 1 2 ma marks the first cooling in the north pacific comparable in intensity to those of the late pleistocene since 1 26 ma high amplitude high frequency variations in the north pacific polar assemblage have continued to occur providing further evidence that the interglacial glacial extremes that characterize the late pleistocene have persisted since this time the frequency distribution of north pacific winter sea surface temperature time series for the late pliocene and early pleistocene includes major concentrations of variance power centered at a frequency corresponding to a 41 ka period with variance dominance of the temperature time series for the last 800 ka shifting to lower frequencies 110 ka period

horizontal shear flow in the mantle beneath the tonga arc in the tonga subduction zone a systematic shear deformation in which deep material moves relatively to the south is superimposed on the mode of down dip shortening which is characteristic of the region we interpret this observation together with a tectonic reconstruction in terms of the interaction between the subducting lithospheric slab and a horizontal shear flow in the mantle

controls of the structure of subducted slabs numerical simulations of subducting slabs are formulated in which the shape and dip of the slab are determined by the dynamics of the flow rather than imposed a priori the dip of slabs is a function of the time since the initiation of subduction slabs fold develop a kink in dip and thicken on entry into a high viscosity lower mantle comparison of the simulations with seismic observations suggest that the lower mantle is at least 10 30 times more viscous than the upper mantle

nonlinear response of pliocene climate to orbital forcing evidence from the eastern equatorial pacific a high resolution radiolarian record from the eastern equatorial pacific is used to examine paleoceanographic events during the pliocene from 1 8 to 4 9 ma these data provide a means to evaluate the response of the equatorial pacific to the onset of northern hemisphere glaciation near 2 47 ma and to the closure of the isthmus of panama near 3 2 ma responses are recorded both as long term trends and as variations within the milankovitch band 104 to 105 years radiolarian based sea surface temperature sst estimates indicate a mean cooling from lower to upper pliocene but no significant changes occur near the onset of northern hemisphere glaciation and no abrupt changes are observed before or after the closure of the isthmus spectral and cross spectral analyses of this sst time series indicate concentrations of variance within the precession 23 kyr and obliquity 41 kyr bands in several but not all time intervals variance is commonly concentrated at lower frequencies cross spectral analysis of sst time series and calculated orbital variations indicates that the sst record is not coherent with and therefore not linearly related to orbital variations detailed time series analysis of the sst and 18o time series indicates the absence of a direct linear relationship but rather the presence of a strongly nonlinear relationship between pliocene climatic variables and earth s orbit high coherence is seen between the amplitude of the 23 kyr component of sst and the amplitude modulation of precession the amplitude modulation of the 41 kyr component of the sst record is dissimilar to the amplitude modulation of obliquity but is extremely similar to and coherent with the amplitude of the 23 kyr component of sst the amplitude modulations of sst and 18o at a range of frequency bands have the same features features that roughly follow the modulation of precession this suggests a strongly nonlinear response of equatorial pacific sst and planktonic 18o to orbital forcing in the pliocene such a response extending over a number of frequency bands has not previously been seen in a paleoclimatic data set but may be consistent with proposed models describing paleoclimatic variability although the reasons for this response are still unclear it is evident that strong nonlinearities in climate are present during this time

uplift of the central andean plateau and bending of the bolivian orocline the topography of the central andes can be considered the primary tectonic signal of late cenozoic mountain building in an arid region where the effects of uplift and magmatism are little obscured by denudation the spatial coverage of the topographic signal is more complete than that for sparsely sampled geological and geophysical data a color coded image of digitized topography between 12 s and 37 s highlights the altiplano puna one of the world s most remarkable plateaus and reveals important physiographic clues about the formation of that major feature the topographic data combined with information on structure magmatism seismicity and paleomagnetism support a simple kinematical model for the late cenozoic evolution of the central andes the model does not require collisional effects or enormous volumes of intrusive additions to the crust but instead calls upon plausible amounts of crustal shortening and lithospheric thinning the model interrelates andean uplift a changing geometry of the subducted nazca plate and a changing outline in map view of the leading edge of the south american plate crustal shortening has accommodated convergence between the chilean peruvian forearc and the south american foreland the altiplano puna plateau can be constructed by a combination of crustal shortening and thickening and lithospheric thinning above a shallow dipping 20 30 subducted plate the seawardly concave bend of western south america the bolivian orocline was enhanced but not completely produced by an along strike variation in the amount of late cenozoic shortening maximum shortening in bolivia both produced the widest part of the plateau and increased the seaward concavity of the bolivian orocline the along strike variations of shortening are hypothesized to result from corresponding along strike variations in the width of a weakened zone in the overriding plate weakening occurs above the wedge of asthenosphere located between the subducted and overriding plates hence the width of the zone of weakening depends upon the dip of the subducted plate two types of shortening are recognized 1 a widespread basin and range laramide like shortening that characterizes modem activity in the sierras pampeanas and late miocene deformation of the altiplano puna and 2 on the eastern side of the cordilleras and plateau an east verging foreland fold thrust belt in which the underthrust foreland compresses and thickens the ductile lower crust and produces a plateau uplift of the upper crust the second type of shortening can be applied to plio quaternary deformations throughout the central andes but with a substantial narrowing of the region of plateau uplift in peru and south of 28 s a proposed monoclinal flexure of the upper crust on the western side of the plateau uplift explains the remarkably simple and regular morphology of the main western slope of the central andes the monocline is located above the tip of the asthenospheric wedge between the converging plates it is postulated to occur above the western limit of lower crustal thickening in the regions of horizontal subduction the monocline can be associated with a late miocene asthenospheric wedge tip

geologic time series from deep sea sediments time scales and distortion by bioturbation abstract in order to evaluate the time and frequency domain characteristics of geologic time series obtained from deep sea sediments two important problems must be considered first an accurate time scale must be established for the sediment records to be studied an accurate time scale is a basic requirement for such techniques as spectral and cross spectral analysis secondly the distortion of a geologic time series by sedimentation processes such as bioturbation must be considered if inferences are to be made about the mechanisms controlling climatic and oceanographic processes reflected in the geologic record recently it has been suggested that the periodic components in paleoclimatic records which can be related to variations in the earth s orbital parameters can be used as a means for providing very detailed control for late pleistocene sediment sections a time scale is estimated by tuning an individual frequency component in a climate record to an orbital parameter the tuning procedure commonly uses the technique of bandpass filtering which removes all frequency components except the one equivalent to an orbital parameter and then adjusting the geologic time scale so that the filtered record and orbital record align since bandpass filtering will produce a periodic record from any data set a question is raised as to the circularity of this tuning approach a series of numerical experiments using synthetic time series are used to determine if the tuning of a geologic record artificially induces improved correlation and increased periodic character in a climatic record these experiments indicate that the tuning does not artificially induce these effects thus it is suggested that the careful use of orbital tuning provides a powerful means of defining a geologic time scale the degree to which bioturbation distorts a paleoclimatic record is studied by comparing the spectra of oxygen isotopic records obtained from cores with sedimentation rates ranging from 1 to 12 cm 1000 yr 1 it is shown that the systematic change in the variance spectra of these records that is predicted by models of bioturbation cannot be observed for sedimentation rates greater than 2 cm 1000 yr 1 this can be explained in part if bioturbation is a function of sedimentation rates if this is the case then bioturbation must proceed to depth of tens to hundreds of centimeters however it seems likely that the degree to which bioturbation alters the climate signal preserved in deep sea sediment cores has been overestimated

a numerical assessment of slab strength during high and low angle subduction and implications for laramide orogenesis a finite element computer program for one dimensional heat flow was modified to apply to two dimensional subduction where conductive heat transport parallel to the upper and lower plates was neglected one dimensional finite element columns were cut displaced and reassembled to model heat advection calculated temperature depth dependent pressure stress levels associated with dislocation creep in olivine at constant strain rate and a 1 gpa deviatoric stress maximum were used to calculate the form of the brittle plastic failure envelope for profiles across subducted slabs the failure envelopes were then used to calculate maximum sustainable compression parallel to the slab model results suggest the following 1 in high angle subduction without frictional heating and horizontal heat flow temperatures along the interplate contact zone remain below about 120 c which is far below temperatures necessary for olivine plasticity without substantial frictional heating or lateral heat flow from the adjacent magmatic arc the subducted slab will gain rather than lose strength before its upper surface contacts asthenosphere in addition viscous coupling due to olivine plasticity and associated ablation of the upper plate will be unimportant slab strength reduction due to heating occurs rapidly after the top of the slab comes into contact with asthenosphere 2 in low angle subduction the top of the slab remains in contact with overlying lithosphere and does not contact asthenosphere slab strength reduction due to conductive heating from above occurs slowly because the upper plate is cooled by the lower plate and is the source of progressively less heat temperatures along the interplate contact zone may be sufficient for olivine plasticity at large distances from the trench and viscous coupling between plates and associated ablation of the upper plate may be important in this tectonic setting 3 low angle subduction for distances of over 1500 km during the laramide orogeny in southwestern north america appears viable at 50 ma slab strength beneath the black hills of south dakota is predicted to have been comparable to or greater than strength at the trench

a quantitative study of lower cretaceous cyclic sequences from the atlantic ocean and the vocontian basin se france milankovitch type periodicity was identified within lower cretaceous cyclic sequences from the western central atlantic ocean and the vocontian basin southeast france through application of the walsh spectral method to detailed stratigraphic sections the repetitious bedding in these sequences appears to be the product of cyclic climatic changes controlled by the orbital forcing obliquity seems to be the most important and persistent orbital element responsible for cyclic sedimentation whereas precession is less important upper berriasian and valanginian strata from the atlantic ocean contain indications of eccentricity induced periodicity and of a nonlinear response of the climatic system variations in orbitally tuned pelagic sedimentation rates in the atlantic ocean correlate with the changes in the character of cycles and with the occurrence of calcareous turbidites therefore this suite of sedimentary features are responding to evolving paleoceanographic conditions and possible relative sea level fluctuations the orbitally tuned durations of the valanginian and hauterivian stages in the vocontian basin are a minimum of 5 9 and 5 3 m y respectively

deep seismicity and modes of deformation in tonga subduction zone the geometry of deep 400 km earthquakes in the subducting western edge of the pacific plate when related to the local distribution of seismic activity reveals a complex pattern of planes of shear failure enabling the modes of deformation by which the descending plate shortens and thickens to be identified the deepest part of the slab is displaced laterally suggesting that it fails to penetrate the discontinuity at 670 km

collision rotation and the initiation of subduction in the evolution of sulawesi indonesia the island of sulawesi indonesia has been shaped and deformed as a result of collision with the sula platform a sliver of continental material from the northern margin of australia new guinea the collision has resulted in rotation of the north volcanic arm of sulawesi and the development of the accretionary wedge of the north sulawesi trench the north sulawesi trench changes laterally from a zone of no active deformation in the eastern part to a wide accretionary wedge in the west early stages of thrusting produce a steep frontal slope 8 16 indicative of relatively high basal shear stress whereas the more advanced western zone of thrusting produces a gentle 2 slope consistent with low basal shear stress reported paleomagnetic data suggest post late eocene counter clockwise rotation of the north arm and the offshore geophysics are explained by a pivot of the north arm with respect to the celebes basin about the eastern end of the arc convergence between the north banda basin and southeast sulawesi is documented by the presence of the tolo thrust its outcrop is strongly arcuate and its accretionary wedge varies in width from a minimum of a few kilometers at each end to a maximum of 30 40 km in the central part the northern end transforms to the leftlateral matano fault with a reported offset of 20 km the southern end of the thrust projects toward the deformed rocks of buton but the structural relations there are not clear the matano fault zone appears to connect westward with the palu fault which forms the western transform of the north sulawesi trench the palu matano fault system acts as a trench trench transform between the north sulawesi trench and the tolo thrust and this system is described by the same rotation pole as that for the sulawesi north arm

convection in three dimensions with surface plates generation of toroidal flow this work presents numerical calculations of mantle convection that incorporate some of the basic observational constraints imposed by plate tectonics the model is three dimensional and includes surface plates it allows plate velocity to change dynamically according to the forces which result from convection we show that plates are an effective means of introducing a toroidal component into the flow field after initial transients the plate motion is nearly parallel to transform faults and in the direction that tends to minimizes the toroidal flow field the toroidal field decays with depth from its value at the surface the poloidal field is relatively constant throughout the layer but falls off slightly at the top and bottom boundaries layered viscosity increasing with depth causes the toroidal field to decay more rapidly effectively confining it to the upper low viscosity layer the effect of viscosity layering on the poloidal field is relatively small which we attribute to its generation by temperature variations distributed throughout the system the generation of toroidal flow by surface plates would seem to account for the observed nearly equal energy of toroidal and poloidal fields of plate motions on the earth a low viscosity region in the upper mantle will cause the toroidal flow to decay significantly before reaching the lower mantle the resulting concentration of toroidal flow in the upper mantle may result in more thorough mixing there and account for some of the geochemical and isotopic differences proposed to exist between the upper and lower mantles

a thin viscous sheet model for continental deformation summary for the purposes of describing its large scale and long term deformation the continental lithosphere is regarded as a continuum obeying a newtonian or a power law rheology the flow of a thin sheet of power law material overlying an inviscid substrate is studied under the assumption that vertical gradients of the horizontal velocity are negligible a numerical model is used to investigate the deformation of such a sheet under conditions approximating those of continent continent collision the material flows in response to forces applied to its boundaries for example the indenting of one continent by another and to forces in its interior arising from gradients in crustal thickness the horizontal divergence of the flow produces changes in the crustal thickness and hence a time dependent form to the flow itself for a given set of boundary conditions the flow depends on the stress exponent in the power law rheology n and on the argand number which is a measure of the ratio between the stress arising from crustal thickness contrasts and the stress required to deform the material at the ambient strain rates when the effective viscosity of the medium is very high ar o crustal thickness variations do not influence the flow if the material is newtonian n l the deformation associated with an influx of material approximating an indenter is of much greater lateral dimension than the width of the indenter whereas when material has a power law rheology n 3 5 are used the deformation is confined to a region of lateral extent comparable to that of the indenter as ar increases the forces arising from crustal thickness contrasts exert more influence on the flow and the maximum crustal thickness that can be sustained by a given influx of material is related by a simple expression to the effective viscosity of the medium at the ambient strain rates in the limit of a very weak medium ar 10 the lithosphere is unable to sustain appreciable crustal elevation contrasts the results of these numerical experiments show that systems in which the effective viscosities are such that the maximum deviatoric stresses are between 1 kbar

geological evidence bearing on the miocene to recent structural evolution of the new hebrides arc abstract the new hebrides archipelago is a complex reversed arc system that can be divided into four major volcanic provinces the western belt is an early to middle miocene extinct volcanic arc that as a result of polarity reversal is now incorporated into the frontal arc of the present day configuration the eastern belt initially received detritus in the early middle miocene from a tholeiitic arc complex but in the mio pliocene became the locus of a more calc alkaline arc volcanism volcanic activity then ceased in the eastern belt but is well represented as a third and largely submerged marginal province through the pliocene into the early pleistocene the present volcanic line the central chain is essentially a continuation of the marginal province volcanism into recent times initial tectonic events in the new hebrides arc were associated with the regional disruption in the middle miocene of an east facing system with consequent termination of western belt arc volcanism the western belt remained as a landmass during the lowermost late miocene but subsided following a late miocene renewal of island arc volcanism to the east this latest phase was coeval with initial expansion of the north fiji basin and marked the advent of the new hebrides as a westwards migrating reversed arc system during arc migration there were apparent hiatuses in island arc volcanism the most notable being a middle pliocene to late pleistocene period of quiescence in the central sector tectonism in the early pleistocene recent raised the fore arc brought about rifting and extension to the rear and concentrated volcanism along the presently active central chain

geothermal exploration of paleozoic formations in central alberta this study explores the distribution of paleozoic formations in the central alberta basin and investigates rock properties with regard to their usability as geothermal reservoirs the study area of this regional scale investigation is about 160 km 200 km in size and located around edmonton where the basin depth ranges between 1 8 and 3 5 km a three dimensional 3d geological model was developed based on stratigraphic tops from about 7000 wells from the database of the alberta geological survey ags spatial distribution and thickness of deep formations were established in the 3d geological model porosity and permeability of four devonian carbonate formations cooking lake leduc and nisku formations and wabamun group were investigated using data from more than 50 000 core analyses average porosity of the devonian strata in the study area ranges from 4 5 nisku to 8 7 wabamun average permeability is between 3 5 10 15 m2 wabamun and 26 10 15 m2 leduc the distribution of both param

a laboratory investigation of effects of trench migration on the descent of subducted slabs a laboratory investigation of viscous slabs subducted from a migrating trench reveals a range of possible behaviours the slab dip in a uniform mantle is found to be steady or oscillatory depending on the rates of descent and trench migration in addition density and viscosity interfaces are used to model the increased resistance to sinking of stiff slabs through the seismic discontinuity at a depth of about 660 km in earth s mantle if the slab is denser than the lower layer and its dip in the upper layer is steady it can continue to descend through the lower layer in a tabular form and in either a steady or oscillatory manner or it can be laid horizontally on the interface for a distance before descending in a chain of diapirs resulting from gravitational instability at the interface if the slab dip in the upper layer is unstable the slab sinks into the lower layer in a chain of large diapirs at a spacing determined by the frequency of oscillations set by instability of the slab within the upper layer the style of penetration depends on the trench migration speed and the ratio of sinking velocities in the two layers estimates for mantle slabs indicate that they may range across the major regime transitions the experimental system provides a gross simplification of mantle conditions but the results are a testimony to the possibility of a range of complex behaviour of subducted lithosphere they also indicate that the relationship between the structural features in the lower mantle revealed by seismic imaging and present day tectonic processes at the surface may not be obvious tomographic images of western pacific subduction zones and data for the migration rates of associated trenches suggest a dependence of slab behaviour on migration rate similar to that seen in the experiments

distribution and origin of igneous rocks from the landward slopes of the mariana trench implications for its structure and evolution the landward slope of the mariana trench is composed largely of igneous rocks serpentinites and serpentinized ultramafic rocks occur at nearly all structural levels on the slope from depths of 8000 to 1200 m seamountlike features on the trench slope break are the surface expression of serpentinite diapirs cumulate and massive gabbros are found several varieties of volcanic rocks are common including boninites altered and metamorphosed basalts andesites and dacites the chemical characteristics of the volcanic rocks indicate that nearly all are products of island arc volcanism together with the gabbros these volcanic rocks represent what is probably a late eocene arc complex these rocks were probably the first volcanic products to result from the subduction of the pacific plate beneath the phillippine sea plate their exposure on the trench slope today implies a significant amount of tectonic erosion of the landward slope since eocene time most of this removal of material appears to have occurred during the early stages of subduction there are isolated occurrences on the landward slope of rock assemblages including alkalic basalts chert hyaloclastites upper cretaceous siliceous sediments and shallow water limestones these assemblages are very similar to rocks dredged from seamounts on the offshore flankmore of the trench and their presence on the landward slope suggests that since the cessation of vigorous tectonic erosion there has been episodic accretion of seamount fragments to the landward slope less

controls of subduction geometry location of magmatic arcs and tectonics of arc and back arc regions most variation in geometry and angle of inclination of subducted oceanic lithosphere is caused by four interdependent factors combinations of 1 rapid absolute upper plate motion toward the trench and active overriding of the subducted plate 2 rapid relative plate convergence and 3 subduction of intraplate island seamount chains aseismic ridges and oceanic plateaus anomalously low density oceanic lithosphere cause low angle subduction under conditions of low angle subduction the upper surface of the subducted plate is in contact with the base of the overlying plate the wedge of low density asthenosphere is replaced by subducted lithosphere and the width of the arc trench gap either is significantly increased or a magmatic arc is not developed within the overlying plate the fourth factor is age of the subducting lithosphere subduction of young lithosphere produces two opposing tendencies 1 low angle subduction and increased arc trench distance owing to its low density and 2 decreased arc trench distance owing to its higher temperature two factors of secondary importance contribute to variation in subduction zone geometry and arc trench distance accretion of sediment in trenches depresses the upper portion of the subducting oceanic plate and causes the trench axis to migrate seaward prolonged subduction thickens the upper plate depresses the isotherms in the subducted plate and may create a broader arc both factors increase the arc trench gap the four primary factors also control development of other tectonic elements such as regional subsidence for example the amazon basin and a portion of the cretaceous interior seaway of western united states intra arc extension for example the basin and range province foreland fold and thrust belts and laramide style tectonics

pliocene variations in the position of the antarctic convergence in the southwest atlantic middle to late pliocene 4 1 1 9 ma variations in the positions of surface water masses and migrations of the polar front in the southwest atlantic are inferred from a factor analysis of radiolarian assemblages from dsdp site 514 faunal results suggest that surface water masses underwent progressive cooling as the polar front zone pfz advanced northward during the latest gilbert to late gauss chron the late gilbert chron 4 12 3 86 ma was the waning phase of a mid gilbert warm interval northern subantarctic surface water prevailed between 4 12 and 3 97 ma indicative of surface water temperatures warmer than those during the succeeding pliocene a long term sea surface cooling trend began between 3 97 and 3 86 ma which moved the pfz northward and closer to site 514 bringing cooler southern subantarctic surface waters to the site this movement northward of cooler surface waters is inferred to be related to the initial growth of sea ice and ice shelves throughout regions of west antarctica which were deglaciated or unglaciated during the prior warming interval of the gilbert chron a major late gilbert early gauss chron increase in the production of antarctic bottom water aabw resulted in an erosional or nondepositional hiatus at site 514 spanning 3 86 to 3 18 ma increased aabw activity at this time is linked to continued growth of expansive ice shelves and sea ice in the west antarctic region in the east antarctic margin and possibly in the marine embayments of east antarctica decreased aabw production led to a 3 2 ma resumption of sedimentation at site 514 and elsewhere in the southern ocean antarctic surface water assemblages above the hiatus 3 18 3 16 ma suggest that the pfz first passed over site 514 during the 3 86 3 18 ma interval of the hiatus following this initial invasion of antarctic surface waters between 3 16 and 2 94 ma the pfz retreated to a position near or south of the site this slight retreat of the pfz with a corresponding decrease in aabw activity is possibly related to a significant reduction in the area of shelf and sea ice as grounded ice expanded in the west antarctic region we interpret the widespread global enrichment of 18o at approximately 3 2 ma to be partially related to a cooling of surface waters the production of cooler bottom waters and an expansion of west antarctic ice volume associated with its grounding phase between 3 5 and 2 7 ma based on surface water mass migrations at site 514 the west antarctic ice volume increased significantly slightly prior to or at approximately 3 2 ma decreased between 3 14 and 2 7 ma and then increased by 2 67 ma to a more stable configuration with a volume at least as great as that at 3 2 ma during the late gauss 2 94 2 47 ma four northward advances of the pfz led to the migration of antarctic surface waters over the site by 2 67 ma and continued dominance of antarctic surface water through the remainder of the pliocene we suggest that this permanent change from subantarctic to antarctic surface water mass dominance was linked indirectly to the initial growth of northern hemisphere ice and a reduction of sea level we propose a model of simultaneous northern and southern hemisphere ice volume growth induced by interhemisphere positive feedback which is responsible for the major late gauss advance of antarctic surface waters and the global enrichment of 18o at 2 5 ma

the development of slabs in the upper mantle insights from numerical and laboratory experiments we have performed numerical and laboratory experiments to model subduction of oceanic lithosphere in the upper mantle from its beginnings as a gravitational instability to the fully developed slab a two dimensional finite element code is applied to model newtonian creep in the numerical experiments scaled analog media are used in the laboratory a sand mixture models the brittle crust silicone putty simulates creep in the lower crust and mantle lithosphere and glucose syrup is the asthenosphere analog both model approaches show similar results and reproduce first order observations of the subduction process in nature based on density and viscosity heterogeneities in a stokes flow model subduction nucleates slowly and a pronounced slab forms only when the viscosity contrast between oceanic plate and mantle is below a threshold we find that the subduction velocity and angle are time dependent and increase roughly exponentially over tens of millions of years before the slab reaches the 670 km discontinuity the style of subduction is controlled by the prescribed velocity of convergence the density contrast between the plates and the viscosity contrast between the oceanic plate and the mantle these factors can be combined in the buoyancy number f which expresses the ratio between driving slab pull and resisting viscous dissipation in the oceanic plate variations in f control the stress in the plates the speed and the dip of subduction and the rate of trench retreat reproducing the contrasting styles of subduction observed in nature the subduction rate is strongly influenced by the work of bending the lithosphere as it subducts

factors affecting seismic moment release rates in subduction zones the amount of shallow seismic activity in subduction zones varies greatly from region to region we quantify this seismicity by calculating seismic moment release rates and seismic slip rates for 24 subduction zones to calculate the moment release we sum the seismic moment for all interplate thrust type events with surface wave magnitude ms 7 0 occurring from 1904 to 1980 we present a time history plot of the seismic moment release for each subduction zone these exhibit the differences in the seismic release patterns for subduction zones where the time window of our study is not representative the total moment release is corrected using information on repeat times the moment release rates are compared with various subduction parameters in order to determine which factors influence the degree of coupling these parameters include the age of the subducting lithosphere absolute velocities of the upper and subducting plates convergence velocity and length maximum depth and dip of the wadati benioff zone the moment release rate decreases as the age of the subducting lithosphere increases when the zones belonging to a single subducting plate are considered this age versus moment release relation is consistent for the zones in which the pacific cocos philippine sea and indian plates are subducting the moment release rates for the subduction zones in which the pacific plate is subducting are much higher than for zones of other plates with similar age the age versus moment release relation holds among the zones which belong to one subducting plate however zones with similar ages but belonging to different plates do not have the same moment release rates this suggests that within a single plate the age is the dominating factor affecting the strength of seismic coupling but that each plate as a whole has a characteristic moment release budget zones with retreating upper plates tend to have lower moment release rates the moment release rate does not increase with convergence velocity no simple relationship was found between these two parameters the moment release rate depends most clearly on the age of the subducting lithosphere and the absolute velocity of the upper plate these are the two independent subduction zone parameters among the variables considered the other variables depend on these two parameters

near surface thermal chemical boundary layer convection at infinite prandtl number two dimensional numerical experiments chemical differentiation and convective removal of internal heat make the earth s lithosphere a thermal and a chemical boundary layer thin layers of chemically light material form near the earth s surface and become embedded within the cold thermal boundary layer associated with interior heat removal the likelihood of near surface thermal and chemical boundary layer interactions influencing the earth s thermotectonic evolution prompts the models presented herein a simplified system consisting of a chemically light layer within the upper thermal boundary layer of a denser thermally convecting layer is explored through a suite of numerical experiments to see how its dynamic behaviour differs from similar well studied thermal boundary layer systems a major cause of differences between the two systems resides in the ability of the deformable near surface chemical layer to alter the effective upper thermal boundary condition imposed on the convectively unstable layer below in thermal equilibrium regions of chemical boundary layer accumulation locally enforce an effectively near constant heat flux condition on the thermally convecting layer due to the finite thermal conductivity of chemical boundary layer material for cases in which chemical accumulations translate laterally above the unstable layer the thermal coupling condition between chemical boundary layer material and the unstable layer below is one of non equilibrium type i e the thermal condition at the top of the convectively unstable layer is time as well as space variable a second major cause of differences is that for the thermal chemical system chemically induced rheologic variations can offset or compete with those due to temperature more specifically the presence of chemically weak material can lubricate convective downwellings allowing for enhanced overturn of an on average strong upper thermal boundary layer both of these factors have low order effects on internal flow structure and heat loss and lead to dynamic behaviour in which chemical boundary layer deformation is not only driven by flow in the thermally convecting interior layer but also feeds back and alters this flow some implications of this in regard to elucidating how near surface chemical boundary layer deformation e g continental tectonics might interact with and influence mantle convection are discussed

magnetic lineations in marginal basins of the western pacific this paper separates the small oceanic basins around the western perimeter of the pacific ocean into marginal basins that have formed through back arc extension and those that apparently have not and reviews our knowledge of magnetic lineation patterns observed in possible and probable back arc basins magnetic lineations in these basins resemble lineations commonly associated with the world s mid oceanic spreading systems indicating that similar processes of crustal accretion occur in both tectonic environments in some back arc basins of the southwestern pacific magnetic lineation and other evidence suggest that back arc basins can evolve through the interaction and growth of multi plate systems because of the small time and space characteristics of back arc basins compared with the world s major oceanic spreading systems tectonic conditions favourable for the generation of back arc basins are either relaxed rapidly or easily interrupted models proposed to account for back arc basin development include a local models where back arc extension is mechanically driven by the downgoing slab and b global plate kinematics models where conditions favourable for back arc extension are governed by the motion of the overriding plate relative to the trench axis

radiolarian assemblage distributions and paleoceanography of the eastern equatorial pacific ocean during the last 127 000 years abstract mathematical treatment q mode factor analysis of radiolarian species abundance data in equatorial pacific core tops yields four oceanographically meaningful assemblages the distribution of these assemblages in surface sediments can be related to the major surface and subsurface characteristics of the equatorial circulation system western tropical mixed layer flow coastal upwelling and the oxygen minimum layer eastern equatorial divergence upwelling and the equatorial undercurrent and the subtropical water masses of both hemispheres changes in the distribution of radiolarian assemblages during the glacial interglacial cycle spanning the last 127 000 yr are observed in two equatorial pacific piston cores these changes reflect fluctuations in the boundaries of equatorial water masses in response to variations in the position and strength of the tradewind system and tend to precede major glacial interglacial transitions by several thousand years a history of atmospheric circulation in the equatorial pacific has been developed for the last 127 000 yr fluctuations in tradewind intensity north south shifts in the latitudinal position of the tradewind maxima and perhaps a change in the zonality of the tradewinds during the last 127 000 yr have been inferred from downcore quartz abundance and radiolarian assemblage distribution patterns oxygen isotopic fluctuations during the last 127 000 yr indicate that continental ice volume was much greater during glacial stage 2 29 000 b p to 11 000 b p than during glacial stage 4 73 000 b p to 61 000 b p large changes in atmospheric and oceanic circulation patterns also are indicated by the distribution of radiolarian assemblages during the two glacial periods during stage 2 oceanic circulation was more intense the southeast tradewind maximum was closer to the equator and wind stress was more meridional than during stage 4

europa s ice tectonics new insights from physical wax experiments with implications for subduction initiation and global resurfacing processes abstract jupiter s moon europa has one of the youngest geological surfaces in our solar system with an age of 40 90 ma implying an intense history of resurfacing the surface of europa indeed shows abundant evidence of tectonic deformation related to extension strike slip and shortening however observed features related to shortening are scarce compared with pervasive extensive extensional features such as dilational bands and do not suffice as the sole mechanism for recycling aging terranes recently evidence for potential plate tectonics associated with subduction zones has been discovered on europa this could be responsible for recycling most of europa s surface however basic physical parameters needed to initiate subduction on europa such as thickness of the brittle layer deformation rates and orientation of pre existing zones of weakness at which subduction could start are not well understood here we aim to better understand the process and the conditions that could lead to initiation of subduction on europa through physical experiments using wax to simulate europa s two layered i e convective icy crust by deforming the wax strain on europa s surface possibly caused by diurnal tides or its nonsynchronous rotation is simulated our results indicate that subduction could initiate over a broad range of surface thicknesses and deformation rates above a minimum conductive layer thickness but is strongly dependent on the orientation of the pre existing zones of weakness very thin conductive layer experiments however result in a previously undescribed process that we term ductile rolldown which creates surface features similar to double ridges observed on europa thus subduction and ductile rolldown represent physically plausible mechanisms that could play a critical role in resurfacing europa throughout its geologic history these results could yield significant implications for europa s thermal history and evolution habitability and future spacecraft missions

change from atlantic type to pacific type carbonate stratigraphy in the middle pliocene equatorial pacific ocean a paleomagnetically dated piston core from the equatorial pacific ocean v28 179 contains levels of highest carbonate percentage in association with lower 18o values in early pliocene sediments this relationship is typical of the atlantic during the quaternary and is thus referred to as an atlantic type carbonate stratigraphy the change from an atlantic type to a pacific type stratigraphy high carbonate glacial episodes appears to have occurred during the earliest gauss magnetic chron 3 18 3 41 ma of middle pliocene age and may be related to the development of northern hemisphere glaciation samples deposited prior to the development of northern hemisphere glaciation in v28 179 have a significant proportion of low frequency 400 ka period variance in both the benthic o18 and carbonate spectra and an atlantic type correlation both the oxygen isotopic and the carbonate spectra of v28 179 have frequency specific correlation coherency with the calculated pliocene record of the eccentricity of earth s orbit thus orbital parameters may have affected equatorial pacific sedimentation in the pliocene the change from atlantic type to pacific type carbonate stratigraphy in the middle pliocene may have been related to a change in the dominant frequency of variance of carbonate oscillations related to orbital parameters and or to intensification of antarctic bottom water circulation during the middle pliocene near the base of the gauss chron

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

the curvature of wadati benioff zones and the torsional rigidity of subducting plates it is an axiom of plate tectonics that lithospheric plates are to a good approximation torsionally rigid that is inextensible in the plane of the plate1 3 although few scientists would question this assumption when it is applied to oceanic lithosphere at the earth s surface there is no consensus as to whether or not oceanic lithosphere maintains its mechanical integrity during subduction com pare refs 4 7 with refs 8 11 i show here that the shapes of many wadati bemoff zones are inconsistent with the notion that they manifest subducted portions of oceanic plates subject only to very small 1 membrane strains

the late quaternary glaciations as the response of a three component feedback system to earth orbital forcing abstract a climatic feedback system previously described consisting of three prognostic nonlinear equations governing the mass of ice sheets the mass of marine and continental marginal ice and the mean ocean temperature is forced by a representation of the effects of external earth orbital variations with reasonable amplitudes for the eccentricity obliquity and precession forcing the free oscillatory solutions of major period near 100 kyr can be modified in a way that substantially agrees with the 18o derived observations of ice mass evolution in particular a proper structure variance spectrum and phase lock of the major variations are obtained over the last 400 kyr an analysis of the sensitivity of these results to variations in the model parameters and to random perturbations shows that the solution is robust for small changes in all but a few of the equation coefficients concomitant variability in the marine ice mass ocean temperature and net radiation at the top of the atmospher

pliocene pleistocene evolution of eastern tropical pacific surface water circulation and thermocline depth we evaluated the response of eastern equatorial pacific upper water column hydrography to the panama seaway closure and the initiation of large scale northern hemisphere glaciation the 18o gradient between globorotalia tumida and globigerinoides sacculifer indicates a general shoaling of the thermocline between 4 2 and 3 0 ma which may be related to a recorded decrease in southeast trade wind strength the 13c gradients suggest that surface water nutrient concentrations have increased during the last 5 myr and that upwelling intermediate waters began to change sources at 3 2 ma changes in the trans pacific east west temperature and nutrient gradients at 4 0 and 1 5 ma are coincident with the seaway closure and a major phase of northern hemisphere glaciation changes in the sensitivity of surface circulation to milankovitch forcing 18o records linearly related to orbital variations only after 3 2 ma are associated with the seaway closure pleistocene initiation of 100 and 41 kyr cycles in the thermocline depth proxy may suggest that the hydrographic response to milankovitch forcing is enhanced by air sea interactions which maintain the relatively steep pleistocene trans pacific gradients the 13c records are dominated by the 41 kyr period and are usually coherent with high latitude climate suggesting that associated changes in the global carbon reservoir dominate the tropical 13c signals

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

evaluating the success of astronomical tuning pitfalls of using coherence as a criterion for assessing pre pleistocene timescales the imprint of orbital variations on the geological record of climatic variability is well documented especially for the plio pleistocene there is considerable interest in developing very high resolution timescales through the cenozoic and into the mesozoic by tuning geological records to the assumed astronomical forcing since the precession signal is so highly amplitude modulated it is widely believed that high coherence between record and assumed forcing in the precession band is an indication that the timescale is probably correct because coherence is supposed to provide a measure of the degree of common amplitude modulation we show that this is misleading even a sinusoidal variation that has been tuned to an insolation record shows highly significant coherence at the 23 kyr and 19 kyr precession frequencies coherence is a good indication that the tuning has generated a consistent phase relationship but complex demodulation is a better tool for assessing the relationship between amplitude modulation in the data and in the hypothetical forcing

a simple global model of plate dynamics and mantle convection cooling and thickening of lithospheric plates with age and subduction result in large scale horizontal density contrasts tending to drive plate motions and mantle flow we quantify the driving forces associated with these density contrasts to determine if they can drive the observed plate motions first two dimensional models are computed to evaluate the effects of assumed rheologies and boundary conditions we are unable to obtain platelike behavior in viscous models with traction free boundary conditions the piecewise uniform velocities distinctive of plate motion can be imposed as boundary conditions and the dynamic consistency of the models evaluated by determining if the net force on each vanishes if the lithosphere has a newtonian viscous rheology the net force on any plate is a strong function of the effective grid spacing used leading to ambiguities in interpretation incorporating a rigid plastic lithosphere which fails at a critical yield stress into the otherwise viscous model removes these ambiguities the model is extended to the actual three dimensional spherical plate geometry the observed velocities of rigid plastic plates are matched to the solution of the viscous stokes equation at the lithosphere asthenosphere boundary body forces from the seismically observed slabs from the thickening of the lithosphere obtained from the actual lithospheric ages and from the differences in structure between continents and oceans are included interior density contrasts such as those resulting from upwellings from a hot bottom boundary layer are assumed to occur on a scale small compared to plate dimensions and are not included the driving forces from the density contrasts within the plates are calculated and compared to resisting forces resulting from viscous drag computed from the three dimensional global return flow and resistance to deformation at converging boundaries the rms residual torque is 30 of the driving torque the density contrasts within the plates themselves can reasonably account for plate motions body forces from convection in the interior may provide only a small net force on the plates at converging boundaries the lithosphere has a yield stress of 100 bars drag at the base of the plates is 5 bars and resists plate motion the net driving forces from subducting slabs and collisional resistance are localized and approximately balance driving forces from lithospheric thickening are distributed over the areas of the plates as is viscous drag the approximate balance of these two forces predicts plate velocities uncorrelated with plate area as observed the model represents a specific case of boundary layer convection the dynamical results are consistent with either upper mantle or mantle wide convection

quantifying nonlinearity and geometry in time series of climate abstract statistical methods provide the opportunity to quantify what are typically qualitative descriptions of climate change a review of linear and nonlinear processes as applied to time series of climate change is presented methods to quantify nonlinear climatic events are given using the late pleistocene record of glacial interglacial oscillations as an example methods to discriminate climate cycles driven by linear processes from nonlinear processes and to determine the strength of nonlinear interactions are described the adoption of a quantitative approach to descriptions of climate variability allows direct comparison of geologic and paleoclimatic data to results of models that seek to describe modes of climate oscillations

decoupling during continental collision and intra plate deformation abstract fundamental aspects of mountain building such as the dominant deformation mechanism strain localisation and topography development in response to the presence of decoupling zones either along boundaries of weak and strong lithospheres or within weak domains of continental lithosphere are studied by analogue lithospheric models scenarios are investigated where a weak part of the lithosphere is bordered on either side by stronger lithospheres geometries of the weak zone vary from box to wedge to graben shaped representing lateral strength variations as a consequence of earlier deformation phases the rheologic and geometric configuration at the onset of shortening is therefore considered as representative for continental collision and intra plate settings the modelling results show that the geometric and topographic evolution of mountain belts through time is sensitive to the presence of decoupling horizons along plate or tectonic boundaries this is because the decoupled boundaries localise strain control the dominant deformation mechanism and hence the locus of uplift and subsidence accordingly decoupled boundaries favour deformation by thrusting and particularly in the case of an initially wedge shaped weak zone the development of mountain belts with overall triangular shape in contrast lithospheres with coupled boundaries dominantly deform by buckling irrespective of the initial geometry of the weak zone as this deformation mechanism affects all model domains it produces a distinguishable pattern of uplift and subsidence controlled by buckling strong decoupling at mid crustal levels within the weak zone favours vertical strain partitioning and the formation of plateau shaped mountains the modelling results suggest that natural collisional systems should respond to an increase in plate coupling through time by widespread uplift of the mountain belt including its foreland basins which were formerly created under decoupled conditions due to the absence of decoupling zones like plate contacts weak domains in intra plate settings are predicted to deform dominantly by buckling

equatorial pacific sea surface temperatures faunal patterns and carbonate burial during the pliocene abstract a high resolution record of radiolarian faunal abundances from the eastern equatorial pacific is compared to records of carbonate and noncarbonate burial to examine the evolution of eastern tropical pacific climate processes during the pliocene these data provide a means to evaluate the sensitivity of the equatorial pacific to the onset of northern hemisphere glaciation around 2 8 2 5 ma to the closure of the isthmus of panama around 4 4 3 2 ma and to orogeny related weathering changes before 4 0 ma radiolarian faunal assemblages and sea surface temperature sst estimates indicate a gradual cooling from early to late pliocene but no significant changes occur near the onset of northern hemisphere glaciation records of carbonate and noncarbonate mass accumulation show a long term decrease from the miocene pliocene boundary to the upper pliocene greater carbonate burial in the early pliocene relative to the late pliocene parallels a gradual cooling from early to late pliocene and may reflect changes related to isthmus closure or widespread orogeny no significant time domain changes are seen in the eastern equatorial pacific that could be related to the onset of northern hemisphere glaciation evolutive spectral analyses of these equatorial pacific climate parameters indicate that variance in sst and seasonality commonly concentrate at frequencies not linearly related to orbital variations furthermore cross spectral comparisons with a high resolution benthic 18 o record indicate that the surface ocean and carbonate flux share little coherent variance with high latitude climate processes during the pliocene given the high degree of chronostratigraphic control in these records these results suggest that milankovitch band surface ocean processes as well as carbonate burial in the equatorial pacific are decoupled from high latitude climate processes during the pliocene

covariance patterns of foraminiferal dgr18o an evaluation of pliocene ice volume changes near 3 2 million years ago oxygen isotope data for a pliocene interval from 3 6 to 2 8 million years ago show a mean increase 0 5 per mil of benthic delta 18 o at about 3 2 million years ago whereas planktic delta 18 o does not increase this lack of covariance indicates that the event at 3 2 million years did not result in a permanent increase in the ice budget of either the northern or the southern hemisphere

a continuum model of continental deformation above subduction zones application to the andes and the aegean large scale continental tectonics of back arc extensional and andean type compressional environments are investigated by using the thin viscous sheet model to calculate the deformation within continental lithosphere that is subjected to horizontal forces on its plate boundaries and to basal drag from the asthenospheric flow beneath the shear tractions acting on the base of a deformable lithosphere are derived from a corner flow model that assumes a rigid subducting plate and a deformable overlying plate because the calculated shear tractions and the deformation within the overlying plate are interdependent the corner flow and the thin viscous sheet models are solved simultaneously we use a perturbation method to obtain analytical solutions for the velocity and strain rate fields within the overlying continental lithosphere the solutions depend on the angle of subduction the dimensionless thickness of the lithosphere and the ratio of asthenospheric to lithospheric viscosities which governs the viscous coupling between the asthenosphere and the lithosphere calculations are compared with observations from the andes and the aegean our results explain some of the features of the deformation in these regions that have heretofore not been explained by other models our model predicts that in a compressional environment a broad region of uplifted topography will tend to develop above a more steeply dipping slab 30o rather than above a shallower slab 10o 15o this is in accord with observations in the various segments of the central andes for an extensional environment the model predicts that a zone of compression can develop near the trench and that extensional strain rate can increase with distance from the trench as is observed in the aegean we also estimate the effective viscosities of 102o pa s for the asthenosphere 2x 1021 pa s for the aegean lithosphere and 1022 pa s for the andean lithosphere

role of preexisting topography and overburden on strain partitioning of oblique doubly vergent convergent wedges 1 scaled sandbox modeling was used to analyze the three dimensional aspects of strain partitioning along obliquely convergent margins that underwent a subduction polarity reversal event in the experiments the first phase of shortening produced a preexisting topography that affected wedge kinematics during the second phase increasing degrees of obliquity during the first shortening phase produced different amounts of inherited topography and overburden giving way to different effects on wedge response during the second phase of shortening during the second phase the models with an inherited heavy overburden show an anomalous orientation for riedel shears predating the development of a margin parallel strike slip fault these shears are oriented at an angle that is higher than kinematically predicted with respect to the subduction slot suggesting a transient rotation of the stress applied by the obliquity of subduction in contrast experiments with a lighter overburden do not show this inconsistency in addition all the models displayed an unexpectedly high accretion rate during the second phase of shortening which we interpret to be dependent on topographical slope breaks inherited after the first phase of shortening and the critical asymmetrical architecture of doubly vergent wedges

deformation of a weak subducted slab and variation of seismicity with depth oceanic lithosphere is assumed to possess platelike properties whether it is lying at the earth s surface or descending deep into the mantle yet the geometries of subducting slabs indicate significant deformation as they descend through the mantle1 4 suggesting that lithospheric plates might be so weakened during subduction as to act not as a rigid solid but as a viscous fluid numerical and laboratory experiments have shown that fluid slabs can indeed take on realistic profiles5 7 however it has not been clear whether a fluid or weak model of slab dynamics can account for deep earthquakes which are usually ascribed to the deformation of a rigid or strong slab here we show using numerical simulations of slab evolution that a weak slab model is in fact consistent with seismic observations assuming that earthquakes occur at a rate proportional to deformation rate we reproduce the observed variation with depth of seismicity rate and focal mechanisms and the cessation of seismicity at 670 km depth provided that sinking material encounters resistance at depth here modelled as a viscosity jump at 670 km the pattern of seismicity can be explained by any mechanism for deep earthquakes in which the rate of seismicity is proportional to strain rate in the slab

seismic constraints on the morphology of deep slabs residual sphere images from deep earthquakes not only detect the presence of slab associated velocity anomalies but also lend insight into the flow and deformation of lithosphere subducted into the lower mantle we have compared travel times from deep events in the kuril and mariana arcs with the seismic velocity anomalies implied by kinematical models that thicken the slab perpendicular to its plane by reducing the vertical velocity of the flow with depth we assume that the details of the deformation whether the slab buckles imbricates fragments etc are averaged out along the ray paths and hence our models constrain the scale not the mode of slab thickening the deep event travel times are best fit by undeformed models but the ability of the residual sphere method to resolve slab thickness is limited by ray bending effects although the mariana times are consistent with advective thickening factors of 5 or more factors larger than 3 are ruled out by the kuril data for all models examined the data require that slab material extends to depths of 900 1000 km global tomographic models and regional studies which delineate high velocity anomalies in the lower mantle beneath zones of cenozoic subduction are consistent with our results as is recent work on pulse distortion by slab gradients comparison of observed and predicted rates of seismic moment release suggests that if substantial advective thickening does occur it is largely aseismic

seismic evidence for continuity of the deep slab beneath central and eastern peru seismograms of peru brazil deep focus earthquakes recorded on the nearby carnegie broadband station at cuzco peru cus are characterized by an anomalous large amplitude p wave arrival that follows about 1 5 s after direct p the anomalous arrival is interpreted to be an underside wide angle reflection from the upper surface of the descending nazca plate somewhere in the depth range 150 400 km and in a region where the slab is wholly aseismic particle motion analysis of the cus seismograms shows that the direct p wave has arrival angles as predicted from simple ray tracing the large impulsive anomalous phase termed pr also arrives as a p wave but is of opposite polarity and about 1 5 times the amplitude of the direct p arrival the back azimuth of pr is about 20 clockwise from the path of direct p and its angle of emergence is significantly larger forward modeling calculations indicate that the seismic observations are consistent with a model in which pr originates as a p wave that travels up inside the slab from the source and is then reflected at wide angle from the steeply dipping upper plate boundary results from modeling lead to the following conclusions 1 the descending nazca plate exists in and is probably continuous through the aseismic region between 150 and 500 km depth 2 the plate is very steeply dipping 70 below depths of about 150 or 200 km and 3 the velocity discontinuity at the reflection point must be relatively sharp probably less than about 10 km in thickness

an experimental study of subduction and slab migration results of fluid dynamical experiments on the behavior of subducted slabs are presented that address two important characteristics of subduction slab penetration through the transition zone and horizontal slab migration cold negatively buoyant molded slabs of concentrated sucrose solution with viscosities of 3 5 106 p were introduced into a more dilute two layered sucrose solution representing the upper and lower mantle the transition zone was modeled by a step increase in both density and viscosity ratios of slab thickness to upper layer depth the relative viscosities in each layer and the effective rayleigh number were close to mantle values the initial configuration consisted of two plates separated by a trench gap with one plate attached to a dipping slab simulating a developing subduction zone with an overriding plate two different boundary conditions were used at the trailing ridge end of the subducting plate 1 fixed end corresponding to zero ridge motion and 2 free end corresponding to zero ridge resistance based on results from 15 experiments with various combinations of upper layer lower layer and slab densities u l s and viscosities we find that slab penetration depth is sensitive to both the normalized slab density anomaly r s l s u and the dip angle three penetration regimes are recognized r 0 5 in which the slab sinks into the lower layer without distortion 0 2 r 0 5 which permits limited penetration the amount dependent on dip angle in this regime the slab develops a thick root just beneath the transition zone and r 0 2 in which the slab is deflected by the transition and little or no penetration occurs retrograde slab motion and trench migration occurred in nearly every case with the overriding plate moving parallel to but faster than the trench causing the gap between subducting and overriding plates to close with time prograde trench motion was never observed and the trench remained stationary only in the extreme case free ridge condition with equal upper and lower layer densities these results suggest that the transient extension found in many back arc basins may be a direct consequence of the tendency for retrograde subduction and the amount of extension may be governed in part by the degree of slab penetration through the transition zone

isotopic paleoceanography of the caribbean and east pacific role of panama uplift in late neogene time comparisons of carbon isotopic data on benthic foraminifera from deep sea drilling project sites 502 western caribbean and 503 eastern pacific indicate that the difference between the atlantic and the pacific in the per mil enrichment in carbon 13 of total dissolved carbon dioxide increased about 6 million years ago and again 3 million years ago when the difference reached the modern level 1 per mil comparisons of planktonic foraminiferal oxygen isotopic data for the caribbean and the pacific suggest that the salinity of caribbean surface waters began increasing 4 million years ago possibly in response to shoaling of the panama isthmus these results suggest that modern circulation patterns in the caribbean and eastern pacific developed by 3 million years ago in concert with changing tectonic climatic and biogeographic patterns

seabeam survey at the southern end of the manila trench transition between subduction and collision processes offshore mindoro island philippines abstract the morphological and structural study conducted at the southern tip of the manila trench reveals that convergence between the south china sea basin and luzon is accommodated differently depending on the nature of the subducted slab when the oceanic crust is subducted a simple accretionary prism fore arc basin pattern is developed conversely where the continental margin of this basin is subducted intraplate deformation is randomly distributed across the major part of the fore arc area which is fragmented into various crustal microblocks results of seabeam mapping and detailed geophysical surveying conducted at this subduction collision transition zone during the pop2 cruise with r v jean charcot are presented and discussed here and allow us a new insight into the mechanism of such a subduction collision transition zone

a critical assessment of viscous models of trench topography and corner flow stresses for newtonian viscous flow in a simple geometry e g corner flow bending flow are obtained in order to study the effect of imposed velocity boundary conditions stress for a delta function velocity boundary condition decays as 1 r 2 for a step function velocity stress goes as 1 r for a discontinuity in curvature the stress singularity is logarithmic for corner flow which has a discontinuity of velocity at a certain point the corresponding stress has a 1 r singularity however for a more realistic circular slab model the stress singularity becomes logarithmic thus the stress distribution is very sensitive to the boundary conditions and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries topography and seismicity data from northern hoshu japan were used to construct a finite element model with flow assumed tangent to the top of the grid for both newtonian and non newtonian flow power law 3 rheology normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies there is poor agreement purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity elasticity and plasticity appear to be important for the subduction process

buoyancy rather than rheology controls the thickness of the overriding mechanical lithosphere at subduction zones the thickness of earth s mechanical lithosphere is poorly defined to investigate whether rheology controls the thickness of the overriding plate s mechanical lithosphere in subduction zones the thermal structure was modelled numerically assuming a temperature dependent mantle viscosity it was found that the overriding lithosphere was ablated such that very high temperatures reached close to the surface near the apex of the wedge corner leading to unrealistically high heat flow since temperature dependent rheology clearly does not control the thickness of the mechanical lithosphere we suggest that it is instead controlled by buoyancy the source of buoyancy we assume is compositional e g buoyant crust two end member models with crustal thickness of 10 and 70 km respectively were then undertaken these had lower heat flow this work supports the assumption of some earlier workers e g plank and langmuir 1988 who equated the mechanical lithosphere with the crust of the overriding plate

radiolarian and silicoflagellate response to oceanographic changes associated with the 1983 el nino strong seasonal signals in fluxes and composition of siliceous microfossils have been recorded in sediment traps from two sites in the equatorial pacific deployed from december 1982 to march 1984 during the early part of the sampling period the 1982 83 el nino event had a profound effect on the radiolaria and silico flagellates within these two areas during the el nino the radio larian trap assemblages at site c 1 n 139 w most resembled faunal assemblages in western equatorial pacific sediments whereas the trap assemblages resembled equatorial divergence sediments in the latter half of the period at site s 11 n 140 w the radiolaria and silicoflagellate species can be correlated with organic carbon fluxes in general silicoflagellate shell fluxes are correlated to total opal fluxes where radiolarian fluxes do not exhibit this relationship approximately 20 25 of the total count of radiolaria in traps are not present in underlying sediments with a significant loss of the weakly silicified forms however the comparison of trap sample assemblages with underlying sediments reported here shows that dissolution does not alter relative abundances of the fossil species used in palaeoclimatic reconstructions a significant difference is observed between the silicoflagellate trap assemblage and the underlying sediment assemblage

polyphase formation and exhumation of high to ultrahigh pressure rocks in continental subduction zone numerical modeling and application to the sulu ultrahigh pressure terrane in eastern china 1 high to ultrahigh pressure hp uhp metamorphic rocks commonly form and exhume during the early continental collision and many questions related to their origin still remain unresolved we focus our study on explaining the poly metamorphic origins of many hp uhp terranes composed of tectonic units having strongly variable ages peak metamorphic conditions and p t paths these features are especially well characterized for the sulu uhp terrane in eastern china which we have chosen therefore as the reference case we conducted 2 d thermomechanical numerical modeling of continental subduction associated with formation and exhumation of hp uhp rocks our experiments suggest existence of several consequent episodes of u hp rocks exhumation related to the inherently cyclic origin of continental crust subduction detachment exhumation processes three major phases of these processes are identified in our reference model for the sulu uhp terrane 1 first and 2 second exhumation episodes of hp rocks originated in the subduction channel at lithospheric depths and 3 exhumation of uhp rocks originated at asthenospheric depths numerical models also suggest that subducted uhp rocks which are positively buoyant compared to the mantle may detach from the slab forming a flattened plume underplating the overriding lithosphere this sublithospheric plume may exist for several million years being heated to 800 900 c by the surrounding hot mantle at the later stage upward extrusion of hot partially molten rocks from the plume may exhume high temperature ht uhp complexes toward the surface

on the evolution of marginal basins this paper presents a critique of the models proposed to explain the formation of back arc basins those marginal basins located behind active or inactive trench systems and whose origin is inferred to be subduction related several marginal basins are not generically back arc but rather are related to the initiation of subduction continental rifting or plate boundary readjustments models of back arc basin formation must be consistent with the observed characteristics of back arc basins which include the following 1 the crustal structure and magnetic lineation patterns of back arc and major ocean basins are similar implying that the processes of crustal accretion in both tectonic environments are also similar 2 arc trench systems without currently spreading back arc basins are in the majority 3 although petrographically similar to and within the compositional range of midocean ridge basalts morb back arc basin basalts babb show consistent geochemical differences from n type morb and are in many respects transitional toward island arc tholeiites iat 4 distinct babb and iat mantle sources may be in extremely close proximity during the initial rifting of an island arc 5 rifting of the volcanic arc may occur on either side of the line of active volcanoes 50 km and may vary from one side to the other along strike 6 a remnant arc is not always developed 7 back arc spreading may be initiated above subducted lithosphere but with time the center of spreading migrates away from the arc and no longer overlies a benioff zone and 8 there is no simple correlation between the timing of global plate reorganizations and the formation of back arc basins the three main classes of models proposed to explain back arc basin formation are mantle diapirism induced aesthenospheric convection and global plate kinematics because the driving mechanism for the first two classes of models is provided by the local subducting slab these models fail to explain the temporal and spatial distribution of back arc basins the third class of models proposes that back arc basins should form whenever global plate interactions require divergence between the overriding plate and the trench line the latter being observed to move only in a seaward direction with respect to a hot spot reference frame however the tectonic setting of several back arc basins suggests that they represent more than just a passive response to kinematic boundary conditions for example some dynamic back arc force seems required to move the new hebrides tonga and south sandwich forearc plates seaward neither the local nor the global models adequately specify the necessary and or sufficient conditions for back arc basin formation a better understanding of the processes associated with the initiation of rifting and of subduction is required arc rifting and back arc spreading are both secondary phenomena to subduction but they are also primary extensional phenomena fundamentally analogous to continental rifting and the initial spreading of major ocean basins

nonlinear response in the global climate system evidence from benthic oxygen isotopic record in core rc13 110 the milankovitch theory of climate change predicts that variations of the climate system should match the dominant frequencies of the orbital forcing in the 41 and 23 kyr 1 frequency bands such a linear theory would predict that the amplitude variations of the climate response in these bands should match amplitude variations in orbital forcing here we compare amplitude variations of the marine oxygen isotope record with orbital forcing in these bands over the last 700 000 years and find systematic changes through time we express these amplitude mismatches as variations in the glacial response time a measure of the climate system s sensitivity to orbitally induced insolation changes variations in the glacial response time occur in all frequencies bands without strong concentration of variance in any given band and have a red spectrum with larger variations at the longer periods the response time is coherent with 18o at periods of 100 and 41 kyr which suggests that the variations in glacial response time in part reflect internal feedback mechanisms of the global climate system the phase relationship between the estimated glacial response time and the 18o ice volume record is very different at these two frequencies which suggests at least two separate feedback mechanisms the first mechanism enhances the 100 000 year climate cycle by increasing rates of change during major glacial terminations candidates for this feedback include lithospheric depression and rebound enhanced ice calving from large marine based ice sheets and possibly others a second set of mechanisms which is detected in the response to the 41 000 year orbital cycle of earth s obliquity accelerates ice growth events and slows glacial melting some models which include feedbacks between ice sheets sea ice and deep ocean temperatures predict early rapid ice growth followed by slower growth and this general feature is consistent with our analysis while we can not at present identify the specific feedbacks leading to asymmetry of growth and decay rates at different frequency bands the finding of this ice growth acceleration mechanism in the 41 000 year frequency band suggests that high latitude processes where insolation varies most strongly at this rhythm may be involved our finding of systematic changes in climate sensitivity has implications for orbitally tuned chronologies in pleistocene sediments instead of a constant phase shift within a frequency band between orbital forcing and glacial response as has been assumed in the past we suggest a variable phase the largest changes in age estimates for isotopic events are at the glacial terminations which in our chronology are as much as 3500 years older that estimated previously

from subduction to collision control of deep processes on the evolution of convergent plate boundary 1 using laboratory experiments we investigate the dynamics of the collisional process that follows the closure of an oceanic basin the evolution of these experiments systematically shows four successive episodes of deformation which correspond to 1 the initiation of oceanic subduction 2 a mature period of oceanic subduction 3 an episode of continental subduction during which the trench absorbs all the convergence and the superficial tectonic regime does not change significantly within the continental plates and 4 continental collision that starts when the trench locks and convergence is absorbed by a system of thrust faults and folds we observe that the amount of continental material that subducts before the onset of collision depends on the slab pull exerted by the subducted oceanic lithosphere the slab pull force in turn depends on the amount of subducted oceanic material on the thickness of the convective mantle and also on the rheology of the slab our experiments indeed suggest that parts of the oceanic slab may separate from the superficial slab to sink rapidly in the mantle decreasing the slab pull level and triggering the rapid onset of collision we observe two possible modes of slab deformation slab break off and development of viscous instabilities we define two dimensionless rheological numbers to characterize the possible occurrence of these modes of deformation in all cases oceanic closure is followed by episode of continental subduction during which the continental lithosphere may reach depths varying between 50 and 450 km prior to the onset of continental collision

a chronology of pliocene sea level fluctuations the u s middle atlantic coastal plain record abstract the sequence of marine and marginal marine silts and sands of the u s middle atlantic coastal plain preserves a relatively complete record of pliocene sea level highstands a general chronology for the deposition of these units has been developed previously from studies of the microfossil biostratigraphy complemented by a limited number of paleomagnetic and radiometric dates within this broad temporal framework it is currently difficult to assign ages to individual transgressive events with a resolution much better than 0 5 ma the age model proposed here attempts to refine the depositional history of the atlantic coastal plain by correlation with the more continuous deep ocean record the model assigns probable ages to the pliocene transgressions onto the atlantic coastal plain based on evidence for sea level highstands inferred from high resolution deep ocean 18 o records at least five major transgressive regressive events can be correlated with presumed eustatic fluctuations the temporal resolution for the timing of the transgressive events is greatly improved over that of the biostratigraphy alone many aspects of the proposed chronology require verification by more detailed field study but it is offered as a testable working hypothesis

evolution of crustal thickening in the central andes bolivia abstract paleoelevation histories from the central andes in bolivia have suggested that the geodynamic evolution of the region has been punctuated by periods of large scale lithospheric removal that drive rapid increases in elevation at the surface here we evaluate viable times and locations of material loss using a map view reconstruction of the bolivian orocline displacement field to forward model predicted crustal thicknesses two volumetric models are presented that test assumed pre deformation crustal thicknesses of 35 km and 40 km both models predict that modern crustal thicknesses were achieved first in the northern eastern cordillera ec by 30 20 ma but remained below modern in the southern ec until 10 ma the altiplano is predicted to have achieved modern crustal thickness after 10 ma but only with a pre deformation thickness of 50 km including 10 km of sediment at the final stage the models predict 8 25 regional excess crustal volume compared to modern thickness largely concentrated in the northern ec the excess predicted volume from 20 to 0 ma can be accounted for by 1 crustal flow to the wc and or peru 2 localized removal of the lower crust or 3 a combination of the two only models with initial crustal thicknesses 35 km predict excess volumes sufficient to account for potential crustal thickness deficits in peru and allow for lower crustal loss however both initial thickness models predict that modern crustal thicknesses were achieved over the same time periods that paleoelevation histories indicate the development of modern elevations localized removal of lower crust is only necessary in the northern ec where crustal thickness exceeds modern by 20 ma prior to paleoelevation estimates of modern elevations by 15 ma in the altiplano crustal thicknesses match modern values at 10 ma and can only exceed modern values by 5 ma post dating when modern elevations were thought to have been established collectively these models predict that the timing of crustal thickening is consistent with paleoelevation data without requiring large scale removal of lower crust and mantle lithosphere

deep structure and tectonics of the burmese arc constraints from earthquake and gravity data abstract active subduction of the indian plate is currently occurring beneath the burmese arc along an east dipping benioff zone which extends to a depth of about 180 km the overriding burma plate has an appearance of an inland seismic slab that is deflected downwards in the vicinity of the benioff zone a crustal seismic zone some 60 80 km east of the benioff zone correlates to backarc activity a triangular aseismic wedge in the top part of the crust outlines the central belt molasse basin east of the burmese foldbdt fault plane solutions show that the burmese benioff zone is characterized by shallow angle thrusting at its upper edge whereas down tip tensional events dominate its lower edge most of the backarc seismicity is accounted for by the sagaing transform or by the activity of the shan scarp normal fault zone at the margin of the asian plate a gravity anomaly pair with amplitude of 175 mgal coincides with the 1100 km long burmese arc lying in a north south direction the gravity anomalies along a profile in central burma and in adjacent areas of the bengal basin are interpreted in terms of plate subduction as well as near surface mass anomalies this suggests that sediments below the central belt may have an average thickness of the order of 10 km but may be as thick as 15 km at the subduction zone the oceanic crust underlying deeper parts of the bengal basin experiences phase transition at about 30 km depth in a benioff zone environment east of the burmese foldbelt several thrust planes are present within the folded and deformed cretaceous tertiary sediments of the fold belt these are often associated with ophiolites and basic to ultrabasic rocks a low density zone at least 60 km wide underlies the andesitic volcanic axis in the overriding plate

late miocene pliocene magnetic epoch 9 gilbert magnetic epoch calcium carbonate stratigraphy of the equatorial pacific ocean calcium carbonate accumulation patterns in the deep ocean are determined by the surface water by dissolution in under saturated deep water and by postdepositional addition or removal by currents arrhenius 1952 berger 1970 broecker 1971 in the eastern equatorial pacific downcore variations in pleistoncene calcium carbonate concentrations have been shown to be disolution cycles arrhenius 1952 hays and other 1969 saito and others 1975 which have been correlated to the glacial interglacial cycles of oxygen isotope stratigraphy emiliani 1955 adelseck and anderson 1978 extension of these dissolution cycles back to the miocene pliocene bounday has been attemped kaneps 1973b but lack of data from miocene samples has precluded further efforts of this type site descriptions the cenop project has undertaken the study of the late miocene ocean about 7 m y before the present at a time when the antarctic was glaciated based on faunal floral and oxygen isotopic evidence this was a time of relatively stable oceanographic conditions shackleton and kennett 1975 douglas and savin 1973 savin and others 1975 savin 1977 the cores used in this study table 1 were chosen because they contain long continuous sediment records that span this relatively stable late miocene interval with good preservation of calcareous and siliceous microfossils during this interval these core sites lay under the pacific high productivity belt associated with the equatorial divergence arrhenius 1952 ewing and others 1968 in which the carbonate compensation depth ccd is depressed due to the excess supply of calcium carbonate heath 1969 van andel and others 1975 because all three piston cores are under the same general surface and deep water condition changes in the factors that control carbonate sedimentation and preservation should be reflected in a similar way in their sediments however the deep sea drilling project dsdp and piston cores used in this study cover a wide range of water depths from a shallow drilled core dsdp site 158 to piston cores that currently lie below the ccd rc12 66 rc12 65 thus the carbonate preservation at different sites may vary due to differences in depth as well as corrosiveness of the bottom waters

sea surface temperature changes in the north pacific during the late miocene abstract twenty seven radiolarian species and species groups are identified and used in this study their stratigraphic ranges span from late miocene to modern times and they appear to have distributions which have changed little over this interval of time census data for these 27 taxa in surface sediment samples are used to define six assemblages q mode factors the distribution of these assemblages are similar to modern water mass distributions and can be statistically related to modern sea surface temperatures census data on these same 27 taxa have been collected in the upper miocene intervals of six pacific sites these faunal data are described in terms of the modern radiolarian assemblages and based on the multiple regression equations relating the modern assemblages to modern temperatures estimates of late miocene temperatures are made the estimated temperatures indicate that there is an overall cooling trend in the later part of the miocene magnetic chron 8 to chron 5 punctuated by several distinct cooling events in august the low latitude sites appear to be slightly cooler and the mid latitude sites warmer than in modern times in february both the low and mid latitude sites appear to be warmer than modern temperatures at the same locations the relatively warm temperatures near the base of the interval studied may be associated with the passage of warm tropical atlantic waters into the eastern pacific via the isthmus of panama the general cooling trend seen in these records is thought to be associated with the gradual closure of these straits and the increased influence of the high southern latitudes

age dating and the orbital theory of the ice ages development of a high resolution 0 to 300 000 year chronostratigraphy using the concept of orbital tuning a continuous high resolution deep sea chronostratigraphy has been developed spanning the last 300 000 yr the chronology is developed using a stacked oxygen isotope stratigraphy and four different orbital tuning approaches each of which is based upon a different assumption concerning the response of the orbital signal recorded in the data each approach yields a separate chronology the error measured by the standard deviation about the average of these four results which represents the best chronology has an average magnitude of only 2500 yr this small value indicates that the chronology produced is insensitive to the specific orbital tuning technique used excellent convergence between chronologies developed using each of five different paleoclimatological indicators from a single core is also obtained the resultant chronology is also insensitive to the specific indicator used the error associated with each tuning approach is estimated independently and propagated through to the average result the resulting error estimate is independent of that associated with the degree of convergence and has an average magnitude of 3500 yr in excellent agreement with the 2500 yr estimate transfer of the final chronology to the stacked record leads to an estimated error of 1500 yr thus the final chronology has an average error of 5000 yr