Mass‐dependent fractionation of nickel isotopes in meteoritic metal

Abstract— We measured nickel isotopes via multicollector inductively coupled plasma mass spectrometry (MC‐ICPMS) in the bulk metal from 36 meteorites, including chondrites, pallasites, and irons (magmatic and non‐magmatic). The Ni isotopes in these meteorites are mass fractionated; the fractionation spans an overall range of ～0.4‰ amu^?1. The ranges of Ni isotopic compositions (relative to the SRM 986 Ni isotopic standard) in metal from iron meteorites (～0.0 to ～0.3‰ amu^?1) and chondrites (～0.0 to ～0.2‰ amu^?1) are similar, whereas the range in pallasite metal (～–0.1 to 0.0‰ amu^?1) appears distinct. The fractionation of Ni isotopes within a suite of fourteen IIIAB irons (～0.0 to ～0.3‰ amu^?1) spans the entire range measured in all magmatic irons. However, the degree of Ni isotopic fractionation in these samples does not correlate with their Ni content, suggesting that core crystallization did not fractionate Ni isotopes in a systematic way. We also measured the Ni and Fe isotopes in adjacent kamacite and taenite from the Toluca IAB iron meteorite. Nickel isotopes show clearly resolvable fractionation between these two phases; kamacite is heavier relative to taenite by ～0.4‰ amu^?1. In contrast, the Fe isotopes do not show a resolvable fractionation between kamacite and taenite. The observed isotopic compositions of kamacite and taenite can be understood in terms of kinetic fractionation due to diffusion of Ni during cooling of the Fe‐Ni alloy and the development of the Widmanstätten pattern.     

Major changes in ice stream dynamics during deglaciation of the north-western margin of the Laurentide Ice Sheet

Victoria Island lies at the north-western limit of the former North American (Laurentide) Ice Sheet in the Canadian Arctic Archipelago and displays numerous cross-cutting glacial lineations. Previous work suggests that several ice streams operated in this region during the last (Wisconsinan) glaciation and played a major role in ice sheet dynamics and the delivery of icebergs into the Arctic Ocean. This paper produces the first detailed synthesis of their behaviour from the Last Glacial Maximum through to deglaciation (～21–9.5 cal ka BP) based on new mapping and a previously published radiocarbon-constrained ice sheet margin chronology. Over 70 discrete ice flow events (flow-sets) are ‘fitted’ to the ice margin configuration to allow identification of several ice streams ranging in size from large and long-lived (thousands of years) to much smaller and short-lived (hundreds of years). The reconstruction depicts major ice streams in M'Clure Strait and Amundsen Gulf which underwent relatively rapid retreat from the continental shelf edge at some time between ～15.2 and 14.1 cal ka BP: a period which encompasses climatic warming and rapid sea level rise (meltwater pulse-1a). Following this, overall retreat was slower and the ice streams exhibited asynchronous behaviour. The Amundsen Gulf Ice Stream continued to operate during ice margin retreat, whereas the M'Clure Strait Ice Stream ceased operating and was replaced by an ice divide within ～1000 years. This ice divide was subsequently obliterated by another short-lived phase of ice streaming in M'Clintock Channel ～13 cal ka BP. The timing of this large ice discharge event coincides with the onset of the Younger Dryas. Subsequently, a minor ice divide developed once again in M'Clintock Channel, before final deglaciation of the island shortly after 9.5 cal ka BP. It is concluded that large ice streams at the NW margin of the Laurentide Ice Sheet, equivalent in size to the Hudson Strait Ice Stream, underwent major changes during deglaciation, resulting in punctuated delivery of icebergs into the Arctic Ocean. Published radiocarbon dates constrain this punctuated delivery, as far as is possible within the limits imposed by their precision, and we note their coincidence with pulses of meltwater delivery inferred from numerical modelling and ocean sediment cores.     

High-performance direct gravitational N-body simulations on graphics processing units

We present the results of gravitational direct N-body simulations using the commercial graphics processing units (GPU) NVIDIA Quadro FX1400 and GeForce 8800GTX, and compare the results with GRAPE-6Af special purpose hardware. The force evaluation of the N-body problem was implemented in Cg using the GPU directly to speed-up the calculations. The integration of the equations of motions were, running on the host computer, implemented in C using the 4th order predictor–corrector Hermite integrator with block time steps. We find that for a large number of particles (N  10⁴) modern graphics processing units offer an attractive low cost alternative to GRAPE special purpose hardware. A modern GPU continues to give a relatively flat scaling with the number of particles, comparable to that of the GRAPE. The GRAPE is designed to reach double precision, whereas the GPU is intrinsically single-precision. For relatively large time steps, the total energy of the N-body system was conserved better than to one in 10⁶ on the GPU, which is impressive given the single-precision nature of the GPU. For the same time steps, the GRAPE gave somewhat more accurate results, by about an order of magnitude. However, smaller time steps allowed more energy accuracy on the grape, around 10⁻¹¹, whereas for the GPU machine precision saturates around 10⁻⁶ For N  10⁶ the GeForce 8800GTX was about 20 times faster than the host computer. Though still about a factor of a few slower than GRAPE, modern GPUs outperform GRAPE in their low cost, long mean time between failure and the much larger onboard memory; the GRAPE-6Af holds at most 256k particles whereas the GeForce 8800GTX can hold 9 million particles in memory.     

