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.

Inference of inhomogeneous clouds in an exoplaent atmosphere

We present new visible and infrared observations of the hot Jupiter Kepler-7b to determine its atmospheric properties. Our analysis allows us to (1) refine Kepler-7b's relatively large geometric albedo of Ag = 0.35 ± 0.02, (2) place upper limits on Kepler-7b thermal emission that remains undetected in both Spitzer bandpasses and (3) report a westward shift in the Kepler optical phase curve. We argue that Kepler-7b's visible flux cannot be due to thermal emission or Rayleigh scattering from H2 molecules. We therefore conclude that high altitude, optically reflective clouds located west from the substellar point are present in its atmosphere. We find that a silicate-based cloud composition is a possible candidate. Kepler-7b exhibits several properties that may make it particularly amenable to cloud formation in its upper atmosphere. These include a hot deep atmosphere that avoids a cloud cold trap, very low surface gravity to suppress cloud sedimentation, and a planetary equilibrium temperature in a range that allows for silicate clouds to potentially form in the visible atmosphere probed by Kepler. Our analysis does not only present evidence of optically thick clouds on Kepler-7b but also yields the first map of clouds in an exoplanet atmosphere.

The isotopic composition of pore-water from Tournemire argillite (France): An inter-comparison study.

The determination of stable isotope contents of pore-water from consolidated argillaceous rocks remains a critical issue. In order to understand the processes involved in techniques developed for acquiring stable isotope compositions of pore-water, a comparative study between different methods was based on core samples of the Tournemire argillite. It concerns two water extraction techniques based on vacuum distillation and two pore-water equilibration techniques (radial diffusion in liquid phase and diffusive exchange in vapor phase). The water-content values obtained from vacuum distillation at 50 °C are always the lowest, on average 8% lower than the values obtained by heating at 105 °C and 17% lower than the values obtained by heating at 150 °C. The amounts of pore-water estimated from vacuum distillation at 105 °C and 150 °C and from radial diffusion method are in good agreement with those determined by heating. On the contrary, the vapor exchange method provides the highest values of water contents. Concerning stable isotope data, a good agreement was found between those obtained by equilibration techniques and those of fracture water, especially for 2H. Vacuum distillation at high temperature (particularly at 150 °C) also provided results consistent with data of fracture fluids. On the other hand, distillation at 50 °C provides a systematic depletion in heavy isotopes (about –20‰ for 2H and –2.7‰ for 18O) that can be modelled by an incomplete Rayleigh-type distillation process.

Benchmarking the simulation of Cr isotope fractionation

A benchmark problem set consisting of four problem levels was developed for the simulation of Cr isotope fractionation in 1D and 2D domains. The benchmark is based on a recent field study where Cr(VI) reduction and accompanying Cr isotope fractionation occurs abiotically by an aqueous reaction with dissolved Fe 2+ (Wanner et al., 2012., Appl. Geochem., 27, 644–662). The problem set includes simulation of the major processes affecting the Cr isotopic composition such as the dissolution of various Cr(VI) bearing minerals, fractionation during abiotic aqueous Cr(VI) reduction, and non-fractionating precipitation of Cr(III) as sparingly soluble Cr-hydroxide. Accuracy of the presented solutions was ensured by running the problems with four well-established reactive transport modeling codes: TOUGHREACT, MIN3P, CRUNCHFLOW, and FLOTRAN. Results were also compared with an analytical Rayleigh-type fractionation model. An additional constraint on the correctness of the results was obtained by comparing output from the problem levels simulating Cr isotope fractionation with the corresponding ones only simulating bulk concentrations. For all problem levels, model to model comparisons showed excellent agreement, suggesting that for the tested geochemical processes any code is capable of accurately simulating the fate of individual Cr isotopes.

Geochemistry of ODP Leg 135 vitrophyric rhyolite samples

In-situ proton-microprobe analyses are presented for glasses, plagioclases, pyroxenes, olivines, and spinels in eleven samples from Sites 834-836, 839, and 841 (vitrophyric rhyolite), plus a Tongan dacite. Elements analyzed are Mn, Ni, Cu, Zn, Ga, Rb, Sr, Y, Zr, Pb, and Sn (in spinels only). The data are used to calculate two sets of partition coefficients, one set based on the ratio of element in mineral/element in coexisting glass. The second set of coefficients, thought to be more robust, is corrected by application of the Rayleigh fractionation equations, which requires additional use of modal data. Data are presented for phenocryst core-rim phases and microphenocryst-groundmass phases from a few samples. Comparison with published coefficients reveals an overall consistency with those presented here, but with some notable anomalies. Examples are relatively high Zr values for pyroxenes and abnormally low Mn values in olivines and clinopyroxenes from Site 839 lavas. Some anomalies may reflect kinetic effects, but interpretation of the coefficients is complicated, especially in olivines from Sites 836 and 839, by possible crystal-liquid disequilibrium resulting from mixing processes.

