Monte Carlo simulation of GCR neutron capture production of cosmogenic nuclides in stony meteorites and lunar surface

Abstract— A purely physical model based on a Monte Carlo simulation of galactic cosmic ray (GCR) particle interaction with meteoroids is used to investigate neutron interactions down to thermal energies. Experimental and/or evaluated excitation functions are used to calculate neutron capture production rates as a function of the size of the meteoroid and the depth below its surface. Presented are the depth profiles of cosmogenic radionuclides ³⁶Cl, ⁴¹Ca, ⁶⁰Co, ⁵⁹Ni, and ¹²⁹I for meteoroid radii from 10 cm up to 500 cm and a 2π irradiation. Effects of bulk chemical composition on n‐capture processes are studied and discussed for various chondritic and lunar compositions. The mean GCR particle flux over the last 300 ka was determined from the comparison of simulations with measured ⁴¹Ca activities in the Apollo 15 drill core. The determined value significantly differs from that obtained using equivalent models of spallation residue production.     

Synthetic SH seismograms in a laterally varying medium by the discrete wavenumber method

Summary. We present a new method to calculate the SH wavefield produced by a seismic source in a half-space with an irregular buried interface. The diffracting interface is represented by a distribution of body forces. The Green's functions needed to solve the boundary conditions are evaluated using the discrete wavenumber method. Our approach relies on the introduction of a periodicity in the source-medium configuration and on the discretization of the interface at regular spacing. The technique developed is applicable to boundaries of arbitrary shapes and is valid at all frequencies. Some examples of calculation in simple configurations are presented showing the capabilities of the method.     

A field study of the coupled effects of aquifer stratification,fluid density,and groundwater fluctuations on dispersivity assessments

A number of experimental studies have tackled the issue of solute transport parameter assessments either in the laboratory or in the field. But yet, the behavior of a plume in the field under density driven forces, is not well known due to possible development of instabilities. Some field tracer tests on the fate of plumes denser than native groundwater such as those encountered under waste disposal facilities, have pointed out the processes of sinking and splitting at the early stage of migration. The process of dispersion was widely investigated, but the range of dispersivity values obtained from either experimental tests, or numerical and theoretical calculations is still very large, even for the same type of aquifers. These discrepancies were considered to be essentially caused by soil heterogeneities and scale effects. In the meantime, studies on the influence of sinking and fingering have remained more scarce. The objective of the work is to analyze how transport parameters such as dispersivities can be affected by unstable conditions, which lead to plume sinking and fingering. A series of tracer tests were carried out to study under natural conditions, the transport of a dense chloride solution injected in a shallow two-layered aquifer. Two types of experiments were performed: in the first type, source injection was such that the plume could travel downward from one layer to the other of higher pore velocity, and in the second one, the migration took place only in the faster layer. The results suggest some new insights in the processes occurring at the early stages of a dense plume migration moving in a stratified aquifer under groundwater fluctuations, which can be summarized through the following points: (i) Above a stability criterion threshold, a fingering process and a multi modal plume transport take place, but local dispersivities can be cautiously derived, using breakthrough curves matching. (ii) When water table is subject to some cycling or rising, the plume can be significantly distorted in the transverse direction, leading to unusual values of the ratio between longitudinal and transverse dispersivities. (iii) Under stable conditions, for example in the case of straightforward injection in the faster aquifer layer, longitudinal dispersivity is greater than the transverse component as usually encountered, and the obtained transport parameters are closed to macro dispersivity values, which reach their asymptotic limit at very short distances. (iv) The classical scale effect about the varying dispersivity at short distances could be a process mainly due to the distance required for a plume stabilization.     

Reefs at the edge: coral community structure around Rapa,southernmost French Polynesia

Rapa (27°36′ S, 144°20′ W) is a small (~40 km²) volcanic island isolated in the Southern Austral Archipelago, where direct anthropogenic stressors are extremely limited. Here, we present the results of the first quantitative survey of coral community structure across habitats and depths around the island. Despite its geographical isolation in the depauperate South Central Pacific, its small size and unfavourable environmental conditions (competition with macroalgae, low sea surface temperatures, reduced reef accretion), the diversity of scleractinian corals at Rapa is particularly high (112 species from 32 genera, including 37 species of Acropora) in comparison to other French Polynesian islands and subtropical Pacific locations. Our results indicate that the abundance (>100 colonies per 10 m² recorded at nine of the 17 sampling stations) and cover (>40% at four stations) of corals are relatively high for a marginal reef location. Strong spatial heterogeneity was found, with high variation in diversity, abundance, cover and community composition among stations. Variation in community composition was related to habitat types, with distinct assemblages among fringing reefs within bays, reef formations at bay entrances, and those on the submerged platform surrounding the island. On the platform, a depth gradient was detected, with generic richness, abundance and cover generally greater at deeper stations (18–20 m depth) compared with medium‐depth (10–12 m) and shallow (1–3 m) stations. A gradient was also recorded along bays, with increasing coral diversity and abundance from the bay heads to the bay entrances. The coral community at Rapa was characterized by the presence of several taxa not found in other French Polynesian archipelagos and the rarity of others that are common and abundant in the Society and the Tuamotu islands. Another distinctive feature of reef communities at Rapa is the high cover and dominance of macroalgae, particularly in the shallower parts of the surrounding platform, which probably explains the lower densities of coral colonies recorded there. These characteristics of the diversity and biogeographical composition of coral assemblages at Rapa provide considerable ecological grounds for its conservation.     

