ALBERT

Description: Can we determine what controls the spatio-temporal distribution of d-excess and 17 O-excess in precipitation using the LMDZ general circulation model? Climate of the Past, 9, 2173-2193, 2013 Author(s): C. Risi, A. Landais, R. Winkler, and F. Vimeux Combined measurements of the H 2 18 O and HDO isotopic ratios in precipitation, leading to second-order parameter D-excess, have provided additional constraints on past climates compared to the H 2 18 O isotopic ratio alone. More recently, measurements of H 2 17 O have led to another second-order parameter: 17 O-excess. Recent studies suggest that 17 O-excess in polar ice may provide information on evaporative conditions at the moisture source. However, the processes controlling the spatio-temporal distribution of 17 O-excess are still far from being fully understood. We use the isotopic general circulation model (GCM) LMDZ to better understand what controls d-excess and 17 O-excess in precipitation at present-day (PD) and during the last glacial maximum (LGM). The simulation of D-excess and 17 O-excess is evaluated against measurements in meteoric water, water vapor and polar ice cores. A set of sensitivity tests and diagnostics are used to quantify the relative effects of evaporative conditions (sea surface temperature and relative humidity), Rayleigh distillation, mixing between vapors from different origins, precipitation re-evaporation and supersaturation during condensation at low temperature. In LMDZ, simulations suggest that in the tropics convective processes and rain re-evaporation are important controls on precipitation D-excess and 17 O-excess. In higher latitudes, the effect of distillation, mixing between vapors from different origins and supersaturation are the most important controls. For example, the lower d-excess and 17 O-excess at LGM simulated at LGM are mainly due to the supersaturation effect. The effect of supersaturation is however very sensitive to a parameter whose tuning would require more measurements and laboratory experiments. Evaporative conditions had previously been suggested to be key controlling factors of d-excess and 17 O-excess, but LMDZ underestimates their role. More generally, some shortcomings in the simulation of 17 O-excess by LMDZ suggest that general circulation models are not yet the perfect tool to quantify with confidence all processes controlling 17 O-excess.

Description: Magnetic field and dynamic pressure ULF fluctuations in coronal-mass-ejection-driven sheath regions Annales Geophysicae, 31, 1559-1567, 2013 Author(s): E. K. J. Kilpua, H. Hietala, H. E. J. Koskinen, D. Fontaine, and L. Turc Compressed sheath regions form ahead of interplanetary coronal mass ejections (ICMEs) that are sufficiently faster than the preceding solar wind. The turbulent sheath regions are important drivers of magnetospheric activity, but due to their complex internal structure, relatively little is known on the distribution of the magnetic field and plasma variations in them. In this paper we investigate ultra low frequency (ULF) fluctuations in the interplanetary magnetic field (IMF) and in dynamic pressure ( P dyn ) using a superposed epoch analysis of 41 sheath regions observed during solar cycle 23. We find strongest fluctuation power near the shock and in the vicinity of the ICME leading edge. The IMF and P dyn ULF power have different profiles within the sheath; the former is enhanced in the leading part of the sheath, while the latter is increased in the trailing part of the sheath. We also find that the ICME properties affect the level and distribution of the ULF power in sheath regions. For example, sheath regions associated with strong or fast ICMEs, or those that are crossed at intermediate distances from the center, have strongest ULF power and large variation in the power throughout the sheath region. The weaker or slower ICMEs, or those that are crossed centrally, have in general considerably weaker ULF power with relatively smooth profiles. The strong and abrupt decrease of the IMF ULF power at the ICME leading edge could be used to distinguish the ICME from the preceding sheath plasma.

Description: Cluster observation of few-hour-scale evolution of structured plasma in the inner magnetosphere Annales Geophysicae, 31, 1569-1578, 2013 Author(s): M. Yamauchi, I. Dandouras, H. Rème, R. Lundin, and L. M. Kistler Using Cluster Ion Spectrometry (CIS) data from the spacecraft-4 perigee traversals during the 2001–2006 period (nearly 500 traversals after removing those that are highly contaminated by radiation belt particles), we statistically examined the local time distribution of structured trapped ions at sub- to few-keV range as well as inbound–outbound differences of these ion signatures in intensities and energy–latitude dispersion directions. Since the Cluster orbit during this period was almost constant and approximately north–south symmetric at nearly constant local time near the perigee, inbound–outbound differences are attributed to temporal developments in a 1–2 h timescale. Three types of structured ions at sub- to few keV range that are commonly found in the inner magnetosphere are examined: – Energy–latitude dispersed structured ions at less than a few keV, – Short-lived dispersionless ion stripes at wide energy range extending 0.1–10 keV, – Short-lived low-energy ion bursts at less than a few hundred eV. The statistics revealed that the wedge-like dispersed ions are most often observed in the dawn sector (60% of traversals), and a large portion of them show significant enhancement during the traversals at all local times. The short-lived ion stripes are predominantly found near midnight, where most stripes are significantly enhanced during the traversals and are associated with substorm activities with geomagnetic AL 〈 −300 nT. The low-energy bursts are observed at all local times and under all geomagnetic conditions, with moderate peak of the occurrence rate in the afternoon sector. A large portion of them again show significant enhancement or decay during the traversals.

Description: Global and regional sea surface temperature trends during Marine Isotope Stage 11 Climate of the Past, 9, 2231-2252, 2013 Author(s): Y. Milker, R. Rachmayani, M. F. G. Weinkauf, M. Prange, M. Raitzsch, M. Schulz, and M. Kučera The Marine Isotope Stage (MIS) 11 (424–374 ka) was characterized by a protracted deglaciation and an unusually long climatic optimum. It remains unclear to what degree the climate development during this interglacial reflects the unusually weak orbital forcing or greenhouse gas trends. Previously, arguments about the duration and timing of the MIS11 climatic optimum and about the pace of the deglacial warming were based on a small number of key records, which appear to show regional differences. In order to obtain a global signal of climate evolution during MIS11, we compiled a database of 78 sea surface temperature (SST) records from 57 sites spanning MIS11, aligned these individually on the basis of benthic ( N = 28) or planktonic ( N = 31) stable oxygen isotope curves to a common time frame and subjected 48 of them to an empirical orthogonal function (EOF) analysis. The analysis revealed a high commonality among all records, with the principal SST trend explaining almost 49% of the variability. This trend indicates that on the global scale, the surface ocean underwent rapid deglacial warming during Termination V, in pace with carbon dioxide rise, followed by a broad SST optimum centered at ~410 kyr. The second EOF, which explained ~18% of the variability, revealed the existence of a different SST trend, characterized by a delayed onset of the temperature optimum during MIS11 at ~398 kyr, followed by a prolonged warm period lasting beyond 380 kyr. This trend is most consistently manifested in the mid-latitude North Atlantic and Mediterranean Sea and is here attributed to the strength of the Atlantic meridional overturning circulation. A sensitivity analysis indicates that these results are robust to record selection and to age-model uncertainties of up to 3–6 kyr, but more sensitive to SST seasonal attribution and SST uncertainties 〉1 °C. In order to validate the CCSM3 (Community Climate System Model, version 3) predictive potential, the annual and seasonal SST anomalies recorded in a total of 74 proxy records were compared with runs for three time slices representing orbital configuration extremes during the peak interglacial of MIS11. The modeled SST anomalies are characterized by a significantly lower variance compared to the reconstructions. Nevertheless, significant correlations between proxy and model data are found in comparisons on the seasonal basis, indicating that the model captures part of the long-term variability induced by astronomical forcing, which appears to have left a detectable signature in SST trends.

Description: Influence of the terrestrial magnetic field geometry on the cutoff rigidity of cosmic ray particles Annales Geophysicae, 31, 1637-1643, 2013 Author(s): K. Herbst, A. Kopp, and B. Heber Studies of the propagation of charged energetic particles in the Earth's magnetic field go back to Carl Størmer. In the end, his investigations finally lead to the definition of the so-called cutoff rigidity R C ; that is, the minimum momentum per charge a particle must have in order to reach a certain geographical location. Employing Monte Carlo simulations with the PLANETOCOSMICS code we investigate the correlation between the geomagnetic field structure and the cutoff rigidity. We show that the geometry of the magnetic field has a considerable influence on the resulting cutoff rigidity distribution. Furthermore, we will present a simple geometry-based parameter, δ B , which is able to reflect the location-dependent cutoff rigidity. We show that this correlation is also visible in the temporal evolution of the Earth's magnetic field, at least over the last 100 yr. Using latitude scans with neutron monitors, changes of the relative counting rates at different positions are calculated, showing small variations for, e.g., Kiel and Moscow, while large ones occur at Mexico City as well as on the British Virgin Islands.

Description: Electron acceleration at Jupiter: input from cyclotron-resonant interaction with whistler-mode chorus waves Annales Geophysicae, 31, 1619-1630, 2013 Author(s): E. E. Woodfield, R. B. Horne, S. A. Glauert, J. D. Menietti, and Y. Y. Shprits Jupiter has the most intense radiation belts of all the outer planets. It is not yet known how electrons can be accelerated to energies of 10 MeV or more. It has been suggested that cyclotron-resonant wave-particle interactions by chorus waves could accelerate electrons to a few MeV near the orbit of Io. Here we use the chorus wave intensities observed by the Galileo spacecraft to calculate the changes in electron flux as a result of pitch angle and energy diffusion. We show that, when the bandwidth of the waves and its variation with L are taken into account, pitch angle and energy diffusion due to chorus waves is a factor of 8 larger at L -shells greater than 10 than previously shown. We have used the latitudinal wave intensity profile from Galileo data to model the time evolution of the electron flux using the British Antarctic Survey Radiation Belt (BAS) model. This profile confines intense chorus waves near the magnetic equator with a peak intensity at ∼5° latitude. Electron fluxes in the BAS model increase by an order of magnitude for energies around 3 MeV. Extending our results to L = 14 shows that cyclotron-resonant interactions with chorus waves are equally important for electron acceleration beyond L = 10. These results suggest that there is significant electron acceleration by cyclotron-resonant interactions at Jupiter contributing to the creation of Jupiter's radiation belts and also increasing the range of L -shells over which this mechanism should be considered.

Description: The coupling between the solar wind and proton fluxes at GEO Annales Geophysicae, 31, 1631-1636, 2013 Author(s): R. J. Boynton, S. A. Billings, O. A. Amariutei, and I. Moiseenko The relationship between the solar wind and the proton flux at geosynchronous Earth orbit (GEO) is investigated using the error reduction ratio (ERR) analysis. The ERR analysis is able to search for the most appropriate inputs that control the evolution of the system. This approach is a black box method and is able to derive a mathematical model of a system from input-output data. This method is used to analyse eight energy ranges of the proton flux at GEO from 80 keV to 14.5 MeV. The inputs to the algorithm were solar wind velocity, density and pressure; the Dst index; the solar energetic proton (SEP) flux; and a function of the interplanetary magnetic field (IMF) tangential magnitude and clock angle. The results show that for lowest five energy channels (80 to 800 keV) the GEO proton fluxes are controlled by the solar wind velocity with a lag of two to three days. However, above 350 keV, the SEP fluxes, accounts for a significant portion of the GEO proton flux variance. For the highest three energy channels (0.74 to 14.5 MeV), the SEPs account for the majority of the ERR. The results also show an anisotropy of protons with gyrocenters inside GEO and outside GEO, where the protons inside GEO are controlled partly by the Dst index and also an IMF-clock angle function.

Description: A statistical study on O + flux in the dayside magnetosheath Annales Geophysicae, 31, 1005-1010, 2013 Author(s): R. Slapak, H. Nilsson, and L. G. Westerberg Studies on terrestrial oxygen ion (O + ) escape into the interplanetary space have considered a number of different escape paths. Recent observations however suggest a yet insufficiently investigated additional escape route for hot O + : along open magnetic field lines in the high altitude cusp and mantle. Here we present a statistical study on O + flux in the high-latitude dayside magnetosheath. The O + is generally seen relatively close to the magnetopause, consistent with observations of O + flowing primarily tangentially to the magnetopause. We estimate the total escape flux in this region to be ~ 7 × 10 24 s −1 , implying this escape route to significantly contribute to the overall total O + loss into interplanetary space.

Description: Simulation of the influence of historical land cover changes on the global climate Annales Geophysicae, 31, 995-1004, 2013 Author(s): Y. Wang, X. Yan, and Z. Wang In order to estimate biogeophysical effects of historical land cover change on climate during last three centuries, a set of experiments with a climate system model of intermediate complexity (MPM-2) is performed. In response to historical deforestation, the model simulates a decrease in annual mean global temperature in the range of 0.07–0.14 °C based on different grassland albedos. The effect of land cover changes is most pronounced in the middle northern latitudes with maximum cooling reaching approximately 0.6 °C during northern summer. The cooling reaches 0.57 °C during northern spring owing to the large effects of land surface albedo. These results suggest that land cover forcing is important for study on historical climate change and that more research is necessary in the assessment of land management options for climate change mitigation.

Description: IMF effect on the polar cap contraction and expansion during a period of substorms Annales Geophysicae, 31, 1021-1034, 2013 Author(s): A. T. Aikio, T. Pitkänen, I. Honkonen, M. Palmroth, and O. Amm The polar cap boundary (PCB) location and motion in the nightside ionosphere has been studied by using measurements from the EISCAT radars and the MIRACLE magnetometers during a period of four substorms on 18 February 2004. The OMNI database has been used for observations of the solar wind and the Geotail satellite for magnetospheric measurements. In addition, the event was modelled by the GUMICS-4 MHD simulation. The simulation of the PCB location was in a rather good agreement with the experimental estimates at the EISCAT longitude. During the first three substorm expansion phases, neither the local observations nor the global simulation showed any poleward motions of the PCB, even though the electrojets intensified. Rapid poleward motions of the PCB took place only in the early recovery phases of the substorms. Hence, in these cases the nightside reconnection rate was locally higher in the recovery phase than in the expansion phase. In addition, we suggest that the IMF B z component correlated with the nightside tail inclination angle and the PCB location with about a 17-min delay from the bow shock. By taking the delay into account, the IMF northward turnings were associated with dipolarizations of the magnetotail and poleward motions of the PCB in the recovery phase. The mechanism behind this effect should be studied further.

Description: A comparison of bow shock models with Cluster observations during low Alfvén Mach number magnetic clouds Annales Geophysicae, 31, 1011-1019, 2013 Author(s): L. Turc, D. Fontaine, P. Savoini, H. Hietala, and E. K. J. Kilpua Magnetic clouds (MCs) are very geoeffective solar wind structures. Their properties in the interplanetary medium have been extensively studied, yet little is known about their characteristics in the Earth's magnetosheath. The Cluster spacecraft offer the opportunity to observe MCs in the magnetosheath, but before MCs reach the magnetosphere, their structure is altered when they interact with the terrestrial bow shock (BS). The physics taking place at the BS strongly depends on Θ Bn , the angle between the shock normal and the interplanetary magnetic field. However, in situ observations of the BS during an MC's crossing are seldom available. In order to relate magnetosheath observations to solar wind conditions, we need to rely on a model to determine the shock's position and normal direction. Yet during MCs, the models tend to be less accurate, because the Alfvén Mach number ( M A ) is often significantly lower than in regular solar wind. On the contrary, the models are generally optimised for high M A conditions. In this study, we compare the predictions of four widely used models available in the literature (Wu et al., 2000; Chapman and Cairns, 2003; Jeřáb et al., 2005; Měrka et al., 2005b) to Cluster's dayside BS crossings observed during five MC events. Our analysis shows that the Θ Bn angle is well predicted by all four models. On the other hand, the Jeřáb et al. (2005) model yields the best estimates of the BS position during low M A MCs. The other models locate the BS either too far from or too close to Earth. The results of this paper can be directly used to estimate the BS parameters in all studies of MC interaction with Earth's magnetosphere.

Description: Skill and reliability of climate model ensembles at the Last Glacial Maximum and mid-Holocene Climate of the Past, 9, 811-823, 2013 Author(s): J. C. Hargreaves, J. D. Annan, R. Ohgaito, A. Paul, and A. Abe-Ouchi Paleoclimate simulations provide us with an opportunity to critically confront and evaluate the performance of climate models in simulating the response of the climate system to changes in radiative forcing and other boundary conditions. Hargreaves et al. (2011) analysed the reliability of the Paleoclimate Modelling Intercomparison Project, PMIP2 model ensemble with respect to the MARGO sea surface temperature data synthesis (MARGO Project Members, 2009) for the Last Glacial Maximum (LGM, 21 ka BP). Here we extend that work to include a new comprehensive collection of land surface data (Bartlein et al., 2011), and introduce a novel analysis of the predictive skill of the models. We include output from the PMIP3 experiments, from the two models for which suitable data are currently available. We also perform the same analyses for the PMIP2 mid-Holocene (6 ka BP) ensembles and available proxy data sets. Our results are predominantly positive for the LGM, suggesting that as well as the global mean change, the models can reproduce the observed pattern of change on the broadest scales, such as the overall land–sea contrast and polar amplification, although the more detailed sub-continental scale patterns of change remains elusive. In contrast, our results for the mid-Holocene are substantially negative, with the models failing to reproduce the observed changes with any degree of skill. One cause of this problem could be that the globally- and annually-averaged forcing anomaly is very weak at the mid-Holocene, and so the results are dominated by the more localised regional patterns in the parts of globe for which data are available. The root cause of the model-data mismatch at these scales is unclear. If the proxy calibration is itself reliable, then representativity error in the data-model comparison, and missing climate feedbacks in the models are other possible sources of error.

Description: Investigation of radiative effects of the optically thick dust layer over the Indian tropical region Annales Geophysicae, 31, 647-663, 2013 Author(s): S. K. Das, J.-P. Chen, M. Venkat Ratnam, and A. Jayaraman Optical and physical properties of aerosols derived from multi-satellite observations (MODIS-Aqua, OMI-Aura, MISR-Terra, CALIOP-CALIPSO) have been used to estimate radiative effects of the dust layer over southern India. The vertical distribution of aerosol radiative forcing and heating rates are calculated with 100 m resolution in the lower atmosphere, using temperature and relative humidity data from balloon-borne radiosonde observations. The present study investigates the optically thick dust layer of optical thickness 0.18 ± 0.06 at an altitude of 2.5 ± 0.7 km over Gadanki, transported from the Thar Desert, producing radiative forcing and heating rate of 11.5 ± 3.3 W m −2 and 0.6 ± 0.26 K day −1 , respectively, with a forcing efficiency of 43 W m −2 and an effective heating rate of 4 K day −1 per unit dust optical depth. Presence of the dust layer increases radiative forcing by 60% and heating rate by 60 times at that altitude compared to non-dusty cloud-free days. Calculation shows that the radiative effects of the dust layer strongly depend on the boundary layer aerosol type and mass loading. An increase of 25% of heating by the dust layer is found over relatively cleaner regions than urban regions in southern India and further 15% of heating increases over the marine region. Such heating differences in free troposphere may have significant consequences in the atmospheric circulation and hydrological cycle over the tropical Indian region.

