propagation of the great salinity anomaly of the 1990s around the northern north atlantic 1 time series of t and s extending through 2001 are used to describe propagation of the great salinity anomaly of the 1990s gsa 90s comparison of the distance time relations for the gsa 70s 80s and 90s reveals a substantial intensification of the large scale circulation in the northern north atlantic especially in the subarctic gyre between newfoundland and the faroes the advection rate of the gsa 70s 80s and 90s between newfoundland and the faroe shetland channel is conservatively estimated to have been 3 5 10 and 10 cm s respectively the circulation intensification apparently occurred within a decade between the gsa 70s and 80s during the next decade the advection rate increased from 10 to 13 cm s between newfoundland and iceland basin the gsa 90s was advected towards the faroe shetland channel by the northern iceland basin s branch of the north atlantic current whereas the contribution of the southern branch via the rockall trough was minimal

regional patterns of sea level change related to interannual variability and multidecadal trends in the atlantic meridional overturning circulation some studies of ocean climate model experiments suggest that regional changes in dynamic sea level could provide a valuable indicator of trends in the strength of the atlantic meridional overturning circulation moc this paper describes the use of a sequence of global ocean ice model experiments to show that the diagnosed patterns of sea surface height ssh anomalies associated with changes in the moc in the north atlantic na depend critically on the time scales of interest model hindcast simulations for 1958 2004 reproduce the observed pattern of ssh variability with extrema occurring along the gulf stream gs and in the subpolar gyre spg but they also show that the pattern is primarily related to the wind driven variability of moc and gyre circulation on interannual time scales it is reflected also in the leading eof of ssh variability over the na ocean as described in previous studies the pattern however is not useful as a fingerprint of longer term changes in the moc as shown with a companion experiment a multidecadal gradual decline in the moc of 5 sv 1 sv 106 m3 s 1 over 5 decades induces a much broader basin scale ssh rise over the mid to high latitude na with amplitudes of 20 cm the detectability of such a trend is low along the gs since low frequency ssh changes are effectively masked here by strong variability on shorter time scales more favorable signal to noise ratios are found in the spg and the eastern na where a moc trend of 0 1 sv yr 1 would leave a significant imprint in ssh already after about 20 years

trace element modelling of pelite derived granites the presence or absence of a vapour phase during incongruent melt reactions of muscovite and biotite together with the composition of the protolith determines the trace element characteristics of the resulting melt provided that equilibrium melting occurs for those phases that host the tracc elements of interest for granitic melts rb sr and ba provide critical constraints on the conditions that prevailed during melting whereas ree are primarily controlled by accessory phase behaviour mass balance constraints for eutectic granites that are formed by the incongruent melting of muscovite in pelites indicate that melting in the presence of a vapour phase will result in a large melt fraction and deplete the restite in feldspar hence the melt will be characterized by low rb sr and high sr ba ratios in contrast vapour absent melting will result in a smaller melt fraction and an increase in the restitic feldspar consequently high rb sr and low sr ba ratios are predicted vapour absent melting will also enhance the negative eu anomaly in the melt granites that result from the incongruent melting of biotite in the source will be characterized by higher rb concentrations than those that result from the incongruent melting of muscovite the himalayan leucogranites provide an example of unfractionated crustally derived eutectic melts that are enriched in rb but depleted in sr and ba relative to their metasedimentary protoliths these compositions may be generated by the incongruent melting of muscovite as a low melt fraction f 0 1 from a pelitic source under vapour absent conditions

causes of decadal changes of the freshwater content in the arctic ocean abstractdecadal changes of the liquid freshwater content in the arctic ocean are studied with a suite of forward and adjoint model simulations adjoint sensitivities show that freshwater volume changes in the norwegian atlantic current north of the lofoten basin and a salinity maximum in the fram strait and in the canadian archipelago lead to an enhanced northward transport of freshwater the dynamical sensitivities indicate that stronger freshwater export from the arctic is related to an enhanced cyclonic circulation around greenland with an enhanced export through the canadian archipelago and a stronger circulation within the fram strait associated with this circulation around greenland is a large scale cyclonic circulation in the arctic cyclonic wind stress anomalies in the arctic ocean as well as over the nordic seas and parts of the subpolar atlantic are optimal to force the freshwater transport changes results from a simulation over the period 1948 2010 corroborate the result that arctic freshwat

response of the global ocean to greenland and antarctic ice melting 1 we investigate the transient response of the global ocean circulation to enhanced freshwater forcing associated with melting of the greenland and antarctic ice sheets increased freshwater runoff from greenland results in a basin wide response of the north atlantic on timescales of a few years communicated via boundary waves equatorial kelvin waves and westward propagating rossby waves in addition modified air sea interaction plays a fundamental role in setting up the basin scale response of the atlantic circulation in its subpolar and subtropical gyres in particular the modified ocean dynamics and thermodynamics lead to a depression in the central north and south atlantic that would not be expected from linear wave dynamics moreover the heat content increases on basin and global scales in response to anomalous freshwater forcing from greenland suggesting that the ocean s response to enhanced freshwater forcing would be a coupled problem other parts of the world ocean experience a much slower adjustment in response to greenland freshwater forcing communicated via planetary waves but also involving advective diffusive processes especially in the southern ocean over the 50 years considered here most of the sea level increase associated with freshwater input from greenland remains in the atlantic ocean in contrast ice melting around antarctica has a much reduced effect on the global ocean in both cases none of the basins came to a stationary state during the 50 year experiment

the transient atmospheric circulation response to north atlantic sst and sea ice anomalies abstract the objective of this study is to investigate the transient evolution of the wintertime atmospheric circulation response to imposed patterns of sst and sea ice extent anomalies in the north atlantic sector using a large ensemble of experiments with the ncar community climate model version 3 ccm3 the initial adjustment of the atmospheric circulation is characterized by an out of phase relationship between geopotential height anomalies in the lower and upper troposphere localized to the vicinity of the forcing this initial baroclinic response reaches a maximum amplitude in 5 10 days and persists for 2 3 weeks diagnostic results with a linear primitive equation model indicate that this initial response is forced by diabatic heating anomalies in the lower troposphere associated with surface heat flux anomalies generated by the imposed thermal forcing following the initial baroclinic stage of adjustment the response becomes progressively more barotropic and increases in both spatial extent an

on the separability of the north atlantic oscillation and arctic oscillation we address the issue of whether the arctic ao and north atlantic oscillations nao are inseparable forming an annular mode in the northern hemisphere atmospheric circulation this annular mode is the leading empirical orthogonal function of hemispheric sea level pressure slp data explaining the largest amount of its variability we examine whether the nao and ao are inseparable spatial modes of the atmospheric circulation using rotated principal component analysis rpca a methodology that identifies simple and unique patterns of spatial dataset variability rpca of the spring summer and autumn slp fields from 1946 1998 reveal nao and ao like patterns occurring as separate regional teleconnections forming the first and second principal components respectively the rpca based nao dipole pattern is like that observed in many prior studies while the ao like pattern exhibits high slp variability over the kara and laptev seas during winter however and in annual analyses a distinct ao like pattern is not obtained and the two patterns may be inseparable using commonly accepted rpca methods the rpca based ao like mode is significantly linked to north central siberian seasonal air temperatures and to the prevailing direction of motion of the underlying arctic ocean in summer suggesting that the non winter ao like pattern as a stand alone teleconnection separate from the nao contributes significantly to high latitude climate and ocean variability the winter nao ao inseparability is discussed as a possible effect of a shared winter storm track between the northeastern atlantic and the arctic