The relationship between dissolved hydrogen and nitrogen fixation in ocean waters

Fixed nitrogen is a key nutrient involved in regulating global marine productivity and hence the global oceanic carbon cycle. Oceanic nitrogen (N₂) fixation is estimated to supply 8×10¹² moles N y^?1 to the ocean, approximately equal to current riverine and the atmospheric inputs of fixed N, and between 50 and 100% of current estimates of oceanic denitrification. However, the spatial and temporal variability of N₂ fixation remains uncertain, mostly because of the normal low resolution sampling for diazotroph distribution and fixation rates. It is well established that N₂ fixation, mediated by the enzyme nitrogenase, is a source of hydrogen (H₂), but the extent to which it leads to supersaturation of H₂ in oceanic waters is unresolved. Here, we present simultaneous measurements of upper ocean dissolved H₂ concentration (nmol L^?1), and rates of N₂ fixation (μmol N m^?3 d^?1), determined using ¹⁵N₂ tracer techniques (at 7 or 15 m), on a transect from Fiji to Hawaii. We find a significant correlation (r=0.98) between dissolved H₂ and rates of N₂ fixation, with the greatest supersaturation of H₂ and highest rates of N₂ fixation being observed in the subtropical gyres at the southern (～18°S) and northern (18°N) reaches of the transect. The lowest H₂ saturation and N₂ fixation were observed in the equatorial region between 8°S and 14°N. We propose that an empirical relationship between H₂ supersaturations and N₂ fixation measurements could be used to guide sampling for ¹⁵N fixation measurements or to aid the spatial interpolation of such measurements.     

On the heat budget of the Arabian Sea

This study analyzes the heat budget of the Arabian Sea using satellite-derived sea-surface temperature (SST) from 1985 to 1995 along with other data sets. For a better understanding of air–sea interaction, canonical average monthly fields representing the spatial and temporal structure of the various components of the heat balance of the Arabian Sea are constructed from up to 30 years of monthly atmospheric and oceanic data. The SST over the Arabian Sea is not uniform and continually evolves with time. Cooling occurs over most of the basin during November through January and May through July, with the greatest cooling in June and July. Warming occurs over most of the basin during the remainder of the year, with the greatest warming occurring in March and September. Results indicate that the sign of the net heat flux is strongly dependent on the location and month. The effects of net heat flux and penetrative solar radiation strongly influence the change in SST during February and are less important during August and September. Horizontal advection acts to cool the sea surface during the northeast monsoon months. During the southwest monsoon horizontal advection of surface waters warms the SST over approximately the southern half of the basin, while the advection of upwelled water from the Somalia and Oman coasts substantially cools the northern basin. The central Arabian Sea during the southwest monsoon is the only area where the change in SST is balanced by the entrainment and turbulent diffusion at the base of the mixed layer. Agreement between the temporal change in the satellite-derived SST and the change calculated from the conservation of heat equation is surprisingly good given the errors in the measured variables and the bulk formula parameters. Throughout the year, monthly results over half of the basin agree within 3°. Considering that the SST changes between 8° and 12° over the year, this means that our results explain from 62% to 75% of the change in SST over 56% of the Arabian Sea. Two major processes contribute to the discrepancy in the change in SST calculated according to the heat budget equation and the change in SST derived from satellite observations. The first is the effect of the horizontal advection term. The position of the major eddies and currents during the southwest monsoon greatly affects the change in SST due to the large gradient in temperature between the cold upwelled waters along the Somali coast to the warm waters in the interior of the basin. The second major process is the thermocline effect. In areas of shallow mixed-layer depth, high insolation and wind speeds of either less than 3 m/s or greater than 15 m/s, the bulk formulae parameterization of the surface heat fluxes is inappropriate.     

Numerical simulation of the circulation within the Perth Submarine Canyon,Western Australia

Surface and sub-surface currents along the ocean boundary of Western Australia were simulated using Regional Ocean Modelling System (ROMS) to examine the circulation within the Perth Canyon. Two major current systems influenced the circulation within the canyon: (1) The Leeuwin current interacted weakly with the canyon as the majority of the canyon was below the depth of the Leeuwin current and (2) Leeuwin undercurrent interacted strongly with the canyon, forming eddies within the canyon at depths of 400–800 m. The results indicated that within the canyon, the current patterns changed continuously although there were some repeated patterns. Recurrent eddies produced regions where upwelling or downwelling dominated during the model runs. Deep upwelling was stronger within the canyon than elsewhere on the shelf, but vertical transport in the upper ocean was strong everywhere when wind forcing was applied. Upwelling alone appeared to be insufficient to transport nutrients to the euphotic zone because the canyon rims were deep. Increased upwelling, combined with entrapment within eddies and strong upwelling-favourable winds, which could assist mixing, may account for the high productivity attributed to the canyon. The Leeuwin current is otherwise a strong barrier to the upwelling of nutrients.     