Description: Stable isotopic evidence of El Niño-like atmospheric circulation in the Pliocene western United States Climate of the Past, 9, 903-912, 2013 Author(s): M. J. Winnick, J. M. Welker, and C. P. Chamberlain Understanding how the hydrologic cycle has responded to warmer global temperatures in the past is especially important today as concentrations of CO 2 in the atmosphere continue to increase due to human activities. The Pliocene offers an ideal window into a climate system that has equilibrated with current atmospheric p CO 2 . During the Pliocene the western United States was wetter than modern, an observation at odds with our current understanding of future warming scenarios, which involve the expansion and poleward migration of the subtropical dry zone. Here we compare Pliocene oxygen isotope profiles of pedogenic carbonates across the western US to modern isotopic anomalies in precipitation between phases of the El Niño–Southern Oscillation (ENSO). We find that when accounting for seasonality of carbonate formation, isotopic changes through the late Pliocene match modern precipitation isotopic anomalies in El Niño years. Furthermore, isotopic shifts through the late Pliocene mirror changes through the early Pleistocene, which likely represents the southward migration of the westerly storm track caused by growth of the Laurentide ice sheet. We propose that the westerly storm track migrated northward through the late Pliocene with the development of the modern cold tongue in the east equatorial Pacific, then returned southward with widespread glaciation in the Northern Hemisphere – a scenario supported by terrestrial climate proxies across the US. Together these data support the proposed existence of background El Niño-like conditions in western North America during the warm Pliocene. If the earth behaves similarly with future warming, this observation has important implications with regard to the amount and distribution of precipitation in western North America.

Description: Mass-movement and flood-induced deposits in Lake Ledro, southern Alps, Italy: implications for Holocene palaeohydrology and natural hazards Climate of the Past, 9, 825-840, 2013 Author(s): A. Simonneau, E. Chapron, B. Vannière, S. B. Wirth, A. Gilli, C. Di Giovanni, F. S. Anselmetti, M. Desmet, and M. Magny High-resolution seismic profiles and sediment cores from Lake Ledro combined with soil and riverbed samples from the lake's catchment area are used to assess the recurrence of natural hazards (earthquakes and flood events) in the southern Italian Alps during the Holocene. Two well-developed deltas and a flat central basin are identified on seismic profiles in Lake Ledro. Lake sediments have been finely laminated in the basin since 9000 cal. yr BP and frequently interrupted by two types of sedimentary events (SEs): light-coloured massive layers and dark-coloured graded beds. Optical analysis (quantitative organic petrography) of the organic matter present in soil, riverbed and lacustrine samples together with lake sediment bulk density and grain-size analysis illustrate that light-coloured layers consist of a mixture of lacustrine sediments and mainly contain algal particles similar to the ones observed in background sediments. Light-coloured layers thicker than 1.5 cm in the main basin of Lake Ledro are synchronous to numerous coeval mass-wasting deposits remoulding the slopes of the basin. They are interpreted as subaquatic mass-movements triggered by historical and pre-historical regional earthquakes dated to AD 2005, AD 1891, AD 1045 and 1260, 2545, 2595, 3350, 3815, 4740, 7190, 9185 and 11 495 cal. yr BP. Dark-coloured SEs develop high-amplitude reflections in front of the deltas and in the deep central basin. These beds are mainly made of terrestrial organic matter (soils and lignocellulosic debris) and are interpreted as resulting from intense hyperpycnal flood event. Mapping and quantifying the amount of soil material accumulated in the Holocene hyperpycnal flood deposits of the sequence allow estimating that the equivalent soil thickness eroded over the catchment area reached up to 5 mm during the largest Holocene flood events. Such significant soil erosion is interpreted as resulting from the combination of heavy rainfall and snowmelt. The recurrence of flash flood events during the Holocene was, however, not high enough to affect pedogenesis processes and highlight several wet regional periods during the Holocene. The Holocene period is divided into four phases of environmental evolution. Over the first half of the Holocene, a progressive stabilization of the soils present through the catchment of Lake Ledro was associated with a progressive reforestation of the area and only interrupted during the wet 8.2 event when the soil destabilization was particularly important. Lower soil erosion was recorded during the mid-Holocene climatic optimum (8000–4200 cal. yr BP) and associated with higher algal production. Between 4200 and 3100 cal. yr BP, both wetter climate and human activities within the drainage basin drastically increased soil erosion rates. Finally, from 3100 cal. yr BP to the present-day, data suggest increasing and changing human land use.

Description: Using data assimilation to investigate the causes of Southern Hemisphere high latitude cooling from 10 to 8 ka BP Climate of the Past, 9, 887-901, 2013 Author(s): P. Mathiot, H. Goosse, X. Crosta, B. Stenni, M. Braida, H. Renssen, C. J. Van Meerbeeck, V. Masson-Delmotte, A. Mairesse, and S. Dubinkina From 10 to 8 ka BP (thousand years before present), paleoclimate records show an atmospheric and oceanic cooling in the high latitudes of the Southern Hemisphere. During this interval, temperatures estimated from proxy data decrease by 0.8 °C over Antarctica and 1.2 °C over the Southern Ocean. In order to study the causes of this cooling, simulations covering the early Holocene have been performed with the climate model of intermediate complexity LOVECLIM constrained to follow the signal recorded in climate proxies using a data assimilation method based on a particle filtering approach. The selected proxies represent oceanic and atmospheric surface temperature in the Southern Hemisphere derived from terrestrial, marine and glaciological records. Two mechanisms previously suggested to explain the 10–8 ka BP cooling pattern are investigated using the data assimilation approach in our model. The first hypothesis is a change in atmospheric circulation, and the second one is a cooling of the sea surface temperature in the Southern Ocean, driven in our experimental setup by the impact of an increased West Antarctic melting rate on ocean circulation. For the atmosphere hypothesis, the climate state obtained by data assimilation produces a modification of the meridional atmospheric circulation leading to a 0.5 °C Antarctic cooling from 10 to 8 ka BP compared to the simulation without data assimilation, without congruent cooling of the atmospheric and sea surface temperature in the Southern Ocean. For the ocean hypothesis, the increased West Antarctic freshwater flux constrainted by data assimilation (+100 mSv from 10 to 8 ka BP) leads to an oceanic cooling of 0.7 °C and a strengthening of Southern Hemisphere westerlies (+6%). Thus, according to our experiments, the observed cooling in Antarctic and the Southern Ocean proxy records can only be reconciled with the reconstructions by the combination of a modified atmospheric circulation and an enhanced freshwater flux.

Description: Low level jet intensification by mineral dust aerosols Annales Geophysicae, 31, 625-632, 2013 Author(s): O. Alizadeh Choobari, P. Zawar-Reza, and A. Sturman Modification of the intensity of a low level jet (LLJ) and near-surface wind speed by mineral dust is important as it has implications for dust emission and its long-range transport. Using the Weather Research and Forecasting with Chemistry (WRF/Chem) regional model, it is shown that direct radiative forcing by mineral dust reduces temperature in the lower atmosphere, but increases it in the layers aloft. The surface cooling is shown to be associated with a reduction of turbulent kinetic energy (TKE) and hence vertical mixing of horizontal momentum. Changes in the vertical profile of temperature over the regions that are under the influence of a LLJ are shown to result in an intensification of the LLJ and near-surface wind speed, but a decrease of winds aloft. These changes in the wind speed profile differ from results of previous research which suggested a decrease of wind speed in the lower atmosphere and its increase in the upper boundary layer.

Description: Parametric validations of analytical lifetime estimates for radiation belt electron diffusion by whistler waves Annales Geophysicae, 31, 599-624, 2013 Author(s): A. V. Artemyev, D. Mourenas, O. V. Agapitov, and V. V. Krasnoselskikh The lifetimes of electrons trapped in Earth's radiation belts can be calculated from quasi-linear pitch-angle diffusion by whistler-mode waves, provided that their frequency spectrum is broad enough and/or their average amplitude is not too large. Extensive comparisons between improved analytical lifetime estimates and full numerical calculations have been performed in a broad parameter range representative of a large part of the magnetosphere from L ~ 2 to 6. The effects of observed very oblique whistler waves are taken into account in both numerical and analytical calculations. Analytical lifetimes (and pitch-angle diffusion coefficients) are found to be in good agreement with full numerical calculations based on CRRES and Cluster hiss and lightning-generated wave measurements inside the plasmasphere and Cluster lower-band chorus waves measurements in the outer belt for electron energies ranging from 100 keV to 5 MeV. Comparisons with lifetimes recently obtained from electron flux measurements on SAMPEX, SCATHA, SAC-C and DEMETER also show reasonable agreement.

Description: Coherent amplitude modulation of electron-beam-driven Langmuir waves Annales Geophysicae, 31, 633-638, 2013 Author(s): K. Baumgärtel A linear approach to the phenomenon of irregular amplitude modulation of beam-driven Langmuir waves, developed in a previous paper, is extended to explain periodic modulation as well. It comes about by beating of the fastest growing mode of the instability with beam-aligned plasma oscillations. They are naturally generated in a uniform domain of beam–plasma interaction prior to the onset of the instability. Particle-in-cell (PIC) simulations support the results of the linear analysis.

