Ionization of the earth’s atmosphere by solar and galactic cosmic rays

A brief review of the research of atmospheric effects of cosmic rays is presented. Numerical models are discussed, that are capable to compute the cosmic ray induced ionization at a given location and time. Intercomparison of the models, as well as comparison with fragmentary direct measurements of the atmospheric ionization, validates their applicability for the entire atmosphere and the whole range of the solar activity level variations. The effect of sporadic solar energetic particle events is shown to be limited on the global scale, even for the most severe event, but can be very strong locally in polar regions, affecting the physical-chemical properties of the upper atmosphere, especially at high altitudes. Thus, a new methodology is presented to study cosmic ray induced ionization of the atmosphere in full detail using realistic numerical models calibrated to direct observations.     

Late slowdown of the Atlantic Meridional Overturning Circulation during the Last Glacial Inception: New constraints from sedimentary (231Pa/230Th)

Our study gives new constraints on the response of Atlantic Meridional Overturning Circulation (AMOC) export to various forcings during the Last Glacial Inception. The decay corrected excess sedimentary (²³¹Pa/²³⁰Th) activity ratio (hereafter referred to as (Pa/Th)) has been measured over that period in two deep cores from the Western (SU90-11, 44°04′N, 40°01′W, 3645 m) and Eastern (MD01-2446, 39°03′N, 12°37′W, 3547 m) basins of the North Atlantic. Both records display significant changes despite the relatively short half-life of ²³¹Pa (～ 32 kyr) compared to the period we investigate. The (Pa/Th) variability does not correlate to changes in local opal flux normalized to ²³⁰Th. Moreover, the (Pa/Th) profiles display a high degree of coherency with indirect proxies of AMOC activity such as the benthic foraminifera δ¹³C and the mid-latitude summer Sea Surface Temperature in nearby reference cores. These additional pieces of evidence support our interpretation of the (Pa/Th) as reflecting AMOC export. The (Pa/Th) repeatedly underwent rapid changes during the Last Glacial Inception associated with the extension of ice rafted detritus in the North Atlantic, highlighting the control of ice-sheet dynamics through freshwater forcing on AMOC export. AMOC export remains large during periods of ice-sheet growth and its decreases lag the Northern Hemisphere summer insolation forcing. AMOC modulation appears driven by ice-sheet dynamics, itself driven by the seasonal insolation gradient between low and high Northern Hemisphere latitudes and the associated intensity of the meridional oceanic and atmospheric circulation.     

Monitoring forest areas from continental to territorial levels using a sample of medium spatial resolution satellite imagery

A global systematic sampling scheme has been developed by the UN FAO and the EC TREES project to estimate rates of deforestation at global or continental levels at intervals of 5 to 10 years. This global scheme can be intensified to produce results at the national level. In this paper, using surrogate observations, we compare the deforestation estimates derived from these two levels of sampling intensities (one, the global, for the Brazilian Amazon the other, national, for French Guiana) to estimates derived from the official inventories. We also report the precisions that are achieved due to sampling errors and, in the case of French Guiana, compare such precision with the official inventory precision.We extract nine sample data sets from the official wall-to-wall deforestation map derived from satellite interpretations produced for the Brazilian Amazon for the year 2002 to 2003. This global sampling scheme estimate gives 2.81 million ha of deforestation (mean from nine simulated replicates) with a standard error of 0.10 million ha. This compares with the full population estimate from the wall-to-wall interpretations of 2.73 million ha deforested, which is within one standard error of our sampling test estimate. The relative difference between the mean estimate from sampling approach and the full population estimate is 3.1%, and the standard error represents 4.0% of the full population estimate.This global sampling is then intensified to a territorial level with a case study over French Guiana to estimate deforestation between the years 1990 and 2006. For the historical reference period, 1990, Landsat-5 Thematic Mapper data were used. A coverage of SPOT-HRV imagery at 20 m × 20 m resolution acquired at the Cayenne receiving station in French Guiana was used for year 2006.Our estimates from the intensified global sampling scheme over French Guiana are compared with those produced by the national authority to report on deforestation rates under the Kyoto protocol rules for its overseas department. The latter estimates come from a sample of nearly 17,000 plots analyzed from same spatial imagery acquired between year 1990 and year 2006. This sampling scheme is derived from the traditional forest inventory methods carried out by IFN (Inventaire Forestier National). Our intensified global sampling scheme leads to an estimate of 96,650 ha deforested between 1990 and 2006, which is within the 95% confidence interval of the IFN sampling scheme, which gives an estimate of 91,722 ha, representing a relative difference from the IFN of 5.4%.These results demonstrate that the intensification of the global sampling scheme can provide forest area change estimates close to those achieved by official forest inventories (<6%), with precisions of between 4% and 7%, although we only estimate errors from sampling, not from the use of surrogate data.Such methods could be used by developing countries to demonstrate that they are fulfilling requirements for reducing emissions from deforestation in the framework of an REDD (Reducing Emissions from Deforestation in Developing Countries) mechanism under discussion within the United Nations Framework Convention on Climate Change (UNFCCC). Monitoring systems at national levels in tropical countries can also benefit from pan-tropical and regional observations, to ensure consistency between different national monitoring systems.     

Unravelling the morphogenesis of coastal terraces at Cape Laundi (Sumba Island,Indonesia): Insights from numerical models

