Gamma-ray binaries and related systems

After initial claims and a long hiatus, it is now established that several binary stars emit high- (0.1–100 GeV) and very high-energy (>100 GeV) gamma rays. A new class has emerged called “gamma-ray binaries”, since most of their radiated power is emitted beyond 1 MeV. Accreting X-ray binaries, novae and a colliding wind binary (η Car) have also been detected—“related systems” that confirm the ubiquity of particle acceleration in astrophysical sources. Do these systems have anything in common? What drives their high-energy emission? How do the processes involved compare to those in other sources of gamma rays: pulsars, active galactic nuclei, supernova remnants? I review the wealth of observational and theoretical work that have followed these detections, with an emphasis on gamma-ray binaries. I present the current evidence that gamma-ray binaries are driven by rotation-powered pulsars. Binaries are laboratories giving access to different vantage points or physical conditions on a regular timescale as the components revolve on their orbit. I explain the basic ingredients that models of gamma-ray binaries use, the challenges that they currently face, and how they can bring insights into the physics of pulsars. I discuss how gamma-ray emission from microquasars provides a window into the connection between accretion–ejection and acceleration, while η Car and novae raise new questions on the physics of these objects—or on the theory of diffusive shock acceleration. Indeed, explaining the gamma-ray emission from binaries strains our theories of high-energy astrophysical processes, by testing them on scales and in environments that were generally not foreseen, and this is how these detections are most valuable.     

Glass-iron-clay interactions in a radioactive waste geological disposal: An integrated laboratory-scale experiment

Glass-iron-clay setups were reacted at 90 °C for 6-18 months to investigate the coupled interactions between glass alteration, Fe corrosion and clay transformation. The reacted interfaces were probed at the microscopic level using complementary characterization methods (scanning electron microscopy coupled with energy-dispersive X-ray analysis, micro-Raman spectroscopy, micro X-ray diffraction, micro X-ray fluorescence spectroscopy, and micro X-ray absorption near-edge structure spectroscopy). The 10-μm thick Fe foil was fully corroded within 10 months, exposing glass to the pore solution. Iron corrosion led to the formation of a layer containing mostly magnetite, siderite and Fe-rich phyllosilicates with one tetrahedral and one octahedral sheet (TO) or two tetrahedral and one octahedral (TOT) sheet per layer. The clay in contact with this corrosion layer was enriched in siderite (FeCO₃). Glass alteration resulted in the formation of a gel layer whose thickness increased with reaction time (from 20 μm after 6 months to 80 μm after 18 months) and a thin layer of secondary precipitates that concentrated lanthanides, P, and Mo. Assuming conservative behavior of Zr, the Si molar concentration in the gel is about 57% that in the glass. Glass dissolution remained at a rate close to the initial dissolution rate r₀. The data are consistent with glass dissolution sustained by the uptake of dissolved Si and charge-compensating cations on secondary (corrosion) products, thus maintaining the gel porosity open and facilitating the leaching of easily soluble elements.     

Automatic correction and simplification of geological maps and cross-sections for numerical simulations

Incorporating prior geological knowledge in geophysical process models often meets practical meshing challenges and raises the question of how much detail is to be included in the geometric model. We introduce a strategy to automatically repair and simplify geological maps, geological cross-sections and the associated meshes while preserving elementary consistency rules. To identify features breaking validity and/or the thin features potentially problematic when generating a mesh, we associate an exclusion zone with each model feature (horizon, fault). When these zones overlap, both the connectivity and the geometry of the geological layers are automatically modified. The output model enforces specific practical quality criteria on the model topology and geometry that facilitates the generation of a mesh with lower bounds on minimum angles and minimum local entity sizes. Our strategy is demonstrated on an invalid geological cross-section from a real-case study in the Lorraine coal basin. We further explore the impacts of the model modifications on wave propagation simulation. We show that the differences on the seismograms due to model simplifications are relatively small if the magnitude of simplifications is adapted to the physical problem parameters.     

Lateral stirring of large-scale tracer fields by altimetry

Ocean surface fronts and filaments have a strong impact on the global ocean circulation and biogeochemistry. Surface Lagrangian advection with time-evolving altimetric geostrophic velocities can be used to simulate the submesoscale front and filament structures in large-scale tracer fields. We study this technique in the Southern Ocean region south of Tasmania, a domain marked by strong meso- to submesoscale features such as the fronts of the Antarctic Circumpolar Current (ACC). Starting with large-scale surface tracer fields that we stir with altimetric velocities, we determine ‘advected’ fields which compare well with high-resolution in situ or satellite tracer data. We find that fine scales are best represented in a statistical sense after an optimal advection time of ～2 weeks, with enhanced signatures of the ACC fronts and better spectral energy. The technique works best in moderate to high EKE regions where lateral advection dominates. This technique may be used to infer the distribution of unresolved small scales in any physical or biogeochemical surface tracer that is dominated by lateral advection. Submesoscale dynamics also impact the subsurface of the ocean, and the Lagrangian advection at depth shows promising results. Finally, we show that climatological tracer fields computed from the advected large-scale fields display improved fine-scale mean features, such as the ACC fronts, which can be useful in the context of ocean modelling.     

