Simulation of the interaction of galactic cosmic‐ray protons with meteoroids: On the production of radionuclides in thick gabbro and iron targets irradiated isotropically with 1.6 GeV protons

Abstract— Thick spherical targets made of gabbro (R = 25 cm) and of steel (R = 10 cm) were irradiated isotropically with 1.6 GeV protons at the Saturne synchrotron at Laboratoire National Saturne (LNS)/CEN Saclay in order to simulate the interaction in space of galactic cosmic‐ray (GCR) protons with stony and iron meteoroids. Proton fluences of 1.32 × 10¹⁴ cm^?2 and 2.45 × 10¹⁴ cm^?2 were received by the gabbro and iron sphere, respectively, which corresponds to cosmic‐ray exposure ages of about 1.6 and 3.0 Ma. Both artificial meteoroids contained large numbers of high‐purity target foils of up to 28 elements at different depths. In these individual target foils, elementary production rates of radionuclides and rare gas isotopes were measured by x‐ and γ‐spectrometry, by low‐level counting, accelerator mass spectrometry (AMS), and by conventional rare gas mass spectrometry. Also samples of the gabbro itself were analyzed. Up to now, for each of the experiments, ～500 target‐product combinations were investigated of which the results for radionuclides are presented here. The experimental production rates show a wide range of depth profiles reflecting the differences between low‐, medium‐, and high‐energy products. The influence of the stony and iron matrices on the production of secondary particles and on particle transport, in general, and consequently on the production rates is clearly exhibited by the phenomenology of the production rates as well as by a detailed theoretical analysis. Theoretical production rates were calculated in an a priori way by folding depth‐dependent spectra of primary and secondary protons and secondary neutrons calculated by Monte Carlo techniques with the excitation functions of the underlying nuclear reactions. Discrepancies of up to a factor of 2 between the experimental and a priori calculated depth profiles are attributed to the poor quality of the mostly theoretical neutron excitation functions. Improved neutron excitation functions were obtained by least‐squares deconvolution techniques from experimental thick‐target production rates of up to five thick‐target experiments in which isotropic irradiations were performed. A posteriori calculations using the adjusted neutron cross sections describe the measured depth profiles of all these simulation experiments within 9%. The thus validated model calculations provide a basis for reliable physical model calculations of the production rates of cosmogenic nuclides in stony and iron meteorites as well as in lunar samples and terrestrial materials.     

Spaceborne visible and thermal infrared lithologic mapping of impact‐exposed subsurface lithologies at the Haughton impact structure,Devon Island,Canadian High Arctic: Applications to Mars

