Semiparametric Bayesian measurement error modeling

This work presents a Bayesian semiparametric approach for dealing with regression models where the covariate is measured with error. Given that (1) the error normality assumption is very restrictive, and (2) assuming a specific elliptical distribution for errors (Student-t for example), may be somewhat presumptuous; there is need for more flexible methods, in terms of assuming only symmetry of errors (admitting unknown kurtosis). In this sense, the main advantage of this extended Bayesian approach is the possibility of considering generalizations of the elliptical family of models by using Dirichlet process priors in dependent and independent situations. Conditional posterior distributions are implemented, allowing the use of Markov Chain Monte Carlo (MCMC), to generate the posterior distributions. An interesting result shown is that the Dirichlet process prior is not updated in the case of the dependent elliptical model. Furthermore, an analysis of a real data set is reported to illustrate the usefulness of our approach, in dealing with outliers. Finally, semiparametric proposed models and parametric normal model are compared, graphically with the posterior distribution density of the coefficients. (C) 2009 Elsevier Inc. All rights reserved.

Geochronological constraints on the age of a Permo-Triassic impact event: U-Pb and Ar-40/Ar-39 results for the 40 km Araguainha structure of central Brazil

Impact cratering has been a fundamental geological process in Earth history with major ramifications for the biosphere. The complexity of shocked and melted rocks within impact structures presents difficulties for accurate and precise radiogenic isotope age determination, hampering the assessment of the effects of an individual event in the geological record. We demonstrate the utility of a multi-chronometer approach in our study of samples from the 40 km diameter Araguainha impact structure of central Brazil. Samples of uplifted basement granite display abundant evidence of shock deformation, but U/Pb ages of shocked zircons and the Ar-40/Ar-39 ages of feldspar from the granite largely preserve the igneous crystallization and cooling history. Mixed results are obtained from in situ Ar-40/Ar-39 spot analyses of shocked igneous biotites in the granite, with deformation along kink-bands resulting in highly localized, partial resetting in these grains. Likewise, spot analyses of perlitic glass from pseudotachylitic breccia samples reflect a combination of argon inheritance from wall rock material, the age of the glass itself, and post-impact devitrification. The timing of crater formation is better assessed using samples of impact-generated melt rock where isotopic resetting is associated with textural evidence of melting and in situ crystallization. Granular aggregates of neocrystallized zircon form a cluster of ten U-Pb ages that yield a "Concordia" age of 247.8 +/- 3.8 Ma. The possibility of Pb loss from this population suggests that this is a minimum age for the impact event. The best evidence for the age of the impact comes from the U-Th-Pb dating of neocrystallized monazite and Ar-40/Ar-39 step heating of three separate populations of post-impact, inclusion-rich quartz grains that are derived from the infill of miarolitic cavities. The Pb-206/U-238 age of 254.5 +/- 3.2 Ma (2 sigma error) and Pb-208/Th-232 age of 255.2 +/- 4.8 Ma (2 sigma error) of monazite, together with the inverse, 18 point isochron age of 254 +/- 10 Ma (MSWD = 0.52) for the inclusion-rich quartz grains yield a weighted mean age of 254.7 +/- 2.5 Ma (0.99%, 2 sigma error) for the impact event. The age of the Araguainha crater overlaps with the timing of the Permo-Triassic boundary, within error, but the calculated energy released by the Araguainha impact is insufficient to be a direct cause of the global mass extinction. However, the regional effects of the Araguainha impact event in the Parana-Karoo Basin may have been substantial. (C) 2012 Elsevier Ltd. All rights reserved.

Evidence of direct detection of interstellar deuterium in the local interstellar medium by IBEX

