Atmospheric Soluble Dust Records from a Tibetan Ice Core: Possible Climate Proxies and Teleconnection With the Pacific Decadal Oscillation

In autumn 2005, a joint expedition between the University of Maine and the Institute of Tibetan Plateau Research recovered three ice cores from Guoqu Glacier (33 degrees 34'37.80 '' N, 91 degrees 10'35.3 '' E, 5720 m above sea level) on the northern side of Mt. Geladaindong, central Tibetan Plateau. Isotopes ( delta(18)O), major soluble ions (Na(+), K(+), Mg(2+), Ca(2+), Cl(-), NO(3)(-), SO(4)(2-)), and radionuclide (beta-activity) measurements from one of the cores revealed a 70-year record (1935-2005). Statistical analysis of major ion time series suggests that atmospheric soluble dust species dominate the chemical signature and that background dust levels conceal marine ion species deposition. The soluble dust time series have interspecies relations and common structure (empirical orthogonal function (EOF) 1), suggesting a similar soluble dust source or transport route. Annual and seasonal correlations between the EOF 1 time series and National Centers for Environmental Prediction/National Center for Atmospheric Research reanalysis climate variables (1948-2004) suggest that the Mt. Geladaindong ice core record provides a proxy for local and regional surface pressure. An approximately threefold decrease of soluble dust concentrations in the middle to late 1970s, accompanied by regional increases in pressure and temperature and decreases in wind velocity, coincides with the major 1976-1977 shift of the Pacific Decadal Oscillation (PDO) from a negative to a positive state. This is the first ice core evidence of a potential teleconnection between central Asian atmospheric soluble dust loading and the PDO. Analysis of temporally longer ice cores from Mt. Geladaindong may enhance understanding of the relationship between the PDO and central Asian atmospheric circulation and subsequent atmospheric soluble dust loading.

Association Between Atmospheric Circulation Patterns and Firn-Ice Core Records from the Inilchek Glacierized Area, Central Tien Shan, Asia

Glacioclimatological research in the central Tien Shan was performed in the summers of 1998 and 1999 on the South Inilchek Glacier at 5100 - 5460 m. A 14.36 m firn-ice core and snow samples were collected and used for stratigraphic, isotopic, and chemical analyses. The firn-ice core and snow records were related to snow pit measurements at an event scale and to meteorological data and synoptic indices of atmospheric circulation at annual and seasonal scales. Linear relationships between the seasonal air temperature and seasonal isotopic composition in accumulated precipitation were established. Changes in the delta(18)O air temperature relationship, in major ion concentration and in the ratios between chemical species, were used to identify different sources of moisture and investigate changes in atmospheric circulation patterns. Precipitation over the central Tien Shan is characterized by the lowest ionic content among the Tien Shan glaciers and indicates its mainly marine origin. In seasons of minimum precipitation, autumn and winter, water vapor was derived from the arid and semiarid regions in central Eurasia and contributed annual maximal solute content to snow accumulation in Tien Shan. The lowest content of major ions was observed in spring and summer layers, which represent maximum seasonal accumulation when moisture originates over the Atlantic Ocean and Mediterranean and Black Seas.

A Record of Atmospheric Co2 During the Last 40,000 Years from the Siple Dome, Antarctica Ice Core

We have measured the CO2 concentration of air occluded during the last 40,000 years in the deep Siple Dome A ( hereafter Siple Dome) ice core, Antarctica. The general trend of CO2 concentration from Siple Dome ice follows the temperature inferred from the isotopic composition of the ice and is mostly in agreement with other Antarctic ice core CO2 records. CO2 rose initially at similar to 17.5 kyr B. P. ( thousand years before 1950), decreased slowly during the Antarctic Cold Reversal, rose during the Younger Dryas, fell to a local minimum at around 8 kyr B. P., and rose continuously since then. The CO2 concentration never reached steady state during the Holocene, as also found in the Taylor Dome and EPICA Dome C ( hereafter Dome C) records. During the last glacial termination, a lag of CO2 versus Siple Dome isotopic temperature is probable. The Siple Dome CO2 concentrations during the last glacial termination and in the Holocene are at certain times greater than in other Antarctic ice cores by up to 20 ppm (mumol CO2/mol air). While in situ production of CO2 is one possible cause of the sporadic elevated levels, the mechanism leading to the enrichment is not yet clear.