from sediment to granite timescales of anatexis in the upper crust abstract granite formation is the culmination of a sequence of events initiated by prograde heating of the protolith and followed by formation of a grain boundary melt melt segregation into a vein network ascent of the melt through the network and finally crystallisation of the melt experimental constraints on the formation of a crustal melt from the incongruent melting of muscovite combined with geochemical studies of anatexis in the high himalaya allow the timescales required for each of these processes to be assessed discordant temperatures determined from monazite and zircon thermometry for himalayan anatectic granites indicate that at least for some intrusives the melt was undersaturated in lree implying that melts have probably been extracted in less than 10 ka experimental studies suggest that some himalayan melts are also undersaturated in zr implying segregation may have occurred within 100 years such short timescales confirm that deformation driven mechanisms are important in extracting these melts from their source the transport distances of himalayan granitic melts of 10 km may be achieved by the ascent of magma through dykes in about 1 day at such rates even the largest granite could theoretically be emplaced in 10 years crystallisation of himalayan melts involves much longer periods if emplaced as thin sheets 100 m wide a timescale of 500 years is required compared with 30 ka for single stage intrusion of the larger laccoliths for composite sheet complexes magma crystallisation rather than melt ascent comprises the rate determining step on the emplacement of the intrusion the overall timescales of melt segregation and emplacement for many orogenic granites are therefore less than 10 ka and possibly less than 1 ka in contrast the timescale required for prograde heating of the protolith is 1 ma since the melt production rate is determined by heat flow into the protolith and not by reaction kinetics for any geologically significant period we conclude that heat flow determined by both the mechanism of heating and the thermal diffusivities of crustal rocks provides the overall rate determining step of crustal anatexis

recent antarctic ice mass loss from radar interferometry and regional climate modelling observed estimates of ice losses in antarctica combined with regional modelling of ice accumulation in the interior suggest that east antarctica is close to a balanced mass budget but large losses of ice occur in the narrow outlet channels of west antarctic glaciers and at the northern tip of the antarctic peninsula

barotropic models of the extratropical response to el nino abstract a series of linear and nonlinear barotropic models are used to interpret the extratropical response to el niflo equatorial surface temperatures as simulated by an atmospheric general circulation model gcm the divergence time mean vorticity tendency due to transients and the zonal mean tlow are specified from the gcm and the deviation of the streamfunction from its zonal mean at an upper tropospheric level is predicted nonlinearsteady state models suggest that the extratropical wave train is primarily forced from the central rather than the western pacific and that subtropical divergence anomalies are of more importance than tropical anomalies these nonlinear solutions can be reproduced with little loss in accuracy by linearizing about the zonally asymmetric climatological flow if one linearizes about the zonally symmetric flow the part of the solution forced from thewestern pacific deteriorates significantly the solution in the tropics and subtropics also deteriorates if advection of v

on the early response of the climate system to a meltwater input from greenland theearlyresponseoftheatmosphere oceansystemtomeltwaterrunofforiginatingfromthegreenlandice sheet is studied using a coupled atmosphere ocean general circulation model aogcm for this purpose aogcm ensemble simulations without and with associated ocean freshening around greenland are compared for freshwater perturbations initiated in northern winter the mean response for the first three months shows the emergence of negative sea surface temperature sst anomalies in the denmark strait in association with enhanced oceanic advection by the east greenland current the response also shows negative sst anomalies in the north atlantic associated with enhanced westerlies at the ocean surface additionally the baroclinic atmospheric cyclonic circulation east of greenland intensifies and anticyclonic circulations with equivalent barotropic structures develop over western europe and the north pacific ocean simulations by the atmospheric component of the aogcm indicate that atmosphere ocean interactions contribute significantly to enhance the response the sensitivity of the coupled system response to the timing of freshwater perturbation is also studied for freshwater perturbations initialized in northern summer the response during the following winter is similar but stronger in magnitude in the northern hemisphere the atmospheric response resembles the arctic oscillation ao mode of variability the association between anomalies in the denmark strait ssts and in the atmosphere east of greenland is consistent with that observedduring previous great salinityanomaly gsa events the results obtainedhighlight the importance ofatmosphere oceaninteractionintheearlyclimateresponsetogreenlandmelting theteleconnectionswith the north pacific and the contribution of gsa events to north atlantic oscillation nao variability

classification seasonality and persistence of low frequency atmospheric circulation patterns abstract orthogonally rotated principle component analysis rpca of northern hemisphere 1 month mean 700 mb heights is used to identify and describe the seasonality and persistence of the major modes of interannual variability the analysis is detailed and comprehensive in that 1 a high resolution approximately equal area 358 point grid is used for the virtually maximum possible 35 year period of record 2 a positive bias in the nmc data base in the early 1950s in the subtropics is largely eliminated for the first time and 3 homogeneous separate analyses of each month of the year are carried out detailing the mouth to month changes in the dominant circulation patterns winter results are similar to those of other recent rpca and teleconnection studies except that some less obvious patterns are identified and further detail of the better known patterns is provided two north south dipole patterns are found over the pacific ocean west pacific oscillation and east pacific pattern and over the atla

the modeled atmospheric response to midlatitude sst anomalies and its dependence on background circulation states abstract the atmospheric response to a midlatitude sst anomaly in the north pacific and its dependence on background flow are examined in a gcm experiments are conducted using the same warm sst anomalies but two different model states perpetual january and perpetual february the atmospheric responses to the sst anomalies are statistically significant in both january and february but are completely different the anomalous circulation in january is characterized by a trough decaying with height immediately downstream of the sst anomalies in february the anomalous circulation is dominated by a downstream ridge growing with height the patterns of the anomalous heights in the two months are nearly orthogonal vorticity and thermodynamic budgets are diagnosed to illustrate how the anomalous circulations are maintained over the sst anomalies low level convergence and ascent are observed in both months in january the anomalous convergence is balanced by a residual due primarily to the forcing by submont

chemical gradient at diffusive interfaces in magma chambers the results of high pressure experiments on diffusion and soret separation in natural silicate melts show that the diffusive behaviour between natural silicic and mafic magmas can be approximately modelled as if the system were a binary mixture of sio2 and other components such as mgo feo cao steady state compositional profiles across a diffusive interface between silicic and mafic magma layers are calculated on the basis of phenomenological relationships for the fluxes of chemical species and heat in the binary mixtures using the experimental data of diffusion and soret coefficients in natural silicate melts the compositional profiles show a curvature with a minimum sio2 value within the interface due to the soret effect and temperature dependence of diffusion coefficient the compositional gradient at the lower half of the diffusive interface is similar to that resulting from the soret separation of a mafic melt regardless of the composition of the silicic magmas these results suggest that picritic magma can be formed in the interfacial region between the mafic and silicic magma layers the compositional gradient explains chemical variation of mafic to picritic inclusions in a mixed andesite of the abu volcano group japan

a reconciled estimate of ice sheet mass balance we combined an ensemble of satellite altimetry interferometry and gravimetry data sets using common geographical regions time intervals and models of surface mass balance and glacial isostatic adjustment to estimate the mass balance of earth s polar ice sheets we find that there is good agreement between different satellite methods especially in greenland and west antarctica and that combining satellite data sets leads to greater certainty between 1992 and 2011 the ice sheets of greenland east antarctica west antarctica and the antarctic peninsula changed in mass by 142 49 14 43 65 26 and 20 14 gigatonnes year 1 respectively since 1992 the polar ice sheets have contributed on average 0 59 0 20 millimeter year 1 to the rate of global sea level rise

