青藏高原南部拉萨地体的变质作用与动力学

拉萨地体位于欧亚板块的最南缘,它在新生代与印度大陆的碰撞形成了青藏高原和喜马拉雅造山带。因此,拉萨地体是揭示青藏高原形成与演化历史的关键之一。拉萨地体中的中、高级变质岩以前被认为是拉萨地体的前寒武纪变质基底。但新近的研究表明,拉萨地体经历了多期和不同类型的变质作用,包括在洋壳俯冲构造体制下发生的新元古代和晚古生代高压变质作用,在陆-陆碰撞环境下发生的早古生代和早中生代中压型变质作用,在洋中脊俯冲过程中发生的晚白垩纪高温/中压变质作用,以及在大陆俯冲带上盘加厚大陆地壳深部发生的两期新生代中压型变质作用。这些变质作用和伴生的岩浆作用表明,拉萨地体经历了从新元古代至新生代的复杂演化过程。(1)北拉萨地体的结晶基底包括新元古代的洋壳岩石,它们很可能是在Rodinia超大陆裂解过程中形成的莫桑比克洋的残余。(2)随着莫桑比克洋的俯冲和东、西冈瓦纳大陆的汇聚,拉萨地体洋壳基底经历了晚新元古代的(～650Ma)的高压变质作用和早古代的(～485Ma)中压型变质作用。这很可能表明北拉萨地体起源于东非造山带的北端。(3)在古特提斯洋向冈瓦纳大陆北缘的俯冲过程中,拉萨地体和羌塘地体经历了中古生代的(～360Ma)岩浆作用。(4)古特提斯洋盆的闭合和南、北拉萨地体的碰撞,导致了晚二叠纪(～260Ma)高压变质带和三叠纪(～220Ma)中压变质带的形成。(5)在新特提斯洋中脊向北的俯冲过程中,拉萨地体经历了晚白垩纪(～90Ma)安第斯型造山作用,形成了高温/中压型变质带和高温的紫苏花岗岩。(6)在早新生代(55～45Ma),印度与欧亚板块的碰撞,导致拉萨地体地壳加厚,形成了中压角闪岩相变质作用和同碰撞岩浆作用。(7)在晚始新世(40～30Ma),随着大陆的继续汇聚,南拉萨地体经历了另一期角闪岩相至麻粒岩相变质作用和深熔作用。拉萨地体的构造演化过程是研究汇聚板块边缘变质作用与动力学的最佳实例。     

拉萨地体东南部林芝杂岩形成与变质演化的锆石U-Pb年代学限定

本文对拉萨地体东南部林芝地区分布的变质岩进行了岩相学和锆石年代学研究。结果表明,林芝杂岩中的变质沉积岩主要由片麻岩和片岩组成,它们经历了中压角闪岩相变质作用。变质岩中的锆石多由继承的碎屑岩浆核和新生的变质边组成。继承锆石核给出了新太古代至晚古生代的年龄范围,其主要年龄峰值在～1560Ma、～1190Ma、～620Ma和～340Ma,而锆石变质边给出了53Ma和27Ma的变质年龄。这一结果表明,林芝杂岩中的变质沉积岩很可能形成在古生代,其物质源区具有Grenville和Pan-Africa期造山作用的构造热事件记录。这一研究和已有的成果进一步证明,拉萨地体起源于Gondwana大陆北缘,在新生代印度与欧亚大陆的碰撞/俯冲过程中,拉萨地体作为俯冲带的上盘经历了多期变质作用改造。本研究为拉萨地体起源与多期构造演化提供了重要信息。     

拉萨地体东部早侏罗纪变质和深熔作用