Meteorites at Meridiani Planum provide evidence for significant amounts of surface and near‐surface water on early Mars

Abstract– Six large iron meteorites have been discovered in the Meridiani Planum region of Mars by the Mars Exploration Rover Opportunity in a nearly 25 km‐long traverse. Herein, we review and synthesize the available data to propose that the discovery and characteristics of the six meteorites could be explained as the result of their impact into a soft and wet surface, sometime during the Noachian or the Hesperian, subsequently to be exposed at the Martian surface through differential erosion. As recorded by its sediments and chemical deposits, Meridiani has been interpreted to have undergone a watery past, including a shallow sea, a playa, an environment of fluctuating ground water, and/or an icy landscape. Meteorites could have been encased upon impact and/or subsequently buried, and kept underground for a long time, shielded from the atmosphere. The meteorites apparently underwent significant chemical weathering due to aqueous alteration, as indicated by cavernous features that suggest differential acidic corrosion removing less resistant material and softer inclusions. During the Amazonian, the almost complete disappearance of surface water and desiccation of the landscape, followed by induration of the sediments and subsequent differential erosion and degradation of Meridiani sediments, including at least 10–80 m of deflation in the last 3–3.5 Gy, would have exposed the buried meteorites. We conclude that the iron meteorites support the hypothesis that Mars once had a denser atmosphere and considerable amounts of water and/or water ice at and/or near the surface.     

Statistical tests of chondrule sorting

Abstract– The variation in sizes of chondrules from one chondrite to the next is thought to be due to some sorting process in the early solar nebula. Hypotheses for the sorting process include chondrule sorting by mass and sorting by some aerodynamic mechanism; one such aerodynamic mechanism is the process of turbulent concentration (TC). We present the results of a series of statistical tests of chondrule data from several different chondrites. The data do not clearly distinguish between various options for the sorting parameter, but we find that the data are inconsistent with being drawn from lognormal or (three‐parameter) Weibull distributions in chondrule radius. We also find that all but one of the chondrule data sets tested are consistent with being drawn from the TC distribution.     

Analysis of Cassini/CIRS limb spectra of Titan acquired during the nominal mission II: Aerosol extinction profiles in the 600–1420 cm spectral range

We have analyzed the continuum emission of limb spectra acquired by the Cassini/CIRS infrared spectrometer in order to derive information on haze extinction in the 3–0.02 mbar range (∼150–350 km). We focused on the 600–1420 cm⁻¹ spectral range and studied nine different limb observations acquired during the Cassini nominal mission at 55°S, 20°S, 5°N, 30°N, 40°N, 45°N, 55°N, 70°N and 80°N. By means of an inversion algorithm solving the radiative transfer equation, we derived the vertical profiles of haze extinction coefficients from 17 spectral ranges of 20-cm⁻¹ wide at each of the nine latitudes. At a given latitude, all extinction vertical profiles retrieved from various spectral intervals between 600 and 1120 cm⁻¹ display similar vertical slopes implying similar spectral characteristics of the material at all altitudes. We calculated a mean vertical extinction profile for each latitude and derived the ratio of the haze scale height (H_haze) to the pressure scale height (H_gas) as a function of altitude. We inferred H_haze/H_gas values varying from 0.8 to 2.4. The aerosol scale height varies with altitude and also with latitude. Overall, the haze extinction does not show strong latitudinal variations but, at 1 mbar, an increase by a factor of 1.5 is observed at the north pole compared to high southern latitudes. The vertical optical depths at 0.5 and 1.7 mbar increase from 55°S to 5°N, remain constant between 5°N and 30°N and display little variation at higher latitudes, except the presence of a slight local maximum at 45°N. The spectral dependence of the haze vertical optical depth is uniform with latitude and displays three main spectral features centered at 630 cm⁻¹, 745 cm⁻¹ and 1390 cm⁻¹, the latter showing a wide tail extending down to ∼1000 cm⁻¹. From 600 to 750 cm⁻¹, the optical depth increases by a factor of 3 in contrast with the absorbance of laboratory tholins, which is generally constant. We derived the mass mixing ratio profiles of haze at the nine latitudes. Below the 0.4-mbar level all mass mixing ratio profiles increase with height. Above this pressure level, the profiles at 40°N, 45°N, 55°N, at the edge of the polar vortex, display a decrease-with-height whereas the other profiles increase. The global increase with height of the haze mass mixing ratio suggest a source at high altitudes and a sink at low altitudes. An enrichment of haze is observed at 0.1 mbar around the equator, which could be due to a more efficient photochemistry because of the strongest insolation there or an accumulation of haze due to a balance between sedimentation and upward vertical drag.