on the mechanisms behind salinity anomaly signals of the northern north atlantic abstract hydrographic time series from the northern north atlantic throughout the 20th century show oscillations in temperature and salinity at more or less regular intervals the great salinity anomalies described during the 1970s dickson r r meincke j malmberg s a lee a j 1988 the great salinity anomaly in the north atlantic 1968 1982 progress in oceanography 20 103 151 during the 1980s belkin i m levitus s antonov j malmberg s a 1998 great salinity anomalies in the north atlantic progress in oceanography 41 1 68 and during the 1990s belkin i m 2004 propagation of the great salinity anomaly of the 1990s around the northern north atlantic geophysical research letters 31 8 l08306 doi 10 1029 2003gl019334 have distinct amplitudes and all three of them were interpreted as low salinity anomalies propagating downstream through the anti clockwise circulation system of the northern north atlantic ocean further inspection of time series from the northeast atlantic and the northwest atlantic over the past century shows however several other distinct negative anomalies of lesser amplitudes additionally a number of high salinity anomalies can be identified the present paper analyses further the propagation of the negative and positive anomalies and links them together it is shown that they have varying speeds of propagation and that the varying speeds are correlated across the north atlantic we propose that varying volume fluxes in and out of the arctic basin is the causal mechanism behind the anomaly signals and that the north atlantic oscillation nao partly has influence on the flux variations described periods of large decadal scale amplitudes of the nao coincide with periods of large decadal scale oscillation in the marine climate

increasing rates of ice mass loss from the greenland and antarctic ice sheets revealed by grace 1 we use monthly measurements of time variable gravity from the grace gravity recovery and climate experiment satellite gravity mission to determine the ice mass loss for the greenland and antarctic ice sheets during the period between april 2002 and february 2009 we find that during this time period the mass loss of the ice sheets is not a constant but accelerating with time i e that the grace observations are better represented by a quadratic trend than by a linear one implying that the ice sheets contribution to sea level becomes larger with time in greenland the mass loss increased from 137 gt yr in 2002 2003 to 286 gt yr in 2007 2009 i e an acceleration of 30 11 gt yr2 in 2002 2009 in antarctica the mass loss increased from 104 gt yr in 2002 2006 to 246 gt yr in 2006 2009 i e an acceleration of 26 14 gt yr2 in 2002 2009 the observed acceleration in ice sheet mass loss helps reconcile grace ice mass estimates obtained for different time periods

timescales of disequilibrium melting in the crust constraints from modelling the distribution of multiple trace elements and a case study from the lesser himalayan rocks of sikkim a numerical code has been developed to track the distribution of trace elements in crustal rocks undergoing melting the model handles diffusion with moving boundaries and accounts for the processes of diffusion dissolution and precipitation in a partially molten system among the various input parameters for modelling source composition i e modal abundance variations diffusion coefficients and partition coefficients are found to exert a significant control on the melt chemistry the other inputs such as melt reaction stoichiometry kinetics of melting and grain size of protolith have lesser influence exploration of the general behaviour indicates that for systems in which disequilibrium melting of the kind considered in this paper occurs trace element concentrations may be used to constrain the composition of the protolith or the timescales of melting depending on the specific circumstances after exploring some general features of melting in a pelitic system the model is applied to calculate trace element distributions in migmatites from the lesser himalayan rocks in sikkim india we focus on the distribution of trace elements during the initial stages of melt formation these partially molten rocks show disequilibrium distribution of trace elements and the numerical code is capable of quantitatively reproducing many of the observed patterns the results of the modelling indicate that the observed melts in this zone were formed within 50 000 years and that segregation of melts into leucosome and restite was complete between 50 000 and 250 000 years these short timescales may point to deformation enhanced melt segregation at least on a hand specimen scale it is important to distinguish between timescales of segregation over these scales and timescales of removal of melt on an outcrop scale to form plutons the latter requiring higher degrees of melting and larger distances of migration take longer

mineral leucosome trace element partitioning in a peraluminous migmatite a laser ablation icp ms study abstract low pressure peraluminous migmatites from the pena negra complex central spain stayed partially molten for 5 10 ma until the latest hercynian deformations allowed the segregation of melt as discordant leucosome veins due to long residence within its source the melt had time enough to equilibrate with residual phases other than accessories included within major refractory minerals we estimated crystal melt partition coefficients as the concentration ratios between leucosome samples appearing to be pure melts and mesosome minerals which were analyzed for trace elements with a laser probe coupled to an icp mass spectrometer our data reveal that when biotite is stable li rb cs tl sc v cr ni nb and ta become strongly compatible the role of biotite in fractionating y th u and the ree is insignificant cordierite strongly fractionates li and be and also has some effect on the hree and u garnet produces extreme fractionation of sc y and the hree the ree partition coefficients for ferromagnesian silicates increase with the atomic number this effect being progressively more important through biotite cordierite and garnet k feldspar strongly fractionates ba and pb but plays a secondary role to biotite for rb and cs k feldspar fractionation does not change the lree hree ratio but plagioclase fractionation produces a significant decrease in the lree hree ratio both feldspars greatly fractionate sr and eu monazite fractionation produces a dramatic depletion in ree th and u as well as a decrease in the lree hree and th u ratios apatite also fractionates the ree although in contrast with monazite it increases the lree hree ratio and does not affect the th u ratio zircon fractionation like apatite and garnet produces a strong depletion in the hree and a concomitant increase in the lree hree ratio in contrast with monazite zircon fractionation causes the th u ratio to increase

revisiting greenland ice sheet mass loss observed by grace in this paper we discuss a new method for determining mass time series for 16 hydrological basins representing the greenland system gs whereby we rely on gravity recovery and climate experiment grace mission data in the same analysis we also considered observed mass changes over ellesmere island baffin island iceland and svalbard ebis the summed contribution of the complete system yields a mass loss rate and acceleration of 252 28 gt yr and 22 4 gt yr2 between march 2003 and february 2010 where the error margins follow from two glacial isostatic adjustment gia models and three processing centers providing grace monthly potential coefficient sets we describe the relation between mass losses in the gs and the ebis region and found that the uncertainties in all areas are correlated the summed contribution of ellesmere island baffin island iceland and svalbard yields a mass loss rate of 51 17 gt yr and an acceleration of 13 3 gt yr2 between march 2003 and february 2010 the new regional basin reconstruction method shows that the mass loss within the southeastern basins in the gs has slowed down since 2007 while mass loss in western basins increased showing a progression to the north of greenland

observed atmospheric response to cold season sea ice variability in the arctic abstractthe relation between weekly arctic sea ice concentrations sics from december to april and sea level pressure slp during 1979 2007 is investigated using maximum covariance analysis mca in the north atlantic sector the interaction between the north atlantic oscillation nao and a sic seesaw between the labrador sea and the greenland barents sea dominates the nao drives the seesaw and in return the seesaw precedes a midwinter spring nao like signal of the opposite polarity but with a strengthened northern lobe thus acting as a negative feedback with maximum squared covariance at a lag of 6 weeks statistical significance decreases when slp is considered in the whole northern hemisphere but it increases when north pacific sic is included in the analysis the maximum squared covariance then occurs after 8 weeks resembling a combination of the nao response to the atlantic sic seesaw and the aleutian icelandic low seesaw like response to in phase sic changes in the bering and okhotsk seas w

granite magma formation transport and emplacement in the earth s crust the origin of granites was once a question solely for petrologists and geochemists but in recent years a consensus has emerged that recognizes the essential role of deformation in the segregation transport and emplacement of silica rich melts in the continental crust accepted petrological models are being questioned either because they require unrealistic rheological behaviours of rocks and magmas or because they do not satisfactorily explain the available structural or geophysical data provided flow is continuous mechanical considerations suggest that far from being geologically sluggish granite magmatism is a rapid dynamic process operating at timescales of or 100 000 years irrespective of tectonic setting