(Table 1) Water column biogenic silica concentrations and Si isotopic composition during Marion Dufresne cruise MD166 off South Africa

Southern Ocean biogeochemical processes have an impact on global marine primary production and global elemental cycling, e.g. by likely controlling glacial-interglacial pCO2 variation. In this context, the natural silicon isotopic composition (d30Si) of sedimentary biogenic silica has been used to reconstruct past Si-consumption:supply ratios in the surface waters. We present a new dataset in the Southern Ocean from a IPY-GEOTRACES transect (Bonus-GoodHope) which includes for the first time summer d30Si signatures of suspended biogenic silica (i) for the whole water column at three stations and (ii) in the mixed layer at seven stations from the subtropical zone up to the Weddell Gyre. In general, the isotopic composition of biogenic opal exported to depth was comparable to the opal leaving the mixed layer and did not seem to be affected by any diagenetic processes during settling, even if an effect of biogenic silica dissolution cannot be ruled out in the northern part of the Weddell Gyre. We develop a mechanistic understanding of the processes involved in the modern Si-isotopic balance, by implementing a mixed layer model. We observe that the accumulated biogenic silica (sensu Rayleigh distillation) should satisfactorily describe the d30Si composition of biogenic silica exported out of the mixed layer, within the limit of the current analytical precision on the d30Si. The failures of previous models (Rayleigh and steady state) become apparent especially at the end of the productive period in the mixed layer, when biogenic silica production and export are low. This results from (1) a higher biogenic silica dissolution:production ratio imposing a lower net fractionation factor and (2) a higher Si-supply:Si-uptake ratio supplying light Si-isotopes into the mixed layer. The latter effect is especially expressed when the summer mixed layer becomes strongly Si-depleted, together with a large vertical silicic acid gradient, e.g. in the Polar Front Zone and at the Polar Front.

In-situ Iron isotope measurements and mineral compositions from ODP Hole 206-1256D and IODP Hole 335-U1256D

In-situ Fe isotope measurements have been carried out to estimate the impact of the hydrothermal metamorphic overprint on the Fe isotopic composition of Fe-Ti-oxides and Fe-sulfides of the different lithologies of the drilled rocks from IODP Hole 1256D (eastern equatorial Pacific; 15 Ma crust formed at the East Pacific Rise). Most igneous rocks normally have a very restricted range in their 56Fe/54Fe ratio. In contrast, Fe isotope compositions of hot fluids (> 300 °C) from mid-ocean-ridge spreading centers define a narrow range that is shifted to lower delta 56Fe values by 0.2 per mil - 0.5 per mil as compared to igneous rocks. Therefore, it is expected that mineral phases that contain large amounts of Fe are especially affected by the interaction with a fluid that fractionates Fe isotopes during exsolution/precipitation of those minerals. We have used a femtosecond UV-Laser ablation system to determine mineral 56Fe/54Fe ratios of selected samples with a precision of < 0.1 per mil (2 sigma level) at micrometer-scale. We have found significant variations of the delta 56Fe (IRMM-014) values in the minerals between different samples as well as within samples and mineral grains. The overall observed scale of delta 56Fe (magnetite) in 1256D rocks ranges from - 0.12 to + 0.64 per mil, and of delta 56Fe (ilmenite) from - 0.77 to + 0.01 per mil. Pyrite in the lowermost sheeted dike section is clearly distinguishable from the other investigated lithological units, having positive delta 56Fe values between + 0.29 and + 0.56 per mil, whereas pyrite in the other samples has generally negative delta 56Fe values from - 1.10 to - 0.59 permil. One key observation is that the temperature dependent inter-mineral fractionations of Fe isotopes between magnetite and ilmenite are systematically shifted towards higher values when compared to theoretically expected values, while synthesized, well equilibrated magnetite-ilmenite pairs are compatible with the theoretical predictions. Theoretical considerations including beta-factors of different aqueous Fe-chlorides and Rayleigh-type fractionations in the presence of a hydrous, chlorine-bearing fluid can explain this observation. The disagreement between observed and theoretical equilibrium fractionation, the fact that magnetite, in contrast to ilmenite shows a slight downhole trend in the delta 56Fe values, and the observation of small scale heterogeneities within single mineral grains imply that a general re-equilibration of the magnetite-ilmenite pairs is overprinted by kinetic fractionation effects, caused by the interaction of magnetite/ilmenite with hydrothermal fluids penetrating the upper oceanic crust during cooling, or incomplete re-equilibration at low temperatures. Furthermore, the observation of significant small-scale variations in the 56Fe/54Fe ratios of single minerals in this study highlights the importance of high spatial-resolution-analyses of stable isotope ratios for further investigations.