We report the first in situ measurements of neutral deuterium originating in the local interstellar medium (LISM) in Earth’s orbit. These measurements were performed with the IBEX-Lo camera on NASA’s interstellar boundary explorer (IBEX) satellite. All data from the spring observation periods of 2009 through 2011 have been analysed. In the three years of the IBEX mission time, the observation geometry and orbit allowed for a total observation time of 115.3 days for the LISM. However, the effects of the spinning spacecraft and the stepping through 8 energy channels mean that we are only observing the interstellar wind for a total time of 1.44 days, in which 2 counts for interstellar deuterium were collected. We report here a conservative number, because a possibility of systematic error or additional noise, though eliminated in our analysis to the best of our knowledge, only supports detection at a 1-sigma level. From these observations, we derive a ratio D/H = (5.8 ± 4.4) × 10-4 at 1 AU. After modelling the transport and loss of D and H from the termination shock to Earth’s orbit, we find that our result of D/HLISM = (1.6 ± 1.2) × 10-5 agrees with D/HLIC = (1.6 ± 0.4) × 10-5 for the local interstellar cloud. This weak interstellar signal is extracted from a strong terrestrial background signal consisting of sputter products from the sensor’s conversion surface. As reference, we accurately measure the terrestrial D/H ratio in these sputtered products and then discriminate this terrestrial background source. Because of the diminishing D and H signal at Earth’s orbit during the rising solar activity due to photoionisation losses and increased photon pressure, our result demonstrates that in situ measurements of interstellar deuterium in the inner heliosphere are only possible during solar minimum conditions.

The global carbon budget 1959–2011

Accurate assessments of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the climate policy process, and project future climate change. Present-day analysis requires the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. Here we describe datasets and a methodology developed by the global carbon cycle science community to quantify all major components of the global carbon budget, including their uncertainties. We discuss changes compared to previous estimates, consistency within and among components, and methodology and data limitations. CO2 emissions from fossil fuel combustion and cement production (EFF) are based on energy statistics, while emissions from Land-Use Change (ELUC), including deforestation, are based on combined evidence from land cover change data, fire activity in regions undergoing deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. Finally, the global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms. For the last decade available (2002–2011), EFF was 8.3 ± 0.4 PgC yr−1, ELUC 1.0 ± 0.5 PgC yr−1, GATM 4.3 ± 0.1 PgC yr−1, SOCEAN 2.5 ± 0.5 PgC yr−1, and SLAND 2.6 ± 0.8 PgC yr−1. For year 2011 alone, EFF was 9.5 ± 0.5 PgC yr−1, 3.0 percent above 2010, reflecting a continued trend in these emissions; ELUC was 0.9 ± 0.5 PgC yr−1, approximately constant throughout the decade; GATM was 3.6 ± 0.2 PgC yr−1, SOCEAN was 2.7 ± 0.5 PgC yr−1, and SLAND was 4.1 ± 0.9 PgC yr−1. GATM was low in 2011 compared to the 2002–2011 average because of a high uptake by the land probably in response to natural climate variability associated to La Niña conditions in the Pacific Ocean. The global atmospheric CO2 concentration reached 391.31 ± 0.13 ppm at the end of year 2011. We estimate that EFF will have increased by 2.6% (1.9–3.5%) in 2012 based on projections of gross world product and recent changes in the carbon intensity of the economy. All uncertainties are reported as ±1 sigma (68% confidence assuming Gaussian error distributions that the real value lies within the given interval), reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. This paper is intended to provide a baseline to keep track of annual carbon budgets in the future.