位于青藏高原南部的拉萨地体不仅记录了中生代的新特提斯洋俯冲及随后新生代的印度-欧亚板块陆陆碰撞造山作用,而且还记录了晚古生代-早中生代南、北拉萨地体的拼合作用。本文对拉萨地体东部东久地区的片岩和脉体进行了岩石学和锆石U-Pb年代学研究,表明片岩经历了峰期高角闪岩相的变质作用和部分熔融,中压角闪岩相退变质过程以及晚期的降温、降压过程。片岩记录了峰期矿物组合蓝晶石+石榴石+黑云母+斜长石+钛铁矿+石英,退变质矿物组合石榴石+夕线石+堇青石+黑云母+斜长石+钛铁矿+石英,晚期退变质矿物组合堇青石+黑云母+白云母+绿泥石+斜长石+钛铁矿+石英。相平衡模拟研究表明,片岩的峰期变质作用温度、压力条件约为720℃、0.9GPa;退变质条件约为670℃、0.59GPa以及480℃、0.12GPa。全岩地球化学研究表明,含石榴石长英质脉体具有显著的Eu元素正异常(δEu=3.57),为斜长石堆晶的产物。锆石U-Pb年代学表明,片岩和脉体在早侏罗纪的181 Ma和195 Ma发生了变质和部分熔融作用。本文结合已发表研究结果表明,东久地区的高级变质岩可划分出不同的构造岩片,在早侏罗纪先后经历了相似温、压条件的变质作用,为南、北拉萨地体碰撞造山作用的产物。     

拉萨地体东南缘始新世早期变质作用及其构造意义

本文对位于青藏高原拉萨地体东南缘林芝杂岩中的片麻岩进行了岩石学和锆石U-Pb年代学研究.所研究的样品包括正片麻岩和副片麻岩,它们经历了中压角闪岩相变质作用.岩石地球化学分析结果表明,所研究的正片麻岩的原岩具有钙碱性岛弧岩浆岩的特征.锆石U-Pb年代学分析结果表明,副片麻岩中的碎屑锆石核部为岩浆成因,它们给出的206Pb/238U年龄范围为3012～ 522Ma,其锆石的增生边给出了～51Ma的变质年龄.在正片麻岩中,黑云母片麻岩给出了～67Ma的原岩结晶年龄和～ 55 Ma的变质年龄;石榴石角闪黑云斜长片麻岩给出了～58Ma的原岩结晶年龄和～54Ma的变质年龄.因此,所研究的林芝杂岩并不能代表拉萨地体中的前寒武纪变质基底,而是古生代的沉积岩和晚白垩纪至早新生代的岩浆岩在始新世早期变质而成.这一时期,表壳岩和侵入岩一起经历的中压角闪岩相变质作用很可能跟新特提斯洋俯冲导致的地壳增生、加厚有关.     

喜马拉雅造山带的高压超高压变质作用与印度-亚洲大陆碰撞

喜马拉雅造山带是印度与亚洲大陆碰撞作用的产物,正在进行造山作用,是研究板块构造的天然实验室.高压和超高压变质岩分布在喜马拉雅造山带的核部.这些变质岩具有不同的形成条件、形成时间和形成过程,为印度与亚洲碰撞带的几何学、运动学和动力学提供了重要的限定.含柯石英的超高压变质岩产出在喜马拉雅造山带的西段,它们形成在古新世与始新世之间(53～46Ma),为印度大陆西北边缘高角度超深俯冲作用的产物,并经历了快速俯冲与快速折返过程.在约5 Myr内,超高压变质岩从＞100km的地幔深度折返到了中地壳深度,且仅仅叠加角闪岩相退变质作用.高压榴辉岩产出在喜马拉雅造山带中段,形成时间约为45Ma,为印度大陆低角度深俯冲作用的产物,经历了至少20Myr的长期折返过程,叠加麻粒岩相退变质作用和部分熔融.高压麻粒岩产出在喜马拉雅造山带的东端,是印度大陆东北缘近平俯冲作用的产物,峰期变质作用时间约为35Ma,经历了约20Myr的长期折返过程,叠加了麻粒岩相和角闪岩相退变质作用,并伴随有多期部分熔融.因此,喜马拉雅造山带的变质作用具有明显的时间与空间变化,显示出大陆深俯冲与折返过程的差异性,以及大陆碰撞造山带形成机制的多样性.     

龟山杂岩的形成时代、构造属性及变质条件:对秦岭-桐柏造山带构造演化的指示

