Improving runoff prediction using agronomical information in a cropped,loess covered catchment

Predicting runoff hot spots and hot‐moments within a headwater crop‐catchment is of the utmost importance to reduce adverse effects on aquatic ecosystems by adapting land use management to control runoff. Reliable predictions of runoff patterns during a crop growing season remain challenging. This is mainly due to the large spatial and temporal variations of topsoil hydraulic properties controlled by complex interactions between weather, growing vegetation, and cropping operations. This interaction can significantly modify runoff patterns and few process‐based models can integrate this evolution of topsoil properties during a crop growing season at the catchment scale. Therefore, the purpose of this study was to better constrain the event‐based hydrological model Limburg Soil Erosion Model by incorporating temporal constraints for input topsoil properties during a crop growing season (LISEM). The results of the temporal constraint strategy (TCS) were compared with a classical event per event calibration strategy (EES) using multi‐scale runoff information (from plot to catchment). The EES and TCS approaches were applied in a loess catchment of 47 ha located 30 km northeast of Strasbourg (Alsace, France). A slight decrease of the Nash–Sutcliffe efficiency criterion on runoff discharge for TCS compared to EES was counterbalanced by a clear improvement of the spatial runoff patterns within the catchment. This study showed that limited agronomical and climatic information added during the calibration step improved the spatial runoff predictions of an event‐based model. Reliable prediction of runoff source, connectivity, and dynamics can then be derived and discussed with stakeholders to identify runoff hot spots and hot‐moments for subsequent land use and crop management modifications.     

Bosumtwi impact structure,Ghana: Evidence for fluidized emplacement of the ejecta

The about 10.5 km diameter Bosumtwi impact crater is one of the youngest large impact structures on Earth. The crater rim is readily noticed on topographic maps or in satellite imagery. It defines a circular basin filled by water (Lake Bosumtwi) and lacustrine sediments. The morphology of this impact structure is also characterized by a circular plateau extending beyond the rim and up to 9–10 km from the center of the crater (about 2 crater radii). This feature comprises a shallow ring depression, also described as an annular moat, and a subdued circular ridge at its outer edge. The origin of this outermost feature could so far not be elucidated based on remote sensing data only. Our approach combines detailed topographic analysis, including roughness mapping, with airborne radiometric surveys (mapping near‐surface K, Th, U concentrations) and field observations. This provides evidence that the moat and outer ring are features inherited from the impact event and represent the partially eroded ejecta layer of the Bosumtwi impact structure. The characteristics of the outer ridge indicate that ejecta emplacement was not purely ballistic but requires ejecta fluidization and surface flow. The setting of Bosumtwi ejecta can therefore be considered as a terrestrial analog for rampart craters, which are common on Mars and Venus, and also found on icy bodies of the outer solar system (e.g., Ganymede, Europa, Dione, Tethys, and Charon). Future studies at Bosumtwi may therefore help to elucidate the mechanism of formation of rampart craters.     

From discrete to continuum modelling of boundary value problems in geomechanics: An integrated FEM-DEM approach

Double-scale numerical methods constitute an effective tool for simultaneously representing the complex nature of geomaterials and treating real-scale engineering problems such as a tunnel excavation or a pressuremetre at a reasonable numerical cost. This paper presents an approach coupling discrete elements (DEM) at the microscale with finite elements (FEM) at the macroscale. In this approach, a DEM-based numerical constitutive law is embedded into a standard FEM formulation. In this regard, an exhaustive discussion is presented on how a 2D/3D granular assembly can be used to generate, step by step along the overall computation process, a consistent Numerically Homogenised Law. The paper also focuses on some recent developments including a comprehensive discussion of the efficiency of Newton-like operators, the introduction of a regularisation technique at the macroscale by means of a second gradient framework, and the development of parallelisation techniques to alleviate the computational cost of the proposed approach. Some real-scale problems taking into account the material spatial variability are illustrated, proving the numerical efficiency of the proposed approach and the benefit of a particle-based strategy.     

Northwest Africa 10414, a pigeonite cumulate shergottite