Abstract— This study serves as a proof‐of‐concept for the technique of using visible‐near infrared (VNIR), short‐wavelength infrared (SWIR), and thermal infrared (TIR) spectroscopic observations to map impact‐exposed subsurface lithologies and stratigraphy on Earth or Mars. The topmost layer, three subsurface layers and undisturbed outcrops of the target sequence exposed just 10 km to the northeast of the 23 km diameter Haughton impact structure (Devon Island, Nunavut, Canada) were mapped as distinct spectral units using Landsat 7 ETM+ (VNIR/SWIR) and ASTER (VNIR/SWIR/TIR) multispectral images. Spectral mapping was accomplished by using standard image contrast‐stretching algorithms. Both spectral matching and deconvolution algorithms were applied to image‐derived ASTER TIR emissivity spectra using spectra from a library of laboratory‐measured spectra of minerals (Arizona State University) and whole‐rocks (Ward's). These identifications were made without the use of a priori knowledge from the field (i.e., a “blind” analysis). The results from this analysis suggest a sequence of dolomitic rock (in the crater rim), limestone (wall), gypsum‐rich carbonate (floor), and limestone again (central uplift). These matched compositions agree with the lithologic units and the pre‐impact stratigraphic sequence as mapped during recent field studies of the Haughton impact structure by Osinski et al. (2005a). Further conformation of the identity of image‐derived spectra was confirmed by matching these spectra with laboratory‐measured spectra of samples collected from Haughton. The results from the “blind” remote sensing methods used here suggest that these techniques can also be used to understand subsurface lithologies on Mars, where ground truth knowledge may not be generally available.     

The Puerto Lápice eucrite

Abstract— Puerto Lápice is a new eucrite fall (Castilla‐La Mancha, Spain, 10 May 2007). In this paper, we report its detailed petrography, magnetic characterization, mineral and bulk chemistry, oxygen and noble gas isotope systematics, and radionuclide data. Study of four thin sections from two different specimens reveal that the meteorite is brecciated in nature, and it contains basaltic and granulitic clasts, as well as recrystallized impact melt and breccia fragments. Shock veins are ubiquitous and show evidence of at least three different shock events. Bulk chemical analyses suggest that Puerto Lápice belongs to the main group of basaltic eucrites, although it has a significantly higher Cr content. Oxygen isotopes also confirm that the meteorite is a normal member of the HED suite. Noble gas abundances show typical patterns, with dominant cosmogenic and radiogenic contributions, and indicate an average exposure age of 19 ± 2 Ma, and a Pu‐fission Xe age well within typical eucrite values. Cosmogenic radionuclides suggest a preatmospheric size of about 20–30 cm in diameter.     

Evidence for interchange reconnection between a coronal hole and an adjacent emerging flux region

Coronal holes are regions of dominantly monopolar magnetic field on the Sun where the field is considered to be ‘open’ towards interplanetary space. Magnetic bipoles emerging in proximity to a coronal hole boundary naturally interact with this surrounding open magnetic field. In the case of oppositely aligned polarities between the active region and the coronal hole, we expect interchange reconnection to take place, driven by the coronal expansion of the emerging bipole as well as occasional eruptive events. Using SOHO/EIT and SOHO/MDI data, we present observational evidence of such interchange reconnection by studying AR 10869 which emerged close to a coronal hole. We find closed loops forming between the active region and the coronal hole leading to the retreat of the hole. At the same time, on the far side of the active region, we see dimming of the corona which we interpret as a signature of field line ‘opening’ there, as a consequence of a topological displacement of the ‘open’ field lines of the coronal hole. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Comparison of Five Numerical Codes for Automated Tracing of Coronal Loops

The three-dimensional (3D) modeling of coronal loops and filaments requires algorithms that automatically trace curvilinear features in solar EUV or soft X-ray images. We compare five existing algorithms that have been developed and customized to trace curvilinear features in solar images: i) the oriented-connectivity method (OCM), which is an extension of the Strous pixel-labeling algorithm (developed by Lee, Newman, and Gary); ii) the dynamic aperture-based loop-segmentation method (developed by Lee, Newman, and Gary); iii) unbiased detection of curvilinear structures (developed by Steger, Raghupathy, and Smith); iv) the oriented-direction method (developed by Aschwanden); and v) ridge detection by automated scaling (developed by Inhester). We test the five existing numerical codes with a TRACE image that shows a bipolar active region and contains over 100 discernable loops. We evaluate the performance of the five codes by comparing the cumulative distribution of loop lengths, the median and maximum loop length, the completeness or detection efficiency, the accuracy, and flux sensitivity. These algorithms are useful for the reconstruction of the 3D geometry of coronal loops from stereoscopic observations with the STEREO spacecraft, or for quantitative comparisons of observed EUV loop geometries with (nonlinear force-free) magnetic field extrapolation models.     

On bridges B(pL, qS) around triangular libration points