秦岭-桐柏-大别造山带是华北板块和扬子板块碰撞拼合的产物.桐柏造山带作为古生代秦岭增生造山带与中生代大别碰撞造山带之间的纽带,记录了两大板块之间从洋壳俯冲、地体增生到最终的碰撞造山过程.本文选取桐柏造山带龟山杂岩为研究对象,对其中变沉积岩、斜长角闪岩和花岗质脉体进行了岩相学、矿物学和锆石U-Pb年代学研究.结果指示变沉积岩碎屑锆石年龄峰值主要为～411 Ma,次要年龄峰值为～490 Ma,最年轻四颗碎屑锆石的加权平均年龄为352±11 Ma.斜长角闪岩原岩结晶时代为498±18 Ma,花岗质脉形成时代为412.8±4.4 Ma,两者中变质锆石加权平均年龄为320±10 Ma,指示它们共同经历了早石炭世变质作用.变质峰期温压条件为T=571～632℃,P=0.59～0.76 GPa,属于中压相系角闪岩相变质.区域年代对比分析显示,龟山杂岩变沉积岩物源来自于北侧的北秦岭地体,属于弧前沉积,沉积时代下限为早石炭世早期.结合前人研究,本文认为在古秦岭洋向北的俯冲过程中,上述沉积岩与斜长角闪岩因构造作用拼贴在一起,于～320 Ma发生中压角闪岩相变质作用.该研究结果也指示了古秦岭洋最终闭合时间可能晚于早石炭世早期.     

西藏拉萨地体北部的前寒武纪高压变质作用及构造意义

高压基性麻粒岩以石榴石-单斜辉石-斜长石-石英共生为特征,是研究大洋或大陆地壳俯冲-碰撞的最好载体之一.西藏拉萨地体北部高压基性麻粒岩以构造块体的形式出现在早古生代沉积岩中.高压麻粒岩的原岩是辉长岩,它经历了四期变质作用,从早到晚分别是:角闪岩相变质(0.9～1.0GPa和710～720℃)、高压麻粒岩相峰期(1.55～1.65GPa和730～740℃)、麻粒岩相(0.82GPa和821℃)和角闪岩相退变质(0.60～0.68GPa和520～540℃)作用.整个变质作用的P-T轨迹是顺时针型的,包括一个近等压降温的早期角闪岩相变质过程,近等温升压的晚期进变质过程,以及加温降压的早期退变质过程和降温降压的晚期退变质过程.这表明,高压基性麻粒岩形成在较高地热梯度条件下,并且经历了加热变质过程.因此,纳木错高压基性麻粒岩并不是形成在典型的洋壳俯冲带构造环境,洋中脊俯冲和地幔柱作用是其成因的可能构造控制因素.     

青藏高原拉萨地体北部的前寒武纪变质作用及构造意义

青藏高原南部拉萨地体中分布的角闪岩相至麻粒岩相变质岩一直被认为是前寒武纪变质基底,但并没有获得可靠的年代学证据。本文运用原位锆石U-Pb定年方法,在拉萨地体北部那果地区的变质岩中获得了约720Ma的变质年龄,从而证明拉萨地体北部在新古元代经历了角闪岩相变质作用和近同期的岩浆作用。基于这一成果和在拉萨地体中、南部高级变质岩中陆续获得的中、新生代变质年龄,对拉萨地体变质作用的时、空变化及其成因进行了初步探讨。     

北秦岭造山带的早古生代多期变质作用

北秦岭造山带的秦岭岩群以高级变质岩石为特征,主要包括少量榴辉岩、高压麻粒岩和区域上广泛分布的麻粒岩-角闪岩相变质岩石。年代学研究显示秦岭岩群中不同岩石记录了多期变质作用。已有的定年资料给出北秦岭官坡地区的榴辉岩的年龄为500Ma左右,代表榴辉岩相的变质时代。结合岩相学资料,对两个高压麻粒岩样品的SHRIMP和LA-ICPMS U-Pb测定分别获得504±7Ma 和506±3Ma的年龄,应代表高压麻粒岩相变质时代。这表明高压麻粒岩和相邻的榴辉岩有相近的变质时代,但形成在造山带中不同的构热造环境中。西峡地区的角闪二辉麻粒岩的U-Pb定年给出两组早古生代年龄,一组为440±2Ma,可能代表了中低压麻粒岩相的变质时代,另一组为426±1Ma,应代表区域角闪岩相的变质时代。桐柏山北部的石榴二辉麻粒岩的U-Pb定年数据给出436±1Ma的年龄,为中压麻粒岩相的变质时代。这些资料表明北秦岭造山带经历了早奥陶世的俯冲和地壳增厚作用,并在晚志留世遭受了广泛的巴罗式区域变质作用。     

造山带中成对出现的高压麻粒岩与榴辉岩及其地球动力学意义