application of proton microprobe data to trace element partitioning in volcanic rocks abstract in situ proton microprobe analyses are used to calculate partition coefficients d for mn fe cu zr ga rb sr y zr ba and pb and more restricted data for ni ge nb and as in pyroxenes olivines biotites hornblendes feldspars and fe ti oxides from 32 volcanic rocks all with glassy matrices these range from leucitite and picritic basalt to high silica and peralkaline rhyolites data for pigeonite aenigmatite allanite and leucite are included where these occur as additional phases coefficients have also been rayleigh fractionation corrected this correction being significant for those elements with high coefficients in specific mineral phases and where pronounced crystallisation has occurred partition coefficients are presented as plots vs nbo t ratios of coexisting glasses melts negative correlations are observed based on log transformed data for mn zn fe and y pyroxenes mn zn and fe olivines mn and zn fe ti oxides nb y sr and zr hornblendes sr plagioclase and ba cu zn and fe feldspars these element partition coefficients are not correlated with bulk mineral compositions zr and ga are not correlated with either glass or mineral compositions except zr in hornblende the correlations support the role of melt polymerisation in controlling element partitioning at least for those elements listed although possible effects of temperature cannot be uniquely isolated in the present data set the coefficients presented here are generally consistent with published values for equivalent rock compositions one difference being that some of the higher published coefficients e g sr in plagioclase are not matched in the present data set from a practical viewpoint the relative regularity of the coefficient data presented does facilitate choice of relevant coefficient values that are required for petrogenetic modelling procedures

oceanic isopycnal mixing by coordinate rotation abstract current numerical models of oceanic circulation differentiate between the eddy diffusion and viscosity transport along the geopotential horizontal and vertical directions only in order to model the effect of anisotropic turbulence as diffusive transport along and across density surfaces the isopycnal mixing tensor has been transformed from a diagonal second rank tensor in the isopycnal coordinate system to a tensor containing off diagonal elements in the geopotential coordinate system

north pacific climate response to freshwater forcing in the subarctic north atlantic oceanic and atmospheric pathways abstract sudden changes of the atlantic meridional overturning circulation amoc are believed to have caused large abrupt climate changes over many parts of the globe during the last glacial and deglacial period this study investigates the mechanisms by which a large freshwater input to the subarctic north atlantic and an attendant rapid weakening of the amoc influence north pacific climate by analyzing four different ocean atmosphere coupled general circulation models gcms under present day or preindustrial boundary conditions when the coupled gcms are forced with a 1 sv sv 106 m3 s 1 freshwater flux anomaly in the subarctic north atlantic the amoc nearly shuts down and the north atlantic cools significantly the south atlantic warms slightly shifting the atlantic intertropical convergence zone southward in addition to this atlantic ocean atmosphere response all of the models exhibit cooling of the north pacific especially along the oceanic frontal zone consistent with paleoclimate recons

geothermobarometry and development of inverted metamorphism in the darjeeling sikkim region of the eastern himalayan abstract the darjeeling sikkim region provides a classic example of inverted himalayan metamorphism the different parageneses of pelitic rocks containing chlorite biotite garnet staurolite kyanite sillimanite plagioclase and k feldspar are documented by a variety of textures resulting from continuous and discontinuous reactions in the different zones microprobe data of coexisting minerals show that xmg varies in the order garnet staurolite biotite chlorite white mica is a solid solution between muscovite and phengite garnet is mostly almandine rich and shows normal growth zoning in the lower part of the main central thrust mct zone and reverse zoning in the upper part of the zone chemographical relations and inferred reactions for different zones are portrayed in afm space in the low grade zones oriented chlorites and micas and rolled garnets grew syntectonically and were succeeded by cross cutting chlorites and micas and garnet rims in the upper zones sillimanite kyanite and staurolite crystallized during a static inter kinematic phase p t contitions of metamorphism estimated through different models of geothermobarometry are estimated to have been 580 c for the garnet zone to a maximum of 770 c for the sillimanite zone the preferred values of pressure range from 5 0 kbar to 7 7 kbar models to explain the inverted metamorphism include overthrusting of a hot high himalayan slab along a c 5 km wide ductile mct zone and the syn or post metamorphic folding of isograds

strontium diffusion kinetics in plagioclase feldspars abstract a pronounced systematic variation by a factor of 104 in strontium tracer diffusion kinetics has been measured in four feldspar specimens that span the range of plagioclase compositions the experiments were carried out in air using isotopically enriched 86sr and analyzed with an ion microprobe the arrhenius relations for strontium diffusion are plagioclase ab d 0 m 2 s q kj g atom sr temperature range c albite 98 1 05e 279 550 1080 oligoclase 73 7 4e 7 261 750 1100 labradorite 36 1 1e 6 295 800 1300 anorthite 4 5 7e 9 267 900 1300 no difference was found within the measurement uncertainty between these results and those for albite under hydrothermal conditions at one kbar 100 mpa strontium diffusion rates parallel to and normal to the albite twin lamellae are the same there is no effect on the strontium diffusion rate in albite between an oxygen partial pressure of 0 2 bar air and 10 19 bar graphite buffer a relation is developed that permits calculation of the strontium diffusion arrhenius parameters for any normal plagioclase composition the large variation in diffusion rate with composition has a major effect on the ease of strontium isotope exchange at high temperatures between plagioclases and coexisting phases in rocks which may prove useful in determining thermal histories

the parameterization of the planetary boundary layer in the ucla general circulation model formulation and results abstract a planetary boundary layer pbl parameterization for general circulation models gcms is presented it uses a mixed layer approach in which the pbl is assumed to be capped by discontinuities in the mean profiles both clear and cloud topped boundary layers are parameterized particular emphasis is placed on the formulation of the coupling between the pbl and both the free atmosphere and cumulus convection for this purpose a modified coordinate is introduced in which the pbl top and the lower boundary are both coordinate surfaces the use of a bulk pbl formulation with this coordinate is extensively discussed results are presented from a july simulation produced by the ucla gcm pbl related variables are shown to illustrate the various regimes the parameterization is capable of simulating

constraints on melting and magma production in the crust major intrusions of granitic rocks are found in several tectonic settings and in all cases crustal melts may contribute to the volumes of granitic magma high grade metamorphism and partial melting of the crust take place predominantly under fluid absent conditions we present a model for calculating the amounts of melt that may be formed by fluid absent breakdown of micas and amphiboles in common crustal rock types pelitic quartzofeldspathic intermediate and mafic melt proportions depend mainly on the kind of source rock the pressure at which melting takes place the temperature and the hydrous mineral h2o content of the source as a consequence of the pressure dependence of water solubility in silicate melts any given source rock will produce more melt by a given fluid absent reaction at lower pressure at a given pressure higher temperature reactions can produce more melt from a given source rock based on a survey of the compositions of common rock types we show that the amounts of melt can vary from 50vol thus crustal rocks vary widely in their fertility as magma sources depending on the types and amounts of hydrous phases they contain in general muscovite breakdown in pelites will yield only small quantities of melt and lead to migmatite formation biotite breakdown in pelites occurs at higher temperature and because most high grade pelites below granulite grade are biotite rich can yield up to about 50 vol of granitoid melt rocks of intermediate composition and hornblende rich mafic rocks are potentially highly fertile magma sources also provided that the high temperatures necessary for biotite and hornblende breakdown are realized pyroxene rich mafic rocks and quartzofeldspathic rocks are much less potentially fertile data suggest that mechanisms exist for the efficient segregation of melt and restite in systems with 30and probably 20vol melt the pressure temperature history of a region can greatly influence crustal source fertility through its control over the occurrence of subsolidus dehydration and melting equilibria

rb and sr diffusion in alkali feldspars with implications for cooling histories of rocks volume diffusion of sr and rb have been measured in orthoclase feldspar as well as sr diffusion in albite using isotope exchange hydrothermal solutions and analysis by ion microprobe for water pressures of 1000 bars the kinetics of diffusion parallel to the c axis are given by the arrhenius relations strontium transport in the or is isotropic to within a factor of four the measurements use a new procedure in which only the elements in question are present in the starting fluid cooling rates of igneous rocks can now be determined from closure temperatures material balance and the rb sr dating method based on the departures of the mineral data from the whole rock isochron for cooling rates less than approximately 100 c my