Description: Model sensitivity to North Atlantic freshwater forcing at 8.2 ka Climate of the Past, 9, 955-968, 2013 Author(s): C. Morrill, A. N. LeGrande, H. Renssen, P. Bakker, and B. L. Otto-Bliesner We compared four simulations of the 8.2 ka event to assess climate model sensitivity and skill in responding to North Atlantic freshwater perturbations. All of the simulations used the same freshwater forcing, 2.5 Sv for one year, applied to either the Hudson Bay (northeastern Canada) or Labrador Sea (between Canada's Labrador coast and Greenland). This freshwater pulse induced a decadal-mean slowdown of 10–25% in the Atlantic Meridional Overturning Circulation (AMOC) of the models and caused a large-scale pattern of climate anomalies that matched proxy evidence for cooling in the Northern Hemisphere and a southward shift of the Intertropical Convergence Zone. The multi-model ensemble generated temperature anomalies that were just half as large as those from quantitative proxy reconstructions, however. Also, the duration of AMOC and climate anomalies in three of the simulations was only several decades, significantly shorter than the duration of ~150 yr in the paleoclimate record. Possible reasons for these discrepancies include incorrect representation of the early Holocene climate and ocean state in the North Atlantic and uncertainties in the freshwater forcing estimates.

Description: Seasonal variation of the ion upflow in the topside ionosphere during SAPS (subauroral polarization stream) periods Annales Geophysicae, 31, 1521-1534, 2013 Author(s): H. Wang and H. Lühr A statistical study has been performed by using two years of DMSP (Defense Meteorological Satellite Program) plasma observations to investigate the seasonal effect of SAPS (subauroral polarization stream) on the ion upflow in the duskside ionosphere of the Northern Hemisphere. There are obvious upflows occurring in the topside ionosphere around the SAPS region, exceeding 200 m s −1 at winter solstice, indicating an important relationship between SAPS and the local plasma upward motion. Both SAPS and ion upward velocities show similar seasonal variations, largest in winter and smallest in summer, irrespective of geomagnetic activity. A good correlation is found and a linear relationship is derived between SAPS and the ion upflow velocities. During December solstice the average upflow flux can reach about 2 × 10 8 cm −2 s −1 for more disturbed periods, which is comparable to the typical upflow flux in the dayside cusp region. The depression of the ion temperatures around the peak SAPS region can be understood in terms of the adiabatic cooling. The hot ion cools down when expanding into the low ion concentration region. The electron temperature elevates around the SAPS region because of the reduced Coulomb cooling in the low ion density region. Both the changes of ion and electron temperatures are larger in winter than in summer, however, for Kp 〈 4 the electron temperatures are almost seasonably independent. The present work highlights the important role of the SAPS-related frictional heating at mid-latitudes on the local formation of the strong upward flow, which might provide a direct ionospheric ion source for the ring current and plasmasphere in the duskside sector.

Description: Magnetic field generation in a jet-sheath plasma via the kinetic Kelvin-Helmholtz instability Annales Geophysicae, 31, 1535-1541, 2013 Author(s): K.-I. Nishikawa, P. Hardee, B. Zhang, I. Duţan, M. Medvedev, E. J. Choi, K. W. Min, J. Niemiec, Y. Mizuno, A. Nordlund, J. T. Frederiksen, H. Sol, M. Pohl, and D. H. Hartmann We have investigated the generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin–Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of the kinetic Kelvin–Helmholtz instability (KKHI) of our jet-sheath configuration is slightly different, even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field E z , perpendicular to the flow boundary, and the magnetic field B y , transverse to the flow direction. After the B y component is excited, an induced electric field E x , parallel to the flow direction, becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios m i / m e = 1836 and m i / m e = 20 are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case (γ j = 1.5) is larger than for a relativistic jet case (γ j = 15).

Description: Ionospheric response to total solar eclipse of 22 July 2009 in different Indian regions Annales Geophysicae, 31, 1549-1558, 2013 Author(s): S. Kumar, A. K. Singh, and R. P. Singh The variability of ionospheric response to the total solar eclipse of 22 July 2009 has been studied analyzing the GPS data recorded at the four Indian low-latitude stations Varanasi (100% obscuration), Kanpur (95% obscuration), Hyderabad (84% obscuration) and Bangalore (72% obscuration). The retrieved ionospheric vertical total electron content (VTEC) shows a significant reduction (reflected by all PRNs (satellites) at all stations) with a maximum of 48% at Varanasi (PRN 14), which decreases to 30% at Bangalore (PRN 14). Data from PRN 31 show a maximum of 54% at Kanpur and 26% at Hyderabad. The maximum decrement in VTEC occurs some time (2–15 min) after the maximum obscuration. The reduction in VTEC compared to the quiet mean VTEC depends on latitude as well as longitude, which also depends on the location of the satellite with respect to the solar eclipse path. The amount of reduction in VTEC decreases as the present obscuration decreases, which is directly related to the electron production by the photoionization process. The analysis of electron density height profile derived from the COSMIC (Constellation Observing System for Meteorology, Ionosphere & Climate) satellite over the Indian region shows significant reduction from 100 km altitude up to 800 km altitude with a maximum of 48% at 360 km altitude. The oscillatory nature in total electron content data at all stations is observed with different wave periods lying between 40 and 120 min, which are attributed to gravity wave effects generated in the lower atmosphere during the total solar eclipse.

Description: A mid-Holocene climate reconstruction for eastern South America Climate of the Past, 9, 2117-2133, 2013 Author(s): L. F. Prado, I. Wainer, C. M. Chiessi, M.-P. Ledru, and B. Turcq The mid-Holocene (6000 calibrated years before present) is a key period in palaeoclimatology because incoming summer insolation was lower than during the late Holocene in the Southern Hemisphere, whereas the opposite happened in the Northern Hemisphere. However, the effects of the decreased austral summer insolation over South American climate have been poorly discussed by palaeodata syntheses. In addition, only a few of the regional studies have characterised the mid-Holocene climate in South America through a multiproxy approach. Here, we present a multiproxy compilation of mid-Holocene palaeoclimate data for eastern South America. We compiled 120 palaeoclimatological datasets, which were published in 84 different papers. The palaeodata analysed here suggest a water deficit scenario in the majority of eastern South America during the mid-Holocene if compared to the late Holocene, with the exception of northeastern Brazil. Low mid-Holocene austral summer insolation caused a reduced land–sea temperature contrast and hence a weakened South American monsoon system circulation. This scenario is represented by a decrease in precipitation over the South Atlantic Convergence Zone area, saltier conditions along the South American continental margin, and lower lake levels.

Description: Glacial fluctuations of the Indian monsoon and their relationship with North Atlantic climate: new data and modelling experiments Climate of the Past, 9, 2135-2151, 2013 Author(s): C. Marzin, N. Kallel, M. Kageyama, J.-C. Duplessy, and P. Braconnot Several paleoclimate records such as from Chinese loess, speleothems or upwelling indicators in marine sediments present large variations of the Asian monsoon system during the last glaciation. Here, we present a new record from the northern Andaman Sea (core MD77-176) which shows the variations of the hydrological cycle of the Bay of Bengal. The high-resolution record of surface water δ 18 O dominantly reflects salinity changes and displays large millennial-scale oscillations over the period 40 000 to 11 000 yr BP. Their timing and sequence suggests that events of high (resp. low) salinity in the Bay of Bengal, i.e. weak (resp. strong) Indian monsoon, correspond to cold (resp. warm) events in the North Atlantic and Arctic, as documented by the Greenland ice core record. We use the IPSL_CM4 Atmosphere-Ocean coupled General Circulation Model to study the processes that could explain the teleconnection between the Indian monsoon and the North Atlantic climate. We first analyse a numerical experiment in which such a rapid event in the North Atlantic is obtained under glacial conditions by increasing the freshwater flux in the North Atlantic, which results in a reduction of the intensity of the Atlantic meridional overturning circulation. This freshwater hosing results in a weakening of the Indian monsoon rainfall and circulation. The changes in the continental runoff and local hydrological cycle are responsible for an increase in salinity in the Bay of Bengal. This therefore compares favourably with the new sea water δ 18 O record presented here and the hypothesis of synchronous cold North Atlantic and weak Indian monsoon events. Additional sensitivity experiments are produced with the LMDZ atmospheric model to analyse the teleconnection mechanisms between the North Atlantic and the Indian monsoon. The changes over the tropical Atlantic are shown to be essential in triggering perturbations of the subtropical jet over Africa and Eurasia, that in turn affect the intensity of the Indian monsoon. These relationships are also found to be valid in additional coupled model simulations in which the Atlantic meridional overturning circulation (AMOC) is forced to resume.

Description: Mid-Holocene ocean and vegetation feedbacks over East Asia Climate of the Past, 9, 2153-2171, 2013 Author(s): Z. Tian and D. Jiang Mid-Holocene ocean and vegetation feedbacks over East Asia are investigated by a set of numerical experiments performed with the version 4 of the Community Climate System Model (CCSM4). With reference to the pre-industrial period, most of the mid-Holocene annual and seasonal surface-air temperature and precipitation changes are found to result from a direct response of the atmosphere to insolation forcing, while dynamic ocean and vegetation modulate regional climate of East Asia to some extent. Because of its thermal inertia, the dynamic ocean induced an additional warming of 0.2 K for the annual mean, 0.5 K in winter (December–February), 0.0003 K in summer (June–August), and 1.0 K in autumn (September–November), but a cooling of 0.6 K in spring (March–May) averaged over China, and it counteracted (amplified) the direct effect of insolation forcing for the annual mean and in winter and autumn (spring) for that period. The dynamic vegetation had an area-average impact of no more than 0.4 K on the mid-Holocene annual and seasonal temperatures over China, with an average cooling of 0.2 K for the annual mean. On the other hand, ocean feedback induced a small increase of precipitation in winter (0.04 mm day −1 ) and autumn (0.05 mm day −1 ), but a reduction for the annual mean (0.14 mm day −1 ) and in spring (0.29 mm day −1 ) and summer (0.34 mm day −1 ) over China, while it also suppressed the East Asian summer monsoon rainfall. The effect of dynamic vegetation on the mid-Holocene annual and seasonal precipitation was comparatively small, ranging from −0.03 mm day −1 to 0.06 mm day −1 averaged over China. In comparison, the CCSM4 simulated annual and winter cooling over China agrees with simulations within the Paleoclimate Modeling Intercomparison Project (PMIP), but the results are contrary to the warming reconstructed from multiple proxy data for the mid-Holocene. Ocean feedback narrows this model–data mismatch, whereas vegetation feedback plays an opposite role but with a level of uncertainty.

Description: Water-soluble organic carbon in snow and ice deposited at Alpine, Greenland, and Antarctic sites: a critical review of available data and their atmospheric relevance Climate of the Past, 9, 2195-2211, 2013 Author(s): M. Legrand, S. Preunkert, B. Jourdain, J. Guilhermet, X. Fa{ï}n, I. Alekhina, and J. R. Petit While it is now recognized that organic matter dominates the present-day atmospheric aerosol load over continents, its sources remain poorly known. The studies of organic species or organic fractions trapped in ice cores may help to overcome this lack of knowledge. Available data on the dissolved (or total) organic carbon (DOC or TOC) content of snow and ice often appear largely inconsistent, and, until now, no critical review has been conducted to understand the causes of these inconsistencies. To draw a more consistent picture of the organic carbon amount present in solid precipitation that accumulates on cold glaciers, we here review available data and, when needed, complete the data set with analyses of selected samples. The different data sets are then discussed by considering the age (modern versus pre-industrial, Holocene versus Last glacial Maximum) and type (surface snow, firn, or ice) of investigated samples, the deployed method, and the applied contamination control. Finally, the OC (DOC or TOC) levels of Antarctic, Greenland, and Alpine ice cores are compared and discussed with respect to natural (biomass burning, vegetation emissions) and anthropogenic sources (fossil fuel combustion) contributing to atmospheric OC aerosol.

Description: Application of nonlinear autoregressive moving average exogenous input models to geospace: advances in understanding and space weather forecasts Annales Geophysicae, 31, 1579-1589, 2013 Author(s): R. J. Boynton, M. A. Balikhin, S. A. Billings, and O. A. Amariutei The nonlinear autoregressive moving average with exogenous inputs (NARMAX) system identification technique is applied to various aspects of the magnetospheres dynamics. It is shown, from an example system, how the inputs to a system can be found from the error reduction ratio (ERR) analysis, a key concept of the NARMAX approach. The application of the NARMAX approach to the Dst (disturbance storm time) index and the electron fluxes at geostationary Earth orbit (GEO) are reviewed, revealing new insight into the physics of the system. The review of studies into the Dst index illustrate how the NARMAX approach is able to find a coupling function for the Dst index from data, which was then analytically justified from first principles. While the review of the electron flux demonstrates how NARMAX is able to reveal new insight into the physics of the acceleration and loss processes within the radiation belt.

Description: Major dust events in Europe during marine isotope stage 5 (130–74 ka): a climatic interpretation of the "markers" Climate of the Past, 9, 2213-2230, 2013 Author(s): D.-D. Rousseau, M. Ghil, G. Kukla, A. Sima, P. Antoine, M. Fuchs, C. Hatté, F. Lagroix, M. Debret, and O. Moine At present, major dust storms are occurring at mid-latitudes in the Middle East and Asia, as well as at low latitudes in Northern Africa and in Australia. Western Europe, though, does not experience such dramatic climate events, except for some African dust reaching it from the Sahara. This modern situation is of particular interest, in the context of future climate projections, since the present interglacial is usually interpreted, in this context, as an analog of the warm Eemian interval. European terrestrial records show, however, major dust events during the penultimate interglacial and early glacial. These events are easily observed in loess records by their whitish-color deposits, which lie above and below dark chernozem paleosols in Central European records of Marine Isotope Stage (MIS) 5 age. We describe here the base of the Dolni Vestonice (DV) loess sequence, Czech Republic, as the reference of such records. The dust is deposited during intervals that are characterized by poor vegetation – manifested by high δ 13 C values and low magnetic susceptibility – while the fine sand and clay in the deposits shows grain sizes that are clearly different from the overlying pleniglacial loess deposits. Some of these dust events have been previously described as "Markers" or Marker Silts (MS) by one of us (G. Kukla), and are dated at about 111–109 ka and 93–92 ka, with a third and last one slightly visible at about 75–73 ka. Other events correspond to the loess material of Kukla's cycles, and are described as eolian silts (ES); they are observed in the same DV sequence and are dated at about 106–105 ka, 88–86 ka, and 78.5–77 ka. These dates are determined by considering the OSL ages with their errors measured on the studied sequence, and the comparison with Greenland ice-core and European speleothem chronologies. The fine eolian deposits mentioned above, MS as well as ES, correspond to short events that lasted about 2 ka; they are synchronous with re-advances of the polar front over the North Atlantic, as observed in marine sediment cores. These deposits also correlate with important changes observed in European vegetation. Some ES and MS events appear to be coeval with significant dust peaks recorded in the Greenland ice cores, while others are not. This decoupling between the European eolian and Greenland dust depositions is of considerable interest, as it differs from the fully glacial situation, in which the Eurasian loess sedimentation mimics the Greenland dust record. Previous field observations supported an interpretation of MS events as caused by continental dust storms. We show here, by a comparison with speleothems of the same age found in the northern Alps, that different atmospheric-circulation modes seem to be responsible for the two categories of dust events, MS vs. ES.

Description: Climate, people, fire and vegetation: new insights into vegetation dynamics in the Eastern Mediterranean since the 1st century AD Climate of the Past, 9, 57-87, 2013 Author(s): J. Bakker, E. Paulissen, D. Kaniewski, J. Poblome, V. De Laet, G. Verstraeten, and M. Waelkens Anatolia forms a bridge between Europe, Africa and Asia and is influenced by all three continents in terms of climate, vegetation and human civilisation. Unfortunately, well-dated palynological records focussing on the period from the end of the classical Roman period until subrecent times are rare for Anatolia and completely absent for southwest Turkey, resulting in a lacuna in knowledge concerning the interactions of climatic change, human impact, and environmental change in this important region. Two well-dated palaeoecological records from the Western Taurus Mountains, Turkey, provide a first relatively detailed record of vegetation dynamics from late Roman times until the present in SW Turkey. Combining pollen, non-pollen palynomorphs, charcoal, sedimentological, archaeological data, and newly developed multivariate numerical analyses allows for the disentangling of climatic and anthropogenic influences on vegetation change. Results show changes in both the regional pollen signal as well as local soil sediment characteristics match shifts in regional climatic conditions. Both climatic as well as anthropogenic change had a strong influence on vegetation dynamics and land use. A moist environmental trend during the late-3rd century caused an increase in marshes and wetlands in the moister valley floors, limiting possibilities for intensive crop cultivation at such locations. A mid-7th century shift to pastoralism coincided with a climatic deterioration as well as the start of Arab incursions into the region, the former driving the way in which the vegetation developed afterwards. Resurgence in agriculture was observed in the study during the mid-10th century AD, coinciding with the Medieval Climate Anomaly. An abrupt mid-12th century decrease in agriculture is linked to socio-political change, rather than the onset of the Little Ice Age. Similarly, gradual deforestation occurring from the 16th century onwards has been linked to changes in land use during Ottoman times. The pollen data reveal that a fast rise in Pinus pollen after the end of the Beyşehir Occupation Phase need not always occur. The notion of high Pinus pollen percentages indicating an open landscape incapable of countering the influx of pine pollen is also deemed unrealistic. While multiple fires occurred in the region through time, extended fire periods, as had occurred during the Bronze Age and Beyşehir Occupation Phase, did not occur, and no signs of local fire activity were observed. Fires were never a major influence on vegetation dynamics. While no complete overview of post-BO Phase fire events can be presented, the available data indicates that fires in the vicinity of Gravgaz may have been linked to anthropogenic activity in the wider surroundings of the marsh. Fires in the vicinity of Bereket appeared to be linked to increased abundance of pine forests. There was no link with specifically wet or dry environmental conditions at either site. While this study reveals much new information concerning the impact of climate change and human occupation on the environment, more studies from SW Turkey are required in order to properly quantify the range of the observed phenomena and the magnitude of their impacts.

Description: Support vector machines for TEC seismo-ionospheric anomalies detection Annales Geophysicae, 31, 173-186, 2013 Author(s): M. Akhoondzadeh Using time series prediction methods, it is possible to pursue the behaviors of earthquake precursors in the future and to announce early warnings when the differences between the predicted value and the observed value exceed the predefined threshold value. Support Vector Machines (SVMs) are widely used due to their many advantages for classification and regression tasks. This study is concerned with investigating the Total Electron Content (TEC) time series by using a SVM to detect seismo-ionospheric anomalous variations induced by the three powerful earthquakes of Tohoku (11 March 2011), Haiti (12 January 2010) and Samoa (29 September 2009). The duration of TEC time series dataset is 49, 46 and 71 days, for Tohoku, Haiti and Samoa earthquakes, respectively, with each at time resolution of 2 h. In the case of Tohoku earthquake, the results show that the difference between the predicted value obtained from the SVM method and the observed value reaches the maximum value (i.e., 129.31 TECU) at earthquake time in a period of high geomagnetic activities. The SVM method detected a considerable number of anomalous occurrences 1 and 2 days prior to the Haiti earthquake and also 1 and 5 days before the Samoa earthquake in a period of low geomagnetic activities. In order to show that the method is acting sensibly with regard to the results extracted during nonevent and event TEC data, i.e., to perform some null-hypothesis tests in which the methods would also be calibrated, the same period of data from the previous year of the Samoa earthquake date has been taken into the account. Further to this, in this study, the detected TEC anomalies using the SVM method were compared to the previous results (Akhoondzadeh and Saradjian, 2011; Akhoondzadeh, 2012) obtained from the mean, median, wavelet and Kalman filter methods. The SVM detected anomalies are similar to those detected using the previous methods. It can be concluded that SVM can be a suitable learning method to detect the novelty changes of a nonlinear time series such as variations of earthquake precursors.

Description: Multi-periodic climate dynamics: spectral analysis of long-term instrumental and proxy temperature records Climate of the Past, 9, 447-452, 2013 Author(s): H.-J. Lüdecke, A. Hempelmann, and C. O. Weiss The longest six instrumental temperature records of monthly means reach back maximally to 1757 AD and were recorded in Europe. All six show a V-shape, with temperature drop in the 19th and rise in the 20th century. Proxy temperature time series of Antarctic ice cores show this same characteristic shape, indicating this pattern as a global phenomenon. We used the mean of the six instrumental records for analysis by discrete Fourier transform (DFT), wavelets, and the detrended fluctuation analysis (DFA). For comparison, a stalagmite record was also analyzed by DFT. The harmonic decomposition of the abovementioned mean shows only six significant frequencies above periods over 30 yr. The Pearson correlation between the mean, smoothed by a 15-yr running average (boxcar) and the reconstruction using the six significant frequencies, yields r = 0.961. This good agreement has a 〉 99.9% confidence level confirmed by Monte Carlo simulations. It shows that the climate dynamics is governed at present by periodic oscillations. We find indications that observed periodicities result from intrinsic dynamics.

Description: Magnetosheath dynamic pressure enhancements: occurrence and typical properties Annales Geophysicae, 31, 319-331, 2013 Author(s): M. O. Archer and T. S. Horbury The first comprehensive statistical study of large-amplitude (〉 100%) transient enhancements of the magnetosheath dynamic pressure reveals events of up to ~ 15 times the ambient dynamic pressure with durations up to 3 min and an average duration of around 30 s, predominantly downstream of the quasi-parallel shock. The dynamic pressure transients are most often dominated by velocity increases along with a small fractional increase in the density, though the velocity is generally only deflected by a few degrees. Superposed wavelet transforms of the magnetic field show that, whilst most enhancements exhibit changes in the magnetosheath magnetic field, the majority are not associated with changes in the Interplanetary Magnetic Field (IMF). However, there is a minority of enhancements that do appear to be associated with solar wind discontinuities which cannot be explained simply by random events. In general, it is found that during periods of magnetosheath dynamic pressure enhancements the IMF is steadier than usual. This suggests that a stable foreshock and hence foreshock structures or processes may be important in the generation of the majority of magnetosheath dynamic pressure enhancements.

Description: The East Asian Summer Monsoon at mid-Holocene: results from PMIP3 simulations Climate of the Past, 9, 453-466, 2013 Author(s): W. Zheng, B. Wu, J. He, and Y. Yu Ten Coupled General Circulation Models (CGCMs) participated in the third phase of Paleoclimate Modelling Intercomparison Project (PMIP3) are assessed for the East Asian Summer Monsoon (EASM) in both the pre-Industrial (PI, 0 ka) and mid-Holocene (MH, 6 ka) simulations. Results show that the PMIP3 model median captures well the large-scale characteristics of the EASM, including the two distinct features of the Meiyu rainbelt and the stepwise meridional displacement of the monsoonal rainbelt. At mid-Holocene, the PMIP3 model median shows significant warming (cooling) during boreal summer (winter) over Eurasia continent that are dominated by the changes of insolation. However, the PMIP3 models fail to simulate a warmer annual mean and winter surface air temperature (TAS) over eastern China as derived from proxy records. The EASM at MH are featured by the changes of large-scale circulation over Eastern China while the changes of precipitation are not significant over its sub-domains of the Southern China and the lower reaches of Yangzi River. The inter-model differences for the monsoon precipitation can be associated with different configurations of the changes in large-scale circulation and the water vapour content, of which the former determines the sign of precipitation changes. The large model spread for the TAS over Tibetan Plateau has a positive relationship with the precipitation in the lower reaches of Yangzi River, yet this relationship does not apply to those PMIP3 models in which the monsoonal precipitation is more sensitive to the changes of large-scale circulation. Except that the PMIP3 model median captured the warming of annual mean TAS over Tibetan Plateau, no significant improvements can be concluded when compared with the PMIP2 models results.

Description: Response of polar mesosphere summer echoes to geomagnetic disturbances in the Southern and Northern Hemispheres: the importance of nitric oxide Annales Geophysicae, 31, 333-347, 2013 Author(s): S. Kirkwood, E. Belova, P. Dalin, M. Mihalikova, D. Mikhaylova, D. Murtagh, H. Nilsson, K. Satheesan, J. Urban, and I. Wolf The relationship between polar mesosphere summer echoes (PMSE) and geomagnetic disturbances (represented by magnetic K indices) is examined. Calibrated PMSE reflectivities for the period May 2006–February 2012 are used from two 52.0/54.5 MHz radars located in Arctic Sweden (68° N, geomagnetic latitude 65°) and at two different sites in Queen Maud Land, Antarctica (73°/72° S, geomagnetic latitudes 62°/63°). In both the Northern Hemisphere (NH) and the Southern Hemisphere (SH) there is a strong increase in mean PMSE reflectivity between quiet and disturbed geomagnetic conditions. Mean volume reflectivities are slightly lower at the SH locations compared to the NH, but the position of the peak in the lognormal distribution of PMSE reflectivities is close to the same at both NH and SH locations, and varies only slightly with magnetic disturbance level. Differences between the sites, and between geomagnetic disturbance levels, are primarily due to differences in the high-reflectivity tail of the distribution. PMSE occurrence rates are essentially the same at both NH and SH locations during most of the PMSE season when a sufficiently low detection threshold is used so that the peak in the lognormal distribution is included. When the local-time dependence of the PMSE response to geomagnetic disturbance level is considered, the response in the NH is found to be immediate at most local times, but delayed by several hours in the afternoon sector and absent in the early evening. At the SH sites, at lower magnetic latitude, there is a delayed response (by several hours) at almost all local times. At the NH (auroral zone) site, the dependence on magnetic disturbance is highest during evening-to-morning hours. At the SH (sub-auroral) sites the response to magnetic disturbance is weaker but persists throughout the day. While the immediate response to magnetic activity can be qualitatively explained by changes in electron density resulting from energetic particle precipitation, the delayed response can largely be explained by changes in nitric oxide concentrations. Observations of nitric oxide concentration at PMSE heights by the Odin satellite support this hypothesis. Sensitivity to geomagnetic disturbances, including nitric oxide produced during these disturbances, can explain previously reported differences between sites in the auroral zone and those at higher or lower magnetic latitudes. The several-day lifetime of nitric oxide can also explain earlier reported discrepancies between high correlations for average conditions (year-by-year PMSE reflectivities and K indices) and low correlations for minute-to-day timescales.

Description: Observations of poleward-propagating large-scale traveling ionospheric disturbances in southern China Annales Geophysicae, 31, 377-385, 2013 Author(s): F. Ding, W. Wan, B. Ning, B. Zhao, Q. Li, Y. Wang, L. Hu, R. Zhang, and B. Xiong We report here on two cases of poleward-propagating large-scale traveling ionospheric disturbances (LSTIDs) in China during a medium-scale storm between 27 May and 1 June 2011. The observations were conducted by making use of the Global Positioning System network and ionosondes in China and Southeast Asia. One northeastward-propagating LSTID occurred on the morning of 30 May, while the other was observed during the nighttime of 1 June. Both poleward-traveling LSTIDs occurred during the storm's recovery phase in southern China's low-latitude region (geomagnetic latitude ~ 7.3–24° N) and experienced severe dissipation during their propagation from south to north. Although the initial relative amplitude of the nighttime LSTID was ~ 60% larger than that of the morning event, the nighttime event dissipated more quickly than the morning event because of a strong nighttime enhancement in background total electronic content (TEC) during storm time, which led to strong ion-drag dissipation during the evening. The poleward-propagating LSTIDs exhibit a narrower latitudinal range, a smaller amplitude, and a slightly higher elevation compared with the equatorward-moving LSTIDs observed in the same region. Given these features, the poleward-propagating LSTIDs were likely excited by some local source near southern China. Excitation of secondary LSTIDs during the dissipation of some primary medium-scale disturbances from the lower atmosphere is a possible mechanism.

Description: Paleohydrology reconstruction and Holocene climate variability in the South Adriatic Sea Climate of the Past, 9, 499-515, 2013 Author(s): G. Siani, M. Magny, M. Paterne, M. Debret, and M. Fontugne Holocene paleohydrology reconstruction is derived combining planktonic and benthic stable oxygen and carbon isotopes, sea surface temperatures (SSTs) and oxygen isotope composition of seawater (δ 18 O w ) from a high sedimentation core collected in the South Adriatic Sea (SAS). Core chronology is based on 10 AMS 14 C measures on planktonic foraminifera and tephra layers. Results reveal two contrasted paleohydrological periods that reflect (i) a marked lowering of δ 18 O w /salinity during the early to mid-Holocene (11.5 ka to 6.3 ka), including the two-step sapropel S1 deposition, followed during the mid- to upper Holocene by (ii) a prevailing period of increased salinity and enhanced arid conditions in the South Adriatic Basin. Superimposed on these trends, short-term centennial-scale hydrological events punctuated the Holocene period in the SAS. During the early to mid-Holocene, two main SST coolings together with prominent δ 18 O w /salinity lowering delineate the sapropel S1 interruption and the post-sapropel phase between 7.3 to 6.3 ka. After 6 ka, centennial-scale δ 18 O w and G. bulloides δ 13 C lowering, mostly centered between 3 to 0.6 ka, reflect short-term hydrological changes related to more intensive runoff of the Po and/or Apennine rivers. These short-term events, even of lesser amplitude compared to the early to mid-Holocene period, may have induced a lowering of sea surface density and consequently reduced and/or inhibited the formation of deep bottom waters in the SAS. Comparison of the emerging centennial- to millennial-scale hydrological record with previous climatic records from the central Mediterranean area and north of the Alps reveal possible synchronicities (within the radiocarbon-dating uncertainty) between phases of lower salinity in the SAS and periods of wetter climatic conditions around the north-central Adriatic Sea. Finally, wavelet analyses provide new clues about the potential origin of climate variability in the SAS, confirming the evidence for a mid-Holocene transition in the central Mediterranean climate and the dominance of a ~1670-yr periodicity after 6 ka, reflecting a plausible connection with the North Atlantic climate system.

Description: Variability of the ocean heat content during the last millennium – an assessment with the ECHO-g Model Climate of the Past, 9, 547-565, 2013 Author(s): P. Ortega, M. Montoya, F. González-Rouco, H. Beltrami, and D. Swingedouw Studies addressing climate variability during the last millennium generally focus on variables with a direct influence on climate variability, like the fast thermal response to varying radiative forcing, or the large-scale changes in atmospheric dynamics (e.g. North Atlantic Oscillation). The ocean responds to these variations by slowly integrating in depth the upper heat flux changes, thus producing a delayed influence on ocean heat content (OHC) that can later impact low frequency SST (sea surface temperature) variability through reemergence processes. In this study, both the externally and internally driven variations of the OHC during the last millennium are investigated using a set of fully coupled simulations with the ECHO-G (coupled climate model ECHAMA4 and ocean model HOPE-G) atmosphere–ocean general circulation model (AOGCM). When compared to observations for the last 55 yr, the model tends to overestimate the global trends and underestimate the decadal OHC variability. Extending the analysis back to the last one thousand years, the main impact of the radiative forcing is an OHC increase at high latitudes, explained to some extent by a reduction in cloud cover and the subsequent increase of short-wave radiation at the surface. This OHC response is dominated by the effect of volcanism in the preindustrial era, and by the fast increase of GHGs during the last 150 yr. Likewise, salient impacts from internal climate variability are observed at regional scales. For instance, upper temperature in the equatorial Pacific is controlled by ENSO (El Niño Southern Oscillation) variability from interannual to multidecadal timescales. Also, both the Pacific Decadal Oscillation (PDO) and the Atlantic Multidecadal Oscillation (AMO) modulate intermittently the interdecadal OHC variability in the North Pacific and Mid Atlantic, respectively. The NAO, through its influence on North Atlantic surface heat fluxes and convection, also plays an important role on the OHC at multiple timescales, leading first to a cooling in the Labrador and Irminger seas, and later on to a North Atlantic warming, associated with a delayed impact on the AMO.

Description: Modulation of Late Cretaceous and Cenozoic climate by variable drawdown of atmospheric p CO 2 from weathering of basaltic provinces on continents drifting through the equatorial humid belt Climate of the Past, 9, 525-546, 2013 Author(s): D. V. Kent and G. Muttoni The small reservoir of carbon dioxide in the atmosphere ( p CO 2 ) that modulates climate through the greenhouse effect reflects a delicate balance between large fluxes of sources and sinks. The major long-term source of CO 2 is global outgassing from sea-floor spreading, subduction, hotspot activity, and metamorphism; the ultimate sink is through weathering of continental silicates and deposition of carbonates. Most carbon cycle models are driven by changes in the source flux scaled to variable rates of ocean floor production, but ocean floor production may not be distinguishable from being steady since 180 Ma. We evaluate potential changes in sources and sinks of CO 2 for the past 120 Ma in a paleogeographic context. Our new calculations show that decarbonation of pelagic sediments by Tethyan subduction contributed only modestly to generally high p CO 2 levels from the Late Cretaceous until the early Eocene, and thus shutdown of this CO 2 source with the collision of India and Asia at the early Eocene climate optimum at around 50 Ma was inadequate to account for the large and prolonged decrease in p CO 2 that eventually allowed the growth of significant Antarctic ice sheets by around 34 Ma. Instead, variation in area of continental basalt terranes in the equatorial humid belt (5° S–5° N) seems to be a dominant factor controlling how much CO 2 is retained in the atmosphere via the silicate weathering feedback. The arrival of the highly weatherable Deccan Traps in the equatorial humid belt at around 50 Ma was decisive in initiating the long-term slide to lower atmospheric p CO 2 , which was pushed further down by the emplacement of the 30 Ma Ethiopian Traps near the equator and the southerly tectonic extrusion of SE Asia, an arc terrane that presently is estimated to account for 1/4 of CO 2 consumption from all basaltic provinces that account for ~1/3 of the total CO 2 consumption by continental silicate weathering (Dessert et al., 2003). A negative climate-feedback mechanism that (usually) inhibits the complete collapse of atmospheric p CO 2 is the accelerating formation of thick cation-deficient soils that retard chemical weathering of the underlying bedrock. Nevertheless, equatorial climate seems to be relatively insensitive to p CO 2 greenhouse forcing and thus with availability of some rejuvenating relief as in arc terranes or thick basaltic provinces, silicate weathering in this venue is not subject to a strong negative feedback, providing an avenue for ice ages. The safety valve that prevents excessive atmospheric p CO 2 levels is the triggering of silicate weathering of continental areas and basaltic provinces in the temperate humid belt. Excess organic carbon burial seems to have played a negligible role in atmospheric p CO 2 over the Late Cretaceous and Cenozoic.

Description: A new method for solving the MHD equations in the magnetosheath Annales Geophysicae, 31, 419-437, 2013 Author(s): C. Nabert, K.-H. Glassmeier, and F. Plaschke We present a new analytical method to derive steady-state magnetohydrodynamic (MHD) solutions of the magnetosheath in different levels of approximation. With this method, we calculate the magnetosheath's density, velocity, and magnetic field distribution as well as its geometry. Thereby, the solution depends on the geomagnetic dipole moment and solar wind conditions only. To simplify the representation, we restrict our model to northward IMF with the solar wind flow along the stagnation streamline. The sheath's geometry, with its boundaries, bow shock and magnetopause, is determined self-consistently. Our model is stationary and time relaxation has not to be considered as in global MHD simulations. Our method uses series expansion to transfer the MHD equations into a new set of ordinary differential equations. The number of equations is related to the level of approximation considered including different physical processes. These equations can be solved numerically; however, an analytical approach for the lowest-order approximation is also presented. This yields explicit expressions, not only for the flow and field variations but also for the magnetosheath thickness, depending on the solar wind parameters. Results are compared to THEMIS data and offer a detailed explanation of, e.g., the pile-up process and the corresponding plasma depletion layer, the bow shock and magnetopause geometry, the magnetosheath thickness, and the flow deceleration.

Description: Investigation of a mesospheric bore event over northern China Annales Geophysicae, 31, 409-418, 2013 Author(s): Q. Li, J. Xu, J. Yue, X. Liu, W. Yuan, B. Ning, S. Guan, and J. P. Younger A mesospheric bore event was observed using an OH all-sky airglow imager (ASAI) at Xinglong (40.2° N, 117.4° E), in northern China, on the night of 8–9 January 2011. Simultaneous observations by a Doppler meteor radar, a broadband sodium lidar, and TIMED/SABER OH intensity and temperature measurements are used to investigate the characteristics and environment of the bore propagation and the possible relations with the Na density perturbations. The bore propagated from northeast to southwest and divided the sky into bright and dark halves. The calculations show that the bore has an average phase velocity of 68 m s −1 . The crests following the bore have a horizontal wavelength of ~ 22 km. These parameters are consistent with the hydraulic jump theory proposed by Dewan and Picard, as well as the previous bore reports. Simultaneous wind measurements from the Doppler meteor radar at Shisanling (40.3° N, 116.2° E) and temperature data from SABER on board the TIMED satellite are used to characterize the propagating environment of the bore. The result shows that a thermal-Doppler duct exists near the OH layer that supports the horizontal propagation of the bore. Simultaneous Na lidar observations at Yanqing (40.4° N, 116.0° E) suggest that there is a downward displacement of Na density during the passage of the mesospheric bore event.

Description: The ultra-fast Kelvin waves in the equatorial ionosphere: observations and modeling Annales Geophysicae, 31, 209-215, 2013 Author(s): A. N. Onohara, I. S. Batista, and H. Takahashi The main purpose of this study is to investigate the vertical coupling between the mesosphere and lower thermosphere (MLT) region and the ionosphere through ultra-fast Kelvin (UFK) waves in the equatorial atmosphere. The effect of UFK waves on the ionospheric parameters was estimated using an ionospheric model which calculates electrostatic potential in the E-region and solves coupled electrodynamics of the equatorial ionosphere in the E- and F-regions. The UFK wave was observed in the South American equatorial region during February–March 2005. The MLT wind data obtained by meteor radar at São João do Cariri (7.5° S, 37.5° W) and ionospheric F-layer bottom height ( h 'F) observed by ionosonde at Fortaleza (3.9° S; 38.4° W) were used in order to calculate the wave characteristics and amplitude of oscillation. The simulation results showed that the combined electrodynamical effect of tides and UFK waves in the MLT region could explain the oscillations observed in the ionospheric parameters.

Description: Payload charging events in the mesosphere and their impact on Langmuir type electric probes Annales Geophysicae, 31, 187-196, 2013 Author(s): T. A. Bekkeng, A. Barjatya, U.-P. Hoppe, A. Pedersen, J. I. Moen, M. Friedrich, and M. Rapp Three sounding rockets were launched from Andøya Rocket Range in the ECOMA campaign in December 2010. The aim was to study the evolution of meteoric smoke particles during a major meteor shower. Of the various instruments onboard the rocket payload, this paper presents the data from a multi-Needle Langmuir Probe (m-NLP) and a charged dust detector. The payload floating potential, as observed using the m-NLP instrument, shows charging events on two of the three flights. These charging events cannot be explained using a simple charging model, and have implications towards the use of fixed bias Langmuir probes on sounding rockets investigating mesospheric altitudes. We show that for a reliable use of a single fixed bias Langmuir probe as a high spatial resolution relative density measurement, each payload should also carry an additional instrument to measure payload floating potential, and an instrument that is immune to spacecraft charging and measures absolute plasma density.

Description: Electric potential differences across auroral generator interfaces Annales Geophysicae, 31, 251-261, 2013 Author(s): J. De Keyser and M. Echim Strong localized high-altitude auroral electric fields, such as those observed by Cluster, are often associated with magnetospheric interfaces. The type of high-altitude electric field profile (monopolar, bipolar, or more complicated) depends on the properties of the plasmas on either side of the interface, as well as on the total electric potential difference across the structure. The present paper explores the role of this cross-field electric potential difference in the situation where the interface is a tangential discontinuity. A self-consistent Vlasov description is used to determine the equilibrium configuration for different values of the transverse potential difference. A major observation is that there exist limits to the potential difference, beyond which no equilibrium configuration of the interface can be sustained. It is further demonstrated how the plasma densities and temperatures affect the type of electric field profile in the transition, with monopolar electric fields appearing primarily when the temperature contrast is large. These findings strongly support the observed association of monopolar fields with the plasma sheet boundary. The role of shear flow tangent to the interface is also examined.

Description: Variability of ionospheric TEC during solar and geomagnetic minima (2008 and 2009): external high speed stream drivers Annales Geophysicae, 31, 263-276, 2013 Author(s): O. P. Verkhoglyadova, B. T. Tsurutani, A. J. Mannucci, M. G. Mlynczak, L. A. Hunt, and T. Runge We study solar wind–ionosphere coupling through the late declining phase/solar minimum and geomagnetic minimum phases during the last solar cycle (SC23) – 2008 and 2009. This interval was characterized by sequences of high-speed solar wind streams (HSSs). The concomitant geomagnetic response was moderate geomagnetic storms and high-intensity, long-duration continuous auroral activity (HILDCAA) events. The JPL Global Ionospheric Map (GIM) software and the GPS total electron content (TEC) database were used to calculate the vertical TEC (VTEC) and estimate daily averaged values in separate latitude and local time ranges. Our results show distinct low- and mid-latitude VTEC responses to HSSs during this interval, with the low-latitude daytime daily averaged values increasing by up to 33 TECU (annual average of ~20 TECU) near local noon (12:00 to 14:00 LT) in 2008. In 2009 during the minimum geomagnetic activity (MGA) interval, the response to HSSs was a maximum of ~30 TECU increases with a slightly lower average value than in 2008. There was a weak nighttime ionospheric response to the HSSs. A well-studied solar cycle declining phase interval, 10–22 October 2003, was analyzed for comparative purposes, with daytime low-latitude VTEC peak values of up to ~58 TECU (event average of ~55 TECU). The ionospheric VTEC changes during 2008–2009 were similar but ~60% less intense on average. There is an evidence of correlations of filtered daily averaged VTEC data with Ap index and solar wind speed. We use the infrared NO and CO 2 emission data obtained with SABER on TIMED as a proxy for the radiation balance of the thermosphere. It is shown that infrared emissions increase during HSS events possibly due to increased energy input into the auroral region associated with HILDCAAs. The 2008–2009 HSS intervals were ~85% less intense than the 2003 early declining phase event, with annual averages of daily infrared NO emission power of ~ 3.3 × 10 10 W and 2.7 × 10 10 W in 2008 and 2009, respectively. The roles of disturbance dynamos caused by high-latitude winds (due to particle precipitation and Joule heating in the auroral zones) and of prompt penetrating electric fields (PPEFs) in the solar wind–ionosphere coupling during these intervals are discussed. A correlation between geoeffective interplanetary electric field components and HSS intervals is shown. Both PPEF and disturbance dynamo mechanisms could play important roles in solar wind–ionosphere coupling during prolonged (up to days) external driving within HILDCAA intervals.

Description: Palynological evidence for gradual vegetation and climate changes during the African Humid Period termination at 13°N from a Mega-Lake Chad sedimentary sequence Climate of the Past, 9, 223-241, 2013 Author(s): P. G. C. Amaral, A. Vincens, J. Guiot, G. Buchet, P. Deschamps, J.-C. Doumnang, and F. Sylvestre Located at the transition between the Saharan and Sahelian zones, at the center of one of the largest endorheic basins, Lake Chad is ideally located to record regional environmental changes that occurred in the past. However, until now, no Holocene archive was directly cored in this lake. In this paper, we present pollen data from the first sedimentary sequence collected in Lake Chad (13° N; 14° E; Sahel region). Dated between ca. 6700 and ca. 5000 cal yr BP, this record is continuous and encompasses part of the termination of the African Humid Period (AHP). Vegetation reconstructions are based on standard analyses of pollen diagrams and are strengthened by quantitative approaches. Potential biomes are reconstructed using the biomization method and mean annual precipitation ( P ann ) is estimated using the modern analogues technique. Results show that, between ca. 6700 and ca. 6050 cal yr BP, a vegetation close to humid woodland or humid savanna, including elements currently found further southward, thrived in the vicinity of the Mega-Lake Chad in place of the modern dry woodland, steppe and desert vegetation. At the same time, montane forest populations extended further southward on the Adamawa Plateau. The high abundance of lowland humid pollen taxa, particularly of Uapaca , is interpreted as the result of a northward migration of the corresponding plants during the AHP. This preferential zonal occurrence of these taxa in Lake Chad Basin (LCB) (rather than extrazonal) is driven by more humid local and regional climate conditions at this latitude, as shown by mean P ann estimated values of ca. 800 (−400/+700) mm during this period. However, we cannot rule out that an increase of the Chari–Logone inputs into the Mega-Lake Chad might have also contributed to control the abundance of these taxa. Changes in the structure and floristic composition of the vegetation towards more open and drier formations occurred after ca. 6050 cal yr BP, following a decrease in mean P ann estimates to approximately 600 (−230/+600) mm. But, the constant significant presence of lowland humid taxa until ca. 5000 cal yr BP, contemporaneous with a slight increase in steppic taxa, demonstrates that at this date, the modern vegetation was not yet established in the vicinity of Lake Chad. Our data indicate that between ca. 6700 and ca. 5000 cal yr BP vegetation and climate changes must have occurred progressively, but that century-scale climate variability was superimposed on this long-term mid-Holocene drying trend as observed around ca. 6300 cal yr BP, where pollen data indicate more humid conditions.

Description: What could have caused pre-industrial biomass burning emissions to exceed current rates? Climate of the Past, 9, 289-306, 2013 Author(s): G. R. van der Werf, W. Peters, T. T. van Leeuwen, and L. Giglio Recent studies based on trace gas mixing ratios in ice cores and charcoal data indicate that biomass burning emissions over the past millennium exceeded contemporary emissions by up to a factor of 4 for certain time periods. This is surprising because various sources of biomass burning are linked with population density, which has increased over the past centuries. We have analysed how emissions from several landscape biomass burning sources could have fluctuated to yield emissions that are in correspondence with recent results based on ice core mixing ratios of carbon monoxide (CO) and its isotopic signature measured at South Pole station (SPO). Based on estimates of contemporary landscape fire emissions and the TM5 chemical transport model driven by present-day atmospheric transport and OH concentrations, we found that CO mixing ratios at SPO are more sensitive to emissions from South America and Australia than from Africa, and are relatively insensitive to emissions from the Northern Hemisphere. We then explored how various landscape biomass burning sources may have varied over the past centuries and what the resulting emissions and corresponding CO mixing ratio at SPO would be, using population density variations to reconstruct sources driven by humans (e.g., fuelwood burning) and a new model to relate savanna emissions to changes in fire return times. We found that to match the observed ice core CO data, all savannas in the Southern Hemisphere had to burn annually, or bi-annually in combination with deforestation and slash and burn agriculture exceeding current levels, despite much lower population densities and lack of machinery to aid the deforestation process. While possible, these scenarios are unlikely and in conflict with current literature. However, we do show the large potential for increased emissions from savannas in a pre-industrial world. This is mainly because in the past, fuel beds were probably less fragmented compared to the current situation; satellite data indicates that the majority of savannas have not burned in the past 10 yr, even in Africa, which is considered "the burning continent". Although we have not considered increased charcoal burning or changes in OH concentrations as potential causes for the elevated CO concentrations found at SPO, it is unlikely they can explain the large increase found in the CO concentrations in ice core data. Confirmation of the CO ice core data would therefore call for radical new thinking about causes of variable global fire rates over recent centuries.

Description: Case study of stratospheric gravity waves of convective origin over Arctic Scandinavia – VHF radar observations and numerical modelling Annales Geophysicae, 31, 239-250, 2013 Author(s): A. Réchou, J. Arnault, P. Dalin, and S. Kirkwood Orography is a well-known source of gravity and inertia-gravity waves in the atmosphere. Other sources, such as convection, are also known to be potentially important but the large amplitude of orographic waves over Scandinavia has generally precluded the possibility to study such other sources experimentally in this region. In order to better understand the origin of stratospheric gravity waves observed by the VHF radar ESRAD (Esrange MST radar) over Kiruna, in Arctic Sweden (67.88° N, 21.10° E), observations have been compared to simulations made using the Weather Research and Forecasting model (WRF) with and without the effects of orography and clouds. This case study concerns gravity waves observed from 00:00 UTC on 18 February to 12:00 UTC on 20 February 2007. We focus on the wave signatures in the static stability field and vertical wind deduced from the simulations and from the observations as these are the parameters which are provided by the observations with the best height coverage. As is common at this site, orographic gravity waves were produced over the Scandinavian mountains and observed by the radar. However, at the same time, southward propagation of fronts in the Barents Sea created short-period waves which propagated into the stratosphere and were transported, embedded in the cyclonic winds, over the radar site.

Description: Long-term trends in the ionospheric F2 region with different solar activity indices Annales Geophysicae, 31, 291-303, 2013 Author(s): J. Mielich and J. Bremer A new comprehensive data collection by Damboldt and Suessmann (2012a) with monthly fo F2 and M(3000)F2 median values is an excellent basis for the derivation of long-term trends in the ionospheric F2 region. Ionospheric trends have been derived only for stations with data series of at least 22 years (124 stations with fo F2 data and 113 stations with M(3000)F2 data) using a twofold regression analysis depending on solar and geomagnetic activity. Three main results have been derived: Firstly, it could be shown that the solar 10.7 cm radio flux F10.7 is a better index for the description of the solar activity than the relative solar sunspot number R as well as the solar EUV proxy E10.7. Secondly, the global mean fo F2 and hm F2 trends derived for the interval between 1948 and 2006 are in surprisingly good agreement with model calculations of an increasing atmospheric greenhouse effect (Rishbeth and Roble, 1992). Thirdly, during the years 2007 until 2009, the hm F2 values and to a smaller amount the fo F2 values strongly decrease. The reason for this effect is a reduction of the thermospheric density and ionization due to a markedly reduced solar EUV irradiation and extremely small geomagnetic activity during the solar cycle 23/24 minimum.

Description: High-temperature thermomagnetic properties of vivianite nodules, Lake El'gygytgyn, Northeast Russia Climate of the Past, 9, 433-446, 2013 Author(s): P. S. Minyuk, T. V. Subbotnikova, L. L. Brown, and K. J. Murdock Vivianite, a hydrated iron phosphate, is abundant in sediments of Lake El'gygytgyn, located in the Anadyr Mountains of central Chukotka, northeastern Russia (67°30′ N, 172°05′ E). Magnetic measurements, including mass-specific low-field AC magnetic susceptibility, field-dependent magnetic susceptibility, hysteresis parameters, temperature dependence of the induced magnetization, as well as susceptibility in different heating media, provide ample information on vivianite nodules. Electron microprobe analyses, electron microscopy and energy dispersive spectroscopy were used to identify diagnostic minerals. Vivianite nodules are abundant in both sediments of cold (anoxic) and warm (oxic) stages. Magnetic susceptibility of the nodules varies from 0.78 × 10 −6 m 3 kg −1 to 1.72 × 10 −6 m 3 kg −1 (average = 1.05 × 10 −6 m 3 kg −1 ) and is higher than the susceptibility of sediments from the cold intervals. Magnetic properties of vivianite are due to the respective product of oxidation as well as sediment and mineral inclusions. Three types of curves for high-temperature dependent susceptibility of vivianite indicate different degrees of oxidation and inclusions in the nodules. Vivianite acts as a reductant and reduces hematite to magnetite and masks the goethite–hematite transition during heating. Heating vivianite and sulfur mixtures stimulates the formation of monoclinic pyrrhotite. An additive of arsenic inhibits the formation of magnetite prior to its Curie temperature. Heating selective vivianite and pyrite mixtures leads to formation of several different minerals – magnetite, monoclinic pyrrhotite, and hexagonal pyrrhotite, and makes it difficult to interpret the thermomagnetic curves.

Description: Large-scale features of Pliocene climate: results from the Pliocene Model Intercomparison Project Climate of the Past, 9, 191-209, 2013 Author(s): A. M. Haywood, D. J. Hill, A. M. Dolan, B. L. Otto-Bliesner, F. Bragg, W.-L. Chan, M. A. Chandler, C. Contoux, H. J. Dowsett, A. Jost, Y. Kamae, G. Lohmann, D. J. Lunt, A. Abe-Ouchi, S. J. Pickering, G. Ramstein, N. A. Rosenbloom, U. Salzmann, L. Sohl, C. Stepanek, H. Ueda, Q. Yan, and Z. Zhang Climate and environments of the mid-Pliocene warm period (3.264 to 3.025 Ma) have been extensively studied. Whilst numerical models have shed light on the nature of climate at the time, uncertainties in their predictions have not been systematically examined. The Pliocene Model Intercomparison Project quantifies uncertainties in model outputs through a coordinated multi-model and multi-model/data intercomparison. Whilst commonalities in model outputs for the Pliocene are clearly evident, we show substantial variation in the sensitivity of models to the implementation of Pliocene boundary conditions. Models appear able to reproduce many regional changes in temperature reconstructed from geological proxies. However, data/model comparison highlights that models potentially underestimate polar amplification. To assert this conclusion with greater confidence, limitations in the time-averaged proxy data currently available must be addressed. Furthermore, sensitivity tests exploring the known unknowns in modelling Pliocene climate specifically relevant to the high latitudes are essential (e.g. palaeogeography, gateways, orbital forcing and trace gasses). Estimates of longer-term sensitivity to CO 2 (also known as Earth System Sensitivity; ESS), support previous work suggesting that ESS is greater than Climate Sensitivity (CS), and suggest that the ratio of ESS to CS is between 1 and 2, with a "best" estimate of 1.5.

Description: Response of methane emissions from wetlands to the Last Glacial Maximum and an idealized Dansgaard–Oeschger climate event: insights from two models of different complexity Climate of the Past, 9, 149-171, 2013 Author(s): B. Ringeval, P. O. Hopcroft, P. J. Valdes, P. Ciais, G. Ramstein, A. J. Dolman, and M. Kageyama The role of different sources and sinks of CH 4 in changes in atmospheric methane ([CH 4 ]) concentration during the last 100 000 yr is still not fully understood. In particular, the magnitude of the change in wetland CH 4 emissions at the Last Glacial Maximum (LGM) relative to the pre-industrial period (PI), as well as during abrupt climatic warming or Dansgaard–Oeschger (D–O) events of the last glacial period, is largely unconstrained. In the present study, we aim to understand the uncertainties related to the parameterization of the wetland CH 4 emission models relevant to these time periods by using two wetland models of different complexity (SDGVM and ORCHIDEE). These models have been forced by identical climate fields from low-resolution coupled atmosphere–ocean general circulation model (FAMOUS) simulations of these time periods. Both emission models simulate a large decrease in emissions during LGM in comparison to PI consistent with ice core observations and previous modelling studies. The global reduction is much larger in ORCHIDEE than in SDGVM (respectively −67 and −46%), and whilst the differences can be partially explained by different model sensitivities to temperature, the major reason for spatial differences between the models is the inclusion of freezing of soil water in ORCHIDEE and the resultant impact on methanogenesis substrate availability in boreal regions. Besides, a sensitivity test performed with ORCHIDEE in which the methanogenesis substrate sensitivity to the precipitations is modified to be more realistic gives a LGM reduction of −36%. The range of the global LGM decrease is still prone to uncertainty, and here we underline its sensitivity to different process parameterizations. Over the course of an idealized D–O warming, the magnitude of the change in wetland CH 4 emissions simulated by the two models at global scale is very similar at around 15 Tg yr −1 , but this is only around 25% of the ice-core measured changes in [CH 4 ]. The two models do show regional differences in emission sensitivity to climate with much larger magnitudes of northern and southern tropical anomalies in ORCHIDEE. However, the simulated northern and southern tropical anomalies partially compensate each other in both models limiting the net flux change. Future work may need to consider the inclusion of more detailed wetland processes (e.g. linked to permafrost or tropical floodplains), other non-wetland CH 4 sources or different patterns of D–O climate change in order to be able to reconcile emission estimates with the ice-core data for rapid CH 4 events.

Description: Tephrostratigraphic studies on a sediment core from Lake Prespa in the Balkans Climate of the Past, 9, 267-287, 2013 Author(s): M. Damaschke, R. Sulpizio, G. Zanchetta, B. Wagner, A. Böhm, N. Nowaczyk, J. Rethemeyer, and A. Hilgers A detailed tephrostratigraphic record, which dates back to Marine Isotope Stage (MIS) 5b (ca. 91 kyr), has been established from a 17.76 m long core (Co1215) from Lake Prespa (Macedonia, Albania and Greece). A total of eleven tephra and cryptotephra layers (PT0915-1 to PT0915-11) were identified, using XRF scanning, magnetic susceptibility measurements, and macro- and microscopic inspection of the sediments. The major element composition of glass shards and/or micro-pumice fragments indicates that the tephras and cryptotephras originate from the explosive volcanism of Italy. Eight tephra and cryptotephra layers were correlated with specific volcanic eruptions: the AD 512 eruption of Somma-Vesuvius (1438 cal yr BP), the Mercato eruption of Somma-Vesuvius (8890 ± 90 cal yr BP), the Tufi Biancastri/LN1-LN2 eruption of the Campi Flegrei (14 749 ± 523 cal yr BP and 15 551 ± 621 cal yr BP), the SMP1-e/Y-3 eruption of the Campi Flegrei (30 000–31 000 cal yr BP), the Campanian Ignimbrite/Y-5 eruption of the Campi Flegrei (39 280 ± 110 cal yr BP), the SMP1-a event of Ischia Island (around 44 000 cal yr BP) and the Green Tuff/Y-6 eruption of Pantelleria Island (around 45 000 cal yr BP). One tephra could be attributed to the volcanic activity of Mount Etna, but probably represents an unknown eruption at ca. 60 000 cal yr BP. Cryptotephras PT0915-6 and PT0915-10 remain unclassified so far, but according to the presented age-depth model these would have been deposited around 35 000 and 48 500 cal yr BP, respectively. Some of the tephras and cryptotephras are recognised for the first time in the Balkan region. The tephrostratigraphic work provides important information about ash dispersal and explosion patterns of source volcanoes and can be used to correlate and date geographically distant paleoenvironmental and archaeological archives in the central Mediterranean region. Moreover, the tephrostratigraphic work in combination with radiocarbon and electron spin resonance (ESR) dating is a precondition for paleoclimatic reconstructions inferred from the sediment succession Co1215.

Description: Tree-ring–based summer mean temperature variations in the Adamello–Presanella Group (Italian Central Alps), 1610–2008 AD Climate of the Past, 9, 211-221, 2013 Author(s): A. Coppola, G. Leonelli, M. C. Salvatore, M. Pelfini, and C. Baroni Climate records from remote mountain sites and for century-long periods are usually lacking for most continents and also for the European Alps. However, detailed reconstructions of climate parameters for pre-instrumental periods in mountain areas, suffering of glacial retreat caused by recent global warming, are needed in the view of a better comprehension of the environmental dynamics. We present here the first annually-resolved reconstruction of summer (JJA) mean temperature for the Adamello–Presanella Group (Central European Alps), one of the most glaciated mountain groups of the Italian Central Alps. The reconstruction has been based on four larch tree-ring width chronologies derived from living trees sampled in four valleys surrounding the Group. The reconstruction spans from 1610 to 2008 and the statistical verification of the reconstruction demonstrates the positive skill of the tree-ring dataset in tracking summer temperature variability also in the recent period.

Description: Simultaneous measurement of OI 557.7 nm, O 2 (0, 1) Atmospheric Band and OH (6, 2) Meinel Band nightglow at Kolhapur (17° N), India Annales Geophysicae, 31, 197-208, 2013 Author(s): N. Parihar, A. Taori, S. Gurubaran, and G. K. Mukherjee Near-simultaneous measurements of OI 557.7 nm, O 2 (0, 1) Atmospheric Band and OH (6, 2) Meinel Band nightglow were carried out at Kolhapur (17° N), India during February–March 2007. Atmospheric temperatures around 87 and 94 km were derived from the knowledge of intensity measurements of spectral features OH (6, 2) Meinel Band and O 2 Atmospheric Band, respectively. An account of the behaviour of derived temperatures has been presented. The nocturnal behaviour of OH and O 2 temperatures is governed by the waves of tidal origin, whereas the signatures of planetary wave-like oscillations is noted in the night-to-night variation of two temperatures. This is probably the first report of planetary waves observed in nightglow temperatures from the Indian subcontinent.

Description: Refilling process in the plasmasphere: a 3-D statistical characterization based on Cluster density observations Annales Geophysicae, 31, 217-237, 2013 Author(s): G. Lointier, F. Darrouzet, P. M. E. Décréau, X. Vallières, S. Kougblénou, J. G. Trotignon, and J.-L. Rauch The Cluster mission offers an excellent opportunity to investigate the evolution of the plasma population in a large part of the inner magnetosphere, explored near its orbit's perigee, over a complete solar cycle. The WHISPER sounder, on board each satellite of the mission, is particularly suitable to study the electron density in this region, between 0.2 and 80 cm −3 . Compiling WHISPER observations during 1339 perigee passes distributed over more than three years of the Cluster mission, we present first results of a statistical analysis dedicated to the study of the electron density morphology and dynamics along and across magnetic field lines between L = 2 and L = 10. In this study, we examine a specific topic: the refilling of the plasmasphere and trough regions during extended periods of quiet magnetic conditions. To do so, we survey the evolution of the ap index during the days preceding each perigee crossing and sort out electron density profiles along the orbit according to three classes, namely after respectively less than 2 days, between 2 and 4 days, and more than 4 days of quiet magnetic conditions (ap ≤ 15 nT) following an active episode (ap 〉 15 nT). This leads to three independent data subsets. Comparisons between density distributions in the 3-D plasmasphere and trough regions at the three stages of quiet magnetosphere provide novel views about the distribution of matter inside the inner magnetosphere during several days of low activity. Clear signatures of a refilling process inside an expended plasmasphere in formation are noted. A plasmapause-like boundary, at L ~ 6 for all MLT sectors, is formed after 3 to 4 days and expends somewhat further after that. In the outer part of the plasmasphere ( L ~ 8), latitudinal profiles of median density values vary essentially according to the MLT sector considered rather than according to the refilling duration. The shape of these density profiles indicates that magnetic flux tubes are not fully replenished after 6 days of quiet conditions. In addition, the outer plasmasphere in the night and dawn sectors (22:00 to 10:00 MLT range) maintains an overall clear deficit of ionospheric population, when compared to the situation in the noon and dusk sectors (10:00 to 22:00 MLT range).

Description: Spatiotemporal variability analysis of diffuse radiation in China during 1981–2010 Annales Geophysicae, 31, 277-289, 2013 Author(s): X. L. Ren, H. L. He, L. Zhang, L. Zhou, G. R. Yu, and J. W. Fan Solar radiation is the primary driver of terrestrial plant photosynthesis and the diffuse component can enhance canopy light use efficiency (LUE), which in turn influences the carbon balance of terrestrial ecosystems. In this study we calculated the spatial data of diffuse radiation in China from 1981 to 2010, using a radiation decomposition model and spatial interpolation method based on observational data. Furthermore, we explored the spatiotemporal characteristics of diffuse radiation using GIS and trend analysis techniques. The results show the following: (1) The spatial patterns of perennial average of annual diffuse radiation during 1981–2010 are complex and inhomogeneous in China, generally lower in the north and higher in the south and west. The perennial average ranges from 1730.20 to 3064.41 MJ m −2 yr −1 across the whole country. (2) There is an increasing trend of annual diffuse radiation in China from 1981 to 2010 on the whole, with mean increasing amplitude of 7.03 MJ m −2 yr −1 per decade. Whereas a significant downtrend was observed in the first 10 years, distinct anomalies in 1982, 1983, 1991 and 1992 occurred due to the eruptions of El Chinchon and Pinatubo. (3) The spatial distribution of the temporal variability of diffuse radiation showed significant regional heterogeneity in addition to the seasonal differences. Northwestern China has the most evident downtrend, with highest decreasing rate of 6% per decade, while the Tibetan Plateau has the most evident uptrend, with highest increasing rate of up to 9% per decade. Such quantitative spatiotemporal characteristics of diffuse radiation are essential in regional scale modeling of terrestrial carbon dynamics.

Description: A new global reconstruction of temperature changes at the Last Glacial Maximum Climate of the Past, 9, 367-376, 2013 Author(s): J. D. Annan and J. C. Hargreaves Some recent compilations of proxy data both on land and ocean (MARGO Project Members, 2009; Bartlein et al., 2011; Shakun et al., 2012), have provided a new opportunity for an improved assessment of the overall climatic state of the Last Glacial Maximum. In this paper, we combine these proxy data with the ensemble of structurally diverse state of the art climate models which participated in the PMIP2 project (Braconnot et al., 2007) to generate a spatially complete reconstruction of surface air (and sea surface) temperatures. We test a variety of approaches, and show that multiple linear regression performs well for this application. Our reconstruction is significantly different to and more accurate than previous approaches and we obtain an estimated global mean cooling of 4.0 ± 0.8 °C (95% CI).

Description: Reconstruction of drip-water δ 18 O based on calcite oxygen and clumped isotopes of speleothems from Bunker Cave (Germany) Climate of the Past, 9, 377-391, 2013 Author(s): T. Kluge, H. P. Affek, T. Marx, W. Aeschbach-Hertig, D. F. C. Riechelmann, D. Scholz, S. Riechelmann, A. Immenhauser, D. K. Richter, J. Fohlmeister, A. Wackerbarth, A. Mangini, and C. Spötl The geochemical signature of many speleothems used for reconstruction of past continental climates is affected by kinetic isotope fractionation. This limits quantitative paleoclimate reconstruction and, in cases where the kinetic fractionation varies with time, also affects relative paleoclimate interpretations. In carbonate archive research, clumped isotope thermometry is typically used as proxy for absolute temperatures. In the case of speleothems, however, clumped isotopes provide a sensitive indicator for disequilibrium effects. The extent of kinetic fractionation co-varies in Δ 47 and δ 18 O so that it can be used to account for disequilibrium in δ 18 O and to extract the past drip-water composition. Here we apply this approach to stalagmites from Bunker Cave (Germany) and calculate drip-water δ 18 O w values for the Eemian, MIS3, and the Holocene, relying on independent temperature estimates and accounting for disequilibrium. Applying the co-variation method to modern calcite precipitates yields drip-water δ 18 O w values in agreement with modern cave drip-water δ 18 O w of −7.9 ± 0.3‰, despite large and variable disequilibrium effects in both calcite δ 18 O c and Δ 47 . Reconstructed paleo-drip-water δ 18 O w values are lower during colder periods (e.g., MIS3: −8.6 ± 0.4‰ and the early Holocene at 11 ka: −9.7 ± 0.2‰) and show higher values during warmer climatic periods (e.g., the Eemian: −7.6 ± 0.2‰ and the Holocene Climatic Optimum: −7.2 ± 0.3‰). This new approach offers a unique possibility for quantitative climate reconstruction including the assessment of past hydrological conditions while accounting for disequilibrium effects.

Description: Proxy benchmarks for intercomparison of 8.2 ka simulations Climate of the Past, 9, 423-432, 2013 Author(s): C. Morrill, D. M. Anderson, B. A. Bauer, R. Buckner, E. P. Gille, W. S. Gross, M. Hartman, and A. Shah The Paleoclimate Modelling Intercomparison Project (PMIP3) now includes the 8.2 ka event as a test of model sensitivity to North Atlantic freshwater forcing. To provide benchmarks for intercomparison, we compiled and analyzed high-resolution records spanning this event. Two previously-described anomaly patterns that emerge are cooling around the North Atlantic and drier conditions in the Northern Hemisphere tropics. Newer to this compilation are more robustly-defined wetter conditions in the Southern Hemisphere tropics and regionally-limited warming in the Southern Hemisphere. Most anomalies around the globe lasted on the order of 100 to 150 yr. More quantitative reconstructions are now available and indicate cooling of ~ 1 °C and a ~ 20% decrease in precipitation in parts of Europe as well as spatial gradients in δ 18 O from the high to low latitudes. Unresolved questions remain about the seasonality of the climate response to freshwater forcing and the extent to which the bipolar seesaw operated in the early Holocene.

Description: Large-scale temperature response to external forcing in simulations and reconstructions of the last millennium Climate of the Past, 9, 393-421, 2013 Author(s): L. Fernández-Donado, J. F. González-Rouco, C. C. Raible, C. M. Ammann, D. Barriopedro, E. García-Bustamante, J. H. Jungclaus, S. J. Lorenz, J. Luterbacher, S. J. Phipps, J. Servonnat, D. Swingedouw, S. F. B. Tett, S. Wagner, P. Yiou, and E. Zorita Understanding natural climate variability and its driving factors is crucial to assessing future climate change. Therefore, comparing proxy-based climate reconstructions with forcing factors as well as comparing these with paleoclimate model simulations is key to gaining insights into the relative roles of internal versus forced variability. A review of the state of modelling of the climate of the last millennium prior to the CMIP5–PMIP3 (Coupled Model Intercomparison Project Phase 5–Paleoclimate Modelling Intercomparison Project Phase 3) coordinated effort is presented and compared to the available temperature reconstructions. Simulations and reconstructions broadly agree on reproducing the major temperature changes and suggest an overall linear response to external forcing on multidecadal or longer timescales. Internal variability is found to have an important influence at hemispheric and global scales. The spatial distribution of simulated temperature changes during the transition from the Medieval Climate Anomaly to the Little Ice Age disagrees with that found in the reconstructions. Thus, either internal variability is a possible major player in shaping temperature changes through the millennium or the model simulations have problems realistically representing the response pattern to external forcing. A last millennium transient climate response (LMTCR) is defined to provide a quantitative framework for analysing the consistency between simulated and reconstructed climate. Beyond an overall agreement between simulated and reconstructed LMTCR ranges, this analysis is able to single out specific discrepancies between some reconstructions and the ensemble of simulations. The disagreement is found in the cases where the reconstructions show reduced covariability with external forcings or when they present high rates of temperature change.

Description: Greenland ice sheet contribution to sea level rise during the last interglacial period: a modelling study driven and constrained by ice core data Climate of the Past, 9, 353-366, 2013 Author(s): A. Quiquet, C. Ritz, H. J. Punge, and D. Salas y Mélia As pointed out by the forth assessment report of the Intergovernmental Panel on Climate Change, IPCC-AR4 (Meehl et al., 2007), the contribution of the two major ice sheets, Antarctica and Greenland, to global sea level rise, is a subject of key importance for the scientific community. By the end of the next century, a 3–5 °C warming is expected in Greenland. Similar temperatures in this region were reached during the last interglacial (LIG) period, 130–115 ka BP, due to a change in orbital configuration rather than to an anthropogenic forcing. Ice core evidence suggests that the Greenland ice sheet (GIS) survived this warm period, but great uncertainties remain about the total Greenland ice reduction during the LIG. Here we perform long-term simulations of the GIS using an improved ice sheet model. Both the methodologies chosen to reconstruct palaeoclimate and to calibrate the model are strongly based on proxy data. We suggest a relatively low contribution to LIG sea level rise from Greenland melting, ranging from 0.7 to 1.5 m of sea level equivalent, contrasting with previous studies. Our results suggest an important contribution of the Antarctic ice sheet to the LIG highstand.

Description: Dansgaard–Oeschger events: bifurcation points in the climate system Climate of the Past, 9, 323-333, 2013 Author(s): A. A. Cimatoribus, S. S. Drijfhout, V. Livina, and G. van der Schrier Dansgaard–Oeschger events are a prominent mode of variability in the records of the last glacial cycle. Various prototype models have been proposed to explain these rapid climate fluctuations, and no agreement has emerged on which may be the more correct for describing the palaeoclimatic signal. In this work, we assess the bimodality of the system, reconstructing the topology of the multi-dimensional attractor over which the climate system evolves. We use high-resolution ice core isotope data to investigate the statistical properties of the climate fluctuations in the period before the onset of the abrupt change. We show that Dansgaard–Oeschger events have weak early warning signals if the ensemble of events is considered. We find that the statistics are consistent with the switches between two different climate equilibrium states in response to a changing external forcing (e.g. solar, ice sheets), either forcing directly the transition or pacing it through stochastic resonance. These findings are most consistent with a model that associates Dansgaard–Oeschger with changing boundary conditions, and with the presence of a bifurcation point.

Description: Reconstruction of northeast Asia spring temperature 1784–1990 Climate of the Past, 9, 261-266, 2013 Author(s): M. Ohyama, H. Yonenobu, J.-N. Choi, W.-K. Park, M. Hanzawa, and M. Suzuki We describe a first attempt of wide-area dendroclimatic reconstruction, based upon seven spring temperature-sensitive chronologies from the ring widths of living trees, in Japan and Korea. Mean March–May temperature derived from a gridded land air temperature dataset (CRUTEM4) between 35–40° N and 125–140° E was reconstructed for the period of AD 1784–1990. Of the seven, two Japanese chronologies were eliminated during the calibration trials. The reconstruction accounted for 19.4% of the temperature variance in the calibration period, and is considered to be skillful for estimating interannual-to-interdecadal variations and not for long-term change. This reconstruction showed remarkably similar fluctuations to regional dendroclimatic reconstructions in Japan and Korea, indicating the past spatial coherency of spring temperatures in the region. The reconstruction was validated against other climate proxies. A fairly good agreement was found with cold periods as estimated from documentary records in southeast China and Japan. The west Japan temperature series recovered from instrumental records also showed a reasonable agreement with the reconstruction. On the other hand, the reconstruction did not show clear abrupt depressions after the Laki and the Tambora eruptions. These comparisons revealed that dendroclimatic spatial reconstruction in this area offers a good potential for reconstructing long-term and large-scale past temperature patterns for northeast Asia.

Description: The Medieval Climate Anomaly and the Little Ice Age in the eastern Ecuadorian Andes Climate of the Past, 9, 307-321, 2013 Author(s): M.-P. Ledru, V. Jomelli, P. Samaniego, M. Vuille, S. Hidalgo, M. Herrera, and C. Ceron To better characterize the climate variability of the last millennium in the high Andes, we analyzed the pollen content of a 1150-yr-old sediment core collected in a bog located at 3800 m a.s.l. in the páramo in the eastern Cordillera in Ecuador. An upslope convective index based on the ratio between cloud transported pollen from the Andean forest to the bog (T) and Poaceae pollen frequencies, related to the edaphic moisture of the páramo (P), was defined. This index was used to distinguish changes in the atmospheric moisture from the soil moisture content of the páramo and their associated patterns of interdecadal El Niño–Southern Oscillation (ENSO) variability and South American summer monsoon (SASM) activity. Results show that between 850 and 1250 AD, the Medieval Climate Anomaly interval was warm and moist with a high transported pollen/Poaceae pollen (T/P) index linked to high ENSO variability and weak SASM activity. Between 1250 and 1550 AD, a dry climate prevailed, characterized by an abrupt decrease in the T/P index and therefore no upslope cloud convection, related to lower ENSO variability and with significant impact on the floristic composition of the páramo. During the Little Ice Age, two phases were observed: first, a wet phase between 1550 and 1750 AD linked to low ENSO variability in the Pacific and warm south equatorial Atlantic sea surface temperatures (SSTs) favored the return of a wet páramo, and then a cold and dry phase between 1750 and 1800 AD associated with low ENSO variability and weak SASM activity resulted in drying of the páramo. The current warm period marks the beginning of a climate characterized by high convective activity – the highest in the last millennium – and weaker SASM activity modifying the water storage of the páramo. Our results show that the páramo is progressively losing its capacity for water storage and that the interdecadal variability of both tropical Pacific and Atlantic SSTs matter for Andean climate patterns, although many teleconnection mechanisms are still poorly understood.

Description: Modern isotope hydrology and controls on δD of plant leaf waxes at Lake El'gygytgyn, NE Russia Climate of the Past, 9, 335-352, 2013 Author(s): K. M. K. Wilkie, B. Chapligin, H. Meyer, S. Burns, S. Petsch, and J. Brigham-Grette Stable isotope data from lipid biomarkers and diatom silica recovered from lake sediment cores hold great promise for paleoclimate and paleohydrological reconstructions. However, these records rely on accurate calibration with modern precipitation and hydrologic processes and only limited data exist on the controls on the δD values for n -alkanoic acids from plant leaf waxes. Here we investigate the stable isotopic composition of modern precipitation, streams, lake water and ice cover, and use these data to constrain isotope systematics of the Lake El'gygytgyn Basin hydrology. Compound-specific hydrogen isotope ratios determined from n -alkanoic acids from modern vegetation are compared with modern precipitation and lake core top sediments. Multi-species net (apparent) fractionation values between source water (precipitation) and modern vegetation (e.g., ϵ wax/precip mean value is −107 ± 12‰) agree with previous results and suggest a consistent offset between source waters and the δD values of alkanoic acids. We conclude that although there may be some bias towards a winter precipitation signal, overall δD values from leaf wax n -alkanoic acids record annual average precipitation within the El'gygytgyn Basin. A net fractionation calculated for 200-yr-integrated lake sediments yields ϵ 30/precip = −96 ± 8‰ and can provide robust net "apparent" fractionation to be used in future paleohydrological reconstructions.

Description: Does Antarctic glaciation cool the world? Climate of the Past, 9, 173-189, 2013 Author(s): A. Goldner, M. Huber, and R. Caballero In this study, we compare the simulated climatic impact of adding an Antarctic ice sheet (AIS) to the "greenhouse world" of the Eocene and removing the AIS from the modern world. The modern global mean surface temperature anomaly (Δ T ) induced by Antarctic Glaciation depends on the background CO 2 levels and ranges from −1.22 to −0.18 K. The Eocene Δ T is nearly constant at ~−0.25 K. We calculate an climate sensitivity parameter S [Antarctica] which we define as Δ T divided by the change in effective radiative forcing (Δ Q Antarctica ) which includes some fast feedbacks imposed by prescribing the glacial properties of Antarctica. The main difference between the modern and Eocene responses is that a negative cloud feedback warms much of the Earth's surface as a large AIS is introduced in the Eocene, whereas this cloud feedback is weakly positive and acts in combination with positive sea-ice feedbacks to enhance cooling introduced by adding an ice sheet in the modern. Because of the importance of cloud feedbacks in determining the final temperature sensitivity of the AIS, our results are likely to be model dependent. Nevertheless, these model results suggest that the effective radiative forcing and feedbacks induced by the AIS did not significantly decrease global mean surface temperature across the Eocene–Oligocene transition (EOT −34.1 to 33.6 Ma) and that other factors like declining atmospheric CO 2 are more important for cooling across the EOT. The results illustrate that the efficacy of AIS forcing in the Eocene is not necessarily close to one and is likely to be model and state dependent. This implies that using EOT paleoclimate proxy data by itself to estimate climate sensitivity for future climate prediction requires climate models and consequently these estimates will have large uncertainty, largely due to uncertainties in modelling low clouds.

Description: Modern sedimentation patterns in Lake El'gygytgyn, NE Russia, derived from surface sediment and inlet streams samples Climate of the Past, 9, 135-148, 2013 Author(s): V. Wennrich, A. Francke, A. Dehnert, O. Juschus, T. Leipe, C. Vogt, J. Brigham-Grette, P. S. Minyuk, M. Melles, and El'gygytgyn Science Party Lake El'gygytgyn/NE Russia holds a continuous 3.58 Ma sediment record, which is regarded as the most long-lasting climate archive of the terrestrial Arctic. Based on multi-proxy geochemical, mineralogical, and granulometric analyses of surface sediment, inlet stream and bedrock samples, supplemented by statistical methods, major processes influencing the modern sedimentation in the lake were investigated. Grain-size parameters and chemical elements linked to the input of feldspars from acidic bedrock indicate a wind-induced two-cell current system as major driver of sediment transport and accumulation processes in Lake El'gygytgyn. The distribution of mafic rock related elements in the sediment on the lake floor can be traced back to the input of weathering products of basaltic rocks in the catchment. Obvious similarities in the spatial variability of manganese and heavy metals indicate sorption or co-precipitation of these elements with Fe and Mn hydroxides and oxides. But the similar distribution of organic matter and clay contents might also point to a fixation to organic components and clay minerals. An enrichment of mercury in the inlet streams might be indicative of neotectonic activity around the lake. The results of this study add to the fundamental knowledge of the modern lake processes of Lake El'gygytgyn and its lake-catchment interactions, and thus, yield crucial insights for the interpretation of paleo-data from this unique archive.

Description: Past climate changes and permafrost depth at the Lake El'gygytgyn site: implications from data and thermal modeling Climate of the Past, 9, 119-133, 2013 Author(s): D. Mottaghy, G. Schwamborn, and V. Rath This study focuses on the temperature field observed in boreholes drilled as part of interdisciplinary scientific campaign targeting the El'gygytgyn Crater Lake in NE Russia. Temperature data are available from two sites: the lake borehole 5011-1 located near the center of the lake reaching 400 m depth, and the land borehole 5011-3 at the rim of the lake, with a depth of 140 m. Constraints on permafrost depth and past climate changes are derived from numerical simulation of the thermal regime associated with the lake-related talik structure. The thermal properties of the subsurface needed for these simulations are based on laboratory measurements of representative cores from the quaternary sediments and the underlying impact-affected rock, complemented by further information from geophysical logs and data from published literature. The temperature observations in the lake borehole 5011-1 are dominated by thermal perturbations related to the drilling process, and thus only give reliable values for the lowermost value in the borehole. Undisturbed temperature data recorded over more than two years are available in the 140 m deep land-based borehole 5011-3. The analysis of these observations allows determination of not only the recent mean annual ground surface temperature, but also the ground surface temperature history, though with large uncertainties. Although the depth of this borehole is by far too insufficient for a complete reconstruction of past temperatures back to the Last Glacial Maximum, it still affects the thermal regime, and thus permafrost depth. This effect is constrained by numerical modeling: assuming that the lake borehole observations are hardly influenced by the past changes in surface air temperature, an estimate of steady-state conditions is possible, leading to a meaningful value of 14 ± 5 K for the post-glacial warming. The strong curvature of the temperature data in shallower depths around 60 m can be explained by a comparatively large amplitude of the Little Ice Age (up to 4 K), with low temperatures prevailing far into the 20th century. Other mechanisms, like varying porosity, may also have an influence on the temperature profile, however, our modeling studies imply a major contribution from recent climate changes.

Description: A biomarker record of Lake El'gygytgyn, Far East Russian Arctic: investigating sources of organic matter and carbon cycling during marine isotope stages 1–3 Climate of the Past, 9, 243-260, 2013 Author(s): A. R. Holland, S. T. Petsch, I. S. Castañeda, K. M. Wilkie, S. J. Burns, and J. Brigham-Grette Arctic paleoenvironmental archives serve as sensitive recorders of past climate change. Lake El'gygytgyn (Far East Russian Arctic) is a high-latitude crater impact lake that contains a continuous sediment record influenced by neither glaciation nor glacial erosion since the time of impact 3.58 Ma ago. Prior research on sediments collected from Lake El'gygytgyn suggest times of permanent ice cover and anoxia corresponding to global glacial intervals, during which the sediments are laminated and are characterized by the co-occurrence of high total organic carbon, microscopic magnetite grains that show etching and dissolution, and negative excursions in bulk sediment organic matter carbon isotope (δ 13 C) values. Here we investigate the abundance and carbon isotopic composition of lipid biomarkers recovered from Lake El'gygytgyn sediments spanning marine isotope stages 1–3 to identify key sources of organic matter (OM) to lake sediments, to establish which OM sources drive the negative δ 13 C excursion exhibited by bulk sediment OM, and to explore if there are molecular and isotopic signatures of anoxia in the lake during glaciation. We find that during marine isotope stages 1–3, direct evidence for water column anoxia is lacking. A ~4‰ negative excursion in bulk sediment δ 13 C values during the Local Last Glacial Maximum (LLGM) is accompanied by more protracted, higher magnitude negative excursions in n -alkanoic acid and n -alkanol δ 13 C values that begin 20 kyr in advance of the LLGM. In contrast, n -alkanes and the C 30 n -alkanoic acid do not exhibit a negative δ 13 C excursion at this time. Our results indicate that the C 24 , C 26 and C 28 n -alkanoic acids do not derive entirely from terrestrial OM sources, while the C 30 n -alkanoic acid at Lake El'gygytgyn is a robust indicator of terrestrial OM contributions. Overall, our results strongly support the presence of a nutrient-poor water column, which is mostly isolated from atmospheric carbon dioxide during glaciation at Lake El'gygytgyn.

Description: Comment on "Storming the Bastille: the effect of electric fields on the ionospheric F-layer" by Rishbeth et al. (2010) Annales Geophysicae, 31, 145-150, 2013 Author(s): B. T. Tsurutani, A. J. Mannuccci, O. P. Verkhoglyadova, and G. S. Lakhina No abstract available.

Description: A first chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core Climate of the Past, 9, 2713-2730, 2013 Author(s): S. O. Rasmussen, P. M. Abbott, T. Blunier, A. J. Bourne, E. Brook, S. L. Buchardt, C. Buizert, J. Chappellaz, H. B. Clausen, E. Cook, D. Dahl-Jensen, S. M. Davies, M. Guillevic, S. Kipfstuhl, T. Laepple, I. K. Seierstad, J. P. Severinghaus, J. P. Steffensen, C. Stowasser, A. Svensson, P. Vallelonga, B. M. Vinther, F. Wilhelms, and M. Winstrup A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) and its model extension (GICC05modelext) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the electrical conductivity measurement (ECM) and dielectrical profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide five additional horizons used for the timescale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard–Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age–gas age difference (Δage) has been reconstructed using a coupled firn densification-heat diffusion model. Temperature and accumulation inputs to the Δage model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from δ 15 N of nitrogen and high-resolution methane data during the abrupt onset of Greenland interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, named GICC05modelext-NEEM-1. Based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed. Together, the timescale and accumulation reconstruction provide the necessary basis for further analysis of the records from NEEM.

Description: The East Asian winter monsoon variability in response to precession during the past 150 000 yr Climate of the Past, 9, 2777-2788, 2013 Author(s): M. Yamamoto, H. Sai, M.-T. Chen, and M. Zhao The response of the East Asian winter monsoon variability to orbital forcing is still unclear, and hypotheses are controversial. We present a 150 000 yr record of sea surface temperature difference (ΔSST) between the South China Sea and other Western Pacific Warm Pool regions as a proxy for the intensity of the Asian winter monsoon, because the winter cooling of the South China Sea is caused by the cooling of surface water at the northern margin and the southward advection of cooled water due to winter monsoon winds. The ΔSST showed dominant precession cycles during the past 150 000 yr. The ΔSST varies at precessional band and supports the hypothesis that monsoon is regulated by insolation changes at low-latitudes (Kutzbach, 1981), but contradicts previous suggestions based on marine and loess records that eccentricity controls variability on glacial–interglacial timescales. Maximum winter monsoon intensity corresponds to the May perihelion at precessional band, which is not fully consistent with the Kutzbach model of maximum winter monsoon at the June perihelion. Variation in the East Asian winter monsoon was anti-phased with the Indian summer monsoon, suggesting a linkage of dynamics between these two monsoon systems on an orbital timescale.

Description: Statistical study of foreshock cavitons Annales Geophysicae, 31, 2163-2178, 2013 Author(s): P. Kajdič, X. Blanco-Cano, N. Omidi, K. Meziane, C. T. Russell, J.-A. Sauvaud, I. Dandouras, and B. Lavraud In this work we perform a statistical analysis of 92 foreshock cavitons observed with the Cluster spacecraft 1 during the period 2001–2006. We analyze time intervals during which the spacecraft was located in the Earth's foreshock with durations longer than 10 min. Together these amount to ~ 50 days. The cavitons are transient structures in the Earth's foreshock. Their main signatures in the data include simultaneous depletions of the magnetic field intensity and plasma density, which are surrounded by a rim of enhanced values of these two quantities. Cavitons form due to nonlinear interaction of transverse and compressive ultra-low frequency (ULF) waves and are therefore always surrounded by intense compressive ULF fluctuations. They are carried by the solar wind towards the bow shock. This work represents the first systematic study of a large sample of foreshock cavitons. We find that cavitons appear for a wide range of solar wind and interplanetary magnetic field conditions and are therefore a common feature upstream of Earth's quasi-parallel bow shock with an average occurrence rate of ~ 2 events per day. We also discuss their observational properties in the context of other known upstream phenomena and show that the cavitons are a distinct structure in the foreshock.

Description: Theoretical constraints on the cross-tail width of bursty bulk flows Annales Geophysicae, 31, 2179-2192, 2013 Author(s): C. X. Chen The characteristic cross-tail width of bursty bulk flows (BBFs) in earth's plasma sheet was investigated at two stages of its life, one at its onset, the other when it is fully developed. Equilibrium domains with gradient of magnetic field are constructed. Interchange instability analysis of such domains yields the most unstable mode with the half wave length comparable with the observed cross-tail width of a flow burst and the inverse of growth rate comparable with its duration. The thickness of the plasma sheet for the most unstable mode is also comparable to the width of BBFs in the north–south direction. We found that viscosity, the dimension of the unstable domain, the thickness of the plasma sheet and gradient of the magnetic field together determine the most unstable mode. The ion Larmor radius plays an important role in viscosity as half effective mean free path. For a fully developed flow, however, velocity-caused pressure difference between the leading and trailing sides of a flow burst also plays a role. The equatorial cross section of flow is reshaped and its cross-tail width is changed as well. Representing the surrounding medium with empirical magnetic field and plasma models, the force balance of the fast flow is analyzed. The cross-section area of flow burst is estimated to be one to several square earth radii, and the cross-tail width of fast flow is estimated to be 1 to 3 earth radii, which is consistent with observations of BBFs.

Description: Ray tracing model of the auroral kilometric radiation generation in the 3-D plasma cavity Annales Geophysicae, 31, 1645-1652, 2013 Author(s): T. M. Burinskaya Propagation and amplification of the auroral kilometric radiation (AKR) in a three-dimensional plasma cavity is investigated using the approximation of the geometrical optics, and taking into account both the slightly relativistic electrons propagating inside a cavity and the background cold electrons. It is shown that the global magnetic field inhomogeneity plays a key role in a wave escape from a thin plasma cavity. The main contribution to the AKR spectrum is made by waves initially generated with the component of group velocity directed to the Earth and with the optimum relationship between the wave vector components, controlling the value of the linear grow rate and duration of the ray lifetime inside a source.

Description: Empirical regional models for the short-term forecast of M3000F2 during not quiet geomagnetic conditions over Europe Annales Geophysicae, 31, 1653-1671, 2013 Author(s): M. Pietrella Twelve empirical local models have been developed for the long-term prediction of the ionospheric characteristic M3000F2 , and then used as starting point for the development of a short-term forecasting empirical regional model of M3000F2 under not quiet geomagnetic conditions. Under the assumption that the monthly median measurements of M3000F2 are linearly correlated to the solar activity, a set of regression coefficients were calculated over 12 months and 24 h for each of 12 ionospheric observatories located in the European area, and then used for the long-term prediction of M3000F2 at each station under consideration. Based on the 12 long-term prediction empirical local models of M3000F2 , an empirical regional model for the prediction of the monthly median field of M3000F2 over Europe (indicated as RM_M3000F2 ) was developed. Thanks to the IFELM_foF2 models, which are able to provide short-term forecasts of the critical frequency of the F2 layer ( fo F2 STF ) up to three hours in advance, it was possible to considerer the Brudley–Dudeney algorithm as a function of fo F2 STF to correct RM_M3000F2 and thus obtain an empirical regional model for the short-term forecasting of M3000F2 (indicated as RM_M3000F2_BD ) up to three hours in advance under not quiet geomagnetic conditions. From the long-term predictions of M3000F2 provided by the IRI model, an empirical regional model for the forecast of the monthly median field of M3000F2 over Europe (indicated as IRI_RM_M3000F2 ) was derived. IRI_RM_M3000F2 predictions were modified with the Bradley–Dudeney correction factor, and another empirical regional model for the short-term forecasting of M3000F2 (indicated as IRI_RM_M3000F2_BD ) up to three hours ahead under not quiet geomagnetic conditions was obtained. The main results achieved comparing the performance of RM_M3000F2 , RM_M3000F2_BD , IRI_RM_M3000F2 , and IRI_RM_M3000F2_BD are (1) in the case of moderate geomagnetic activity, the Bradley–Dudeney correction factor does not improve significantly the predictions; (2) under disturbed geomagnetic conditions, the Bradley–Dudeney formula improves the predictions of RM_M3000F2 in the entire European area; (3) in the case of very disturbed geomagnetic conditions, the Bradley–Dudeney algorithm is very effective in improving the performance of IRI_RM_M3000F2 ; (4) under moderate geomagnetic conditions, the long-term prediction maps of M3000F2 generated by RM_M3000F2 can be considered as short-term forecasting maps providing very satisfactory results because quiet geomagnetic conditions are not so diverse from moderate geomagnetic conditions; (5) the forecasting maps originated by RM_M3000F2 , RM_M3000F2_BD , and IRI_RM_M3000F2_BD show some regions where the forecasts are not satisfactory, but also wide sectors where the M3000F2 forecasts quite faithfully match the M3000F2 observations, and therefore RM_M3000F2 , RM_M3000F2_BD , and IRI_RM_M3000F2_BD could be exploited to produce short-term forecasting maps of M3000F2 over Europe up to 3 h in advance.

Description: On the supply of heavy planetary material to the magnetotail of Mercury Annales Geophysicae, 31, 1673-1679, 2013 Author(s): D. C. Delcourt We examine the transport of low-energy heavy ions of planetary origin (O + , Na + , Ca + ) in the magnetosphere of Mercury. We show that, in contrast to Earth, these ions are abruptly energized after ejection into the magnetosphere due to enhanced curvature-related parallel acceleration. Regardless of their mass-to-charge ratio, the parallel speed of these ions is rapidly raised up to ~ 2 V E × B (denoting by V E × B the magnitude of the local E × B drift speed), in a like manner to Fermi-type acceleration by a moving magnetic mirror. This parallel energization is such that ions with very low initial energies (a few tenths of eVs) can overcome gravity and, regardless of species or convection rate, are transported over comparable distances into the nightside magnetosphere. The region of space where these ions reach the magnetotail is found to extend over altitudes similar to those where enhanced densities are noticeable in the MESSENGER data, viz., from ~ 1000 km up to ~ 6000 km in the pre-midnight sector. The observed density enhancements may thus follow from E × B related focusing of planetary material of dayside origin into the magnetotail. Due to the planetary magnetic field offset, an asymmetry is found between drift paths anchored in the Northern and Southern hemispheres, which puts forward a predominant role of heavy material originating in the Northern Hemisphere in populating the innermost region of Mercury's magnetotail.

Description: Evaluating the applicability of the finite element method for modelling of geoelectric fields Annales Geophysicae, 31, 1689-1698, 2013 Author(s): B. Dong, D. W. Danskin, R. J. Pirjola, D. H. Boteler, and Z. Z. Wang Geomagnetically induced currents in power systems are due to space weather events which create geomagnetic disturbances accompanied by electric fields at the surface of the Earth. The purpose of this paper is to evaluate the use of the finite element method (FEM) to calculate the magnetic and electric fields to which long transmission lines of power systems on the Earth are exposed. The well-known technique of FEM is used for the first time to simulate magnetic and electric fields applicable to power systems. Several test cases are modelled and compared with known solutions. It is shown that FEM is an effective modelling technique that can be applied to determine the electric fields which affect power systems. FEM enables an increased capability beyond the traditional methods for modelling electric and magnetic fields with layered earth conductivity structures, as spatially more complex structures can be considered using FEM. As an example results are presented for induction, due to a line current source, in adjacent regions with different layered conductivity structures. The results show the electric field away from the interface is the same as calculated for a single region; however near the interface the electric field is influenced by both regions.

Description: Stereoscopic determination of all-sky altitude map of aurora using two ground-based Nikon DSLR cameras Annales Geophysicae, 31, 1543-1548, 2013 Author(s): R. Kataoka, Y. Miyoshi, K. Shigematsu, D. Hampton, Y. Mori, T. Kubo, A. Yamashita, M. Tanaka, T. Takahei, T. Nakai, H. Miyahara, and K. Shiokawa A new stereoscopic measurement technique is developed to obtain an all-sky altitude map of aurora using two ground-based digital single-lens reflex (DSLR) cameras. Two identical full-color all-sky cameras were set with an 8 km separation across the Chatanika area in Alaska (Poker Flat Research Range and Aurora Borealis Lodge) to find localized emission height with the maximum correlation of the apparent patterns in the localized pixels applying a method of the geographical coordinate transform. It is found that a typical ray structure of discrete aurora shows the broad altitude distribution above 100 km, while a typical patchy structure of pulsating aurora shows the narrow altitude distribution of less than 100 km. Because of its portability and low cost of the DSLR camera systems, the new technique may open a unique opportunity not only for scientists but also for night-sky photographers to complementarily attend the aurora science to potentially form a dense observation network.

Description: Mid- and late Holocene dust deposition in western Europe: the Misten peat bog (Hautes Fagnes – Belgium) Climate of the Past, 9, 2285-2298, 2013 Author(s): M. Allan, G. Le Roux, N. Piotrowska, J. Beghin, E. Javaux, M. Court-Picon, N. Mattielli, S. Verheyden, and N. Fagel Dust deposition in southern Belgium is estimated from the geochemical signature of an ombrotrophic peatland. The rare earth elements (REE) and lithogenic elements concentrations, as well as Nd isotopes, were determined by HR-ICP-MS and MC-ICP-MS, respectively, along an ~6 m peat section covering 5300 yr, from 2000 to 7300 cal BP, dated by the 14 C method. Changes in REE concentration in the peat correlate with those of Ti, Al, Sc and Zr that are lithogenic conservative elements, suggesting that REE are immobile in the studied peat bogs and can be used as tracers of dust deposition. Peat humification and testate amoebae were used to evaluate hydroclimatic conditions. The range of dust deposition varied from 0.03 to 4.0 g m −2 yr −1 . The highest dust fluxes were observed from 2750 to 2550 cal BP and from 5150 to 4750 cal BP, and correspond to cold periods. The ϵNd values show a large variability from −13 to −5, identifying three major sources of dusts: local soils, distal volcanic and desert particles.

Description: Inferred changes in El Niño–Southern Oscillation variance over the past six centuries Climate of the Past, 9, 2269-2284, 2013 Author(s): S. McGregor, A. Timmermann, M. H. England, O. Elison Timm, and A. T. Wittenberg It is vital to understand how the El Niño–Southern Oscillation (ENSO) has responded to past changes in natural and anthropogenic forcings, in order to better understand and predict its response to future greenhouse warming. To date, however, the instrumental record is too brief to fully characterize natural ENSO variability, while large discrepancies exist amongst paleo-proxy reconstructions of ENSO. These paleo-proxy reconstructions have typically attempted to reconstruct ENSO's temporal evolution, rather than the variance of these temporal changes. Here a new approach is developed that synthesizes the variance changes from various proxy data sets to provide a unified and updated estimate of past ENSO variance. The method is tested using surrogate data from two coupled general circulation model (CGCM) simulations. It is shown that in the presence of dating uncertainties, synthesizing variance information provides a more robust estimate of ENSO variance than synthesizing the raw data and then identifying its running variance. We also examine whether good temporal correspondence between proxy data and instrumental ENSO records implies a good representation of ENSO variance. In the climate modeling framework we show that a significant improvement in reconstructing ENSO variance changes is found when combining information from diverse ENSO-teleconnected source regions, rather than by relying on a single well-correlated location. This suggests that ENSO variance estimates derived from a single site should be viewed with caution. Finally, synthesizing existing ENSO reconstructions to arrive at a better estimate of past ENSO variance changes, we find robust evidence that the ENSO variance for any 30 yr period during the interval 1590–1880 was considerably lower than that observed during 1979–2009.

Description: Causes of Greenland temperature variability over the past 4000 yr: implications for northern hemispheric temperature changes Climate of the Past, 9, 2299-2317, 2013 Author(s): T. Kobashi, K. Goto-Azuma, J. E. Box, C.-C. Gao, and T. Nakaegawa Precise understanding of Greenland temperature variability is important in two ways. First, Greenland ice sheet melting associated with rising temperature is a major global sea level forcing, potentially affecting large populations in coming centuries. Second, Greenland temperatures are highly affected by North Atlantic Oscillation/Arctic Oscillation (NAO/AO) and Atlantic multidecadal oscillation (AMO). In our earlier study, we found that Greenland temperature deviated negatively (positively) from northern hemispheric (NH) temperature trend during stronger (weaker) solar activity owing to changes in atmospheric/oceanic changes (e.g. NAO/AO) over the past 800 yr (Kobashi et al., 2013). Therefore, a precise Greenland temperature record can provide important constraints on the past atmospheric/oceanic circulation in the region and beyond. Here, we investigated Greenland temperature variability over the past 4000 yr reconstructed from argon and nitrogen isotopes from trapped air in a GISP2 ice core, using a one-dimensional energy balance model with orbital, solar, volcanic, greenhouse gas, and aerosol forcings. The modelled northern Northern Hemisphere (NH) temperature exhibits a cooling trend over the past 4000 yr as observed for the reconstructed Greenland temperature through decreasing annual average insolation. With consideration of the negative influence of solar variability, the modelled and observed Greenland temperatures agree with correlation coefficients of r = 0.34–0.36 ( p = 0.1–0.04) in 21 yr running means (RMs) and r = 0.38–0.45 ( p = 0.1–0.05) on a centennial timescale (101 yr RMs). Thus, the model can explain 14 to 20% of variance of the observed Greenland temperature in multidecadal to centennial timescales with a 90–96% confidence interval, suggesting that a weak but persistent negative solar influence on Greenland temperature continued over the past 4000 yr. Then, we estimated the distribution of multidecadal NH and northern high-latitude temperatures over the past 4000 yr constrained by the climate model and Greenland temperatures. Estimated northern NH temperature and NH average temperature from the model and the Greenland temperature agree with published multi-proxy temperature records with r = 0.35–0.60 in a 92–99% confidence interval over the past 2000 yr. We found that greenhouse gases played two important roles over the past 4000 yr for the rapid warming during the 20th century and slightly cooler temperature during the early period of the past 4000 yr. Lastly, our analysis indicated that the current average temperature (1990–2010) or higher temperatures occurred at a frequency of 1.3 times per 1000 yr for northern high latitudes and 0.36 times per 4000 yr for NH temperatures, respectively, indicating that the current multidecadal NH temperature (1990–2010) is more likely unprecedented than not ( p = 0.36) for the past 4000 yr.

Description: Different ocean states and transient characteristics in Last Glacial Maximum simulations and implications for deglaciation Climate of the Past, 9, 2319-2333, 2013 Author(s): X. Zhang, G. Lohmann, G. Knorr, and X. Xu The last deglaciation is one of the best constrained global-scale climate changes documented by climate archives. Nevertheless, understanding of the underlying dynamics is still limited, especially with respect to abrupt climate shifts and associated changes in the Atlantic meridional overturning circulation (AMOC) during glacial and deglacial periods. A fundamental issue is how to obtain an appropriate climate state at the Last Glacial Maximum (LGM, 21 000 yr before present, 21 ka BP) that can be used as an initial condition for deglaciation. With the aid of a comprehensive climate model, we found that initial ocean states play an important role on the equilibrium timescale of the simulated glacial ocean. Independent of the initialization, the climatological surface characteristics are similar and quasi-stationary, even when trends in the deep ocean are still significant, which provides an explanation for the large spread of simulated LGM ocean states among the Paleoclimate Modeling Intercomparison Project phase 2 (PMIP2) models. Accordingly, we emphasize that caution must be taken when alleged quasi-stationary states, inferred on the basis of surface properties, are used as a reference for both model inter-comparison and data model comparison. The simulated ocean state with the most realistic AMOC is characterized by a pronounced vertical stratification, in line with reconstructions. Hosing experiments further suggest that the response of the glacial ocean is dependent on the ocean background state, i.e. only the state with robust stratification shows an overshoot behavior in the North Atlantic. We propose that the salinity stratification represents a key control on the AMOC pattern and its transient response to perturbations. Furthermore, additional experiments suggest that the stratified deep ocean formed prior to the LGM during a time of minimum obliquity (~ 27 ka BP). This indicates that changes in the glacial deep ocean already occur before the last deglaciation. In combination, these findings represent a new paradigm for the LGM and the last deglaciation, which challenges the conventional evaluation of glacial and deglacial AMOC changes based on an ocean state derived from 21 ka BP boundary conditions.

Description: Cluster as current sheet surveyor in the magnetotail Annales Geophysicae, 31, 1605-1610, 2013 Author(s): Y. Narita, R. Nakamura, and W. Baumjohann A novel analysis technique is presented to estimate the current sheet thickness unambiguously and directly, without associating time series data with spatial structure. The technique is a combination of eigenvalue analysis and minimum variance estimator adapted to Harris current sheet geometry, and needs one-time, four-point magnetic field data as provided by the Cluster spacecraft. Two current sheet parameters, thickness and distance to the spacecraft, can be determined at each time step of the magnetic field measurements. An example is shown from a Cluster magnetotail crossing under quiet magnetospheric conditions, yielding the result that the current sheet thickness is on the scale of the proton gyroradius. The analysis technique can also be used to track the dynamical evolution of the current sheet structure in three dimensions.

Description: Changes in the ultra-low frequency wave field during the precursor phase to the Sichuan earthquake: DEMETER observations Annales Geophysicae, 31, 1597-1603, 2013 Author(s): S. N. Walker, V. Kadirkamanathan, and O. A. Pokhotelov Electromagnetic phenomena observed in association with increases in seismic activity have been studied for several decades. These phenomena are generated during the precursory phases of an earthquake as well as during the main event. Their occurrence during the precursory phases may be used in short-term prediction of a large earthquake. In this paper, we examine ultra-low frequency (ULF) electric field data from the DEMETER satellite during the period leading up to the Sichuan earthquake. It is shown that there is an increase in ULF wave activity observed as DEMETER passes in the vicinity of the earthquake epicentre. This increase is most obvious at lower frequencies. Examination of the ULF spectra shows the possible occurrence of geomagnetic pearl pulsations, resulting from the passage of atmospheric gravity waves generated in the vicinity of the earthquake epicentre.

Description: Determination of wave vectors using the phase differencing method Annales Geophysicae, 31, 1611-1617, 2013 Author(s): S. N. Walker and I. Moiseenko Due to the collisionless nature of space plasmas, plasma waves play an important role in the redistribution of energy between the various particle populations in many regions of geospace. In order to fully comprehend such mechanisms it is necessary to characterise the nature of the waves present. This involves the determination of properties such as wave vector k . There are a number of methods used to calculate k based on the multipoint measurements that are now available. These methods rely on the fact that the same wave packet is simultaneously observed at two or more locations whose separation is small in comparison to the correlation length of the wave packet. This limitation restricts the analysis to low frequency (MHD) waves. In this paper we propose an extension to the phase differencing method to enable the correlation of measurements that were not made simultaneously but differ temporally by a number of wave periods. The method is illustrated using measurements of magnetosonic waves from the Cluster STAFF search coil magnetometer. It is shown that it is possible to identify wave packets whose coherence length is much less than the separation between the measurement locations. The resulting dispersion is found to agree with theoretical results.

Description: Effective CO 2 lifetime and future CO 2 levels based on fit function Annales Geophysicae, 31, 1591-1596, 2013 Author(s): G. R. Sonnemann and M. Grygalashvyly The estimated global CO 2 emission rates and the measured atmospheric CO 2 concentrations show that only a certain share of the emitted CO 2 accumulates in the atmosphere. For given atmospheric emissions of CO 2 , the effective lifetime determines its accumulation in the atmosphere and, consequently, its impact on the future global warming. We found that on average the inferred effective lifetime of CO 2 decreases as its atmospheric concentration increases, reducing the rate of its accumulation in the atmosphere. We derived a power function that fits the varying lifetimes. Based on this fitting function, we calculated the increase of CO 2 for different scenarios of future global emission rates.

Description: Why could ice ages be unpredictable? Climate of the Past, 9, 2253-2267, 2013 Author(s): M. Crucifix It is commonly accepted that the variations of Earth's orbit and obliquity control the timing of Pleistocene glacial–interglacial cycles. Evidence comes from power spectrum analysis of palaeoclimate records and from inspection of the timing of glacial and deglacial transitions. However, we do not know how tight this control is. Is it, for example, conceivable that random climatic fluctuations could cause a delay in deglaciation, bad enough to skip a full precession or obliquity cycle and subsequently modify the sequence of ice ages? To address this question, seven previously published conceptual models of ice ages are analysed by reference to the notion of generalised synchronisation. Insight is being gained by comparing the effects of the astronomical forcing with idealised forcings composed of only one or two periodic components. In general, the richness of the astronomical forcing allows for synchronisation over a wider range of parameters, compared to periodic forcing. Hence, glacial cycles may conceivably have remained paced by the astronomical forcing throughout the Pleistocene. However, all the models examined here show regimes of strong structural dependence on parameters. This means that small variations in parameters or random fluctuations may cause significant shifts in the succession of ice ages. Whether the actual system actually resides in such a regime depends on the amplitude of the effects associated with the astronomical forcing, which significantly differ across the different models studied here. The possibility of synchronisation on eccentricity is also discussed and it is shown that a high Rayleigh number on eccentricity, as recently found in observations, is no guarantee of reliable synchronisation.

Description: Enhanced EISCAT UHF backscatter during high-energy auroral electron precipitation Annales Geophysicae, 31, 1681-1687, 2013 Author(s): N. M. Schlatter, N. Ivchenko, T. Sergienko, B. Gustavsson, and B. U. E. Brändström Natural enhancements in the backscattered power of incoherent scatter radars up to 5 orders of magnitudes above the thermal backscatter are sometimes observed at high latitudes. Recently observations of enhancements in the backscattered power including a feature at zero Doppler shift have been reported. These enhancements are limited in altitude to tens of kilometers. The zero Doppler shift feature has been interpreted as a signature of electron density cavitation. Enhanced plasma lines during these observations have also been reported. We report on the first EISCAT UHF observations of enhanced backscattered radar power including a zero Doppler shift feature. The enhancements originated from two distinct and intermittent layers at about 200 km altitude. The altitude extent of the enhancements, observed during auroral high-energy electron precipitation, was 〈 2 km.

Description: Seasonal dependence of the longitudinal variations of nighttime ionospheric electron density and equivalent winds at southern midlatitudes Annales Geophysicae, 31, 1699-1708, 2013 Author(s): X. Luan and X. Dou It has been indicated that the observed Weddell Sea anomaly (WSA) appeared to be an extreme manifestation of the longitudinal variations in the Southern Hemisphere, since the WSA is characterized by greater evening electron density than the daytime density in the region near the Weddell Sea. In the present study, the longitudinal variations of the nighttime F2-layer peak electron density at southern midlatitudes are analyzed using the observations of the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) satellites between 2006 and 2008. It is found that significant longitudinal difference (〉 150%) relative to the minimum density at each local time prevails in all seasons, although the WSA phenomenon is only evident in summer under this solar minimum condition. Another interesting feature is that in summer, the maximum longitudinal differences occur around midnight (~ 23:00–00:00 LT) rather than in the evening (19:00–21:00 LT) in the evening, when the most prominent electron density enhancement occurs for the WSA phenomenon. Thus the seasonal–local time patterns of the electron density longitudinal variations during nighttime at southern midlatitudes cannot be simply explained in terms of the WSA. Meanwhile, the variations of the geomagnetic configuration and the equivalent magnetic meridional winds/upward plasma drifts are analyzed to explore their contributions to the longitudinal variations of the nighttime electron density. The maximum longitudinal differences are associated with the strongest wind-induced vertical plasma drifts after 21:00 LT in the Western Hemisphere. Besides the magnetic declination–zonal wind effects, the geographic meridional winds and the magnetic inclination also have significant effects on the upward plasma drifts and the resultant electron density.