在一些典型碰撞造山带中,高压麻粒岩与榴辉岩在空间和时间上密切相关,它们之间的关系对揭示碰撞造山带的造山过程和造山机制具有重要意义.本文以中国西部的南阿尔金、柴北缘及中部的北秦岭造山带为例,详细陈述了这3个地区榴辉岩和相关的高压麻粒岩的野外关系、变质演化和形成时代,目的是要建立大陆碰撞造山带中榴辉岩和相关高压麻粒岩形成的地球动力学背景模式.南阿尔金榴辉岩呈近东西向分布在江尕勒萨依,玉石矿沟一带,与含夕线石副片麻岩、花岗质片麻岩和少量大理岩构成榴辉岩一片麻岩单元,榴辉岩中含有柯石英假象,其峰期变质条件为P=2.8～3.0GPa,T=730～850℃,并在抬升过程中经历了角闪岩-麻粒岩相的叠加;大量年代学研究显示其峰期变质时代为485～500Ma.南阿尔金高压麻粒岩分布在巴什瓦克地区,包括高压基性麻粒岩和高压长英质麻粒岩,它们与超基性岩构成了一个大约5km宽的构造岩石单元,与周围角闪岩相的片麻岩为韧性剪切带接触.长英质麻粒岩和基性麻粒岩的峰期组合均具有蓝晶石和三元长石(已变成条纹长石),形成的温压条件为T=930～1020℃,P=1.8～2.5GPa,并在退变质过程中经历了中压麻粒岩相变质作用叠加.锆石SHRIMP测定显示巴什瓦克高压麻粒岩的峰期变质时代为493～497Ma.都兰地区的榴辉岩分布柴北缘HP-UHP变质带的东端,在榴辉岩和围岩副片麻岩中均发现有柯石英保存,形成的峰期温压条件为T=670～730℃和P=2.7～3.25GPa,退变质阶段经过了角闪岩相的叠加;榴辉岩相变质时代为420～450Mao都兰地区的高压麻粒岩分布在阿尔茨托山西部,高压麻粒岩包括基性麻粒岩长英质麻粒岩,基性麻粒岩的峰期矿物组合为Grt+Cpx+Pl±Ky±Zo+Rt±Qtz,长英质麻粒岩的峰期矿物组合为：Grt+Kf+Ky+Pl+Qtz.峰期变质条件为T=800～925℃,P=1.4～1.85GPa,退变质阶段经历了角闪岩-绿片岩的改造,高压麻粒岩的变质时代为420～450Ma.北秦岭榴辉岩分布在官坡-双槐树一带,榴辉岩的峰期变质组合为Grt+Omp±Phe+Qtz+Rt,所计算的峰期温压条件为T=680～770℃和P=2.25～2.65GPa,年代学数据显示榴辉岩的变质时代为500Ma左右.北秦岭高压麻粒岩分布在含榴辉岩单元的南侧松树沟一带,包括高压基性麻粒岩和高压长英质麻粒岩,与超基性岩在空间上密切伴生,高压麻粒岩的峰期温压条件为T=850～925℃,P=1.45～1.80GPa,锆石U-Pb年代学研究显示其峰期变质时代为485～507Ma.以上三个实例显示,出现在同一造山带、在空间上伴生的高压麻粒岩和榴辉岩有各自不同的变质演化历史,但榴辉岩中的榴辉岩相变质时代和相邻的高压麻粒岩中的高压麻粒岩相变质作用时代相同或相近,这种成对出现的榴辉岩和高压麻粒岩代表了它们同时形成在造山带中不同的构造环境中,即榴辉岩的形成于大陆俯冲带中,而高压麻粒岩可能形成在俯冲带之上增厚的大陆地壳根部.     

摩擦桩基桩土间极限摩阻力取值问题探讨

通过对广珠东线高速公路横沥大桥的试桩及土体的工程地质条件分析,总结出影响摩擦桩基桩土间极限摩阻力取值的一般问题以及解决问题的方法和措施。     

从榴辉岩与围岩的关系论苏鲁榴辉岩的形成与折返