The morphology of coastal sequences provides fundamental observations to unravel past sea level (SL) variations. For that purpose, converting morphometric observations into a SL datum requires understanding their morphogenesis. The long-lasting sequence of coral reef terraces (CRTs) at Cape Laundi (Sumba Island, Indonesia) could serve as a benchmark. Yet, it epitomizes a pitfall that challenges the ultimate goal: the overall chronology of its development remains poorly constrained. The polycyclic nature of the terraces, involving marine erosion and reoccupation of old coral colonies by more recent ones hinders any clear assignment of Marine Isotope Stages (MIS) to specific terraces, in particular the reference datum corresponding to the last Interglacial maximum (i.e., MIS 5e). Thus, to overcome these obstacles, we numerically model the genesis of the sequence, testing a range of eustatic SL (ESL) reconstructions and uplift rates, as well as exploring the parameter space to address reef growth, erosion and sedimentation. A total of 625 model runs allowed us to improve the morpho-chronological constraints of the coastal sequence and, more particularly, to explain the morphogenesis of the several CRTs associated with MIS 5e. Our results suggest that the lowermost main terrace was first constructed during the marine transgression of MIS 5e and was later reshaped during the marine regression of MIS 5e, as well as during the MIS 5c and MIS 5a highstands. Finally, we discuss the general morphology of the sequence and the implications it may have on SL reconstructions. At Cape Laundi, as elsewhere, we emphasize the necessity of addressing the development of CRT sequences with a dynamic approach, that is, considering that a CRT is a landform built continuously throughout the history of SL oscillations, and not simply during a singular SL maximum.     

Nanoscale resetting of the Th/Pb system in an isotopically-closed monazite grain: A combined atom probe and transmission electron microscopy study

Understanding the mechanisms of parent-daughter isotopic mobility at the nanoscale is key to rigorous interpretation of Ue The Pb data and associated dating. Until now, all nanoscale geochronological studies on geological samples have relied on either Transmission Electron Microscope(TEM) or Atom Probe Microscopy(APM) characterizations alone, thus suffering from the respective weaknesses of each technique. Here we focus on monazite crystals from a ~1 Ga, ultrahigh temperature granulite from Rogaland(Norway). This sample has recorded concordant UeP b dates(measured by LA-ICP-MS) that range over 100 My, with the three domains yielding distinct isotopic Ue Pb ages of 1034 ± 6 Ma(D1; Srich core), 1005 ± 7 Ma(D2), and 935 ± 7 Ma(D3), respectively. Combined APM and TEM characterization of these monazite crystals reveal phase separation that led to the isolation of two different radiogenic Pb(Pb*) reservoirs at the nanoscale. The S-rich core of these monazite crystals contains Cae Srich clusters, 5 -10 nm in size, homogenously distributed within the monazite matrix with a mean interparticle distance of 40 -60 nm. The clusters acted as a sink for radiogenic Pb(Pb*) produced in the monazite matrix, which was reset at the nanoscale via Pb diffusion while the grain remained closed at the micro-scale. Compared to the concordant ages given by conventional micro-scale dating of the grain,the apparent nano-scale age of the monazite matrix in between clusters is about 100 Myr younger, which compares remarkably well to the duration of the metamorphic event. This study highlights the capabilities of combined APM-TEM nano-structural and nano-isotopic characterizations in dating and timing of geological events, allowing the detection of processes untraceable with conventional dating methods.     

Threshold of motion and settling velocities of mollusc shell debris: Influence of faunal composition

Bioclastic particles derived from mollusc shell debris can represent a significant fraction of sandy to gravelly sediments in temperate and cool‐water regions with high carbonate productivity. Their reworking and subsequent transport and deposition by waves and currents is highly dependent on the shape and density of the particles. In this study, the hydrodynamic behaviour of shell debris produced by eight mollusc species is investigated for several grain sizes in terms of settling velocity (measurements in a settling tube) and threshold of motion under unidirectional current (flume experiments using an acoustic profiler). Consistent interspecific differences in settling velocity and critical bed shear stress are found, related to differences in shell density, shell structure imaged by scanning electron microscopy and grain shape. Drag coefficients are proposed for each mollusc species, based on an interpolation of settling velocity data. Depending on the shell species, the critical bed shear stress values obtained for bioclastic particles fall within or slightly below empirical envelopes established for siliciclastic particles, despite very low settling velocity values. The results suggest that settling velocity, often used to describe the entrainment of sediment particles through the equivalent diameter, is not a suitable parameter to predict the initiation of motion of shell debris. The influence of the flat shape of bioclastic particles on the initiation of motion under oscillatory flows and during bedload and saltation transport is yet to be elucidated.     

Automating conceptual models to easily assess trap integrity and oil preservation risks associated with fault reactivation

Extensive Neogene fault reactivation and leakage in the Timor Sea, Australian North-West Shelf, has long been identified as the likely cause for traps in the region being dry or underfilled, despite widespread evidence for hydrocarbon charge and larger palaeo-hydrocarbon columns than those preserved today. A structural model has been previously proposed (Gartrell et al., 2005) to explain the distribution of palaeo- and live hydrocarbon columns. This basic geometrical model is used to define the volume of closure that is protected from fault reactivation, thereby providing the potential for preserving hydrocarbons. Although the model is retrospectively successful in identifying which traps contain hydrocarbons, application to other areas is subjective. To assess the application of the Gartrell et al. (2005) concept a computational model has been generated and applied to a regional 3250 km² seismic interpretation, allowing sensitivity testing of the critical input parameters.The automation has allowed us to demonstrate that the model successfully identifies dry traps and derives reasonable approximations of hydrocarbon-water contacts in commercial accumulations. Therefore, the basic geometrical model is deemed valid to apply to other hydrocarbon provinces, both frontier and mature, that are adversely impacted by fault reactivation after hydrocarbon charge.