Geochemistry and soil characteristics of sediment core GeoB4234

Stable isotopes of sedimentary nitrogen and organic carbon are widely used as proxy variables for biogeochemical parameters and processes in the water column. In order to investigate alterations of the primary isotopic signal by sedimentary diagenetic processes, we determined concentrations and isotopic compositions of inorganic nitrogen (IN), organic nitrogen (ON), total nitrogen (TN), and total organic carbon (TOC) on one short core recovered from sediments of the eastern subtropical Atlantic, between the Canary Islands and the Moroccan coast. Changes with depth in concentration and isotopic composition of the different fractions were related to early diagenetic conditions indicated by pore water concentrations of oxygen, nitrate, and ammonium. Additionally, the nature of the organic matter was investigated by Rock-Eval pyrolysis and microscopic analysis. A decrease in ON during aerobic organic matter degradation is accompanied by an increase of the 15N/14N ratio. Changes in the isotopic composition of ON can be described by Rayleigh fractionation kinetics which are probably related to microbial metabolism. The influence of IN depleted in 15N on the bulk sedimentary (TN) isotope signal increases due to organic matter degradation, compensating partly the isotopic changes in ON. In anoxic sediments, fixation of ammonium between clay lattices results in a decrease of stable nitrogen isotope ratio of IN and TN. Changes in the carbon isotopic composition of TOC have to be explained by Rayleigh fractionation in combination with different remineralization kinetics of organic compounds with different isotopic composition. We have found no evidence for preferential preservation of terrestrial organic carbon. Instead, both TOC and refractory organic carbon are dominated by marine organic matter. Refractory organic carbon is depleted in 13C compared to TOC.

(Table 2) Annual mass fluxes, 230Th composition and organic geochemistry of sediment trap samples in the South Altantic

The scavenging of 231Pa and 230Th was investigated in the Atlantic Sector of the Southern Ocean by combining results from sediment trap and in situ filtration studies. We present the first high-resolution profile of dissolved 230Th and 231Pa in surface waters across the ACC, showing a dramatic southward increase of both radionuclides around the southern ACC Front at 533S. High dissolved 231Pa/230Th ratios combined with low 230Th/231Pa fractionation factors (F) in these surface waters result in extremely high 231Pa94/230Th94 ratios of material collected in the shallow traps. Particulate 231Pa94/230Th94 ratios in a shallow trap near Bouvet Island increase continuously during the productive period in austral summer, and drop back in the low flux period. This behavior, following the Rayleigh fractionation principle, is interpreted to be due to an increase in the dissolved 231Pa/230Th ratio in the euphotic zone resulting from preferential scavenging of 230Th relative to 231Pa, even in opal-dominated regions. In the post-bloom stage, the depleted radionuclide concentrations are replenished by upwelling of Circumpolar Deep Water. The high particulate 231Pa94/230Th94 signal is weakened during downward transport of the bloom particles in the water column by incorporation of deep suspended particles, which have a lower 231Pa94/230Th94 ratio. It is shown that under the special hydrographic conditions in the Southern Ocean scavenging from the upper water column significantly influences the budgets of 230Th and 231Pa in the sediment. Nevertheless, the budgets are still made up primarily by scavenging from the large standing stock of deep suspended particles.

A laboratory-scale experiment and a numerical simulation of unusual spiral plumes in a High-Prandtl-number fluid, links supplementary material

We experimentally and numerically investigated the generation of plumes from a local heat source (LHS) and studied the interaction of these plumes with cellular convective motion (CCM) in a rectangular cavity filled with silicon oil at a Prandtl number (Pr) of approximately two thousand. The LHS is generated using a 0.2-W green laser beam. A roll-type CCM is generated by vertically heating one side of the cavity. The CCM may lead to the formation of an unusual spiral convective plume that resembles a vertical Archimedes spiral. A similar plume is obtained in a direct numerical simulation. We discuss the physical mechanism for the formation of a spiral plume and the application of the results to mantle convection problems. We also estimate the Reynolds (Re) and Rayleigh (Ra) numbers and apply self-similarity theory to convection in the Earth's mantle. Spiral plumes can be used to interpret mantle tomography results over the last decade.

