A comparison of PMIP2 model simulations and the MARGO proxy reconstruction for tropical sea surface temperatures at last glacial maximum

Results from multiple model simulations are used to understand the tropical sea surface temperature (SST) response to the reduced greenhouse gas concentrations and large continental ice sheets of the last glacial maximum (LGM). We present LGM simulations from the Paleoclimate Modelling Intercomparison Project, Phase 2 (PMIP2) and compare these simulations to proxy data collated and harmonized within the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface Project (MARGO). Five atmosphere–ocean coupled climate models (AOGCMs) and one coupled model of intermediate complexity have PMIP2 ocean results available for LGM. The models give a range of tropical (defined for this paper as 15°S–15°N) SST cooling of 1.0–2.4°C, comparable to the MARGO estimate of annual cooling of 1.7 ± 1°C. The models simulate greater SST cooling in the tropical Atlantic than tropical Pacific, but interbasin and intrabasin variations of cooling are much smaller than those found in the MARGO reconstruction. The simulated tropical coolings are relatively insensitive to season, a feature also present in the MARGO transferred-based estimates calculated from planktonic foraminiferal assemblages for the Indian and Pacific Oceans. These assemblages indicate seasonality in cooling in the Atlantic basin, with greater cooling in northern summer than northern winter, not captured by the model simulations. Biases in the simulations of the tropical upwelling and thermocline found in the preindustrial control simulations remain for the LGM simulations and are partly responsible for the more homogeneous spatial and temporal LGM tropical cooling simulated by the models. The PMIP2 LGM simulations give estimates for the climate sensitivity parameter of 0.67°–0.83°C per Wm⁻², which translates to equilibrium climate sensitivity for doubling of atmospheric CO₂ of 2.6–3.1°C.     

ITCZ rather than ENSO signature for abrupt climate changes across the tropical Pacific?

Latitudinal movements of the Intertropical Convergence Zone (ITCZ), analogous to its present-day seasonal shifts, and El Niño Southern Oscillation (ENSO)-type variability both potentially impacted rainfall changes at the millennial timescale during the last glacial period. In this study we compare tropical Pacific sedimentary records of paleoprecipitation to decipher which climate mechanism was responsible for the past rainfall changes. We find that latitudinal movements of the ITCZ are consistent with the observed rainfall patterns, challenging the ENSO hypothesis for explaining the rapid rainfall changes at low latitudes. The ITCZ-related mechanism appears to reflect large-scale atmospheric rearrangements over the tropical belt, with a pronounced Heinrich-Dansgaard/Oeschger signature. This observation is coherent with the simulated tropical rainfall anomalies induced by a weakening of the Atlantic thermohaline circulation in modeling experiments.     

Evaluation of site effects by means of Joint Analysis of Sonogram and Standard Spectral Ratio (JASSSR)

In order to evaluate thefrequency dependence of site effects, thejoint analysis of the standard spectralratio technique with a sonogram approach isproposed (Joint Analysis of Sonogram andStandard Spectral Ratio). Two tools aresuggested: the enhanced sonograms andthe differential sonograms. Inparticular, the local geology effects areinvestigated by considering both theamplification of the seismic signals andtheir prolongation. Tests carried out onsynthetic signals show the capability ofthe sonogram approach to evaluate thefrequency-dependent lengthening of theseismic signal. The JASSSR is applied todata recorded by a dense local networkdeployed during the Umbria-Marche seismicsequence of 1997 and to data recorded inthe Volvi Test site (Euroseistest, Greece).The results show that using the JASSSR gives amore complete overview of the effects oflocal geology on seismic ground motion than maybe obtained by using the StandardSpectral Ratio alone. In particular, theeffect of secondary arrivals, both on theamplitude and on the lengthening of theseismic ground motion, can be pointed out.In addition, the characteristic frequencyof these arrivals can be evaluated.     

Seismic Microzonation: A Comparison between Geotechnical and Seismological Approaches in Pointe-à-Pitre (French West Indies)