Beach recovery from extreme storm activity during the 2013–14 winter along the Atlantic coast of Europe

The storm sequence of the 2013–14 winter left many beaches along the Atlantic coast of Europe in their most eroded state for decades. Understanding how beaches recover from such extreme events is essential for coastal managers, especially in light of potential regional increases in storminess due to climate change. Here we analyse a unique dataset of decadal beach morphological changes along the west coast of Europe to investigate the post-2013–14 winter recovery. We show that the recovery signature is site specific and multi-annual, with one studied beach fully recovered after 2 years, and the others only partially recovered after 4 years. During the recovery phase, winter waves primarily control the timescales of beach recovery, as energetic winter conditions stall the recovery process whereas moderate winter conditions accelerate it. This inter-annual variability is well correlated with climate indices. On exposed beaches, an equilibrium model showed significant skill in reproducing the post-storm recovery and thus can be used to investigate the recovery process in more detail. © 2018 John Wiley & Sons, Ltd.     

Signification des sables tertiaires de l’Ardenne méridionale : du nouveau sur la capture de la Meuse lorraine

Résumé

La question de la capture de la Meuse lorraine est abordée par le biais de l’étude des sables tertiaires qui jalonnent les plateaux de l’Ardenne méridionale. Toujours polygéniques, ces sables ont une histoire plus ou moins longue qui débute à l’Éocène où l’altération poussée du socle produit un épais manteau sableux. Certains gisements in situ fossilisés sont restés quasiment intacts (Régniowez) tandis que d’autres (Doisches) portent la marque de remaniements successifs (marins et/ou continentaux) qui se suivent jusqu’au Miocène sans que leur filiation avec le matériau d’origine, attestée par la nature de leurs barylites, ne se perde. Comme les sables fluviatiles miocènes (in fine) du plateau ardennais ne portent pas la signature minéralogique d’une Meuse lorraine, à l’inverse de la Traînée mosane (datée du Miocène moyen), la capture se trouve à la fois confirmée et datée. L’ancien réseau de la Meuse lorraine, quant à lui, n’a pu se mettre en place qu’après l’Oligocène, époque à laquelle un glacis incliné de l’Ardenne vers le Bassin parisien fonctionnait encore. © 2000 Éditions scientifiques et médicales Elsevier SAS     

Atmospheric response to the North Atlantic Ocean variability on seasonal to decadal time scales

The NCEP twentieth century reanalyis and a 500-year control simulation with the IPSL-CM5 climate model are used to assess the influence of ocean-atmosphere coupling in the North Atlantic region at seasonal to decadal time scales. At the seasonal scale, the air-sea interaction patterns are similar in the model and observations. In both, a statistically significant summer sea surface temperature (SST) anomaly with a horseshoe shape leads an atmospheric signal that resembles the North Atlantic Oscillation (NAO) during the winter. The air-sea interactions in the model thus seem realistic, although the amplitude of the atmospheric signal is half that observed, and it is detected throughout the cold season, while it is significant only in late fall and early winter in the observations. In both model and observations, the North Atlantic horseshoe SST anomaly pattern is in part generated by the spring and summer internal atmospheric variability. In the model, the influence of the ocean dynamics can be assessed and is found to contribute to the SST anomaly, in particular at the decadal scale. Indeed, the North Atlantic SST anomalies that follow an intensification of the Atlantic meridional overturning circulation (AMOC) by about 9 years, or an intensification of a clockwise intergyre gyre in the Atlantic Ocean by 6 years, resemble the horseshoe pattern, and are also similar to the model Atlantic Multidecadal Oscillation (AMO). As the AMOC is shown to have a significant impact on the winter NAO, most strongly when it leads by 9 years, the decadal interactions in the model are consistent with the seasonal analysis. In the observations, there is also a strong correlation between the AMO and the SST horseshoe pattern that influences the NAO. The analogy with the coupled model suggests that the natural variability of the AMOC and the gyre circulation might influence the climate of the North Atlantic region at the decadal scale.     

The Multiscale TROPIcal CatchmentS critical zone observatory M-TROPICS dataset II: Land use,hydrology and sediment production monitoring in Houay Pano,northern Lao PDR