A 700-Year Record of Atmospheric Circulation Developed from the Law Dome Ice Core, East Antarctica

A 700-year, high-resolution, multivariate ice core record from Dome Summit South (DSS) (66degrees46'S, 112degrees48'E; 1370 m), Law Dome, is used to investigate sea level pressure (SLP) variability in the region of East Antarctica. Empirical orthogonal function (EOF) analysis reveals that the first EOF (LDEOF1) of the combined glaciochemical, oxygen isotope ratio, and accumulation rate record from DSS represents most of the variability in sea salt seen in the record. LDEOF1 is positively correlated (at least 95% confidence level) to instrumental June mean SLP across most of East Antarctica. Over the last 700 years, LDEOF1 levels at Law Dome were the highest during the nineteenth century, suggesting an increase in intensification of winter circulation during this period. The Law Dome DSS oxygen isotope ratio series also indicates that the nineteenth century had the coldest winters of any century in the record. In contrast, LDEOF1 levels were the lowest at Law Dome during the eighteenth century, suggesting a significant shift in the patterns and/or intensity of East Antarctic atmospheric circulation between the eighteenth and the nineteenth centuries. The LDEOF1 sea salt record is characterized by significant decadal-scale variability with a strong 25-year periodic structure.

Ice Layers as an Indicator of Summer Warmth and Atmospheric Blocking in Alaska

Samples were collected from a snow pit and shallow urn core near Kahiltna Pass (2970 m a.s.l.), Denali National Park, Alaska, USA, in May 2008. The record spans autumn 2003 to spring 2008 and reveals clusters of ice layers interpreted as summertime intervals of above-freezing temperatures. High correlation coefficients (0.75-1.00) between annual ice-layer thickness and regional summertime station temperatures for 4 years (n=4) indicate ice-layer thickness is a good proxy for mean and extreme summertime temperatures across Alaska, at least over the short period of record. A Rex-block (aka high-over-low) pattern, a downstream trough over Hudson Bay, Canada, and an upstream trough over eastern Siberia occurred during the three melting events that lasted at least 2 weeks. About half of all shorter melting events were associated with a cut-off low traversing the Gulf of Alaska. We hypothesize that a surface-to-bedrock core extracted from this location would provide a high-quality record of summer temperature and atmospheric blocking variability for the last several hundred years.

Influence of Equatorial Diatom Processes on Si Deposition and Atmospheric CO(2) Cycles at Glacial/Interglacial Timescales

The causes of the glacial cycle remain unknown, although the primary driver is changes in atmospheric CO(2), likely controlled by the biological pump and biogeochemical cycles. The two most important regions of the ocean for exchange of CO(2) with the atmosphere are the equatorial Pacific and the Southern Ocean ( SO), the former a net source and the latter a net sink under present conditions. The equatorial Pacific has been shown to be a Si(OH)(4)-limited ecosystem, a consequence of the low source Si(OH)(4) concentrations in upwelled water that has its origin in the SO. This teleconnection for nutrients between the two regions suggests an oscillatory relationship that may influence or control glacial cycles. Opal mass accumulation rate (MAR) data and delta(15)N measurements in equatorial cores are interpreted with predictions from a one- dimensional Si(OH)(4)-limited ecosystem model (CoSINE) for the equatorial Pacific. The results suggest that equatorial Pacific surface CO(2) processes are in opposite phase to that of the global atmosphere, providing a negative feedback to the glacial cycle. This negative feedback is implemented through the effect of the SO on the equatorial Si(OH)(4) supply. An alternative hypothesis, that the whole ocean becomes Si(OH)(4) poor during cooling periods, is suggested by low opal MAR in cores from both equatorial and Antarctic regions, perhaps as a result of low river input. terminations in this scenario would result from blooms of coccolithophorids triggered by low Si(OH)(4) concentrations.

Diurnal variations in middle-atmospheric water vapor by ground-based microwave radiometry