Description: Extracting gravity wave parameters during the September 2002 Southern Hemisphere major sudden stratospheric warming using a SANAE imaging riometer Annales Geophysicae, 31, 1709-1719, 2013 Author(s): N. Mbatha, V. Sivakumar, H. Bencherif, and S. Malinga Using absorption data measured by imaging riometer for ionospheric studies (IRIS) located at the South Africa National Antarctic Expedition (SANAE), Antarctica (72° S, 3° W), we extracted the parameters of gravity waves (GW) of periods between 40 and 50 min during late winter/spring of the year 2002, a period of the unprecedented major sudden stratospheric warming (SSW) in the Southern Hemisphere middle atmosphere. During this period, an unprecedented substantial increase of temperature by about 25–30 K throughout the stratosphere was observed. During the period of the occurrence of the major stratospheric warming, there was a reduction of both the GW horizontal phase speeds and the horizontal wavelengths at 90 km. The GW phase speeds and horizontal wavelengths were observed to reach minimum values of about 7 m s −1 and 19 km, respectively, while during the quiet period the average value of the phase speed and horizontal wavelength was approximately 23 m s −1 and 62 km, respectively. The observed event is discussed in terms of momentum flux and also a potential interaction of gravity waves, planetary waves and mean circulation.

Description: Using ice-flow models to evaluate potential sites of million year-old ice in Antarctica Climate of the Past, 9, 2335-2345, 2013 Author(s): B. Van Liefferinge and F. Pattyn Finding suitable potential sites for an undisturbed record of million-year old ice in Antarctica requires slow-moving ice (preferably an ice divide) and basal conditions that are not disturbed by large topographic variations. Furthermore, ice should be thick and cold basal conditions should prevail, since basal melting would destroy the bottom layers. However, thick ice (needed to resolve the signal at sufficient high resolution) increases basal temperatures, which is a conflicting condition for finding a suitable drill site. In addition, slow moving areas in the center of ice sheets are also low-accumulation areas, and low accumulation reduces potential cooling of the ice through vertical advection. While boundary conditions such as ice thickness and accumulation rates are relatively well constrained, the major uncertainty in determining basal thermal conditions resides in the geothermal heat flow (GHF) underneath the ice sheet. We explore uncertainties in existing GHF data sets and their effect on basal temperatures of the Antarctic Ice Sheet, and propose an updated method based on Pattyn (2010) to improve existing GHF data sets in agreement with known basal temperatures and their gradients to reduce this uncertainty. Both complementary methods lead to a better comprehension of basal temperature sensitivity and a characterization of potential ice coring sites within these uncertainties. The combination of both modeling approaches show that the most likely oldest ice sites are situated near the divide areas (close to existing deep drilling sites, but in areas of smaller ice thickness) and across the Gamburtsev Subglacial Mountains.