The global carbon budget 1959–2011

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the climate policy process, and project future climate change. Present-day analysis requires the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. Here we describe datasets and a methodology developed by the global carbon cycle science community to quantify all major components of the global carbon budget, including their uncertainties. We discuss changes compared to previous estimates, consistency within and among components, and methodology and data limitations. Based on energy statistics, we estimate that the global emissions of CO2 from fossil fuel combustion and cement production were 9.5 ± 0.5 PgC yr−1 in 2011, 3.0 percent above 2010 levels. We project these emissions will increase by 2.6% (1.9–3.5%) in 2012 based on projections of Gross World Product and recent changes in the carbon intensity of the economy. Global net CO2 emissions from Land-Use Change, including deforestation, are more difficult to update annually because of data availability, but combined evidence from land cover change data, fire activity in regions undergoing deforestation and models suggests those net emissions were 0.9 ± 0.5 PgC yr−1 in 2011. The global atmospheric CO2 concentration is measured directly and reached 391.38 ± 0.13 ppm at the end of year 2011, increasing 1.70 ± 0.09 ppm yr−1 or 3.6 ± 0.2 PgC yr−1 in 2011. Estimates from four ocean models suggest that the ocean CO2 sink was 2.6 ± 0.5 PgC yr−1 in 2011, implying a global residual terrestrial CO2 sink of 4.1 ± 0.9 PgC yr−1. All uncertainties are reported as ±1 sigma (68% confidence assuming Gaussian error distributions that the real value lies within the given interval), reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. This paper is intended to provide a baseline to keep track of annual carbon budgets in the future.

Model - data comparison of western African rainfall evolutions during the Last Interglacial (130-115 ka)

This dataset characterizes the evolution of western African precipitation indicated by marine sediment geochemical records in comparison to transient simulations using CCSM3 global climate model throughout the Last Interglacial (130-115 ka). It contains (1) defined tie-points (age models), newly published stable isotopes of benthic foraminifera and Al/Si log-ratios of eight marine sediment cores from the western African margin and (2) annual and seasonal rainfall anomalies (relative to pre-industrial values) for six characteristic latitudinal bands in western Africa simulated by CCSM3 (two transient simulations: one non-accelerated and one accelerated experiment).

Radiocarbon ages and 230Th data of five sediment cores from the Björn and Gardar Drift, North Atlantic

In locations of rapid sediment accumulation receiving substantial amounts of laterally transported material the timescales of transport and accurate quantification of the transported material are at the focus of intense research. Here we present radiocarbon data obtained on co-occurring planktic foraminifera, marine haptophyte biomarkers (alkenones) and total organic carbon (TOC) coupled with excess Thorium-230 (230Thxs) measurements on four sediment cores retrieved in 1649-2879 m water depth from two such high accumulation drift deposits in the Northeast Atlantic, Björn and Gardar Drifts. While 230Thxs inventories imply strong sediment focussing, no age offsets are observed between planktic foraminifera and alkenones, suggesting that redistribution of sediments is rapid and occurs soon after formation of marine organic matter, or that transported material contains negligible amounts of alkenones. An isotopic mass balance calculation based on radiocarbon concentrations of co-occurring sediment components leads us to estimate that transported sediment components contain up to 12% of fossil organic matter that is free of or very poor in alkenones, but nevertheless appears to consist of a mixture of fresh and eroded fossil material. Considering all available constraints to characterize transported material, our results show that although focussing factors calculated from bulk sediment 230Thxs inventories may allow useful approximations of bulk redeposition, they do not provide a unique estimate of the amount of each laterally transported sediment component. Furthermore, our findings provide evidence that the occurrence of lateral sediment redistribution alone does not always hinder the use of multiple proxies but that individual sediment fractions are affected to variable extents by sediment focussing.

Micro-facies and µXRF element scanning data of Dead Sea sediment cores DSEn and 5017-1