Abstract. In this paper, we compare the diurnal variations in middle-atmospheric water vapor as measured by two ground-based microwave radiometers in the Alpine region near Bern, Switzerland. The observational data set is also compared to data from the chemistry–climate model WACCM. Due to the small diurnal variations of usually less than 1%, averages over extended time periods are required. Therefore, two time periods of five months each, December to April and June to October, were taken for the comparison. The diurnal variations from the observational data agree well with each other in amplitude and phase. The linear correlation coefficients range from 0.8 in the upper stratosphere to 0.5 in the upper mesosphere. The observed diurnal variability is significant at all pressure levels within the sensitivity of the instruments. Comparing our observations with WACCM, we find that the agreement of the phase of the diurnal cycle between observations and model is better from December to April than from June to October. The amplitudes of the diurnal variations for both time periods increase with altitude in WACCM, but remain approximately constant at 0.05 ppm in the observations. The WACCM data are used to separate the processes that lead to diurnal variations in middle-atmospheric water vapor above Bern. The dominating processes were found to be meridional advection below 0.1 hPa, vertical advection between 0.1 and 0.02 hPa and (photo-)chemistry above 0.02 hPa. The contribution of zonal advection is small. The highest diurnal variations in water vapor as seen in the WACCM data are found in the mesopause region during the time period from June to October with diurnal amplitudes of 0.2 ppm (approximately 5% in relative units).

Influence of microphysical schemes on atmospheric water in the Weather Research and Forecasting model

This study examines how different microphysical parameterization schemes influence orographically induced precipitation and the distributions of hydrometeors and water vapour for midlatitude summer conditions in the Weather Research and Forecasting (WRF) model. A high-resolution two-dimensional idealized simulation is used to assess the differences between the schemes in which a moist air flow is interacting with a bell-shaped 2 km high mountain. Periodic lateral boundary conditions are chosen to recirculate atmospheric water in the domain. It is found that the 13 selected microphysical schemes conserve the water in the model domain. The gain or loss of water is less than 0.81% over a simulation time interval of 61 days. The differences of the microphysical schemes in terms of the distributions of water vapour, hydrometeors and accumulated precipitation are presented and discussed. The Kessler scheme, the only scheme without ice-phase processes, shows final values of cloud liquid water 14 times greater than the other schemes. The differences among the other schemes are not as extreme, but still they differ up to 79% in water vapour, up to 10 times in hydrometeors and up to 64% in accumulated precipitation at the end of the simulation. The microphysical schemes also differ in the surface evaporation rate. The WRF single-moment 3-class scheme has the highest surface evaporation rate compensated by the highest precipitation rate. The different distributions of hydrometeors and water vapour of the microphysical schemes induce differences up to 49 W m−2 in the downwelling shortwave radiation and up to 33 W m−2 in the downwelling longwave radiation.

Impact of geomagnetic excursions on atmospheric chemistry and dynamics

Geomagnetic excursions, i.e. short periods in time with much weaker geomagnetic fields and substantial changes in the position of the geomagnetic pole, occurred repeatedly in the Earth's history, e.g. the Laschamp event about 41 kyr ago. Although the next such excursion is certain to come, little is known about the timing and possible consequences for the state of the atmosphere and the ecosystems. Here we use the global chemistry climate model SOCOL-MPIOM to simulate the effects of geomagnetic excursions on atmospheric ionization, chemistry and dynamics. Our simulations show significantly increased concentrations of nitrogen oxides (NOx) in the entire stratosphere, especially over Antarctica (+15%), due to enhanced ionization by galactic cosmic rays. Hydrogen oxides (HOx) are also produced in greater amounts (up to +40%) in the tropical and subtropical lower stratosphere, while their destruction by reactions with enhanced NOx prevails over the poles and in high altitudes (by −5%). Stratospheric ozone concentrations decrease globally above 20 km by 1–2% and at the northern hemispheric tropopause by up to 5% owing to the accelerated NOx-induced destruction. A 5% increase is found in the southern lower stratosphere and troposphere. In response to these changes in ozone and the concomitant changes in atmospheric heating rates, the Arctic vortex intensifies in boreal winter, while the Antarctic vortex weakens in austral winter and spring. Surface wind anomalies show significant intensification of the southern westerlies at their poleward edge during austral winter and a pronounced northward shift in spring. Major impacts on the global climate seem unlikely.

ANALYTICAL MODELS OF EXOPLANETARY ATMOSPHERES. I. ATMOSPHERIC DYNAMICS VIA THE SHALLOW WATER SYSTEM

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical, and spherical), rotation, magnetic tension, and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag, and magnetic drag), and magnetic tension are included. The global atmospheric structure is largely controlled by a single key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations, and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to develop decisively physical intuition and as a reference point for three-dimensional magnetohydrodynamic simulations of atmospheric circulation.

Atmospheric Dynamics of Hot Exoplanets