Description: Pre-LGM Northern Hemisphere ice sheet topography Climate of the Past, 9, 2365-2378, 2013 Author(s): J. Kleman, J. Fastook, K. Ebert, J. Nilsson, and R. Caballero We here reconstruct the paleotopography of Northern Hemisphere ice sheets during the glacial maxima of marine isotope stages (MIS) 5b and 4.We employ a combined approach, blending geologically based reconstruction and numerical modeling, to arrive at probable ice sheet extents and topographies for each of these two time slices. For a physically based 3-D calculation based on geologically derived 2-D constraints, we use the University of Maine Ice Sheet Model (UMISM) to calculate ice sheet thickness and topography. The approach and ice sheet modeling strategy is designed to provide robust data sets of sufficient resolution for atmospheric circulation experiments for these previously elusive time periods. Two tunable parameters, a temperature scaling function applied to a spliced Vostok–GRIP record, and spatial adjustment of the climatic pole position, were employed iteratively to achieve a good fit to geological constraints where such were available. The model credibly reproduces the first-order pattern of size and location of geologically indicated ice sheets during marine isotope stages (MIS) 5b (86.2 kyr model age) and 4 (64 kyr model age). From the interglacial state of two north–south obstacles to atmospheric circulation (Rocky Mountains and Greenland), by MIS 5b the emergence of combined Quebec–central Arctic and Scandinavian–Barents-Kara ice sheets had increased the number of such highland obstacles to four. The number of major ice sheets remained constant through MIS 4, but the merging of the Cordilleran and the proto-Laurentide Ice Sheet produced a single continent-wide North American ice sheet at the LGM.