Laminated lake sediments from the Dead Sea basin provide high-resolution records of climatic variability in the eastern Mediterranean region, which is especially sensitive to changing climatic conditions. In this study, we aim on detailed reconstruction of climatic fluctuations and related changes in the frequency of flood and dust deposition events at ca. 3300 and especially at 2800 cal. yr BP from high-resolution sediment records of the Dead Sea basin. A ca. 4-m-thick, mostly varved sediment section from the western margin of the Dead Sea (DSEn - Ein Gedi profile) was analysed and correlated to the new International Continental Scientific Drilling Program (ICDP) Dead Sea Deep Drilling Project core 5017-1 from the deep basin. To detect even single event layers, we applied a multi-proxy approach of high-resolution microscopic thin section analyses, micro-X-ray fluorescence (µ-XRF) element scanning and magnetic susceptibility measurements, supported by grain size data and palynological analyses. Based on radiocarbon and varve dating, two pronounced dry periods were detected at ~3500-3300 and ~3000-2400 cal. yr BP which are differently expressed in the sediment records. In the shallow-water core (DSEn), the older dry period is characterised by a thick sand deposit, whereas the sedimentological change at 2800 cal. yr BP is less pronounced and characterised mainly by an enhanced frequency of coarse detrital layers interpreted as erosion events. In the 5017-1 deep-basin core, both dry periods are depicted by halite deposits. The onset of the younger dry period coincides with the Homeric Grand Solar Minimum at ca. 2800 cal. yr BP. Our results suggest that during this period, the Dead Sea region experienced an overall dry climate, superimposed by an increased occurrence of flash floods caused by a change in synoptic weather patterns.

Geochemistry and paleontology of sediment core GeoB3359-3

A geochemical and paleontological reconstruction of paleoproductivity, upwelling intensity and sea surface temperature (SST) off central Chile at 35°S (GeoB3359-3) reveals marked changes from the Last Glacial Maximum (LGM) through the Early Holocene. Surface-water productivity was determined by the interaction between the atmospheric (the Southern Westerlies) and oceanographic (the Antarctic Circumpolar Current, ACC) systems from the LGM through early Termination I (TI). The northward shift of the climate zones during the LGM brought the ACC, as the main macronutrient source, closer to the GeoB3359-3, SST lowered, and surface water productivity and accumulation rates of biogenic components enhanced. With the poleward return of the Southern Westerlies and the ACC, the subtropical high-pressure system became the dominant atmospheric component southward till 35°S during the late TI and Early Holocene and caused surface water productivity to increase through enhanced upwelling.

Geochemistry on the Kekesai composite intrusion

The late Carboniferous to Permian is a critical period for final amalgamation of the Central Asian Orogenic Belt (CAOB), which is characterized by voluminous igneous rocks, particularly granitoids. The Kekesai composite granitoid porphyry intrusion, situated in the Chinese western Tianshan (southwest part of CAOB) includes two intrusive phases, a monzogranite phase, intruded by a granodiorite phase. LA-ICPMS U-Pb zircon analyses suggest that the monzogranitic rocks formed at 305.5±1.1 Ma, with a wide age range of inherited zircons (358-488 Ma and 1208-1391 Ma), whereas the granodioritic rocks formed at 288.7±1.5 Ma. The monzogranitic and granodioritic phases have similar geochemical features and Sr-Nd-Hf isotopic compositions. They exhibit high and variable SiO2 (66-71 wt.%) and MgO (0.41-2.14 wt.%) contents with some arc-like geochemical characteristics (e.g., enrichment of large ion lithophile elements and negative anomalies of Nb, Ta and Ti) and relatively high initial 87Sr/86Sr ratios (ISr=0.7055-0.7059), low positive eNd(t) (+0.84 to +1.03) as well as a large variation in Hf isotopic compositions with eHf(t) between +3.43 to +14.8, implying both of them were derived from similar source materials. These geochemical characteristics suggest that they might be mainly derived from the partial melting of arc-derived Mesoproterozoic mafic lower crust with involvement of a mantle-derived component in variable proportions by mantle-derived magma underplating. The presence of late-Ordovician to earliest early Carboniferous inherited zircons and the Hf isotopic compositions in the monzogranitic sample, similar to that of the widespread juvenile arc rocks, indicates some crust contamination during magma emplacement. Our new data, combined with previous studies, imply that extensive post-collisional magmatism due to underplating of mantle-derived magma, could plausibly be explained by slab break-off regime.

Fission track analyses of ODP Hole 129-800A samples