位于华北和扬子两板块碰撞带中的苏鲁榴辉岩形成的温压条件不但是超高压，而且是高温。榴辉岩的PTt轨迹表明其为陆-陆磁撞俯冲带的产物。榴辉岩的区域性围岩花岗质片麻岩为新元古代同碰撞期花岗岩，榴辉岩及其他直接围岩皆呈包体存在于其中，并见新元古代花岗岩呈脉状侵入榴辉岩包体中。区域性围岩新元古代花岗岩的锆石中发现有柯石英、绿辉石等包裹体，表明新元古代花岗岩的组成物质也经受过超高压变质作用，且榴辉岩与围岩新元古代花岗岩的锆石U-Pb体系同位素年龄基本相同。但新元古代花岗岩所记录的变质作用和变形作用期次（或阶段）却少于榴辉岩。椐上述可得如下推断：超高压榴辉岩与新元古代花岗岩岩浆是同时在碰撞带底部（俯冲板块前部）形成的；榴辉岩的第一折返阶段是由新元古代花岗岩岩浆携带上升的，其第二折返阶段是和新元古代花岗岩一起由逆冲及区域性隆起而上升，遭受剥蚀。     

The Importance of the Precipitation and the Susceptibility of the Slopes for the Triggering of Landslides Along the Roads

In order to characterise the influence of the heavyrains on the observed landslides during the 1996–1997hydrological cycle, rainfall records for the last 100years are analysed from 104 stations in easternAndalusia. Regarding the amounts of rain recordedbetween October 1996 and March 1997 in the 104stations studied, 31 presented new all-time records;15 presented values that were 80–100% of thepre-1995 record; 49 stations, 80–50%; and 9stations, < 50%. A map has been devised of thesusceptibility of the materials through which thesouth-eastern Andalusian road network crosses,together with an inventory of the damage caused byinstability phenomena on banks and cuttings of theroad network during the winter of 1996–1997. Therelationships between the rainfall during the studyperiod, the damage caused to the road network and thesusceptibility of the materials affected are analysed.The results indicate that there is a clearcorrespondence between the rainfall recorded and thesusceptibility of the materials with the inventorieddamage. It is concluded that the widespread seriousdamage caused in early 1997 to the roads andsurrounding areas in the Alpujarra region and thecoast of the Province of Granada was mainly caused bythe extraordinarily heavy rains. However, considerablyless damage was observed where the susceptibility ofthe terrain is low, thus highlighting the extremeusefulness of terrain-susceptibility maps for riskprevention and territorial planning.     

某高速公路下伏煤矿采空区稳定性分析

在论述某高速公路下伏砦脖煤矿采空区地质、采矿和工程地质特征的基础上, 进行了稳定性数值模拟分析, 定性与定量地分析与评价了该煤矿采空区的地表变形特征及稳定性。研究结果表明: 该煤矿采空区的变形尚未完成, 对拟建的高速公路将产生很大的危害, 必须采取相应的工程治理措施。      

混凝剂处理钻井废泥浆液的研究

通过烧杯搅拌实验对混凝剂处理钻井废泥浆液进行研究,从混凝剂适应p H值范围、混凝效果、沉降速度三方面研究比较,找出一种较为理想的混凝剂,并分析了影响混凝剂性能的主要因素,确定了混凝剂的最佳投放剂量。      

黄河源区水环境变化及黄河出现冬季断流的原因

自1954年有水文观测资料以来，黄河曾在青海省玛多县黄河沿水文站发生过3次断流。本文在分析黄河源区水环境特征及其影响因素的基础上指出，鄂陵湖、扎陵湖的环湖融区调节能力低，当遇到连续干旱、冬季其调节水量不足以维系黄河径流时便会发生断流，这是断流的主因。湖水位降低、开采沙金、过度放牧等自然和人为因素也会对黄河发生断流产生影响。鄂陵湖口附近黄河上修建的水电站开始蓄水，提高了两湖及环湖融区的调节能力，今后黄河冬季出现断流的可能性将大为降低。     

Unity of the national and the international in the socio-cultural development of the peoples of the USSR

International unity is becoming ever stronger in this country owing to an increasing similarity in the development of the cultural environment. This comprises the provision of all the country's republics with a sufficient number of schools, theatres, and other institutions and cultural information media in accordance with the needs of the population. An important part is played by the rise in ‘the general educational level, as well as the level of professional qualifications and skills. Among all the Soviet nations and nationalities, this rise being more rapid among formerly backward peoples. Prominent among the factors of internationalization is the progressive development of the nationalities’ cultural resources, while professional culture is being increasingly brought within the reach of the masses.The implementation of the nationalities policy promotes the all-round development of all Soviet nations and nationalities, their drawing together, the upsurge of the individual capabilities of every Soviet citizen.     

