IceChrono1: a probabilistic model to compute a common and optimal chronology for several ice cores

Polar ice cores provide exceptional archives of past environmental conditions. The dating of ice cores and the estimation of the age-scale uncertainty are essential to interpret the climate and environmental records that they contain. It is, however, a complex problem which involves different methods. Here, we present IceChrono1, a new probabilistic model integrating various sources of chronological information to produce a common and optimized chronology for several ice cores, as well as its uncertainty. IceChrono1 is based on the inversion of three quantities: the surface accumulation rate, the lock-in depth (LID) of air bubbles and the thinning function. The chronological information integrated into the model are models of the sedimentation process (accumulation of snow, densification of snow into ice and air trapping, ice flow), ice- and air-dated horizons, ice and air depth intervals with known durations, Δdepth observations (depth shift between synchronous events recorded in the ice and in the air) and finally air and ice stratigraphic links in between ice cores. The optimization is formulated as a least squares problem, implying that all densities of probabilities are assumed to be Gaussian. It is numerically solved using the Levenberg–Marquardt algorithm and a numerical evaluation of the model's Jacobian. IceChrono follows an approach similar to that of the Datice model which was recently used to produce the AICC2012 (Antarctic ice core chronology) for four Antarctic ice cores and one Greenland ice core. IceChrono1 provides improvements and simplifications with respect to Datice from the mathematical, numerical and programming point of views. The capabilities of IceChrono1 are demonstrated on a case study similar to the AICC2012 dating experiment. We find results similar to those of Datice, within a few centuries, which is a confirmation of both IceChrono1 and Datice codes. We also test new functionalities with respect to the original version of Datice: observations as ice intervals with known durations, correlated observations, observations as air intervals with known durations and observations as mixed ice–air stratigraphic links. IceChrono1 is freely available under the General Public License v3 open source license

The Antarctic ice core chronology (AICC2012): an optimized multi-parameter and multi-site dating approach for the last 120 thousand years

The deep polar ice cores provide reference records commonly employed in global correlation of past climate events. However, temporal divergences reaching up to several thousand years (ka) exist between ice cores over the last climatic cycle. In this context, we are hereby introducing the Antarctic Ice Core Chronology 2012 (AICC2012), a new and coherent timescale developed for four Antarctic ice cores, namely Vostok, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML) and Talos Dome (TALDICE), alongside the Greenlandic NGRIP record. The AICC2012 timescale has been constructed using the Bayesian tool Datice (Lemieux-Dudon et al., 2010) that combines glaciological inputs and data constraints, including a wide range of relative and absolute gas and ice stratigraphic markers. We focus here on the last 120 ka, whereas the companion paper by Bazin et al. (2013) focuses on the interval 120–800 ka. Compared to previous timescales, AICC2012 presents an improved timing for the last glacial inception, respecting the glaciological constraints of all analyzed records. Moreover, with the addition of numerous new stratigraphic markers and improved calculation of the lock-in depth (LID) based on δ15N data employed as the Datice background scenario, the AICC2012 presents a slightly improved timing for the bipolar sequence of events over Marine Isotope Stage 3 associated with the seesaw mechanism, with maximum differences of about 600 yr with respect to the previous Datice-derived chronology of Lemieux-Dudon et al. (2010), hereafter denoted LD2010. Our improved scenario confirms the regional differences for the millennial scale variability over the last glacial period: while the EDC isotopic record (events of triangular shape) displays peaks roughly at the same time as the NGRIP abrupt isotopic increases, the EDML isotopic record (events characterized by broader peaks or even extended periods of high isotope values) reached the isotopic maximum several centuries before. It is expected that the future contribution of both other long ice core records and other types of chronological constraints to the Datice tool will lead to further refinements in the ice core chronologies beyond the AICC2012 chronology. For the time being however, we recommend that AICC2012 be used as the preferred chronology for the Vostok, EDC, EDML and TALDICE ice core records, both over the last glacial cycle (this study), and beyond (following Bazin et al., 2013). The ages for NGRIP in AICC2012 are virtually identical to those of GICC05 for the last 60.2 ka, whereas the ages beyond are independent of those in GICC05modelext (as in the construction of AICC2012, the GICC05modelext was included only via the background scenarios and not as age markers). As such, where issues of phasing between Antarctic records included in AICC2012 and NGRIP are involved, the NGRIP ages in AICC2012 should therefore be taken to avoid introducing false offsets. However for issues involving only Greenland ice cores, there is not yet a strong basis to recommend superseding GICC05modelext as the recommended age scale for Greenland ice cores