Description: GNSS troposphere tomography based on two-step reconstructions using GPS observations and COSMIC profiles Annales Geophysicae, 31, 1805-1815, 2013 Author(s): P. Xia, C. Cai, and Z. Liu Traditionally, balloon-based radiosonde soundings are used to study the spatial distribution of atmospheric water vapour. However, this approach cannot be frequently employed due to its high cost. In contrast, GPS tomography technique can obtain water vapour in a high temporal resolution. In the tomography technique, an iterative or non-iterative reconstruction algorithm is usually utilised to overcome rank deficiency of observation equations for water vapour inversion. However, the single iterative or non-iterative reconstruction algorithm has their limitations. For instance, the iterative reconstruction algorithm requires accurate initial values of water vapour while the non-iterative reconstruction algorithm needs proper constraint conditions. To overcome these drawbacks, we present a combined iterative and non-iterative reconstruction approach for the three-dimensional (3-D) water vapour inversion using GPS observations and COSMIC profiles. In this approach, the non-iterative reconstruction algorithm is first used to estimate water vapour density based on a priori water vapour information derived from COSMIC radio occultation data. The estimates are then employed as initial values in the iterative reconstruction algorithm. The largest advantage of this approach is that precise initial values of water vapour density that are essential in the iterative reconstruction algorithm can be obtained. This combined reconstruction algorithm (CRA) is evaluated using 10-day GPS observations in Hong Kong and COSMIC profiles. The test results indicate that the water vapor accuracy from CRA is 16 and 14% higher than that of iterative and non-iterative reconstruction approaches, respectively. In addition, the tomography results obtained from the CRA are further validated using radiosonde data. Results indicate that water vapour densities derived from the CRA agree with radiosonde results very well at altitudes above 2.5 km. The average RMS value of their differences above 2.5 km is 0.44 g m −3 .