Pointe-à-Pitre, the main city of Guadeloupe in the French West Indies, has on several occasions been partially destroyed by major historical earthquakes. Moreover, a post-seismic assessment of the damage from the 1985 Montserrat earthquake indicates that the town is prone to site effects. Consequently, from 1996 to 1998, BRGM conducted a seismic microzonation study based on geotechnical and geological data. At the same time, three seismological studies were being conducted – two based on earthquake recordings using a time-series analysis and the classical spectral ratio (CSR) method (CETE/LCPC and BRGM), and the third based on noise measurement at 400 points using the horizontal-to-vertical noise ratio (HVNR) method (CETE/LCPC). The objective of this paper is not to carry out a new microzonation study by taking into account all the results, but rather to show in what respects the results of these different methods are in agreement or not. A comparison of the results of the seismological studies with the geotechnical microzonation shows that they are in fairly good agreement, albeit with some discrepancies. The results indicate that the seismological methods and the geotechnical data are highly complementary and should be used together in compiling seismic transfer-function microzonation maps. This revised version was published online in July 2006 with corrections to the Cover Date.     

Sea-level estimates during the last deglaciation based on δO and accelerator mass spectrometry C ages measured in Globigerina bulloides

Coupled measurements of δ¹⁸O and accelerator mass spectrometry (AMS) ¹⁴C in a particular species of planktonic foraminifera may be used to calculate sea-level estimates for the last deglaciation. Of critical importance for this type of study is a knowledge of the seasonality of foraminiferal growth, which can be provided by δ¹⁸O measurements of modern shells (core tops, plankton tows). Isotopic (δ¹⁸O, AMS-¹⁴C dating) and faunal records (transfer function sea surface temperature) were obtained from two cores in the North Atlantic at about 37°N. The locations were chosen to obtain high sedimentation rate records removed from the major ice-melt discharge areas of the last deglaciation. Based upon Globigerina bulloides data, four δ¹⁸O-based sea-level estimates were calculated: −67 ± 7 m at 12,200 yr B.P. and −24 ± 8 m at about 8200 yr B.P. for core SU 81-18; −83 ± 10 m at 12,200 yr B.P. and −13 ± 11 m at about 8500 yr B.P. for core SU 81-14. Using a second working hypothesis concerning the seasonability of G. bulloides growth, it is suggested that the sea-level rose by about 40 m during the millennium which followed 14,500 yr B.P.     

Using a maximum simplicity paleoclimate model to simulate millennial variability during the last four glacial periods

Many studies have documented the existence of millennial-scale variability in the Earth system during the last glacial period. An increasing number of studies document the occurrence of similar millennial variability during glacial periods previous to the last one. Here we use the simplest possible thermal-bipolar seesaw model to consider this variability for the last four glacial periods. We invert this model and use the high-pass filtered Vostok stable isotope records to make a first, tentative, attempt to estimate high-latitude N. Hemisphere temperature variability over the last four glacial periods, beyond the reach of Greenland ice-core records. The model result is compared against the Vostok methane record, which shows rapid variations in parallel to Greenland temperature records during the last glacial period. A further comparison is carried out against the planktonic oxygen isotope of north Atlantic core ODP 980. There is agreement between the records on the existence of similar millennial-scale variability during the last three glacial periods with very similar characteristics to the variability during the last glacial cycle.     

Experimental study of the ground motion on a large scale topographic hill at Kitherion (Greece)

We present the results of an experiment about the effects on the seismic ground motion of a hill of large dimensions (6km long, 3km large and 700m high) near Corinth (Greece). We installed 7 seismometers across the hill, and analyzed the ground motion with different methods: the classical spectral ratios (CSR) and the horizontal to vertical spectral ratios calculated both on noise (HVNR) and earthquake data (RF). The amplification at the resonant frequency (around 0.7Hz) is low, as expected, but one of the two stations installed at the top exhibits large amplification (up to a factor of 10) at 3Hz. The comparison between the H/V results and the CSR ones suggests that the former method is able to point out the fundamental frequencies of a hill. Amplifications and spatial localization predicted by numerical 3D modeling are consistent with the experimental data only at frequencies lower than 1Hz.     

Blind deconvolution methodology for site-response evaluation exclusively from ground-surface seismic recordings