Mountain regions of the humid tropics are characterized by steep slopes and heavy rains. These regions are thus prone to both high surface runoff and soil erosion. In Southeast Asia, uplands are also subject to rapid land-use change, predominantly as a result of increased population pressure and market forces. Since 1998, the Houay Pano site, located in northern Lao PDR (19.85°N 102.17°E) within the Mekong basin, aims at assessing the long-term impact of the conversion of traditional slash-and-burn cultivation systems to commercial perennial monocultures such as teak tree plantations, on the catchment hydrological response and sediment yield. The instrumented site monitors hydro-meteorological and soil loss parameters at both microplot (1 m²) and small catchment (0.6 km²) scales. The monitored catchment is part of the network of critical zone observatories named Multiscale TROPIcal CatchmentS (M-TROPICS). The data shared by M-TROPICS in Houay Pano are (1) rainfall, (2) air temperature, air relative humidity, wind speed, and global radiation, (3) catchment land use, (4) stream water level, suspended particulate matter, bed particulate matter and stones, (5) soil surface features, and (6) soil surface runoff and soil detachment. The dataset has already been used to interpret suspended particulate matter and bed particulate matter sources and dynamics, to assess the impact of land-use change on catchment hydrology, soil erosion, and sediment yields, to understand bacteria fate and weed seed transport across the catchment, and to build catchment-scale models focused on hydrology and water quality issues. The dataset may be further used to, for example, assess the role of headwater catchments in large tropical river basin hydrology, support the interpretation of new variables measured in the catchment (e.g., contaminants other than faecal bacteria), and assess the relative impacts of both climate and land-use change on the catchment.     

A semi-empirical wind set-up forecasting model for Lake Champlain

The precision of Lake Champlain's water level estimation is a key component in the flood forecasting process for the Richelieu River. Hydrological models do not typically take into consideration the effects of the wind on the water level (also known as the wind set-up). The objective of this study is to create an empirical wind set-up forecast model for Lake Champlain during high wind events. The proposed model uses wind speed and direction across the Lake, as well as wind gusts as inputs. The model is calibrated to a subset of observations and evaluated on an independent sample, considering four wind speed bins. It is tested and compared to a variant of the Zuider Zee equation on 20 wind set-up events that occurred between 2017 and 2019 using hindcast data from five different numerical weather prediction systems (GDPS, RDPS, HRDPS, NOAA and ECMWF). A quantile mapping-based forecast calibration scheme is implemented for each of the forecast products to correct their biases. Results show that events are successfully predicted by the proposed model at least 72 h in advance. These results are better than the other comparative models found in the literature and tested herein. Overall, significant improvements are obtained by including wind speed and wind gusts from different weather stations.     

Contribution of magmatism,partial melting buffering and localized crustal thinning on the late Variscan thermal structure of the Agly massif (French Pyrenees)

Low pressure-high temperature (LPHT) metamorphism, with geothermal gradients in the order of 50–100°C/km, is a common feature of the late evolution of collisional orogens. These abnormal thermal conditions may be the results of complex interactions between magmatism, metamorphism and deformation. The Agly massif, in the French Pyrenees, preserves the metamorphic footprints of the late Variscan thermal structure of an almost continuous section from the upper and middle continental crust. The upper crust is characterized by a very high geothermal gradient of ~55°C/km, evolving from greenschist to amphibolite facies, while the middle crust, exposed in a gneissic core, exhibits granulite facies conditions with a near isothermal geothermal gradient (<8°C/km) between 740 and 790°C. The abnormal and discontinuous crustal geothermal gradient, dated at c. 305 Ma on syn-granulitic monazite by LA-ICP-MS, is interpreted to be the result of magmatic intrusions at different structural levels in the crust: the Ansignan charnockite (c. 305 Ma) in the deepest part of the gneissic core, the Tournefort granodiorite (c. 308 Ma) at the interface between the gneissic core and the upper crust and the Saint-Arnac granite (c. 304 Ma) in the upper section of the massif. The heat input from these magmas combined with the thermal buffering effect of the biotite dehydration-melting reaction resulted in the near isothermal geothermal gradient in the gneissic core (melt-enhanced geotherm). The higher geothermal gradient (>50°C/km) in the upper crust is only due to conduction between the hot middle crust and the Earth's surface. The estimated maximum finite pressure range suggests that ~10 to 12 km of crust are exposed in the Agly massif while the present-day thickness does not exceed 5–6 km. This pressure/depth gap is consistent with the presence of several normal mylonitic shear zones that could have contributed to the subtraction of ~5 km of the rock pile. Monazite U–Th–Pb ages carried out on monazite overgrowths from a highly mylonitized sample suggest that this vertical thinning of the massif occurred at c. 296–300 Ma. This later Variscan extension might have slightly perturbed the 305 Ma geothermal gradient, resulting in an apparent higher conductive geothermal gradient in the upper crust. Although the Agly massif has been affected by Cretaceous extension and Eocene Alpine compression, we suggest that most of the present-day thickness of the column rock was acquired by the end of the Palaeozoic.