Description: Corrigendum to "Estimates of aerosol absorption over India using multi-satellite retrieval" published in Ann. Geophys., 31, 1773–1778, 2013 Annales Geophysicae, 31, 1791-1791, 2013 Author(s): D. Narasimhan and S. K. Satheesh No abstract available.

Description: Eurasian Arctic climate over the past millennium as recorded in the Akademii Nauk ice core (Severnaya Zemlya) Climate of the Past, 9, 2379-2389, 2013 Author(s): T. Opel, D. Fritzsche, and H. Meyer Understanding recent Arctic climate change requires detailed information on past changes, in particular on a regional scale. The extension of the depth–age relation of the Akademii Nauk (AN) ice core from Severnaya Zemlya (SZ) to the last 1100 yr provides new perspectives on past climate fluctuations in the Barents and Kara seas region. Here, we present the easternmost high-resolution ice-core climate proxy records (δ 18 O and sodium) from the Arctic. Multi-annual AN δ 18 O data as near-surface air-temperature proxies reveal major temperature changes over the last millennium, including the absolute minimum around 1800 and the unprecedented warming to a double-peak maximum in the early 20th century. The long-term cooling trend in δ 18 O is related to a decline in summer insolation but also to the growth of the AN ice cap as indicated by decreasing sodium concentrations. Neither a pronounced Medieval Climate Anomaly nor a Little Ice Age are detectable in the AN δ 18 O record. In contrast, there is evidence of several abrupt warming and cooling events, such as in the 15th and 16th centuries, partly accompanied by corresponding changes in sodium concentrations. These abrupt changes are assumed to be related to sea-ice cover variability in the Barents and Kara seas region, which might be caused by shifts in atmospheric circulation patterns. Our results indicate a significant impact of internal climate variability on Arctic climate change in the last millennium.