Northwest Africa (NWA) 10414 is an unusual shergottite with a cumulate texture. It contains 73% coarse prismatic pigeonite, plus 18% interstitial maskelynite, 2% Si‐rich mesostasis, 2% merrillite, and minor chromite‐ulvöspinel. It contains no olivine, and only ~3% augite. Phase compositions are pigeonite (En_68‐43Fs_27‐48Wo_5‐15) and maskelynite An_~54‐36, more sodic than most maskelynite in shergottites. Chromite‐ulvöspinel composition plots between the earliest and most fractionated spinel‐group minerals in olivine‐phyric shergottites. NWA 10414 mineralogically resembles the contact facies between Elephant Moraine 79001 lithologic units A and B, with abundant pigeonite phenocrysts, though it is coarser grained. Its most Mg‐rich pigeonite also has a similar composition to the earliest crystallized pyroxenes in several other shergottites, including Shergotty. The Shergotty intercumulus liquid composition crystallizes pigeonite with a similar composition range to NWA 10414 pigeonite, using PETROLOG. Olivine‐phyric shergottite NWA 6234, with a pure magma composition, produces an even better match to this pigeonite composition range, after olivine crystallization. These observations suggest that after the accumulation of olivine from an olivine‐phyric shergottite magma, the daughter liquid could precipitate pigeonite locally to form this pigeonite cumulate, before the crystallization of overlying liquid as a normal basaltic shergottite.     

In situ LA-ICP-MS U–Pb titanite dating of Na–Ca metasomatism in orogenic belts: the North Pyrenean example

In the Pyrenees, in association with the rotation of the Iberian plate around Europe during the Mid-Cretaceous, a Na–Ca metasomatism is recognized as a complementary record of the hydrothermal activity that led to Na-metasomatism (albitization) and talc–chlorite mineralization. It affected metasedimentary rocks as well as Hercynian granites. In situ laser ablation ICP-MS U–Pb analyses of titanite grains formed in albitites during metasomatism date the Na–Ca metasomatism between 110 and 92 Ma. The temperature of the Na–Ca metasomatism is estimated to be approximately 550 °C. Both the time constraints and temperature estimates suggest that the Na–Ca metasomatism is related to the low-P high-T North Pyrenean metamorphism.     

An elastic–viscous–plastic sea ice model formulated on Arakawa B and C grids

The elastic–viscous–plastic (EVP) sea ice rheology has been introduced in the large-scale Louvain-la-Neuve sea-Ice Model, version 2 (LIM2), and its performance has been evaluated. Centred difference versions of the rheology have been implemented on both an Arakawa B grid and a C grid, and their performance have been intercompared in coupled simulations with the Nucleus for European Modelling of the Ocean (NEMO) model. Integrations with both implementations lead to fairly similar results which compare well with observations and with previous LIM simulations. The C grid version, however, offers a number of advantages: (a) easier ice coupling with NEMO, which is itself defined on a C grid; (b) possibility of representing ice transport across one-cell-wide straits and passages; (c) better representation of inertial-plastic compressive waves. For these reasons, we recommend the use of the C grid EVP formulation of the ice dynamics in future LIM applications.     

Using the seismology of non-magnetic chemically peculiar stars as a probe of dynamical processes in stellar interiors

Chemical composition is a good tracer of the hydrodynamical processes that occur in stars as they often lead to mixing and particle transport. By comparing abundances predicted by models and those observed in stars we can infer some constraints on those mixing processes. As pulsations in the stars are often very sensitive to chemical composition, we can use asteroseismology to probe the internal chemical composition of stars where no direct observations are possible. In this paper I focus on main sequence stars Am, λ Bootis, and HgMn stars and discuss what we can learn of mixing processes in these stars from seismology.     

Predicting the impacts of climate change on genetic diversity in an endangered lizard species

Many endangered species persist as a series of isolated populations, with some populations more genetically diverse than others. If climate change disproportionately threatens the most diverse populations, the species’ ability to adapt (and hence its long-term viability) may be affected more severely than would be apparent by its numerical reduction. In the present study, we combine genetic data with modelling of species distributions under climate change to document this situation in an endangered lizard (Eulamprus leuraensis) from montane southeastern Australia. The species is known from only about 40 isolated swamps. Genetic diversity of lizard populations is greater in some sites than others, presumably reflecting consistently high habitat suitability over evolutionary time. Species distribution modelling suggests that the most genetically diverse populations are the ones most at risk from climate change, so that global warming will erode the species’ genetic variability faster than it curtails the species’ geographic distribution.