The European coastal zone: characterization and first assessment of ecosystem metabolism

The geomorphic, oceanographic, terrestrial and anthropogenic attributes of the European coastal zone are described and published data on ecosystem function (primary production and respiration) are reviewed. Four regions are considered: the Baltic Sea, Mediterranean Sea, Black Sea and the European Atlantic coast including the North Sea. The metabolic database (194 papers) suffers from a non-homogeneous geographical coverage with no usable data for the Black Sea which was therefore excluded from this part of our study. Pelagic gross primary production in European open shelves is, by far, the most documented parameter with an estimated mean of 41 mmol C m⁻² d⁻¹, the lowest value is reported in the Mediterranean Sea (21 mmol C m⁻² d⁻¹) and the highest one in the Atlantic/North Sea area (51 mmol C m⁻² d⁻¹). Microphytobenthic primary production, mostly measured in shallow areas, is extrapolated to the entire 0–200 m depth range. Its contribution to total primary production is low in all regions (mean: 1.5 mmol C m⁻² d⁻¹). Although macrophyte beds are very productive, a regional production estimate is not provided in this study because their geographical distribution along the European coastline remains unknown. Measurements of pelagic community respiration are clearly too sparse, especially below the euphotic zone, to yield an accurate picture of the fate of organic matter produced in the water column. With a mean value of 17 mmol C m⁻² d⁻¹, benthic community respiration consumes approximately 40% of the pelagic organic matter production. Estuaries generally exhibit high metabolic rates and a large range of variation in all parameters, except microphytobenthic primary production. Finally, the problem of eutrophication in Europe is discussed and the metabolic data obtained in the framework of the Land–Ocean Interactions in the Coastal Zone (LOICZ) project are compared with available direct measurements of net ecosystem production.     

Comparaison de catalogues de crateres lunaires

Résumé La partie systématique des différences entre les positions des cratères de plusieurs catalogues est développée en harmoniques sphériques. L'application du test en ² détermine l'ordre du développement. La méthode a été appliquée pour comparer le système de référence de Arthur avec les catalogues de Schrutka-Rechtenstamm et Gavrilov. On obtient les surfacesf (l, b) des différences systématiques de ces catalogues. La zone centrale de la face visible de la lune (± 40° en latitude et longitude) est bien déterminée, mais la zone marginale présente des différences systématiques qui peuvent atteindre le kilomètre.

---

The systematic part of the differences between positions in crater catalogues is expanded into spherical harmonics. The application of the ² test determines the highest order of the expansion. The method was applied to compare the Arthur System with the Schrutka-Rechtenstamm and Gavrilov catalogues. The surfacesf (l, b) of systematic differences from these catalogues were obtained. The central zone of the near side of the Moon (± 40° in latitude and longitude) is well determined but the marginal zone presents systematic differences that can be as large as one kilometer.

     

Numerical simulation of the propagation of P waves in fractured media

We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) ⁻¹ in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.     

Graphical tools based on Turc-Budyko plots to detect changes in catchment behaviour

Abstract

This paper aims at revisiting the use of non-dimensional representations for catchment change and model robustness analysis. As well as being helpful for catchment classification according to hydroclimatic conditions, the Turc-Budyko representation enables visualization of temporal changes in these conditions and/or modification of the hydrological behaviour of the catchment. This brings a new perspective to hydrological analysis, different from the classical time series plots, which helps when interpreting catchment functioning with respect to hydroclimatic constraints. These tools do not require statistical analyses of observed series and are therefore very simple to implement. Four case studies are considered here to illustrate the use of the Q/P = f(P/E_p) graph. When combined with the inter-annual Turc-Mezentsev formula, this visual framework enables anticipation of potential difficulties for models based on the same hypotheses, solely using the analysis of observed data.     

A critical zone observatory dedicated to suspended sediment transport: The meso-scale Galabre catchment (southern French Alps)

The 20 km² Galabre catchment belongs to the French network of critical zone observatories (OZCAR; Gaillardet et al., Vadose Zone Journal, 2018, 17(1), 1–24). It is representative of the sedimentary lithology and meteorological forcing found in Mediterranean and mountainous areas. Due to the presence of highly erodible and sloping badlands on various lithologies, the site was instrumented in 2007 to understand the dynamics of suspended sediments (SS) in such areas. Two meteorological stations including measurements of air temperature, wind speed and direction, air moisture, rainfall intensity, raindrop size and velocity distribution were installed both in the upper and lower part of the catchment. At the catchment outlet, a gauging station records the water level, temperature and turbidity (10 min time-step). Stream water samples are collected automatically to estimate SS concentration-turbidity relationships, allowing quantification of SS fluxes with known uncertainty. The sediment samples are further characterized by measuring their particle size distributions and by applying a low-cost sediment fingerprinting approach using spectrocolorimetric tracers. Thus, the contributions of badlands located on different lithologies to total SS flux are quantified at a high temporal resolution, providing the opportunity to better analyse the links between meteorological forcing variability and watershed hydrosedimentary response. The set of measurements was extended to the dissolved phase in 2017. Both stream water electrical conductivity and major ion concentrations are measured each week and every 3 h during storm events. This extension of measurements to the dissolved phase will allow progress in understanding both the origin of the water during the events and the partitioning between particulate and dissolved fluxes of solutes in the critical zone. All data sets are available at https://doi.osug.fr/public/DRAIXBLEONE_GAL/index.html .