An optimized multi-proxy, multi-site Antarctic ice and gas orbital chronology (AICC2012): 120-800 ka

An accurate and coherent chronological framework is essential for the interpretation of climatic and environmental records obtained from deep polar ice cores. Until now, one common ice core age scale had been developed based on an inverse dating method (Datice), combining glaciological modelling with absolute and stratigraphic markers between 4 ice cores covering the last 50 ka (thousands of years before present) (Lemieux-Dudon et al., 2010). In this paper, together with the companion paper of Veres et al. (2013), we present an extension of this work back to 800 ka for the NGRIP, TALDICE, EDML, Vostok and EDC ice cores using an improved version of the Datice tool. The AICC2012 (Antarctic Ice Core Chronology 2012) chronology includes numerous new gas and ice stratigraphic links as well as improved evaluation of background and associated variance scenarios. This paper concentrates on the long timescales between 120–800 ka. In this framework, new measurements of δ18Oatm over Marine Isotope Stage (MIS) 11–12 on EDC and a complete δ18Oatm record of the TALDICE ice cores permit us to derive additional orbital gas age constraints. The coherency of the different orbitally deduced ages (from δ18Oatm, δO2/N2 and air content) has been verified before implementation in AICC2012. The new chronology is now independent of other archives and shows only small differences, most of the time within the original uncertainty range calculated by Datice, when compared with the previous ice core reference age scale EDC3, the Dome F chronology, or using a comparison between speleothems and methane. For instance, the largest deviation between AICC2012 and EDC3 (5.4 ka) is obtained around MIS 12. Despite significant modifications of the chronological constraints around MIS 5, now independent of speleothem records in AICC2012, the date of Termination II is very close to the EDC3 one

Black sea observing system

The ultimate goal of modern operational oceanography are end user oriented products with high scientific quality. Beneficiaries are the governmental services, coast and offshore based enterprises and research institutions that make use of the products generated by operational oceanography. Direct users are coastal managers, shipping, search and rescue, oil spill combat, offshore industry, ports, fishing, tourism, and recreation industry. Indirect beneficiaries, through climate forecasting based on ocean observations, are food, energy, water and medical suppliers. Availability of updated information on the actual state as well as forecast of marine environment is essential for the success and safety of maritime operations in the offshore industry. Various systems for the collection and presentation of marine data for the needs of different users have been developed and put in operation in the Black Sea. The systems are located both along the coast and in the open sea and the information they provide is used by both the maritime industry and the widest range of users. The Black Sea Monitoring and Forecasting Center in the frame of the Copernicus Marine Service is providing regular and systematic information about the physical state of the ocean, marine ecosystem and wave conditions in the Black Sea area, assimilating observations, keeping efficient operations, advanced technology and high quality modeling products. Combining and optimizing in situ, remote sensing, modeling and forecasting into a Black Sea observing system is a task that has to be solved, and that will allow to get a more complete and comprehensive picture of the state of the marine environment as well as to forecast future changes of physical and biogeochemical state of the Black Sea and the Black Sea ecosystem

Abrupt Ice Age Shifts in Southern Westerlies and Antarctic Climate Forced from the North