Microsedimentological characterization as indicators of sedimentary processes and climate changes during Lateglacial at Laguna Potrok Aike, Santa Cruz, Argentina

Palaeoclimatic and paleoenvironmental high latitude records in the Southern Hemisphere are scarce compared to the northern counterpart. However, understanding global evolution of environmental systems during sudden climate changes is inseparable from an equivalent knowledge of both Hemispheres. In this context, a high-resolution study of lacustrine sediments from Laguna Potrok Aike, Santa Cruz province, Patagonia, Argentina, was conducted for the Lateglacial period using concurrent X-Ray Fluorescence (XRF) and Scanning electron microscope analyses. Peaks of Ca/Si and Mn, and occurrences of the green alga Phacotus lenticularis have been interpreted as variations in ventilation of the water column from 13.6 to 11.1 ka cal. BP. During this interval, mild climate conditions during the Younger Dryas are characterized by relatively weak westerlies favouring the formation of a stratified water body as indicated by preserved manganese and Ca/Si peaks and high Total Organic Carbon (TOC) values. In this environment, water in the epilimnion can reach sufficiently high temperature to allow P. lenticularis to grow. Colder conditions are marked by peaks in Ca without P. lenticularis and occur during the Antarctic Cold Reversal (ACR). In this Lateglacial interval, micropumices were also detected in large amount. Image analysis of thin sections allowed the counting and size measurement of detrital particles and micropumices separately. Micropumices significantly influence the iron and titanium content, hence preventing to use them as proxies of detrital input in this interval.

Chemical and mineral compositions of basalts and volcanic glasses from DSDP Sites 65-483 and 65-485

Major and rare earth element (REE) data for basalts from Holes 483, 483B, and 485A of DSDP Leg 65, East Pacific Rise, mouth of the Gulf of California, support a simple fractional crystallization model for the genesis of rocks from this suite. The petrography and mineral chemistry (presented in detail elsewhere) provide no evidence for magma mixing, but rather a simple multistage cooling process. Based on its lowest TiO2 content (0.88%), FeO*/MgO ratio (0.95 with total Fe as FeO), and Mg# (100 Mg/Mg + Fe" = 70), sample 483-17-2-(78-83) has been selected as the most primitive primary magma of the samples analyzed. This is supported by the REE data which show this sample has the lowest total REE content, a La/Sm_cn (chondrite-normalized) = 0.36, and Eu/Sm_cn = 1.05. Because other samples analyzed have higher SiO2, lower Mg#, and a negative Eu anomaly (Eu/Sm_cn as low as 0.89), they are most likely derivative magmas. Wright-Doherty and trace element modelling support fractional crystallization of 14.1% plagioclase (An88), 6.7% olivine (Fo86), and 4.7% clinopyroxene (Wo41En49Fs10) from 483-17-2-(78-83) to form the least differentiated sample with Mg# = 63. The La/Sm_cn of this derivative magma is almost identical to the parent magma (0.35 to 0.36), but the other samples have higher La/Sm_cn (0.45 to 0.51), more total REE, and lower Mg# (60 to 56). Both Wright-Doherty and trace element modelling indicate that the primary magma chosen cannot produce these more evolved samples. For the major elements, the TiO2 and P2O5 are too low in the calculated versus the observed (1.38 to 1.90; 0.11 to 0.17, respectively, for example). Rayleigh fractionation calculates a lower La/Sm_cn and requires about 60% crystal removal versus 40% for the Wright-Doherty. These more evolved samples must be derived from a parent magma different from the one selected here and, unfortunately, not sampled in this study. A magma formed by a smaller degree of partial melting with slightly more residual clinopyroxene left in the mantle than for sample 483-17-2-(78-83) is required.

Seawater carbonate chemistry and biological processes of foraminifera, Globigerinoides ruber and Globigerinella siphonifera during experiments, 2011