barotropic wave propagation and instability and atmospheric teleconnection patterns abstract a global barotropic model linearized about the 300 mb climatological mean january flow is perturbed by applying a series of localized forcings distributed throughout the tropics and subtropics structures which resemble the observed pacific north american and east atlantic teleconnection patterns noted by wallace and gutzler 1981 tend to recur in the responses similar patterns are found to result from the dispersion of isolated initial perturbations placed at a variety of locations in the tropics and midlatitudes it is shown that these structures are related to the most rapidly growing mode associated with barotropic instability of the zonally varying climatological basic state in the absence of damping this mode has an e folding time of about a week and a period close to 50 days in localized regions the instantaneous growth rates can be competitive with those of baroclinic instability these episodes of rapid local barotropic growth are interspersed with intervals in which the local

magmatic residence times of zoned phenocrysts introduction and application of the binary element diffusion modelling bedm technique this paper describes a general technique binary element diffusion modelling bedm for determining single crystal residence times in magmas that relies on modelling the diffusion of two or more elements in the crystal bedm has the advantage over other diffusion based models in that it does not need a precisely defined initial compositional profile for the crystal at zero time and instead requires that the concentrations of the two elements are correlated during crystallisation any differences subsequently observed between the two elements are caused by intracrystalline diffusion during residence in hot magma these differences are removed by artificially ageing the slower diffusing of the two elements and the amount of time taken to undo the difference between the elements is simply related to the crystal residence time decoupling time at high temperatures the bedm principle is demonstrated using artificial data and is then applied to literature data for sr and ba in a zoned sanidine crystal from the bishop tuff anderson et al in j petrol 41 3 449 473 2000 for this crystal the method gives a residence time estimate of 114 ka at 800 c which is then compared with estimates from other methods in theory the method can be further expanded for use as a geothermometer as well as geochronometer however this is not easily possible with the diffusivity data currently available

introduction to special section mechanisms and consequences of melt segregation from crustal protoliths much of petrological effort has been concentrated on attempts to discern define and interpret regularly recurring patterns of igneous and metamorphic petrogenesis i do not wish to deny the existence of discernable order in petrological phenomena but every geological event is unique if we concentrate too heavily on discernment of order and pattern may we perhaps overlook or underrate the unique quality of each igneous or metamorphic episode possibly uniqueness may have as great significance as conformity to pattern when we attempt to fit the phenomena of petrology into the broad framework of geology francis j turner presidential address to the mineralogical society of america november 11 1969

loss of metals from pelites during regional metamorphism in aluminous metapelites the ratio h2o k2o decreases with increasing metamorphic grade and degree of reaction this ratio is a very practical indicator for the progress of the mineral reconstitution during progressive metamorphism with decreasing values of the ratio h2o k2o the cu concentration and the following element ratios also decrease either continuously or in stepwise fashion tl k2o ba k2o pb k2o bi k2o hg k2o sr na2o zn fe2 mg cd fe2 mg rb k2o remains approximately constant in the aluminous metapelites of the damara orogen in namibia the following losses occur between the biotite isograd and anatexis 61 cu 20 tl 34 ba 59 pb 86 bi 46 hg 30 sr 25 zn 31 cd thus the potential of regional metamorphism to form hydrothermal deposits in the low grade environment should not be neglected

arctic oscillation or north atlantic oscillation abstract the definition and interpretation of the arctic oscillation ao are examined and compared with those of the north atlantic oscillation nao it is shown that the nao reflects the correlations between the surface pressure variability at its centers of action whereas this is not the case for the ao the nao pattern can be identified in a physically consistent way in principal component analysis applied to various fields in the euro atlantic region a similar identification is found in the pacific region for the pacific north american pna pattern but no such identification is found here for the ao the ao does reflect the tendency for the zonal winds at 35 and 55 n to anticorrelate in both the atlantic and pacific regions associated with the nao and pna because climatological features in the two ocean basins are at different latitudes the zonally symmetric nature of the ao does not mean that it represents a simple modulation of the circumpolar flow an increase in the ao or nao implies stro

geochemistry of granitic melts produced during the incongruent melting of muscovite implications for the extraction of himalayan leucogranite magmas granites formed from the incongruent melting of muscovite preserve evidence of chemical disequilibrium both between restite and melt and within individual phases during crystallization miocene leucogranites from the himalayas exemplify melt formation from a pelitic source through muscovite breakdown trace elements rb sr and ba reside almost entirely in micas and feldspars in the protolith and therefore are appropriate for quantitative modeling equilibrium between melt and restite was established for sr under fluid phase absent conditions but some disequilibrium is suggested by rb abundances consistent with known diffusion rates for these elements abundances of zr and ce in the granites which provide essential structural constituents of the accessory phases zircon and monazite indicate that the himalayan granite melts were saturated in these elements if formed at temperatures of 670 700 c but slightly undersaturated if formed at somewhat higher temperatures ion microprobe traverses of y and heavy rare earth element abundances across single garnets from himalayan leucogranites provide evidence for disequilibrium between melt and phenocryst during garnet growth thus constraining the evolution of magmatic compositions during melt formation and crystallization in some cases progressive melting of the protolith and crystallization of garnet from the melt appear to have been concomitant indicating that minimum solidus temperatures were maintained throughout the melting event this suggests that thermal buffering rather than melt extraction provided the rate determining step in magma evolution

intercrystalline stable isotope diffusion a fast grain boundary model we formulated a numerical model for stable isotope interdiffusion which predicts the temperatures recorded between two or more minerals and the intragranular distribution of stable isotopes in each mineral as functions of mineral grain sizes and shapes diffusivities modes equilibrium isotopic fractionations and the cooling rate of a rock one of the principal assumptions of the model is that grain boundaries are regions of rapid transport of stable isotopes this fast grain boundary fgb model describes interdiffusion between any number of mineral grains assuming that local equilibrium and mass balance restrictions apply on the grain boundaries throughout the volume modeled the model can be used for a rock containing any number of minerals and number of grain sizes of each mineral several grain shapes and any thermal history or domain size desired previous models describing stable isotope interdiffusion upon cooling have been based on dodson s equation or an equivalent numerical analogue the closure temperature of dodson is the average bulk temperature recorded between a mineral and an infinite reservoir by using dodson s equation these models have treated the closure temperature as an innate characteristic of a given mineral independent of the amounts and diffusion rates of other minerals such models do not accurately describe the mass balance of many stable isotope interdiffusion problems existing models for cation interdiffusion could be applied to stable isotopes with some modifications but only describe exchange between two minerals under specific conditions the results of fgb calculations differ considerably from the predictions of dodson s equation in many rock types of interest actual calculations using the fgb model indicate that closure temperature and diffusion profiles are as strongly functions of modal abundance and relative differences in diffusion coefficient as they are functions of grain size and cooling rate closure temperatures recorded between two minerals which exchanged stable isotopes by diffusion are a function of modal abundance and differences in diffusion coefficient and may differ from that predicted by dodson s equation by hundreds of degrees c either or both of two minerals may preserve detectable zonation which may in some instances be larger in the faster diffusing mineral rocks containing three or more minerals can record a large span of fractionations resulting from closed system processes alone the results of fgb diffusion modeling indicate that the effects of diffusive exchange must be evaluated before interpreting mineral fractionations concordant or discordant recorded within any rock in which diffusion could have acted over observable scales the predictions of this model are applicable to thermometry evaluation of open or closed system retrogression and determination of cooling rates or diffusion coefficients

rapid barotropic sea level rise from ice sheet melting 1 sea level rise associated with idealized greenland and antarctic ice sheet melting events is examined using a global coupled ocean sea ice model that has a free surface formulation and thus can simulate fast barotropic motions the perturbation experiments follow the coordinated ocean ice reference experiment core version iii all regions of the global ocean experience a sea level rise within 7 8 days of the initialization of a polar meltwater input of 0 1 sv 1 sv 106 m3 s 1 the fast adjustment contrasts sharply with the slower adjustment associated with the smaller steric sea level evolution that is also connected with melt events the global mean sea level rises by 9 mm yr 1 when this forcing is applied either from greenland or antarctica nevertheless horizontal inter basin gradients in sea level remain for climate adaption in low lying coastal and island regions it is critical that the barotropic sea level signal associated with melt events is taken into consideration as it leads to a fast sea level rise from melting ice sheets for the bulk of the global ocean a linear relation between sea level rise and global meltwater input is further supported by experiments in which idealized melting occurs only in a region east or west of the antarctic peninsula and when melting rates are varied between 0 01 sv and 1 0 sv the results indicate that in ocean models that do not explicitly represent the barotropic signal the barotropic component of sea level rise can be added off line to the simulated steric signal