Description: Re-evaluation of the age model for North Atlantic Ocean Site 982 – arguments for a return to the original chronology Climate of the Past, 9, 2391-2397, 2013 Author(s): K. T. Lawrence, I. Bailey, and M. E. Raymo Recently, the veracity of the published chronology for the Pliocene section of North Atlantic Ocean Drilling Program Site 982 was called into question. Here, we examine the robustness of the original age model as well as the proposed age model revision. The proposed revision is predicated on an apparent misidentification of the depth to the Gauss–Matuyama (G/M) polarity chronozone reversal boundary (2.581 Ma) based on preliminary shipboard paleomagnetic data, and offers a new chronology that includes a hiatus between ~3.2 and 3 Ma. However, an even more accurate shore-based, u-channel-derived polarity chronozone stratigraphy for the past ~2.7 Ma supports the shipboard composite stratigraphy and demonstrates that the original estimate of the depth of the G/M reversal in the Site 982 record is correct. Thus, the main justification forwarded to support the revised chronology no longer exists. We demonstrate that the proposed revision results in a pronounced anomaly in sedimentation rates proximal to the proposed hiatus, erroneous assignment of marine-isotope stages in the Site 982 Pliocene benthic stable oxygen isotope stratigraphy, and a markedly worse correlation of proxy records between this site and other regional paleoclimate data. We conclude that the original chronology for Site 982 is a far more accurate age model than that which arises from the published revision. We strongly recommend the use of the original chronology for all future work at Site 982.

Description: Frequency variations of gravity waves interacting with a time-varying tide Annales Geophysicae, 31, 1731-1743, 2013 Author(s): C. M. Huang, S. D. Zhang, F. Yi, K. M. Huang, Y. H. Zhang, Q. Gan, and Y. Gong Using a nonlinear, 2-D time-dependent numerical model, we simulate the propagation of gravity waves (GWs) in a time-varying tide. Our simulations show that when a GW packet propagates in a time-varying tidal-wind environment, not only its intrinsic frequency but also its ground-based frequency would change significantly. The tidal horizontal-wind acceleration dominates the GW frequency variation. Positive (negative) accelerations induce frequency increases (decreases) with time. More interestingly, tidal-wind acceleration near the critical layers always causes the GW frequency to increase, which may partially explain the observations that high-frequency GW components are more dominant in the middle and upper atmosphere than in the lower atmosphere. The combination of the increased ground-based frequency of propagating GWs in a time-varying tidal-wind field and the transient nature of the critical layer induced by a time-varying tidal zonal wind creates favorable conditions for GWs to penetrate their originally expected critical layers. Consequently, GWs have an impact on the background atmosphere at much higher altitudes than expected, which indicates that the dynamical effects of tidal–GW interactions are more complicated than usually taken into account by GW parameterizations in global models.

Description: Holocene climate variability in the winter rainfall zone of South Africa Climate of the Past, 9, 2347-2364, 2013 Author(s): S. Weldeab, J.-B. W. Stuut, R. R. Schneider, and W. Siebel We established a multi-proxy time series comprising analyses of major elements in bulk sediments, Sr and Nd isotopes, grain size of terrigenous fraction, and δ 18 O and δ 13 C in tests of Neogloboquadrina pachyderma (sinistral) from a marine sediment sequence recovered off the Orange River. The records reveal coherent patterns of variability that reflect changes in wind strength, precipitation over the river catchments, and upwelling of cold and nutrient-rich coastal waters off western South Africa. The wettest episode of the Holocene in the winter rainfall zone (WRZ) of South Africa occurred during the "Little Ice Age" (700–100 cal years BP) most likely in response to a northward shift of the austral westerlies. Wet phases and strengthened coastal water upwellings are companied by a decrease of Agulhas water leakage into the South Atlantic and a reduced dust incursion over Antarctica, as indicated in previous studies. A continuous aridification trend in the WRZ and a weakening of the southern Benguela Upwelling System (BUS) between 9000 and 5500 cal years BP parallel with increase of dust deposition over Antarctica and an enhanced leakage of warm Agulhas water into the eastern South Atlantic. The temporal relationship between precipitation changes in the WRZ, the thermal state of the coastal surface water, and leakage of warm water in the South Atlantic, and variation of dust incursion over Antarctica suggests a causal link that most likely was related to latitudinal shifts of the Southern Hemisphere westerlies. Our results of the mid-Holocene time interval may serve as an analogue to a possible long-term consequence of the current and future southward shift of the westerlies. Furthermore, warming of the coastal surface water as a result of warm Agulhas water incursion into the southern BUS may affect coastal fog formation.