The Southern Hemisphere (SH) mid-latitude westerly winds play a central role in the global climate system via Southern Ocean upwelling, carbon exchange with the deep ocean, Agulhas Leakage, and Antarctic ice sheet stability. Meridional shifts in the SH westerlies have been hypothesized in response to abrupt North Atlantic Dansgaard-Oeschger (DO) climatic events of the last ice age, in parallel with the well-documented shifts of the intertropical convergence zone. Shifting moisture pathways to West Antarctica are consistent with this view, but may represent a Pacific teleconnection pattern. The full SH atmospheric-circulation response to the DO cycle, as well as its impact on Antarctic temperature, have so far remained unclear. Here we use five volcanically-synchronized ice cores to show that the Antarctic temperature response to the DO cycle can be understood as the superposition of two modes: a spatially homogeneous oceanic “bipolar seesaw” mode that lags Northern Hemisphere (NH) climate by about 200 years, and a spatially heterogeneous atmospheric mode that is synchronous with NH abrupt events. Temperature anomalies of the atmospheric mode are similar to those associated with present-day Southern Annular Mode (SAM) variability, rather than the Pacific South America (PSA) pattern. Moreover, deuterium excess records suggest a zonally coherent migration of the SH westerlies over all ocean basins in phase with NH climate. Our work provides a simple conceptual framework for understanding the circum-Antarctic temperature response to abrupt NH climate change. We provide observational evidence for abrupt shifts in the SH westerlies, with ramifications for global ocean circulation and atmospheric CO₂. These coupled changes highlight the necessity of a global, rather than a purely North Atlantic, perspective on the DO cycle

Copernicus Marine Service Ocean State Report

This is the final version. Available from Taylor & Francis via the DOI in this record

Формування та розвиток загальної теорії стійкості (середина XVIII ст. — 30-і рр. ХХ ст.)

У статті розглянуто історію вивчення стійкості (середина XVIII — початок XX ст., світовий контекст). Досліджено внесок А. Пуанкаре та О.М. Ляпунова в розвиток загальної теорії стійкості. Показано розвиток їх ідей у працях російських та українських учених.В статье рассмотрена история изучения устойчивости (середина XVIII — начало XX в., мировой контекст). Исследован вклад французского ученого А. Пуанкаре и русского ученого А.М. Ляпунова в развитие общей теории устойчивости. Показано дальнейшее развитие их идей в трудах русских и украинских ученых.The history of basic research in stability is given. Contributions from H.Poincaré, a French mathematician, mechanic and physicist, and O. Lapunov, a soviet mathematician and mechanic (working in the Kharkiv university) to development of the general theory of stability are shown. In 1892—1902, O. Lyapunov constructed an original robust mathematical apparatus to study stability of motion. Development of ideas and methods of H.Poincar of H.Poincar³e and O. Lapunov in works of later Ukrainian and Russian scientists is shown

Extending an oceanographic variational scheme to allow for affordable hybrid and four-dimensional data assimilation

The traditional formulation of three-dimensional variational (3DVAR) data assimilation schemes for oceanographic applications neglects the temporal evolution of background errors within and across assimilation temporal windows. Such a simplification may be limiting for many climate (e.g. reanalyses) and operational (e.g. medium-range forecast) applications. This work explores possible extensions of the OceanVar data assimilation code aiming at overcoming these limitations. General formulations are proposed and implemented in order to extend the 3DVAR scheme of OceanVar into a simplified hybrid (ensemble-variational) four-dimensional variational (4DVAR) assimilation scheme, where (i) background-error covariances combine stationary and flow-dependent components through an augmented control vector and (ii) a simplified tangent-linear and adjoint model, which assumes that only temperature and salinity are independent variables. These extensions are shown to allow the background-error covariances to follow the time-varying structure typical of climate modes like ENSO, and to shape the analysis increments in agreement with the underlying ocean circulation, respectively. The two extensions are cross-compared in terms of computational time cost and accuracy and further combined together into a hybrid 4DVAR scheme. The hybrid formulation provides in general largely positive impact at short forecast ranges, while 4DVAR at long ones. The hybrid 4DVAR scheme improves the verification skill scores in most cases