gneiss domes and orogeny many gneiss domes record positive feedback between decompression and partial melting of orogenic middle crust exhumed orogens are riddled with gneiss domes cored by migmatites that underwent dehydration melting during decompression the decreasing buoyancy associated with increasing melt fraction drives further decompression at near isothermal conditions as the partially molten crust rises diapirically this combination of processes may explain the generation and retention of large volumes of crustally derived melt recorded in many deep seated migmatite terranes and inferred for active orogens in exhumed orogens the signature of the rapid ascent of partially molten crust is a gneiss dome cored by migmatite granite the large volume of material involved in the vertical transfer of partially molten crust indicates that the formation of gneiss domes is an efficient mechanism for heat advection during orogenesis

decadal variability of hydrography in the upper northern north atlantic in 1948 1990 we investigate the variability of the north atlantic subarctic gyre in recent decades from time series of station temperature and salinity decadal variability stronger at the surface is identified which exhibits vertical coherence over a layer deeper than the late winter mixed layer in the northwestern atlantic it corresponds to the layer with a component of water from the arctic ocean or from the canadian arctic the spatial coherence of the signal is investigated an empirical orthogonal function decomposition of lagged time series indicates that a single pattern explains 70 of the variance in upper ocean salt content corresponding to a propagating signal from the west to the northeast in the subarctic gyre the most likely interpretation is that the salinity signal originates in the slope currents of the labrador sea and is diffused advected eastward of the grand banks over the near western atlantic in the northwestern atlantic temperature fluctuations are strongly correlated to salinity fluctuations and are aligned along the average t s characteristics this signal suggests large variations in the outflow of fresh cold water in the slope current and is strongly correlated with ice cover a basin scale atmospheric circulation of weakened westerlies at 55 n weaker northwesterlies west of greenland and weaker southerlies over the central and eastern north atlantic is associated with the high salinity and warm water phase of the first principal component this circulation pattern leads fluctuations in the northeast atlantic and lags those in the northwestern part of the basin the wind indices also suggest that the fluctuations of the fresh water outflow occur during intervals of anomalously northerly winds either east of greenland 1965 1968 1969 or off the canadian archipelago 1983 1984

great salinity anomalies in the north atlantic abstract we revisited the great salinity anomaly of the 1970s gsa 70s dickson et al 1988 and documented the newly identified great salinity anomaly of the 1980s hence termed gsa 80s both propagated around the north atlantic in a similar fashion the advective mechanism initially proposed to explain the observed sequence of low salinity low temperature events during the gsa 70s apparently holds also for the gsa 80s the latter was successively observed in the west greenland current 1982 the labrador current 1983 the flemish pass 1984 near the charlie gibbs fracture zone 1984 1985 in the rockall channel 1985 south of iceland 1985 1988 in the north sea 1986 1987 norwegian sea 1987 1988 barents sea 1988 1989 and iceland sea 1989 1990 the advection speed of the gsa 80s seems to be greater than the one of the gsa 70s the 1980s anomaly reached the barents sea 6 to 7 years after peaking in the west greenland current while the 1970s anomaly traveled the same route in 8 to 10 years these anomalies however seem to be of different origin the gsa 70s was apparently boosted remotely by a freshwater sea ice pulse from the arctic via fram strait consequently the gsa 70s was accompanied by a large sea ice extent anomaly in the greenland and iceland seas which propagated into the labrador sea in contrast the gsa 80s was likely formed locally in the labrador sea baffin bay mainly because of the extremely severe winters of the early 1980s but supplemented with a possible contribution of the arctic freshwater outflow via the canadian archipelago facilitated by strong northerly winds which would have enhanced stability and ice formation this anomaly was also associated with a positive sea ice extent anomaly in the labrador sea baffin bay which however had no upstream precursor in the greenland sea thus the gsas are not necessarily caused solely by an increased export of freshwater and sea ice from the arctic via fram strait these results are corroborated by the early 1990s data when a new fresh cold anomaly was formed in the labrador sea and accompanied by a large positive sea ice extent anomaly the harsh winters of the early 1990s were however confined to the labrador sea baffin bay area while the atmospheric and oceanic conditions in the greenland iceland and irminger seas were normal the labrador sea baffin bay area appears therefore to play a key role in formation of gsas as well as in propagation of the gsas formed upstream a likely contribution of the enhanced canadian archipelago freshwater outflow to the gsa formation also seems to be significant two major modes of the gsa origin are thus identified remote generated by an enhanced arctic ocean freshwater export via either fram strait or the canadian archipelago and local resulting from severe winters in the labrador sea baffin bay both modes should be taken into account to model decadal variability of the coupled ocean ice atmosphere system in the north atlantic arctic basin

isopycnal mixing in ocean circulation models abstract a subgrid scale form for mesoscale eddy mixing on isopycnal surfaces is proposed for use in non eddy resolving ocean circulation models the mixing is applied in isopycnal coordinates to isopycnal layer thickness or inverse density gradient as well as to passive scalars temperature and salinity the transformation of these mixing forms to physical coordinates is also presented

long term coordinated changes in the convective activity of the north atlantic the north atlantic is a peculiarly convective ocean the convective renewal of intermediate and deep waters in the labrador sea and greenland iceland sea both contribute significantly to the production and export of north atlantic deep water thus helping to drive the global thermohaline circulation while the formation and spreading of 18 degree water at shallow to intermediate depths off the us eastern seaboard is a major element in the circulation and hydrographic character of the west atlantic for as long as time series of adequate precision have been available to us it has been apparent that the intensity of convection at each of these sites and the hydrographic character of their products have been subject to major interannual change as shown by aagaard 1968 clarke et al 1990 and meincke et al 1992 for the greenland sea in the ows bravo record from the labrador sea eg lazier 1980 et seq and at the panulirus hydrostation s site in the northern sargasso off bermuda eg jenkins 1982 talley and raymer 1982 this paper reviews the recent history of these changes showing that the major convective centres of the greenland and labrador seas are currently at opposite convective extrema in our postwar record with vertical exchange at the former site limited to 1000 m or so but with labrador sea convection reaching deeper than previously observed to over 2300 m as a result greenland sea deep water has become progressively warmer and more saline since the early 70 s due to increased horizontal exchange with the arctic ocean through fram strait while the labrador sea water has become progressively colder and fresher over the same period through increased vertical exchange most recently convection has become deep enough there to reach into the more saline nadw which underlies it so that cooler but now saltier and denser lsw has resulted

calculation of river discharge and prediction of lake height from satellite radar altimetry example for the lake chad basin 1 the application of satellite radar altimetry to the determination of lake and river elevations has been used in numerous projects and is well validated here we show that with the aid of ground based information this technique can be extended to determine river discharge and predict downstream lake and marsh height the lake chad basin provides an ideal case study due to its well known hydrology and complex lake and marsh morphology and because prediction of lake and marsh height has been identified as potentially useful to people living in the region altimetric stage measurements from the topex poseidon satellite at the chari ouham confluence estimate river discharge about 500 km downstream at n djamena 10 days in advance r2 0 9562 via simple linear correlation methods the stage measurements successfully estimate the height of the permanent waters of the lake 600 km downstream 39 days in advance r2 0 9297 predicting the water height on the western marshes of the lake bed is poorer r2 0 7958 due to a change in response time of the local stage to the seasonal floods coincident with an observed increase in mean water level in the latter half of the 1990s before 1997 a 96 day phase lag results in the best fit r2 0 6463 after 1997 the best fit is obtained with a 66 day phase lag r2 0 8139 the excellent river discharge and lake height predictions show that altimetry is a useful tool where ground based data are difficult to obtain and where rapid water resource assessment is desirable