Four models of fission track annealing in apatite are compared with measured fission track lengths in samples from Site 800 in the East Mariana Basin, Ocean Drilling Program Leg 129, given an independently determined temperature history. The temperature history of Site 800 was calculated using a one-dimensional, compactive, conductive heat flow model assuming two end-member thermal cases: one for cooling of Jurassic ocean crust that has experienced no subsequent heating, and one for cooling of Cretaceous ocean crust. Because the samples analyzed were only shallowly buried and because the tectonic history of the area since sample deposition is simple, resolution of the temperature history is high. The maximum temperature experienced by the sampled bed is between 16°-21°C and occurs at 96 Ma; temperatures since the Cretaceous have dropped in spite of continued pelagic sediment deposition because heat flow has continued to decay exponentially and bottom-water temperatures have dropped. Fission tracks observed within apatite grains from the sampled bed are 14.6 +/- 0.1 µm (1 sigma) long. Given the proposed temperature history of the samples, one unpublished and three published models of fission track annealing predict mean track lengths from 14.8 to 15.9 µm. These models require temperatures as much as 40°C higher than the calculated paleotemperature maximum of the sampled bed to produce the same degree of track annealing. Measured and predicted values are different because annealing models are based on extrapolation of high temperature laboratory data to geologic times. The model that makes the closest prediction is based on the greatest number of experiments performed at low temperature and on an apatite having composition closest to that of the core samples.

Geochemistry of phosphatic crusts from the shelf off Peru

Authigenic phosphorite crusts from the shelf off Peru (9°40°S to 13°30°S) consist of a facies with phosphatic coated grains covered by younger phosphatic laminite. The crusts are composed of carbonate fluorapatite, which probably formed via an amorphous precursor close to the sediment water interface as indicated by low F/P2O5 ratios, Sr and Ca isotopes, as well as rare earth element patterns agreeing with seawater-dominated fluids. Small negative Ce anomalies and U enrichment in the laminite suggest suboxic conditions close to the sediment-water interface during its formation. Increased contents of chalcophilic elements and abundant sulfide minerals in the facies with phosphatic coated grains as well as in the laminite denote sulfate reduction and, consequently, point to episodical development of anoxic conditions during phosphogenesis. The Peruvian phosphorites formed episodically over an extended period of time lasting from Middle Miocene to Pleistocene. Individual phosphatic coated grains show a succession of phosphatic layers with varying contents of organic matter and sulfide-rich phosphatic layers. Coated grains supposedly formed as a result of episodic suspension caused by high turbulence and shifting redox conditions. Episodic anoxia in the pore water induced pyritization in the outermost carbonate fluorapatite layer. Phosphatic coated grains were later transported to the place of crust formation, where subsequent laminite formation was favored under lower energy conditions. A similar succession of phosphatic layers with varying contents of organic matter and sulfide-rich layers in the laminite suggests a formation mechanism analogous to that of individual coated grains.

Iridium contents and microprobe analyses of DSDP Leg 77 samples

Restudy of Deep Sea Drilling Project Sites 536 and 540 in the southeast Gulf of Mexico gives evidence for a giant wave at Cretaceous-Tertiary boundary time. Five units are recognized: (1) Cenomanian limestone underlies a hiatus in which the five highest Cretaceous stages are missing, possibly because of catastrophic K-T erosion. (2) Pebbly mudstone, 45 m thick, represents a submarine landslide possibly of K-T age. (3) Current-bedded sandstone, more than 2.5 m thick, contains anomalous iridium, tektite glass, and shocked quartz; it is interpreted as ejecta from a nearby impact crater, reworked on the deep-sea floor by the resulting tsunami. (4) A 50-cm interval of calcareous mudstone containing small Cretaceous planktic foraminifera and the Ir peak is interpreted as the silt-size fraction of the Cretaceous material suspended by the impact-generated wave. (5) Calcareous mudstone with basal Tertiary forams and the uppermost tail of the Ir anomaly overlies the disturbed interval, dating the impact and wave event as K-T boundary age. Like Beloc in Haiti and Mimbral in Mexico, Sites 536 and 540 are consistent with a large K-T age impact at the nearby Chicxulub crater.