The characterization of exoplanetary atmospheres has come of age in the last decade, as astronomical techniques now allow for albedos, chemical abundances, temperature profiles and maps, rotation periods and even wind speeds to be measured. Atmospheric dynamics sets the background state of density, temperature and velocity that determines or influences the spectral and temporal appearance of an exoplanetary atmosphere. Hot exoplanets are most amenable to these characterization techniques; in the present review, we focus on highly-irradiated, large exoplanets (the "hot Jupiters"), as astronomical data begin to confront theoretical questions. We summarize the basic atmospheric quantities inferred from the astronomical observations. We review the state of the art by addressing a series of current questions and look towards the future by considering a separate set of exploratory questions. Attaining the next level of understanding will require a concerted effort of constructing multi-faceted, multi-wavelength datasets for benchmark objects. Understanding clouds presents a formidable obstacle, as they introduce degeneracies into the interpretation of spectra, yet their properties and existence are directly influenced by atmospheric dynamics. Confronting general circulation models with these multi-faceted, multi-wavelength datasets will help us understand these and other degeneracies. The coming decade will witness a decisive confrontation of theory and simulation by the next generation of astronomical data.

CONSTRAINING THE ATMOSPHERIC COMPOSITION OF THE DAY-NIGHT TERMINATORS OF HD 189733b: ATMOSPHERIC RETRIEVAL WITH AEROSOLS

A number of observations have shown that Rayleigh scattering by aerosols dominates the transmission spectrum of HD 189733b at wavelengths shortward of 1 μm. In this study, we retrieve a range of aerosol distributions consistent with transmission spectroscopy between 0.3-24 μm that were recently re-analyzed by Pont et al. To constrain the particle size and the optical depth of the aerosol layer, we investigate the degeneracies between aerosol composition, temperature, planetary radius, and molecular abundances that prevent unique solutions for transit spectroscopy. Assuming that the aerosol is composed of MgSiO3, we suggest that a vertically uniform aerosol layer over all pressures with a monodisperse particle size smaller than about 0.1 μm and an optical depth in the range 0.002-0.02 at 1 μm provides statistically meaningful solutions for the day/night terminator regions of HD 189733b. Generally, we find that a uniform aerosol layer provide adequate fits to the data if the optical depth is less than 0.1 and the particle size is smaller than 0.1 μm, irrespective of the atmospheric temperature, planetary radius, aerosol composition, and gaseous molecules. Strong constraints on the aerosol properties are provided by spectra at wavelengths shortward of 1 μm as well as longward of 8 μm, if the aerosol material has absorption features in this region. We show that these are the optimal wavelengths for quantifying the effects of aerosols, which may guide the design of future space observations. The present investigation indicates that the current data offer sufficient information to constrain some of the aerosol properties of HD189733b, but the chemistry in the terminator regions remains uncertain.

Erratum: "Atmospheric retrieval analysis of the directly imaged exoplanet HR 8799 b"(2013, ApJ, 778, 97)

This is an Erratum for the article 2013 ApJ 778 97

ATMOSPHERIC RETRIEVAL ANALYSIS OF THE DIRECTLY IMAGED EXOPLANET HR 8799b

Directly imaged exoplanets are unexplored laboratories for the application of the spectral and temperature retrieval method, where the chemistry and composition of their atmospheres are inferred from inverse modeling of the available data. As a pilot study, we focus on the extrasolar gas giant HR 8799b, for which more than 50 data points are available. We upgrade our non-linear optimal estimation retrieval method to include a phenomenological model of clouds that requires the cloud optical depth and monodisperse particle size to be specified. Previous studies have focused on forward models with assumed values of the exoplanetary properties; there is no consensus on the best-fit values of the radius, mass, surface gravity, and effective temperature of HR 8799b. We show that cloud-free models produce reasonable fits to the data if the atmosphere is of super-solar metallicity and non-solar elemental abundances. Intermediate cloudy models with moderate values of the cloud optical depth and micron-sized particles provide an equally reasonable fit to the data and require a lower mean molecular weight. We report our best-fit values for the radius, mass, surface gravity, and effective temperature of HR 8799b. The mean molecular weight is about 3.8, while the carbon-to-oxygen ratio is about unity due to the prevalence of carbon monoxide. Our study emphasizes the need for robust claims about the nature of an exoplanetary atmosphere to be based on analyses involving both photometry and spectroscopy and inferred from beyond a few photometric data points, such as are typically reported for hot Jupiters.