the great salinity anomaly in the northern north atlantic 1968 1982 the widespread freshening of the upper 500 800m layer of the northern north atlantic which this paper describes represents one of the most persistent and extreme variations in global ocean climate yet observed in this century though a range of explanations have been advanced to explain this event including in situ changes in the surface moisture flux this paper describes the great salinity anomaly as largely though not entirely an advective event traceable around the atlantic subpolar gyre for over 14 years from its origins north of iceland in the mid to late 1960s until its return to the greenland sea in 1981 1982 the overall propagation speed around this subpolar gyre is estimated at about 3cm s 1 of the total salt deficit associated with the anomaly as it passed south along the labrador coast in 1971 1973 about 72 109 tonnes a deficit equivalent to about two thirds of this figure 47 109 tonnes ultimately passed through the faroe shetland channel to the barents sea arctic ocean and greenland sea during the mid 1970s possible effects on deep water formation and on the representativeness of historical section data are discussed

volcanic eruptions large scale modes in the northern hemisphere and the el nino southern oscillation abstract the author analyzes the impact of 13 major stratospheric aerosol producing volcanic eruptions since 1870 on the large scale variability modes of sea level pressure in the northern hemisphere winter the paper focuses on the arctic oscillation ao and the north atlantic oscillation nao to address the question about the physical nature of these modes the hypothesis that the phase of the el nino southern oscillation enso may control the geographical extent of the dominant mode in the northern hemisphere is also investigated as well as the related possibility that the impact of the eruptions may be different according to the phase of enso the author finds that both the ao and the nao are excited in the first winter after the eruptions with statistical significance at the 95 level both the signal and the significance are larger for the nao than for the ao the excitation of the ao and the nao is connected with the excitation of a secondary mode which resembles an augmented pacific north ame

are the north atlantic oscillation and the northern annular mode distinguishable this study addresses the question of whether persistent events of the north atlantic oscillation nao and the northern annular mode nam teleconnection patterns are distinguishable from each other standard daily index time series are used to specify the amplitude of the nao and nam patterns the above question is examined with composites of sea level pressure and 300 and 40 hpa streamfunction along with tests of field significance a null hypothesis is specified that the nao and nam persistent events are indistinguishable this null hypothesis is evaluated by calculating the difference between time averaged nao and nam composites it is found that the null hypothesis cannot be rejected even at the 80 confidence level the wave breaking characteristics during the nam life cycle are also examined both the positive and negative nam phases yield the same wave breaking properties as those for the nao the results suggest that not only are the nao and nam persistent events indistinguishable but that the nao nam events are neither confined to the north atlantic nor are they annular

exhumation history of a section of the sikkim himalayas india records in the metamorphic mineral equilibria and compositional zoning of garnet the exhumation history of pelitic migmatite samples from the high himalayan crystalline complex hhc near the south tibetian detachment system in the sikkim darjeeling section has been determined on the basis of thermo barometric analyses retrograde breakdown reactions and compositional zoning of garnet the peak metamorphic condition is estimated to be 10 4 kbar 800 c from thermo barometric and phase equilibrium constraints the observed retrograde breakdown of garnet to spinel and cordierite requires near isothermal and hence extremely rapid 15 mm yr exhumation up to the depth of 15 km numerical modeling suggests that the initial rapid exhumation must have been followed by a much slower process 2 mm yr up to at least 5 km depth to lead to the development of the observed compositional zoning of garnet the dramatic change of the exhumation velocity vz might reflect a process of tectonic thinning followed by erosion and or horizontal flow at shallow depth assuming that the vz did not change significantly below 5 km depth these results suggest that the hhc samples studied in this work exhumed from a depth of 34 km within 8 ma

ree fractionation and nd isotope disequilibrium during crustal anatexis constraints from himalayan leucogranites abstract the processes controlling rare earth element ree behaviour during crustal anatexis have profound implications not only for ree abundances of crustally derived granites but also for their model nd ages for typical protoliths of anatectic granites such as are exposed in metapelites from the high himalayan crystalline series hhcs a high proportion 90 of bulk rock lree and mree are sited in the accessory phases apatite and monazite and of hree in garnet and apatite mass balance constraints indicate that observed ree abundances of himalayan leucogranites cannot be generated from ree deficient reactant phases in their metasedimentary protoliths muscovite biotite plagioclase and quartz but must result from the dissolution of accessory phases ree and th abundances in anatectic melts are modelled by combining empirical expressions for monazite apatite and zircon solubilities as a function of temperature melt composition and melt water content with analysed accessory phase ree abundances this approach which allows the contribution to the melt from each accessory phase to be identified successfully predicts himalayan leucogranite trace element abundances the lree and th concentrations of a metapelite derived melt are buffered by monazite stability although apatite dissolution contributes up to 10 of lree to the melt together with 50 of the mree and 90 of the hree zircon contributes a negligible proportion of ree to the melt during anatexis dissolution of monazite and apatite in the melt results in 1 increasing ree abundances and lree hree ratios with increasing temperature of melting 2 a negative eu anomaly in the melt 3 high sm nd ratios in the melt relative to its metasedimentary protolith and 4 nd isotope disequilibrium between the melt and its protolith all of these features are observed in the ree geochemistry of himalayan leucogranites the combined effects of sm nd fractionation and nd isotope disequilibrium during anatexis imply that model nd ages calculated for anatectic granites derived from apatite and monazite bearing protoliths should be treated with considerable caution

the generation of global rotational flow by steady idealized tropical divergence abstract tropical convective heating is balanced on the large scale by the adiabatic cooling of ascent the horizontal divergence of the wind above this heating may be viewed as driving the upper tropospheric rotational wind field a vorticity equation model is used to diagnose this relationship it is shown that because of the advection of vorticity by the divergent component of the flow the rossby wave source can be very different from the simple fd source often used in particular an equatorial region of divergence situated in easterly winds can lead to a rossby wave source in the subtropical westerlies where it is extremely effective this part of the source can be relatively insensitive to the longitudinal position of the equatorial divergence a divergence field which is asymmetric about the equator can lead to a quite symmetric rossby wave source for a steady frictionless flow the rossby wave source averaged over regions within closed streamfunction or absolute vorticity contours is under cert

reconsideration of the true versus apparent arctic oscillation abstract the physical reality of the arctic oscillation ao or northern annular mode is considered the data used are mainly the monthly mean sea level pressure slp a schematic figure is first presented to illustrate the relationship between the north atlantic oscillation nao pacific north american oscillation pna system and the ao negative correlation mode between the atlantic and the pacific ao ncm system although the nao pna apparent ao ncm and true ao ncm systems give rise to the same eofs the probability density functions for the time coefficients of the two leading modes are different therefore the discrimination of the two systems is possible several pieces of evidence indicate that in the real world the nao pna and the ao ncm are located on almost the same plane in phase space this means that the nao pna and ao ncm systems have the same variations on the plane in common implying that when the nao pna system is real the ao ncm is unlikely to be real simple independent compon