Stable-isotope records from Dronning Maud Land, Antarctica

The European Project for Ice Coring in Antarctica (EPICA) includes a comprehensive pre-site survey on the inland ice plateau of Dronning Maud Land. This paper focuses on the investigation of the 18O content of shallow firn and ice cores. These cores were dated by profiles derived from dielectric-profiling and continuous flow analysis measurements. The individual records were stacked in order to obtain composite chronologies of 18O contents and accumulation rates with enhanced signal-to-noise variance ratios.These chronologies document variations in the last 200 and 1000 years.The 18O contents and accumulation rates decreased in the 19th century and increased during the 20th century.Using the empirical relationships between stable isotopes, accumulation rates and the 10m firn temperature, the variation of both parameters can be explained by the same temperature history.But other causes for these variations, such as the build-up of the snow cover, cannot be excluded. A marked feature of the 1000 year chronology occurs during the period AD 1180-1530 when the 18O contents remains below the long-term mean. Cross-correlation analyses between five cores from the Weddell Sea region and Dronning Maud Land show that 18O records can in some periods be positively correlated and in others negatively correlated, indicating a complex climatic history in time and space.

Nd, Sr and Pb isotopic composition of eolian dust deposited in the central North Pacific

The silicate fractions of recent pelagic sediments in the central north Pacific Ocean are dominated by eolian dust derived from central Asia. An 11 Myr sedimentary record at ODP Sites 885/886 at 44.7°N, 168.3°W allows the evaluation of how such dust and its sources have changed in response to late Cenozoic climate and tectonics. The extracted eolian fraction contains variable amounts (>70%) of clay minerals with subordinate quartz and plagioclase. Uniform Nd isotopic compositions (epsilon-Nd =38.6 to 310.5) and Sm/Nd ratios (0.170-0.192) for most of the 11 Myr record demonstrate a well-mixed provenance in the basins north of the Tibetan Plateau and the Gobi Desert that was a source of dust long before the oldest preserved Asian loess formed. epsilon-Nd values of up to 36.5 for samples 62.9 Ma indicate <=35 wt% admixture of a young, Kamchatka-like volcanic arc component. The coherence of Pb and Nd in the erosional cycle allows us to constrain the Pb isotopic composition of Asian loess devoid of anthropogenic contamination to 206Pb/204Pb =18.97 +/- 0.06, 207Pb/204Pb =15.67 +/- 0.02, 208Pb/204Pb =39.19 +/- 0.11. 87Sr/86Sr (0.711-0.721) and Rb/Sr ratios (0.39-1.1) vary with dust mineralogy and provide an age indication of ~250 Ma. 40Ar/39Ar ages of six dust samples are uniform around 200 Ma and match the K-Ar ages of modern dust deposited on Hawaii. These data reflect the weighted age average of illite formation. Changes from illite- smectite with significant kaolinite to illite- and chlorite-rich, kaolinite-free assemblages since the late Pliocene document changes in the intensity of chemical weathering in the source region. Such weathering evidently did not disturb the K-Ar systematics, and only induced scatter in the Rb-Sr data. We propose that when smectite forms at the expense of illite, K and Ar are quantitatively lost from what becomes smectite, but are quantitatively retained in adjacent illite layers. 40Ar/39Ar age data, therefore, are insensitive to smectite formation during chemical weathering but date the diagenetic growth of illite, the major K-bearing phase in the dust. Over the past 12 Myr, the dust flux to the north Pacific increased by more than an order of magnitude, documenting a substantial drying of central Asia. This climatic change, however, did not alter the ultimate source of the dust, and neoformational products of chemical weathering always remained subordinate to assemblages reworked by mechanical erosion in dust deposited in eastern Asia and the Pacific Ocean.