Artemia in the Salt Lakes of Russia:the Productivity of Populations,the Reserves of the Cysts and the Fisheries

正Artemia cysts are an extremely important component of aquaculture diets.It is well established that the cultivation of fish and shellfish derive substantial health and growth advantages when Artemia are included in the diets of the     

天山南坡科其喀尔冰川表碛区小气候特征研究

利用天山南坡科其喀尔冰川3号观测站2009年全年的气象观测资料,分析研究了科其喀尔冰川表碛区的小气候特征. 结果表明：总辐射和净辐射夏秋季较高、冬春季较低;反射辐射和地表反照率反之. 与其他地区不同,该区主要受积雪物理性质和下垫面状况的影响,冬春季地表反照率日变化表现为由大到小的变化过程,夏秋季表现为倒U型. 温度年变化表现为夏秋季高、冬春季低,最高月均值出现在8月,为9.4℃,最低月均值出现在1月,为-9.6℃. 受山谷风和冰川风的影响,全年的风向以西北风和西北偏西风为主,风向的日变化以11：00为界发生转向. 受降水和冰川消融等的影响,比湿夏秋季月均值较大,冬春季月均值较小.     

Climate: Is the past the key to the future?

 The climate of the Holocene is not well suited to be the baseline for the climate of the planet. It is an interglacial, a state typical of only 10% of the past few million years. It is a time of relative sea-level stability after a rapid 130-m rise from the lowstand during the last glacial maximum. Physical geologic processes are operating at unusual rates and much of the geochemical system is not in a steady state. During most of the Phanerozoic there have been no continental ice sheets on the earth, and the planet’s meridional temperature gradient has been much less than it is presently. Major factors influencing climate are insolation, greenhouse gases, paleogeography, and vegetation; the first two are discussed in this paper. Changes in the earth’s orbital parameters affect the amount of radiation received from the sun at different latitudes over the course of the year. During the last climate cycle, the waxing and waning of the northern hemisphere continental ice sheets closely followed the changes in summer insolation at the latitude of the northern hemisphere polar circle. The overall intensity of insolation in the northern hemisphere is governed by the precession of the earth’s axis of rotation, and the precession and ellipticity of the earth’s orbit. At the polar circle a meridional minimum of summer insolation becomes alternately more and less pronounced as the obliquity of the earth’s axis of rotation changes. Feedback processes amplify the insolation signal. Greenhouse gases (H₂O, CO₂, CH₄, CFCs) modulate the insolation-driven climate. The atmospheric content of CO₂ during the last glacial maximum was approximately 30% less than during the present interglacial. A variety of possible causes for this change have been postulated. The present burning of fossil fuels, deforestation, and cement manufacture since the beginning of the industrial revolution have added CO₂ to the atmosphere when its content due to glacial-interglacial variation was already at a maximum. Anthropogenic activity has increased the CO₂ content of the atmosphere to 130% of its previous Holocene level, probably higher than at any time during the past few million years. During the Late Cretaceous the atmospheric CO₂ content was probably about four times that of the present, the level to which it may rise at the end of the next century. The results of a Campanian (80 Ma) climate simulation suggest that the positive feedback between CO₂ and another important greenhouse gas, H₂O, raised the earth’s temperature to a level where latent heat transport became much more significant than it is presently, and operated efficiently at all latitudes. Atmospheric high- and low-pressure systems were as much the result of variations in the vapor content of the air as of temperature differences. In our present state of knowledge, future climate change is unpredictable because by adding CO₂ to the atmosphere we are forcing the climate toward a “greenhouse” mode when it is accustomed to moving between the glacial–interglacial “icehouse” states that reflect the waxing and waning of ice sheets. At the same time we are replacing freely transpiring C3 plants with water-conserving C4 plants, producing a global vegetation complex that has no past analog. The past climates of the earth cannot be used as a direct guide to what may occur in the future. To understand what may happen in the future we must learn about the first principles of physics and chemistry related to the earth’s system. The fundamental mechanisms of the climate system are best explored in simulations of the earth’s ancient extreme climates. Received: 7 November 1996/Accepted: 23 January 1997