styles of post subduction collisional orogeny influence of convergence velocity crustal rheology and radiogenic heat production abstract a number of models of continental collision already exist but the role of foregoing oceanic plate subduction or variable coupling between plates still remains to be explored in addition heat generation by radioactive decay may be quite variable for example our geochemical data from low to high grade metamorphic rocks of the lesser and higher himalayan sequences reveal measured heat production values on the order 4 5 w m 3 which is considerably higher than normal inferred upper crustal values 3 therefore we performed 2d numerical experiments of oceanic continental subduction followed by continental collision for varying continental lower crustal strength convergence velocity and radiogenic heating rate oceanic subduction produces thinning and upwarping of the overriding lithospheric wedge such as is observed in andean type collisional settings cascadia and central andes after collision two main orogenic styles develop depending on the values of key parameters either a slow convergence rates weak shear heating and consequently strong rheological coupling between plates and or a weak lower crust generate a two sided thick and narrow symmetrical collisional zone or b high convergence rates strong shear heating and therefore weak rheological coupling between plates and a strong lower crust produce a one sided thin and broad asymmetrical collisional zone in this case a vertical lithospheric body divides two upper crustal wedges that develop distinct deformation patterns in both cases enhanced radiogenic heating favors the formation of midcrustal partially melted channels that propagate toward the surface and the pro foreland partial melting triggers lateral extrusion of metamorphic rocks from different depths 15 23 kbar 700 c on top of coherent lp lt units the geometry of metamorphic isograds inverted shapes of p t paths at different grades and exhumation rates are similar to those observed in parts of the himalayan chain indicating that the models capture the essence of real systems the model results demonstrate that internal heat production rate is a crucial parameter in determining the timing of slab breakoff polarity reversal and exhumation of the inverted metamorphic sequence with respect to the timing of partial melting late stages of some numerical experiments showed interesting geodynamical evolution of the collisional zone with the appearance of polarity reversal roll over of both plates and roll back of the subducting slab irrespective of these details hundreds of km of continental lithosphere are recycled into the mantle in all models supporting the idea that continental crustal consumption is indeed a common and considerable process during collision this affects inferences about extents of convergence and or shortening during orogenesis moreover crustal recycling into the deep mantle during collisional orogeny needs to be considered in models of crustal growth and evolution

the arctic oscillation signature in the wintertime geopotential height and temperature fields the leading empirical orthogonal function of the wintertime sea level pressure field is more strongly coupled to surface air temperature fluctuations over the eurasian continent than the north atlantic oscillation nao it resembles the nao in many respects but its primary center of action covers more of the arctic giving it a more zonally symmetric appearance coupled to strong fluctuations at the 50 hpa level on the intraseasonal interannual and interdecadal time scales this arctic oscillation ao can be interpreted as the surface signature of modulations in the strength of the polar vortex aloft it is proposed that the zonally asymmetric surface air temperature and mid tropospheric circulation anomalies observed in association with the ao may be secondary baroclinic features induced by the land sea contrasts the same modal structure is mirrored in the pronounced trends in winter and springtime surface air temperature sea level pressure and 50 hpa height over the past 30 years parts of eurasia have warmed by as much as several k sea level pressure over parts of the arctic has fallen by 4 hpa and the core of the lower stratospheric polar vortex has cooled by several k these trends can be interpreted as the development of a systematic bias in one of the atmosphere s dominant naturally occurring modes of variability

rare earth diffusion kinetics in garnet experimental studies and applications abstract we determined the diffusion coefficient of sm in almandine garnet as function of temperature at 1 bar and fo 2 corresponding to that of wustite iron buffer and to a limited extent that of a few other selected rare earth elements in almandine and pyrope garnets both garnets were demonstrated to have metastably survived the diffusion annealing at conditions beyond their stability fields the experimental diffusion profiles were analyzed by secondary ion mass spectrometry and in addition by rutherford back scattering spectroscopy for two samples transmission electron microscopic study of an almandine crystal that was diffusion annealed did not reveal any near surface fast diffusion path using reasonable approximations and theoretical analysis of vacancy diffusion the experimental data were used to develop an expression of rare earth element ree diffusion coefficient in garnet as a function of temperature pressure fo 2 ionic radius and matrix composition calculation of the closure temperature for the sm nd decay system in almandine garnet in a metamorphic terrain shows very good agreement with that constrained independently modeling of the ree evolution in melt and residual garnet suggests that for dry melting condition the ree pattern in the melt should commonly conform closely to that expected for equilibrium melting however for much lower solidus temperatures that would prevail in the presence of a h 2 o co 2 fluid the concentration of light ree in the melt could be significantly lower than that under equilibrium melting condition a reported core and rim differences in the ree content of a garnet crystal in a mantle xenolith in kimberlite have been reproduced by assuming that the ree zoning was a consequence of entrapment in a magma derived from an external source for 32 000 yr before the eruption

the role of melt in the uplift and exhumation of orogenic belts abstract in some orogenic belts rapid pulses of uplift occur when melt is present detailed work on the relation of melt to structural and metamorphic features in several localities shows that hot metamorphic rocks and especially anatectic migmatite can be physically transported from the lower to middle crust on decollements the decollements form in the lower crust because even though the lower crust is structurally weak throughout due to thermally activated crystal plastic deformation mechanisms the crust is substantially weakened where melt forms or is intruded melt in the lower crust is present long enough to accomodate strain whereas in the upper crust it crystallizes too fast to be a factor in strain localization following crustal loading due to thrusting exhumation results in the formation of additional melt in the migmatites due to decompression melting support of this sequence of events includes 1 evidence for near isothermal decompression of high temperature metamorphic rocks with early high pressure minerals having formed pre or synchronous with penetrative deformation and with lower pressure minerals having formed after the deformation 2 the occurrence of magma sheets that occur along decollements between high temperature metamorphic rocks and lower temperature metamorphic rocks in either thrust or normal configurations 3 nearly concordant k ar cooling dates for biotite and hornblende suggesting rapid cooling rates and 4 post tectonic intrusions of leucogranite or leucotonalite that may represent the products of decompression melting the crystalline slabs of the high himalaya the compressional surges of the coast orogen of british columbia and alaska the cordillera of southern british columbia and the southern belt of the hercynian chain are examples where melt assisted rapid pulses of uplift appear to have occurred during the terminal phases of uplift and following uplift these orogenic belts were exhumed during extension while melt was still present

the influences of the atlantic multidecadal oscillation on the mean strength of the north pacific subtropical high during boreal winter abstractthe atlantic multidecadal oscillation amo has been shown to be capable of exerting significant influences on the pacific climate in this study the authors analyze reanalysis datasets and conduct forced and coupled experiments with an atmospheric general circulation model agcm to explain why the winter north pacific subtropical high strengthens and expands northwestward during the positive phase of the amo the results show that the tropical atlantic warming associated with the positive amo phase leads to a westward displacement of the pacific walker circulation and a cooling of the tropical pacific ocean thereby inducing anomalous descending motion over the central tropical pacific the descending motion then excites a stationary rossby wave pattern that extends northward to produce a nearly barotropic anticyclone over the north pacific a diagnosis based on the quasigeostrophic vertical velocity equation reveals that the stationary wave pattern also results in enhanced subsidence over the

accelerated antarctic ice loss from satellite gravity measurements accurate quantification of antarctic ice sheet mass balance and its contribution to global sea level rise remains challenging gravity recovery and climate experiment data spanning the period april 2002 to january 2009 confirm earlier estimates of ice loss for antarctica and indicate that east antarctica started losing mass in about 2006

diffusion coupling between trace and major elements and a model for calculation of magma residence times using plagioclase residence times of plagioclase crystals in magma reservoirs can be determined by modeling the compositional zoning of trace elements in these crystals we present a formulation to model diffusion of trace elements in plagioclase paying special attention to certain thermodynamic and kinetic aspects in particular we account for the compositional dependence on anorthite content of the chemical potential and diffusion coefficients of trace elements e g mg the choice of suitable boundary conditions and potential effects of diffusion in more than one dimension we show that contrary to intuition diffusive fluxes of trace elements may be coupled to major element concentration gradients and ignoring such coupling can lead to incorrect estimates of timescales we illustrate application of the model using plagioclase crystals from a suite of gabbroic xenoliths from a holocene dacitic lava flow of volcan san pedro chilean andes 36 s the inferred timescale for metasomatism of the xenoliths by evolved liquids is on the order of 100 30 to 148 yr and serves to illustrate how trace element zoning in plagioclase provides a window into timescales of magmatic processes inaccessible by isotopic or other methods