A novel blind deconvolution methodology for identification of the local site characteristics based on two seismograms recorded on the free surface of a sediment site is presented. The approach does not require recordings at depth nor at a nearby rock outcrop, and eliminates the need for any prior parameterization of source and site characteristics. It considers that the surface recordings are the result of the convolution of the ‘input motion at depth' with transfer functions (channels) representing the characteristics of the transmission path of the waves from the input location to each recording station. The input motion at depth is considered to be the common component in the seismograms (same input in a statistical sense). The channel characteristics are considered to be the part in the seismograms that is non-common, since the travel path of the waves from the input motion location at depth to each recording station is different, due to spatially variable site effects. By means of blind deconvolution, the algorithm eliminates what is common in the seismograms, namely the input motion at depth, and retains what is different, namely the transfer functions of the site from the input location to each recording station. It estimates the site response in both frequency and time domains, and identifies the duration of the site's transfer functions. The methodology is applied herein to synthetic data at realistic sites for performance validation. The blindly estimated results are in almost perfect agreement with the actual site characteristics, indicating that the approach is a promising new tool for seismic site-response identification from recorded data.     

Quantification of soil non-linearity based on simulation

Recent earthquake studies have demonstrated that non-linear behavior of soft soil can be significant compared to other effects affecting seismic motion. Therefore, the question is to know when sediment non-linearity is a first-order effect and when it is not. In this study, we propose a method for quantifying non-linear effects based on simulations. An elasto-plastic model is used to simulate the behavior of four materials (sands and clays). For each computation, the non-linearity is quantified by the use of a ‘non-linear parameter’ and compared to four non-linearity indicators. These computations suggest that the efficiency of an indicator depends on the nature of the soil underlying the recorder and that the most efficient indicator should be based on the high frequencies content rather than on the resonant frequency changes.     

岩石圈伸展的壳/幔拆离模型(Parallel Extension Tectonics):华北克拉通东部早白垩世岩石圈减薄与破坏机理

华北克拉通岩石圈减薄和破坏机理长期以来存在争议,基于岩石学、岩石地球的化学分析研究突出强调深部过程的重要性。前人提出了两种重要模式：包括以拆沉作用为代表的top-down tectonics模型和以热-机械侵蚀与化学侵蚀,或地幔置换、交代作用的bottom-up tectonics模型。然而,对于这两种模式而言尚存在许多无法合理解释的问题,比如在此深部过程中,区域性岩石圈伸展有多大的贡献？地壳伸展构造是作为深部过程的响应,还是同为岩石圈伸展的产物？本文基于早白垩世东亚地区（尤其是华北克拉通东部地区）伸展构造与岩浆活动的综合分析,揭示出华北克拉通东部不同地区伸展构造变形与岩浆活动之间的时、空和成因关系有一定的差异。但整体上看,岩石圈伸展起着主导作用,控制着岩浆上侵和就位,在拆离断层下盘侵入形成各种规模的花岗质为主的侵入体,或于上盘喷发形成火山-沉积岩盆地。在伸展构造发育的不同阶段,可以有伸展早期、伸展期及伸展期后的岩浆活动。岩浆活动的强度及岩浆源区特点有显著的时空变化。一方面,在同一地区不同演化阶段其源区有很大的差异。表现为主体上是早期以古老下地壳源为主,随着壳/幔伸展作用演化,逐渐向混合源或独立幔源的演化。同时,不同地区岩浆源区的变化规律也显著不同。以胶辽地区为例,胶东整体上是壳幔混合源区对于岩浆演化有重要贡献;而辽东地区具有显著的源区演化特点：从剪切早期古老下地壳源区为主,并伴有幔源物质加入,剪切期古老下地壳为主,到剪切晚期和剪切期后以新生下地壳为主。本文认为岩石圈伸展的壳/幔拆离模型（Parallel Extension Tectonics）,可以合理地解释华北克拉通及邻区早白垩世构造-岩浆活动性。在该模型中,遭受伸展的华北克拉通岩石圈发生壳-幔拆离作用。在岩石圈伸展作用期间,地壳层次的拆离作用与岩石圈地幔层次上的拆离作用可以是耦合的或者是解耦的,从而导致华北克拉通岩石圈减薄过程中在地壳尺度上的拆离作用与变质核杂岩的剥露有三种不同的类型：同岩浆活动型伸展（C型：Co-magmatism mode extension）、无岩浆活动型伸展（A型：Amagmatism mode extension）和多阶段混合型（M型：Multi-mode extension）。