When μ is smaller than Routh’s critical value μ ₁ = 0.03852 . . . , two planar Lyapunov families around triangular libration points exist, with the names of long and short period families. There are periodic families which we call bridges connecting these two Lyapunov families. With μ increasing from 0 to 1, how these bridges evolve was studied. The interval (0,1) was divided into six subintervals (0, μ ₅), (μ ₅, μ ₄), (μ ₄, μ ₃), (μ ₃, μ ₂), (μ ₂, μ ₁), (μ ₁, 1), and in each subinterval the families B(pL, qS) were studied, along with the families B(qS, qS′). Especially in the interval (μ ₂, μ ₁), the conclusion that the bridges B(qS, qS′) do not exist was obtained. Connections between the short period family and the bridges B(kS, (k + 1)S) were also studied. With these studies, the structure of the web of periodic families around triangular libration points was enriched.     

Automated probabilistic classification of transients and variables

There is an increasing number of large, digital, synoptic sky surveys, in which repeated observations are obtained over large areas of the sky in multiple epochs. Likewise, there is a growth in the number of (often automated or robotic) follow‐up facilities with varied capabilities in terms of instruments, depth, cadence, wavelengths, etc., most of which are geared toward some specific astrophysical phenomenon. As the number of detected transient events grows, an automated, probabilistic classification of the detected variables and transients becomes increasingly important, so that an optimal use can be made of follow‐up facilities, without unnecessary duplication of effort. We describe a methodology now under development for a prototype event classification system; it involves Bayesian and Machine Learning classifiers, automated incorporation of feedback from follow‐up observations, and discriminated or directed follow‐up requests. This type of methodology may be essential for the massive synoptic sky surveys in the future. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A bulk cloud parameterization in a Venus General Circulation Model

A condensing cloud parameterization is included in a super-rotating Venus General Circulation Model. A parameterization including condensation, evaporation and sedimentation of mono-modal sulfuric acid cloud particles is described. Saturation vapor pressure of sulfuric acid vapor is used to determine cloud formation through instantaneous condensation and destruction through evaporation, while pressure dependent viscosity of a carbon dioxide atmosphere is used to determine sedimentation rates assuming particles fall at their terminal Stokes velocity. Modifications are described to account for the large range of the Reynolds number seen in the Venus atmosphere.Two GCM experiments initialized with 10 ppm-equivalent of sulfuric acid are integrated for 30 Earth years and the results are discussed with reference to “Y” shaped cloud structures observed on Venus. The GCM is able to produce an analog of the “Y” shaped cloud structure through dynamical processes alone, with contributions from the mean westward wind, the equatorial Kelvin wave, and the mid-latitude/polar Mixed Rossby/Gravity waves. The cloud top height in the GCM decreases from equator to pole and latitudinal gradients of cloud top height are comparable to those observed by Pioneer Venus and Venus Express, and those produced in more complex microphysical models of the sulfur cycle on Venus. Differences between the modeled cloud structures and observations are described and dynamical explanations are suggested for the most prominent differences.     

The roughness of the dark side of Iapetus from the 2004 to 2005 flyby

The roughness of a planetary surface offers clues to its past geologic history. We apply a surface roughness model developed by Buratti and Veverka (Buratti, B.J., Veverka, J. [1985]. Icarus 64, 320-328) to Cassini ISS data from the January 1st, 2005 flyby of Iapetus. This model uses the observed scattering behavior to provide a depth to radius factor q quantifying the size of idealized craters on the surface. Our findings indicate that the surface on the dark side is significantly smoother than the surfaces of other icy low-albedo saturnian satellites. We have found that the average depth to radius on the leading (dark) side is 0.084, corresponding to a Hapke mean slope angle of 6°. As compared to the 13-33° Hapke mean slope angle of other icy satellites (Buratti, B.J., and 10 colleagues [2008]. Icarus 193, 309-322), our results present a clearly different picture for the leading surface of Iapetus, suggesting that the dark deposit contributes to the decrease in macroscopic surface roughness of the leading side. Attempts were made to obtain an average depth to radius value for the trailing (bright) side; however the scans of the bright side from this flyby exhibited large variations in albedo, resulting in results that were physically unrealistic.