Specimens of two species of planktic foraminifera, Globigerinoides ruber and Globigerinella siphonifera, were grown under controlled laboratory conditions at a range of temperatures (18-31 °C), salinities (32-44 psu) and pH levels (7.9-8.4). The shells were examined for their calcium isotope compositions (d44/40Ca) and strontium to calcium ratios (Sr/Ca) using Thermal Ionization Mass Spectrometry and Inductively Coupled Plasma Mass Spectrometry. Although the total variation in d44/40Ca (~0.3 per mill) in the studied species is on the same order as the external reproducibility, the data set reveals some apparent trends that are controlled by more than one environmental parameter. There is a well-defined inverse linear relationship between d44/40Ca and Sr/Ca in all experiments, suggesting similar controls on these proxies in foraminiferal calcite independent of species. Analogous to recent results from inorganically precipitated calcite, we suggest that Ca isotope fractionation and Sr partitioning in planktic foraminifera are mainly controlled by precipitation kinetics. This postulation provides us with a unique tool to calculate precipitation rates and draws support from the observation that Sr/Ca ratios are positively correlated with average growth rates. At 25 °C water temperature, precipitation rates in G. siphonifera and G. ruber are calculated to be on the order of 2000 and 3000 µmol/m**2/h, respectively. The lower d44/40Ca observed at 29 °C in both species is consistent with increased precipitation rates at high water temperatures. Salinity response of d44/40Ca (and Sr/Ca) in G. siphonifera implies that this species has the highest precipitation rates at the salinity of its natural habitat, whereas increasing salinities appear to trigger higher precipitation rates in G. ruber. Isotope effects that cannot be explained by precipitation rate in planktic foraminifera can be explained by a biological control, related to a vacuolar pathway for supply of ions during biomineralization and a pH regulation mechanism in these vacuoles. In case of an additional pathway via cross-membrane transport, supplying light Ca for calcification, the d44/40Ca of the reservoir is constrained as -0.2 per mill relative to seawater. Using a Rayleigh distillation model, we calculate that calcification occurs in a semi-open system, where less than half of the Ca supplied by vacuolization is utilized for calcite precipitation. Our findings are relevant for interpreting paleo-proxy data on d44/40Ca and Sr/Ca in foraminifera as well as understanding their biomineralization processes.

Lithospheric thickness of Pangea

Surface wave tomography, using the fundamental Rayleigh wave velocities and those of higher modes between 1 and 4 and periods between 50 and 160 s, is used to image structures with a horizontal resolution of ~250 km and a vertical resolution of ~50 km to depths of ~300 km in the mantle. A new model, PM_v2_2012, obtained from 3×10**6 seismograms, agrees well with earlier lower resolution models. It is combined with temperature estimates from oceanic plate models and with pressure and temperature estimates from the mineral compositions of garnet peridotite nodules to generate a number of estimates of SV(P,T) based on geophysical and petrological observations alone. These are then used to estimate the unrelaxed shear modulus and its derivatives with respect to pressure and temperature, which agree reasonably with values from laboratory experiments. At high temperatures relaxation occurs, causing the shear wave velocity to depend on frequency. This behaviour is parameterised using a viscosity to obtain a Maxwell relaxation time. The relaxation behaviour is described using a dimensionless frequency, which depends on an activation energy E and volume Va. The values of E and Va obtained from the geophysical models agree with those from laboratory experiments on high temperature creep. The resulting expressions are then used to determine the lithospheric thickness from the shear wave velocity variations. The resolution is improved by about a factor of two with respect to earlier models, and clearly resolves the thick lithosphere beneath active intracontinental belts that are now being shortened. The same expressions allow the three dimensional variations of the shear wave attenuation and viscosity to be estimated.

Sulphate concentrations and sulphur isotopic composition in pore fluids from ODP Leg 168 sites

A numerical model of sulfate reduction and isotopic fractionation has been applied to pore fluid SO4**2- and d34S data from four sites drilled during Ocean Drilling Program (ODP) Leg 168 in the Cascadia Basin at 48°N, where basement temperatures reach up to 62°C. There is a source of sulfate both at the top and the bottom of the sediment column due to the presence of basement fluid flow, which promotes bacterial sulfate reduction below the sulfate minimum zone at elevated temperatures. Pore fluid d34S data show the highest values (135 per mil) yet found in the marine environment. The bacterial sulfur isotopic fractionation factor, a, is severely underestimated if the pore fluids of anoxic marine sediments are assumed to be closed systems and Rayleigh fractionation plots yield erroneous values for a by as much as 15 per mil in diffusive and advective pore fluid regimes. Model results are consistent with a = 1.077+/-0.007 with no temperature effect over the range 1.8 to 62°C and no effect of sulfate reduction rate over the range 2 to 10 pmol/ccm/day. The reason for this large isotopic fractionation is unknown, but one difference with previous studies is the very low sulfate reduction rates recorded, about two orders of magnitude lower than literature values that are in the range of µmol/ccm/day to tens of nmol/ccm/day. In general, the greatest 34S depletions are associated with the lowest sulfate reduction rates and vice versa, and it is possible that such extreme fractionation is a characteristic of open systems with low sulfate reduction rates.